Некорректная привязка и/или трансформация

[wowka1319]

Есть растр в PNG и набор lat/lon по эллипсоиду GRS80. Нужно получить GeoTiff в нужной проекции (например WGS84/Меркатор).
Сначала я привязываю растр так:

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

gdal_translate -gcp X Y lon lat -gcp X Y lon lat {и еще много точек} src.png img_with_latlon.tif

Далее получаю проекцию:

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

gdalwarp -s_srs "EPSG:4019" -t_srs "EPSG:3395" img_with_latlon.tif result_epsg3395.tif

-s_srs "EPSG:4019" - тут я указал, что исходный эллипсоид GRS80 (не уверен, что это верный датум)
-t_srs "EPSG:3395" - это меркатор/WGS84 - как я понимаю это уже не очень важно ибо подойдет любая проекция, которая известна целевой ГИС.

В итоге получаю результат в котором присутствует сильное смещение.
Что может быть не так?

Причем не понятно поведение gdalwarp: от параметра почему-то -t_srs ничего не меняется

[wowka1319]

Подозреваю, что gdalwarp кривой (по крайней мере v1.11). Похоже параметр -t_srs не работает, так как изменение его значения ничего не меняет. Задавая те же параметры через код (работаю в C++), все ок.
Задача решена.