Геопривязка с использованием GDAL?

16

У меня есть изображение .tiff, которое я пытаюсь наложить на Google Планета Земля. Для этого мне нужно сначала привязать его. Проекция изображения - UTM с датой NAD83. Я знаю широту / длину четырех углов изображения.

Я делаю следующим образом:

  1. gdalinfo, чтобы получить четыре координаты.
  2. gdal_translate для указания широты / долготы
  3. gdal2tiles.py для мозаичного изображения.

Единственное, что меня беспокоит, так это то, что я никогда не указывал UTM и NAD83, тогда почему он может наложить мой образ? Кто-нибудь знает, как указать эту информацию проекции?

gdal georeferencing ДАК
источник

Ответы:

23

следуйте по этому пути:

  1. Найдите свой пространственный код ссылки из здесь

  2. Узнайте вашу левую верхнюю координату и нижнюю правую координаты.

  3. используйте эту команду для согласования:

    • gdal_translate -of GTiff -a_ullr ullon ullat lrlon lrlat -a_srs EPSG:4269 input.tif output.tif

  4. используйте эту команду для mercator: (epsg: 3857)

    • gdalwarp -of GTiff -t_srs EPSG:3857 input.tif output.tif

  5. затем используйте gdal2tiles.py для мозаики.

Арагон
источник
1
Спасибо за ответ. Но мне нужно наложить мою проекцию UTM на Google Earth, которая использует WGS84 в качестве базы данных. Так не должен ли я установить целевую систему отсчета на ESPG: 4326?
Акд
1
Gdalwarp @Akd Арагона в EPSG: 3857 хорошо для Google. Смотрите wiki.openstreetmap.org/wiki/EPSG:3857 для получения дополнительной информации.
Дэйв Икс
1
В соответствии с документацией и тестом: -a_ullr ulx uly lrx lryили с лат / лон терминами -a_ullr ullon ullat lrlon lrlatэто правильный способ использования команды.
Зета
2
Я отредактировал пост для правильного порядка х и у в -a_ullr.
Никос Александрис
2
Я также пытаюсь привязать географическую привязку, но мне интересно, как мы можем узнать верхнюю левую координату и нижнюю правую координату? Не могли бы вы помочь мне @Aragon?
Manimalis