Форумы GIS-Lab.info

При перепроецировании слоев глобального охвата, во многих случаях возникают артефакты (см. приложение). Кто-нибудь знает, как это можно побороть?
В приложении слой береговой линии из NaturalEarth (ГСК) перепроецирован в World_Van_der_Grinten_I (EPSG:54029)

Предварительно разрежьте объекты по меридиану, на 180° отстоящему от осевого. По-хорошему, надо бы после этого куски, вылезающие за 180° от осевого меридиана, перенести на 360°.
Это решает проблему частично, поскольку один из краёв разреза будет восприниматься как чужой (скажем, точки на долготе +180° программа будет переносить на долготу -180°). Для окончательного решения придётся чуть-чуть сместить этот край внутрь.
Кроме того, полезно на края разреза добавить дополнительные точки, чтобы сохранить линию меридиана.

А можно ли описанный процесс хоть частично автоматизировать, например, с помощью OGR? Можно ли с пом. OGR разрезать слой по заданному меридиану? Гуглом нашел, что вроде как нужно использовать функцию -wrapdateline, но с пом. нее я результата добиться не смог. Артефактные линии лишь изменяют свое положение:

На всякий случай приложу файл для экспериментов. Есть подозрение, что там данные выходят за диапазон -180, 180.

Если резать ogr2ogr, то не так. Однако в данном случае всё уже порезано, как надо. Пришлось лишь один отрезочек подправить, чтобы не залезал за 180.
QGIS-1.8.0, проекция Всё отображается зер гут.
Возможны проблемы при смене датума. Тогда пересчитанные долготы могут вылезать за границу 180 градусов.

ErnieBoyd писал(а):Пришлось лишь один отрезочек подправить, чтобы не залезал за 180.

Можно поподробнее? Как вы этот отрезок нашли и что именно с ним сделали?

Подтверждаю то, что написал ErnieBoyd. Я немного подтянул на запад несколько вершин в самой восточной части (около 180 град.) — стало перепроецироваться "на лету" без артефактов.

bolotoved писал(а):Можно поподробнее? Как вы этот отрезок нашли и что именно с ним сделали?

Сохранил слой как файл MapInfo с расширением ".mif". В текстовом редакторе искал строку "180", нашёл "180.00000044181", исправил на "179.99998". Были ещё долготы 180 ровно, но я их не трогал.

Спасибо, все понял. Значит проблема была в том, что данные выходили за 180 меридиан.

bolotoved писал(а):Значит проблема была в том, что данные выходили за 180 меридиан.

Это "во-вторых". Во-первых, проблема в том, что QGIS сворачивает узлы объектов, вылезающие за 180 градусов от осевого меридиана.
Странная идеология. Могли бы оставить картографу право решать, где и как резать свою карту. Заодно таких конфузов не возникало бы.

Странно всё это. У меня такие же проблемы, всегда, бесят. Нет ли чего-то типа фильтра для отображения, чтобы округлять с отбрасыванием дробной части до +-180 градусов?

Странник писал(а):Нет ли чего-то типа фильтра для отображения, чтобы округлять с отбрасыванием дробной части до +-180 градусов?

Помогает точная обрезка с пом. ogr2ogr: