Message posté par : SANTANNA
----------------------------------------
Hello,
Tu veux bien essayer ce code? Pour la partie intersection, je te laisse réajuster par
rapport à ton scénario, je l'ai fait vis-à-vis de l'atlas direct. Le plus
intéressant est la deuxième partie, je pense.
-----------------
Code :
with_variable(
'geomInt' -- on crée une variable pour récupérer la géométrie en superposition
, collect_geometries(
intersection( $geometry , geometry( @atlas_feature ) )
)
, with_variable( 'nbrParts' -- une autre variable qui stocke le nombre de parts
dans la geom en superposition
, generate_series( 1, num_geometries(@geomInt))
, with_variable( 'geomCibleListe' -- on convertit la geometrie
multipolygone en array, en itérant sur le nbre de parts
, array_foreach(
@nbrParts
, geometry_n( @geomInt, @element ) )
, collect_geometries( -- on regénère une géométrie
multipolygone dont on filtrera les parties indésirées
array_filter( @geomCibleListe
, area(@element) >
100000 -- on indique la surface seuil
)
)
)
)
)
-----------------
----------------------------------------
Le message est situé
https://georezo.net/forum/viewtopic.php?pid=343095#p343095
Pour y répondre : qgis_fr(a)ml.georezo.net ou reply de votre messagerie
Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement
--
Association GeoRezo - le portail géomatique
https://georezo.net