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

перепроецирование Vmap0 в EPSG:3388 с помощью ogr2ogr

Добавлено: 01 фев 2011, 14:17
ponimas
подскажите, как правильно выполнить репроекцию vmap0 в EPSG:3388 Pulkovo 1942 / Caspian Sea Mercator?
Сейчас делаю так

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

ogr2ogr -f "ESRI Shapefile" -t_srs "EPSG:3388" -skipfailures -s_srs "EPSG:4326" new.shp old.shp
в итоге имею вот такую картину http://funkyimg.com/u2/513/819/2011-02- ... _scrot.png
полигоны вместо того чтоб быть разрезанными, теперь растянуты по карте.
как это исправить?

Re: перепроецирование Vmap0 в EPSG:3388 с помощью ogr2ogr

Добавлено: 01 фев 2011, 15:32
bms
EPSG:3388 - это проекция Меркатора, EPSG:4326 - географическая проекция. Так что после перевода из одной проекции в другую карта "деформируется" http://gis-lab.info/qa/proj-sk-faq.html#20

Re: перепроецирование Vmap0 в EPSG:3388 с помощью ogr2ogr

Добавлено: 01 фев 2011, 15:35
ponimas
EPSG:3388 - это проекция Меркатора, EPSG:4326 - географическая проекция. Так что после перевода из одной проекции в другую карта "деформируется"
я в курсе.
взгляните на мой скриншот хотя бы для начала, что б увидеть, что меня не это интерисует.

Re: перепроецирование Vmap0 в EPSG:3388 с помощью ogr2ogr

Добавлено: 01 фев 2011, 15:37
KolesovDmitry
ponimas писал(а):
взгляните на мой скриншот
Лучше вложите его сюда, а то ваша ссылка ведет на какую-то подозрительную страницу (во всяком случае, картинки там нет)

Re: перепроецирование Vmap0 в EPSG:3388 с помощью ogr2ogr

Добавлено: 01 фев 2011, 15:45
ponimas
по-моему ссылка ведёт напрямую на пнг-файл.
присоединять файлы к сообщениям возможности не вижу, поэтому положи файл на другой хостинг.
http://ompldr.org/vNzhycA/2011-02-01-15 ... _scrot.png

Re: перепроецирование Vmap0 в EPSG:3388 с помощью ogr2ogr

Добавлено: 01 фев 2011, 18:31
KolesovDmitry
Скажите, а какой у вас охват территории? Смотрю на картинку -- у меня такого рода ужасы получались, когда я перепроецировал данные, пересекающие 180-й мередиан?

Re: перепроецирование Vmap0 в EPSG:3388 с помощью ogr2ogr

Добавлено: 01 фев 2011, 18:42
ponimas
разобрался.
нужно было сначала просто шейп разрезать тем же ogr2ogr и только потом препроецировать.
когда я перепроецировал данные, пересекающие 180-й мередиан?
дело в том, что данная проекция использует lon_0=51

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

+proj=merc +lon_0=51 +k=1 +x_0=0 +y_0=0 +ellps=krass +units=m +no_defs 
соответственно идёт смещение края отображаемой територии.
и такая каша происходит с данными пересекающими не 180, а -129 меридиан, если я правильно всё понимаю.
после того, как нарезал - всё стало прекрасно отображаться.
http://picthost.ru/images/j1vonn2zgf1nu91hwbyp.png