Смешения тайликов относительно векторной информации

[Denis Rykov]

В растре информация о СК не представлена, вытаскиваем её из шейп-файла:

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

gdalsrsinfo -v -o proj4 forum.prj
Validate Succeeds
+proj=tmerc +lat_0=0 +lon_0=29.5 +k=1 +x_0=3300000 +y_0=-9214.688 +ellps=krass +units=m +no_defs

И используя полученное описание проекции, конвертируем растр:

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

gdalwarp -t_srs "EPSG:4326" -s_srs "+proj=tmerc +lat_0=0 +lon_0=29.5 +k=1 +x_0=3300000 +y_0=-9214.688 +ellps=krass +units=m +no_defs" X-39-55-D-d-4.tif X-39-55-D-d-4_4326.tif

Затем осуществляем разбивку на тайлы:

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

gdal2tiles X-39-55-D-d-4_4326.tif ./tiles

Если вы на Windows, то проще всего установить GDAL, используя OSGeo4W, запустить OSGeo4W.bat, который импортирует переменные окружения, и уже после этого все гдаловские утилиты можно запускать как приведено выше.

Ещё, я сделал небольшой патч к gdal2tiles, который генерирует OpenLayers проект на базе созданных тайлов, используя последнюю версию 2.12, а не 2.7 как сейчас. Вроде обещают через пару недель включить в транк. Прилагаю пропатченную версию, если у кого есть желание - можете протестировать.

[out2002]

Ух ты, спасибо за патчик и информацию. Буду тестить.

[Denis Rykov]

Можно и без предварительной перегонки в 4326 потайлить так:

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

gdal2tiles.py -s "+proj=tmerc +lat_0=0 +lon_0=29.5 +k=1 +x_0=3300000 +y_0=-9214.688 +ellps=krass +units=m +no_defs" X-39-55-D-d-4.tif ./tiles