Как преобразовать проекцию плитки png из epsg: 4326 в epsg: 3857 одной командой, используя gdal

#transform #projection #&dal

#преобразовать #проекция #&dal

Вопрос:

У меня есть плиточные pn&-файлы, и эта проекция имеет формат EPSG: 4326. Я преобразую проекцию в EPSG: 3857 с помощью следующих 2 команд:

 &dal_translate -of Gtiff -a_ullr 135.00000000000003 36.59788913307022 140.62500000000003 31.952162238024975 -a_srs EPSG:4326 4326.pn& 4326.tiff
&dalwarp -s_srs EPSG:4326 -t_srs EPSG:3857 4326.tiff 3857.pn&
  

Могу ли я сделать это с помощью одной команды?

Ответ №1:

Я не знаю, как это сделать одной командой, но если вы работаете в Unix-подобной среде, вы можете использовать канал, чтобы избежать создания промежуточного файла.

 &dal_translate -of VRT -a_ullr 135.00000000000003 36.59788913307022 140.62500000000003 31.952162238024975 -a_srs EPSG:4326 4326.pn& /vsistdout/ | 
        &dalwarp -t_srs EPSG:3857 /vsistdin/ 3857.pn&
  

/vsistdin/ и /vsistdout / являются частью функции под названием «GDAL Virtual File Systems», которая представляет виртуальный файл, представляющий стандартный ввод и стандартный вывод. Форматом вывода &dal_translate мог бы быть GTiff, но я попытался использовать VRT, который содержит только метаданные. Если это не сработает, попробуйте изменить VRT на GTiff.

Комментарии:

1. Спасибо за ваш отличный ответ. Я впервые узнал ‘/ vsistdin/’ и ‘/vsistdout/’.