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

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

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

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

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

gdal georeferencing

— ДАК
источник

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

Найдите свой пространственный код ссылки из здесь
Узнайте вашу левую верхнюю координату и нижнюю правую координаты.
используйте эту команду для согласования:
- gdal_translate -of GTiff -a_ullr ullon ullat lrlon lrlat -a_srs EPSG:4269 input.tif output.tif
используйте эту команду для mercator: (epsg: 3857)
- gdalwarp -of GTiff -t_srs EPSG:3857 input.tif output.tif
затем используйте gdal2tiles.py для мозаики.

— Арагон
источник

Спасибо за ответ. Но мне нужно наложить мою проекцию UTM на Google Earth, которая использует WGS84 в качестве базы данных. Так не должен ли я установить целевую систему отсчета на ESPG: 4326?

— Акд

Gdalwarp @Akd Арагона в EPSG: 3857 хорошо для Google. Смотрите wiki.openstreetmap.org/wiki/EPSG:3857 для получения дополнительной информации.

— Дэйв Икс

В соответствии с документацией и тестом: -a_ullr ulx uly lrx lryили с лат / лон терминами -a_ullr ullon ullat lrlon lrlatэто правильный способ использования команды.

— Зета

Я отредактировал пост для правильного порядка х и у в -a_ullr.

— Никос Александрис

Я также пытаюсь привязать географическую привязку, но мне интересно, как мы можем узнать верхнюю левую координату и нижнюю правую координату? Не могли бы вы помочь мне @Aragon?

— Manimalis