Le portail francophone de la géomatique

Bonjour,

Sous QGIS, qqch à essayer :

    Convertir le polygone "pol" en polyligne : Vecteur > Géométrie > Polygone vers ligne => "polLigne"
    Générer des points le long de la polyligne avec le plugin QChainage => "polLignePts"
    Générer un diagramme de Voronoï à partir de ces points : Vecteur > Géométrie > Polygones de Voronoï => "polLignePtsVor"
    Transformer ce diagramme en ligne : Vecteur > Géométrie > Polygone vers ligne => "polLignePtsVorligne"
    Découper ce diagramme de Voronoï par "pol" : vecteur > géotraitement > découper => "polLignePtsVorligneDcp"
    Exploser les lignes de ce diagramme : boîte à outils traitement > Explode Lines => "polLignePtsVorligneDcpExpl"
    Générer un tampon autour de "polLigne" => "polLigneTampon"
    Sélectionner les lignes de "polLignePtsVorligneDcpExpl" intersectent "polLigneTampon" avec vecteur > requête spatiale > intersecte.
    Manuellement, désélectionner les lignes de départ et d'arrivée.
    Supprimer les lignes sélectionnées
    Fusionner tout l'objet puis scinder les géométries : vecteur > géométrie > scinder => "polLignePtsVorligneDcpExplScind"
    Cibler les autres artefacts, éventuellement sur des critères de longueur (variable $length) et les enlever.

L'outil r.thin de Grass, accessible dans la boîte à outils Traitements de QGIS, fournit également des résultats intéressants. Il faut l'utiliser sur un raster, donc avoir préalablement rasterisé votre couche de polygones (menu Raster, Conversion, Rastériser).