Le portail francophone de la géomatique

Bonjour,

Sylvain, c'est le moteur de recherche qui renvoie de la doc 3.4 (lequel et quels termes utilisés) ou le lien était dans les favoris/l'historique?

Peut-être avec la moitié de ce champ calculé plutôt, si Hippo veut que les zones tampons se touchent juste. Par contre, cela veut dire qu'hormis un contexte de grille régulière, les entités seraient réunies deux à deux. Peut-être pas ce que veut Hippo.
Hippo, quelle est la finalité de l'opération au juste?

Bonjour,

Je me suis penché sur ce sujet que je trouve super intéressant, et j'ai trouvé un moyen (avec Qgis natif, un plugin et aucun script ni SQL) de découper un polygone par une couche de points, les lignes de découpe étant équidistantes aux points proches. Voici ma démarche :

disposer au préalable d'une couche de points (appelée par la suite 'couche_points_ini') et d'une couche de polygone(s) ('couche_polygone_ini')

1- "Matrice des distances" de 'couche_points_ini' vers elle-même (la couche doit avoir un ID unique) => le résultat est une couche multipoints composées des noeuds de départ et d'arrivée pour chaque distance
2- "De morceaux multiples à morceaux uniques" avec le résultat précédent, puis donner un ID unique à chaque point
3- "Points vers lignes" (champ de tri sans importance, champ de regroupement : id_unique précédent)
4- "Supprimer les géométries dupliquées" car le traitement précédent va relier chaque points aux autres dans un sens et dans l'autre
5- "Centroïdes" des lignes précédentes
6- installer le plugin "Water Network Tools" et utiliser le traitement "Split lines at points" avec les centroïdes pour les points et les lignes nettoyées de ses doublons pour les lignes
7- sélectionner et supprimer un demi-chemin sur deux (par exemple ceux qui partent d'un noeud initial) :
    aggregate(
    'id_de_la_couche_points_ini',
    'count',
    "champ_quelconque",
    intersects(buffer($geometry,0.01),start_point(geometry(@parent)))
    )>0

7- "Transect" (500m, les deux côtés) : ce traitement crée des perpendiculaires aux extrémités des lignes
8- sélectionner et supprimer les transect qui intersectent pile (avec un buffer de 0.01 pour que l'intersection fonctionne) les points de couche_points_ini :
    aggregate(
    'id_de_la_couche_points_ini',
    'count',
    "champ_quelconque",
    intersects($geometry,buffer(geometry(@parent),0.01))
    )
    >0

7- "Couper avec des lignes" : découper 'couche_polygone_ini' avec les lignes transect restantes
8- sélectionner et supprimer les polygones découpés ayant une aire < 0.1 ou nulle :
     $area is null or $area<0.1
9- "Accrocher les géométries à la couche" : de la couche résultat précédente vers elle-même, tolérance : 0.1m
10- créer un ID unique pour chaque polygone
11- "Distance au plus proche centre (points)" avec en couche source la couche de polygones obtenue précédemment et en destination 'couche_points_ini' (attribut nom de la couche des centres l'id unique créé précédemment)
12- "Joindre les attributs par valeur de champ" : en entrée la couche de polygones et en entrée 2 la couche de centres créée précédemment. Indiquer l'id unique pour les deux et choisir uniquement l'atttribut "HubName" à ajouter à la couche source
13- "Regrouper" en utilisant le champ "HubName" comme champ de regroupement
14- "Géométrie d'emprise minimale" :  choisir "HubName" comme champ de regroupement et enveloppe convexe comme type de géométrie. ça permet de supprimer les artéfacts liés au regroupement de plein de petits polygones

Résultat : voir le fichier png en PJ

Dernière modification par aheurion (Thu 04 November 2021 08:47)