Перепроецирование с помощью GDAL

[Gerhald]

Все таки дело, видимо, в этом файле привязке.

[alexandr cherepanov]

Если уверены что это 3395 тогда должно быть примерно так

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

2445.9849051256103515625
0.0
0.0
-2440.78067014107666015625
-20036285.35033643719482421875
19993654.859460629461669921875

Если все таки 3857, то

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

2445.9849051256103515625
0
0
-2445.9849051256103515625
-20036285.35033643719482421875
20036285.35033643719482421875

[Gerhald]

В обоих случаях выдает ошибку:

ERROR 1: Unable to compute a GEOTRANSFORM based transformation between pixel/lin
e and georeferenced coordinates for map.tif.

Может нужно этот файл привязки как-то указать в GDAL? Как он его распознает вообще?

[Gerhald]

А вот, посмотрел какая проекция была.
Я качал карту с SAS.Planet и там по умолчанию стояла такая проекция:

Projection of map - Mercator / Google maps (Sphere Radius 6378137) /EPSG:3785.

Обманул вас, наверное. Извиняюсь.

[alexandr cherepanov]

Создаете файл с именем map.tfw в той же папке, что файл map.tif
Вставляете туда параметры привязки для 3857.

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

2445.9849051256103515625
0
0
-2445.9849051256103515625
-20036285.35033643719482421875
20036285.35033643719482421875

Упрощаем команду до минимума
gdalwarp -s_srs EPSG:3857 -t_srs EPSG:4326 -r lanczos -co COMPRESS=LZW map.tif map_4326.tif

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

После создания файла привязки map.tfw нужно выполнить команду

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

gdal_edit.py -a_srs EPSG:3395 map.tif

После этого привязка и система координат окажутся внутри map.tif, и не нужно никаких экспериментов.

[Gerhald]

То же самое...

ERROR 1: The transformation is already "north up" or a transformation between pi
xel/line and georeferenced coordinates cannot be computed for map.tif. There is
no affine transformation and no GCPs. Specify transformation option SRC_METHOD=N
O_GEOTRANSFORM to bypass this check.

Все таки, какую СК писать то?

Попробовал так:
1) gdal_edit.py -a_srs EPSG:3395 map.tif
gdalwarp -s_srs EPSG:3395 -t_srs EPSG:4326 -r lanczos -co COMPRESS=LZW map.tif map_4326.tif
Не работает.
2) gdal_edit.py -a_srs EPSG:3857 map.tif
gdalwarp -s_srs EPSG:3857 -t_srs EPSG:4326 -r lanczos -co COMPRESS=LZW map.tif map_4326.tif
Не работает

[Gerhald]

Причем эти команды никаких ошибок не выдают:
gdal_edit.py -a_srs EPSG:3395 map.tif
gdal_edit.py -a_srs EPSG:3857 map.tif

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

Каюсь, виноват, был не в себе.
В файле привязки map.tfw последнее число должно быть положительным:

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

2445.98490512564
0
0
-2445.98490512564
-20036285.350336737
20036285.350336737

Регистрируем:

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

gdal_edit.py -a_srs EPSG:3395 map.tif

Перепроецируем:

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

 gdalwarp -t_srs EPSG:4326 map.tif map_4326.tif

Всё должно быть OK.

[Gerhald]

Все равно не делает...

ERROR 1: The transformation is already "north up" or a transformation between pi
xel/line and georeferenced coordinates cannot be computed for map.tif. There is
no affine transformation and no GCPs. Specify transformation option SRC_METHOD=N
O_GEOTRANSFORM to bypass this check.

Хоть ты застрелись...

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

Зачем же так мрачно! Выложите, наконец, картинку.

[Gerhald]

Да обычная картинка. Вот уменьшенный вариант, полный много весит...

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

Это просто. Превратим JPEG в GeoTIFF:

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

gdal_translate -of GTiff -a_srs EPSG:3395 -gcp 0 0 -20037508.34 20037508.34 -gcp 3000 0 20037508.34 20037508.34 -gcp 0 3000 -20037508.34 -20037508.34 map_resize.jpg map_resize.tif

Трансформируем в географическую WGS 84:

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

gdalwarp -t_srs EPSG:4326 map_resize.tif map_resize_4326.tif

[Gerhald]

О, здорово!
Теперь преобразует, но неправильно. Видимо из-за того, что у меня размеры оригинальной карты 16384x16384 пикселя. Как правильно теперь прописать координаты для такого изображения?
Объясните в двух словах, как распознавать эти цифры:

-gcp 0 0 -20037508.34 20037508.34 -gcp 3000 0 20037508.34 20037508.34 -gcp 0 3000 -20037508.34 -20037508.34

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

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

gdal_translate -of GTiff -a_srs EPSG:3395 -gcp 0 0 -20037508.34 20037508.34 -gcp 16384 0 20037508.34 20037508.34 -gcp 0 16384 -20037508.34 -20037508.34 map_resize.jpg map_resize.tif

Опция -gcp задаёт точку привязки: пара координат в пикселях растра и пара координат на глобусе.
RTFM