Изменение проекции изображения с помощью Python - PullRequest
Новая
       53

Изменение проекции изображения с помощью Python

1 голос
/ 31 мая 2011

Координаты угла изображения известны:

  • Вверху слева: (70.20315, 20.92749)
  • Вверху справа: (70.20315, 31.45240)
  • Внизу слева: (59,70571, 20,92749)
  • Внизу справа: (59,70571, 31,45240)

(в основном карта Финляндии и немного Швеции / Норвегии / России)

Я не смог найти способ изменить проекцию изображения на оверлейное изображение карт Google, желательно с использованием Python.

Это можно сделать с помощью pyproj и PIL , но я не смог понять, как это сделать. Я также нашел GDAL , а точнее gdalwarp , но не смог найти варианты для правильного изменения проекции.

Example image

1 Ответ

0 голосов
/ 04 июня 2011

(Копировать-вставить из gis answer)

Обширный метод проб и ошибок показал, что он работает нормально, без указания исходной проекции.То же самое может быть достигнуто с помощью только GDAL Python API.Шаги для воспроизведения:

Добавление контрольных точек к TIF (конечно, сначала convert mapimage.png mapimage.tif): 

gdal_translate mapimage.tif mapimage_georeffed.tif -gcp 2300 2300 31.45240 59.70571 -gcp 0 0 20.92749 70.20315 -gcp 2300 0 31.45240 70.20315 -gcp 0 2300 20.92749 59.70571 -a_srs +proj=longlat

Затем, привязки Python GDAL включают gdal2tiles.py, что автоматически генерирует плитки на правой проекциидля карт, а также образцы страниц карты.

Для справки, я не совсем уверен, что проекция является совершенно правильной, но в этой области она, по крайней мере, достаточно хороша.

...