Страница 2 из 2

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

Добавлено: 16 июл 2012, 06:51
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 как сейчас. Вроде обещают через пару недель включить в транк. Прилагаю пропатченную версию, если у кого есть желание - можете протестировать.

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

Добавлено: 16 июл 2012, 13:13
out2002
Ух ты, спасибо за патчик и информацию. Буду тестить.

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

Добавлено: 20 июл 2012, 05:22
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