Блог

Я автоматически привязываю многие исторические изображения к географической привязке, используя подход CV. У меня есть список GCPS и файл TIF без ссылок, который я копирую и применяю к GCPs с помощью функции GDALs SetGCPs (). Это TIF-файл с привязкой к географической привязке, созданный SetGCPs():

Но если я сохраню те же GCPS в текстовый файл и вручную импортирую/применю их к изображению без ссылок, используя полиномиальное преобразование 2-го порядка, то повторное изображение будет: (вы можете увидеть предыдущее изображение немного ниже). Мой вопрос в том, как я могу изменить метод преобразования, используемый функцией SetGCPs ()? Я бы хотел, чтобы это всегда было полиномиальное преобразование 2-го порядка.

Пример кода:

 from osgeo import gdal, osr

sr = osr.SpatialReference()
sr.ImportFromEPSG(26986)

ds = gdal.Open(output_filename, gdal.GA_Update)
ds.SetGCPs(gcps, sr.ExportToWkt())
 

Где gcps-это список объектов gdal GCP.

Спасибо!