Не правильно проходит перепроецирование

[Syllion]

Мне нужно положить снимок сделанный в проекции 4326 на эллипсе на карту. Для карт-основы использую leaflet (который как я понимаю данные берет с OSM). Исходя из этого на выходе мне проекция нужна EPSG:3857

Проблема в том, что у меня почему-то перепроецированные данные отображаются на карте со смещением т.е. вот так:


Использую следующие команды:
gdalwarp -t_srs EPSG:3857 -srcnodata 0 <tif1> <tif2>
gdal_translate -of PNG <tif2> <png>

В чем может быть проблема?

В QGIS результирующий png почему-то ложится нормально:

[SergeyRyzhkov]

tif1 и tif2.
В первой команде не перепутан порядок?

[ Сообщение с мобильного устройства ]

[Syllion]

Вроде бы нет. Я там добавил пример того как те же данные qgis отображает.

[SergeyRyzhkov]

Вроде бы нет. Я там добавил пример того как те же данные qgis отображает.

У вас все правильно. Просто где-было что первым идет target первым, а вторым source
Только поможет пример растра. Растр в студию!

[Syllion]

Прикрепил файл.

Вот полный список используемых команд:
1) Перевод из с 16 в 8 бит с рескейлом изображения
gdal_translate -ot byte -scale <img> <tif1>

2) Перепроецирование изображения с заданием значения NoData
gdalwarp -t_srs EPSG:3857 -srcnodata 0 <tif1> <tif2>

3) Перевод в PNG
gdal_translate -of PNG <tif2> <png>

[SergeyRyzhkov]

Что-то в архиве уже результат, а не исходники (либо я не понял, ссори).
Также желательно сразу показать JS для Leaflet, который Вы используете для визуализации растра

[Syllion]

В архиве данные в проекции 4326. Нужно их трансформировать так, чтобы они легли на карту. Трансформировать нужно в 3857. Пока в северных широтах вылезает вот такой сдвиг, как я на скриншоте привел.

[trir]

Ну если QGIS показывает правильно, то может цепляете неправильно?

[Syllion]

trir, ой вы правы. А не подскажите как мне координаты из 4326 в 3857 трансформировать? Именно вот есть строка долготы и широты. Нужна точно такая же но трансформированная для 3857

[Игорь Белов]

А не подскажите как мне координаты из 4326 в 3857 трансформировать?

gdaltransform

Вы уверены, что именно этого хотите? Растр в Вашем архиве в проекции UTM.

[Syllion]

Раст я перевожу в PNG вроде как правильно, но у меня есть отдельный специальный файл с метаданными, которые записаны в UTM. Их то и надо пересчитать.

По поводу gdaltransform. После запуска примера вот что получается:
> gdaltransform -s_srs EPSG:28992 -t_srs EPSG:31370 177502 311865
ERROR 4: `177502' does not exist in the file system,
and is not recognised as a supported dataset name.

что ему не нравится то? Это ж вроде тупо копи-паста.

[Игорь Белов]

Вот копипаста:

Код: Выделить всё

> gdaltransform -s_srs EPSG:28992 -t_srs EPSG:31370
> 177502 311865

А в Вашем варианте программа считает, что "177502" и "311865" - названия входного и выходного растров.