Extracción de la línea central de un polígono complejo en PostGIS / Python

¿Alguna idea de cómo extraer la línea central del siguiente polígono (Fig1) usando las funciones de PostGIS o los módulos de Python? Estaba probando la función Voronoi, pero eso no es útil para este tipo de forma. Quiero algo como Fig2 blackline.

Fig1: Forma de polígono Fig2: Línea central

python postgis voronoi-thiessen

— Zia
fuente

Es posible que desee ver esta respuesta: gis.stackexchange.com/questions/29863/… Sé que utiliza ArcGIS en su lugar, pero el enfoque puede ser útil.

— WhiteboxDev

¿Qué tal: postgis.net/docs/ST_StraightSkeleton.html ?

— Jakub Kania

@Jakub Estoy tratando de configurar SFCGAL (necesario para la función Esqueleto) para mi base de datos PostGIS ya existente. ¿Aunque he instalado el libSFCGAL.soarchivo en el /usr/local/libdirectorio pero no sé cómo traer ST_StriaghtSkeleton()el db? ¿Alguna idea?

— Zia

PostGIS debe compilarse con soporte SFCGAL, pero esa es una pregunta diferente y no es algo con lo que esté demasiado familiarizado.

— Jakub Kania

@Jakub Gracias. He resuelto de alguna manera esta configuración SFCGAL por mí mismo aquí

— Zia

Respuestas:

Esta puede ser una herramienta incorrecta para usted, pero puede hacerlo simplemente con OpenJUMP y el complemento Skeletonizer. El complemento se puede descargar desde http://sourceforge.net/projects/jump-pilot/files/OpenJUMP_plugins/More%20Plugins/JUMP%20Skeletonizer%20Plugin/

El archivo zip contiene un buen manual. El complemento está destinado a corrientes y cuerpos de agua, pero su polígono me imita una corriente.

Este es su polígono como lo digitalicé. Observe los vértices en las líneas centrales de las entradas / salidas.

ingrese la descripción de la imagen aquí

Puntos de entrada digitalizados en una capa separada. Los puntos deben ajustarse a los vértices en el límite del polígono.

ingrese la descripción de la imagen aquí

Use las herramientas de Skeletonizer y complete los campos.

ingrese la descripción de la imagen aquí

Obtén el resultado

ingrese la descripción de la imagen aquí

— usuario30184
fuente

Lamentablemente, el archivo jar no contiene ningún manual. ¿Podría guiarme cómo agregar este complemento a mi OpenJump-GUI?

— Zia

Mire dentro del zip , en skeletonizer-1.0 \ doc \ debería encontrar el archivo "SkeletonizerUserGuide.doc". La instalación se realiza colocando los archivos jar del skeletonizer-1.0 \ lib \ en el directorio lib \ ext en OpenJUMP.

— user30184

Sin embargo, parece que la mayoría de los frascos o versiones más nuevas de esos ya están incluidos en OpenJUMP, ya sea en \ lib o en \ lib \ ext. Quizás sea mejor copiar solo aquellos que parecen faltar. Es necesario reiniciar después de agregar nuevos frascos.

— user30184

Encontré este frasco aquí . Pero recibo este error

java.lang.Exception: Could not determine shapefile type - data is either all GeometryCollections or empty (Illegal State Exception)

al ejecutar este complemento Skeletonizer

— Zia

Si su polígono es ráster, digitalícelo.

Una vez que su polígono es realmente un polígono :), una primera solución radica, como dijo @WhiteboxDev, es el esqueleto recto, pero tiene el problema de que, además de la línea central que desea, también devuelve una gran cantidad de segmentos que se irradian desde esa línea a El borde de su polígono. Pero si tiene PostGIS 2.2+, puede usarlo ST_ApproximateMedialAxis(), que es lo que desea.

— Marcos Dione
fuente