Extraction de la ligne centrale d'un polygone complexe dans PostGIS / Python

Une idée comment extraire la ligne médiane du polygone suivant (Fig1) en utilisant des fonctions PostGIS ou des modules Python? J'essayais la fonction Voronoi mais ce n'est pas utile pour ce type de forme. Je veux quelque chose comme Fig2 blackline.

Fig1: Forme de polygone Fig2: Ligne centrale

python postgis voronoi-thiessen

— Zia
source

Vous voudrez peut-être regarder cette réponse: gis.stackexchange.com/questions/29863/… Je sais qu'il utilise ArcGIS à la place, mais l'approche peut être utile.

— WhiteboxDev

Que diriez-vous de: postgis.net/docs/ST_StraightSkeleton.html ?

— Jakub Kania

@Jakub J'essaie de configurer SFCGAL (nécessaire pour la fonction Skeleton) sur ma base de données PostGIS déjà existante. Bien que j'aie installé le libSFCGAL.sofichier dans le /usr/local/librépertoire, mais je ne sais pas comment importer ST_StriaghtSkeleton()la base de données? Une idée?

— Zia

PostGIS doit être compilé avec le support SFCGAL, mais c'est une question différente et pas quelque chose que je connais trop bien.

— Jakub Kania

@Jakub Merci. J'ai en quelque sorte résolu cette configuration SFCGAL par moi-même ici

— Zia

Réponses:

Cela peut être un mauvais outil pour vous, mais vous pouvez le faire simplement avec OpenJUMP et le plugin Skeletonizer. Le plugin peut être téléchargé à partir de http://sourceforge.net/projects/jump-pilot/files/OpenJUMP_plugins/More%20Plugins/JUMP%20Skeletonizer%20Plugin/

Le fichier zip contient un bon manuel. Le plugin est destiné aux cours d'eau et aux plans d'eau, mais votre polygone imite un flux pour moi.

Ceci est votre polygone tel que je l'ai numérisé. Remarquez les sommets dans les lignes centrales des entrées / sorties.

entrez la description de l'image ici

Points d'entrée numérisés sur une couche séparée. Les points doivent s'accrocher aux sommets sur la limite du polygone.

entrez la description de l'image ici

Utilisez les outils Skeletonizer et remplissez les champs

entrez la description de l'image ici

Obtenez le résultat

entrez la description de l'image ici

— user30184
source

Malheureusement, le fichier jar ne contient aucun manuel. Pourriez-vous s'il vous plaît me guider comment ajouter ce plugin à mon interface graphique OpenJump?

— Zia

Regardez à l'intérieur du zip , dans skeletonizer-1.0 \ doc \ vous devriez trouver le fichier "SkeletonizerUserGuide.doc". L'installation se fait en déposant les fichiers jar de skeletonizer-1.0 \ lib \ dans le répertoire lib \ ext dans OpenJUMP.

— user30184

Cependant, il semble que la plupart des jars ou des versions plus récentes de ceux-ci soient déjà inclus dans OpenJUMP soit dans \ lib soit dans \ lib \ ext. Il est peut-être préférable de copier uniquement ceux qui semblent manquer. Un redémarrage est nécessaire après l'ajout de nouveaux pots.

— user30184

J'ai trouvé ce pot ici . Mais je reçois cette erreur

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

lors de l'exécution de ce plugin Skeletonizer

— Zia

Si votre polygone est raster, numérisez-le.

Une fois que votre polygone est vraiment un polygone :), une première solution se trouve, comme @WhiteboxDev l'a dit, c'est le squelette droit, mais il a le problème que, outre la ligne médiane que vous voulez, il renvoie également beaucoup de segments rayonnant de cette ligne vers la bordure de votre polygone. Mais si vous avez PostGIS 2.2+, vous pouvez utiliser à la ST_ApproximateMedialAxis()place, ce que vous voulez.

— Marcos Dione
source