Le portail francophone de la géomatique

Bonjour Paul, merci de votre réponse !

je viens d'intégrer la formule, mais la touche OK reste grisée.. ai-je mal rentré la formule ?

La formule a en effet fonctionné, mais cela n'est pas ce que je souhaite, j'ai du mal m'exprimer.
En voyant la photo (de gauche) ci-jointe, on remarque 3 parcelles cadastrales (D423, D424, D425).
En voyant la photo (de droite) ci-jointe, on voit 4 parcelles forestières, dont la parcelle 10c, par exemple, qui surplombe les 3 parcelles cadastrales.
Je souhaiterais pouvoir générer un tableau m'indiquant précisément, quelle est la surface de la 10d dans chacune de ces 3 parcelles (et ceci de manière générale sur toute la propriété).

Est-ce possible selon vous ?

Je vous remercie encore pour votre aide précieuse !

Bien à vous,

Corentin

Alors,
c'est possible, et il faut réaliser 3 actions distinctes
Étape 1 : (si ce n'est pas déjà fait) calculer l'aire des polygones de la couche 'peuplement' --> calculatrice de champs, nouveau champ, formule : $area

Étape 2 : Lancer l'algo qgis "Joindre les attributs par localisation" ( et non "Joindre les attributs par localisation (résumé)")
Couche de base --> parcellesst-martin, joindre la couche --> peuplement, champs à ajouter --> le champ avec les noms des parcelles forestière et les aires calculés en étape 1
l'algo va vous sortir une nouvelle couche --> ça sera la meme que "parcellesst-martin" mais avec des parcelles dupliqués (autant de doublon que d'intersections entre les deux couches)

Étape 3 : sur la nouvelle couche, créer un champs "surface_intersect" (c'est le champs qui vous donnera la surface intersectée)
remplir avec la formule suivante :

area(intersection($geometry,geometry(get_feature('peuplement','id_peuplement', "id_peuplement" ))))

remplacer juste "id_peuplement" par le nom du champ contenant les numéros des parcelles forestières

la table attributaire sera alors un tableau qui vous permettra de faire la correspondance entre chaque parcelles cadastrale et les surfaces intersectées sur les parcelles forestières

Dernière modification par shirosaki (Wed 22 December 2021 14:50)

Ah oui je connaissais pas cet algo, ça permet de simplifier le process assez radicalement !

Oui vous n'avez pas remplacé 'id_peuplement' par son vrai nom, 'id_peuplement' c'était juste un exemple

sinon est ce que vous pouvez nous transmettre un extrait de vos données ? comme ça je vous donne les formules avec les bon noms de champs (les tables peuvent être vides, ou inventées, il me faut juste les noms de champs)

Dernière modification par shirosaki (Wed 22 December 2021 15:52)

je vous remercie
petite question résiduelle : Dans la couche peuplement, est ce que le champ "Parcelle" est unique (est ce que les valeurs peuvent apparaitre plus d'une fois ? deux entités 10a par exemple ? si non on va devoir créer un champ unique préalablement)

On reprend :

Étape 1 : déjà fait d'après votre capture d'écran

Étape 2 : Lancer l'algo qgis "Joindre les attributs par localisation"
Couche de base --> parcellesst-martin, joindre la couche --> peuplement, champs à ajouter --> "Parcelle"
l'algo va vous sortir une nouvelle couche --> ça sera la meme que "parcellesst-martin" mais avec des parcelles dupliqués (autant de doublon que d'intersections entre les deux couches)

Étape 3 : sur la nouvelle couche, créer un champs "surface_intersect" (c'est le champs qui vous donnera la surface intersectée)
remplir avec la formule suivante :

area(intersection($geometry,geometry(get_feature('peuplements','Parcelle', "Parcelle" ))))

Dernière modification par shirosaki (Wed 22 December 2021 16:28)

Ah cool !
Je suis content d'avoir pu vous aider,

Bonnes fêtes également