Как обобщить SpatialPolygons и сохранить топологию в R?

18

Я хотел бы обобщить SpatialPolygonsDataFrame в R, чтобы сделать построение графиков более эффективным. Полигоны, которые я имею, приводят к огромным pdf-графикам, которые медленно загружаются и трудно интегрируются в документ.

Я попытался использовать gSimplifyфункцию из rgeosпакета, но, к сожалению, она не сохраняет топологию полигонов.

Есть ли альтернатива, которую я могу использовать?

r topology generalization

— yellowcap
источник

не совсем решение, но вы могли бы вместо этого построить растр (jpg / png)? Этот вопрос stackoverflow, который обсуждает уменьшение размера pdfs, может быть

— полезен

Построение карт в виде растровых изображений является опцией, но по возможности я бы хотел использовать векторную графику. Качество изображения лучше при использовании векторной графики, особенно при просмотре цифровых версий документов.

— желтая шапка

gSimplify()имеет необязательный параметр, topologyPreserveкоторый по умолчанию равен FALSE. Вы пытались установить это TRUE?

— krlmlr

7

Об этом идет дискуссия на r-sig-geo . Для окончательного ответа вы должны спросить там, потому что есть люди, которые знают идеи пространственного R.

Но вы также можете сделать это в настольных ГИС-приложениях (экспортируйте фигуру с помощью команды writeOGR из rgdal или writePolyShape () из maptools), например, QuantumGIS, GRASS или SAGA.

Для QuantumGIS используйте Vector / Geometry Tools / Упростить геометрию (я протестировал и не сохраняю топологию, но применительно к шейп-файлу admin1 в Румынии, отлично смотрится с 2000 узлами).

Для GRASS используйте v.generalize (см. Руководство для получения информации об алгоритмах, некоторые из них).

Для SAGA вы должны сначала преобразовать полигон в линию (Shapes - Lines / Convert Polygons to Lines), затем упростить линии (Shapes Lines / Line Simplification) и, наконец, преобразовать линии в полигоны (Shapes - Polygons / Convert Lines в Polygons). Я проверил это, и результаты сохранили топологию.

Ура и удачи,

— Никулита Михай
источник

Тем временем я также нашел онлайн-сервис, который выполняет упрощение при сохранении топологий: mapshaper.org . Я использовал его как альтернативу, так как кажется, что это невозможно в R напрямую (QGis на самом деле не сохраняет топологию, и я не знаю других программных пакетов)

— yellowcap

1

Я проверил ваши наблюдения с помощью Saga 2.0.8, и результат не сохранил топологию!

— Фабиан Цейндль

2

Теперь есть замечательный rmapshaperпакет от Andy Teucher, который включает в себя функцию упрощения, которая «выполняет топологически-упрощенное многоугольное упрощение».

Из их репозитория github, пример использования:

states_simp <- ms_simplify(states_sp)

где states_spявляется spatialPolygons*объектом.

См. Пакет README для получения дополнительной информации: https://github.com/ateucher/rmapshaper

— Фил
источник

0

это должно помочь:

gSimplify(..., topologyPreserve=T)из GEOSпакета

— petermeissner
источник

2

К сожалению, опция сохранения топологии в gSimplifyдействительности не сохраняет топологию для нескольких полигонов. Я думаю, что это больше относится к топологии одиночных или мультиполигонов. Для шейп-файла со многими полигонами топология не сохраняется.

— yellowcap