Использование GDAL для привязки растровых материалов

Библиотека GDAL - очень полезный, бесплатно распространяемый набор инструментов и библиотек с открытым исходным кодом. Использование этих инструментов это один из немногих способов осуществить географическую привязку растрового изображения используя исключительно бесплатное (и легальное) ПО.

Привязка растровых изображений, которые кроме файловой (пиксельной) другой системы координат не имеют (например отсканированных изображений) проводится в два шага. Если растр уже имеет географическую или спроектированную систему координат и ее нужно изменить на другую то следует переходить сразу к шагу 2.

Назначение контрольных точек:

На этом этапе, информация о том, какие конечные координаты (Lat, Long) должен получить определенный пиксел (имеющий пиксельные координаты X,Y). Конечные координаты не обязательно должны быть широтой долготой, они также могут быть иметь и спроектированную прямоугольную систему координат:

Как правило, прописывается несколько точек, для трансформации на втором этапе необходимо минимум 2 точки, следующим образом (все контрольные точки прописываются одной командой):

Проверить, правильно ли прописались контрольные точки можно с помощью listgeo:

В случае если операция осуществлена корректно, в тэгах GeoTIFF должна появиться информация о контрольных точках:

Трансформация растра

Обычно на втором этапе производится трансформация растра, если растр уже имеет географическую или спроектированную систему координат и ее нужно изменить на другую, то первый шаг не нужен:

-s_rsr - исходная система координат, в формате строки описания системы координат PROJ. Исходная система координат должна соответствовать системе координат введенных на первом этапе контрольных точек.

-t_srs - конечная система координат, если необходимо только трансформировать растр в систему координат контрольных точек, то она должна быть равна исходной.

либо можно ее вообще не указывать, опять же, в этом случае будет просто произведена аффинная или полиномиальная трансформация растра (подробнее), которые к проективным преобразованиям отношения не имеют:

Примечания:

• Названия файлов не рекомендуется начинать с цифр
• По умолчанию для трансформации используется полиномиальное преобразование, причем его порядок определяется количеством точек, что не всегда удобно, так как если у вас много контрольных точек, то выбрано будет преобразование высокого порядка, которое сильно деформирует растр, что в большинстве случаев неприемлемо. Чтобы установить определенный порядок полиномиального преобразования, нужно использовать ключ -order n, где n - порядок полинома.
• GDAL также позволяет использовать преобразование thin-plate spline, для этого нужно использовать ключ -tps.

Благодарим за идею Григория Морозова

Обсудить в форуме (Комментариев - 0)

Последнее обновление: March 29 2008
(Наверх)