Le portail francophone de la géomatique

Bonsoir,

\1ere méthode : un aggregat par la calculatrice permet cela et plus encore :

aggregate('DEPARTEMENT','max',"INSEE_DEP",intersects(centroid(geometry(@parent)),$geometry))

ici (totalement inutile mais à trnasposer), je suis sur une couche commune en màj sur un champ pour récupérer le code dep de ma couche DEPARTEMENT type admin express.
Il faut voir la couche commune, en mise à jour, comme ta couche de points (ici=centroid(geometry(@parent)), donc pour ton cas seulement geometry(@parent) car couche de points : ok pour intersects avec $geometry de 'DEPARTEMENT' (=ta couche de polygone "contenants"))
le $geometry dans le aggregate concerne la couche cible 'DEPARTEMENT' ; Je te renvoie sur la doc dans l'interface graphique pour les autres méthodes d'aggrégat ('concatenate','collect',...), mais ça devrait faire le taf, à adapter selon ton champs à récupérer et ta couche cible!
aggregate('couche_cible','methode_aggr',"champ_cible",expression filtre,..)

\2eme méthode : extension Refunction à installer, puis dans la calculatrice de champs dans onglet/groupe Reference, de nouvelles fonctionnalités  permettent de croiser deux couches et de récupérer les infos de l'une dans l'autre

Ou jointure spatiale en effet avec génération d'une nouvelle couche

Bonjour à tous,
Je relance le sujet car j'ai un cas à peu près identique mais qui ne fonctionne qu'à moitié !
J'ai donc, sous format shape, une couche de points et une couche de communes.
Je souhaite renseigner dans la couche points, l'attribut "nom" avec le nom de la commune qui les contiennent.
J'ai donc utilisé la formule suivante :

 aggregate( 'communes','concatenate', "nom", intersects($geometry, geometry(@parent)))

où "communes" est le nom de ma couche de communes et "nom", l'attribut que je cherche à renseigner.
J'ai d'abord fait le test sur un extrait de données et cela a très bien fonctionné.
Or quand je veux le faire sur ma couche de points entière, je n'ai rien (même l'aperçu du résultat dans la calculatrice de champs ne renvoie que ", alors qu'avec l'extrait de la couche, je voyais bien apparaître un nom de commune).
J'en déduis qu'il me manque quelque chose, mais je ne vois pas quoi...
Merci de vos éclaircissements.
Cordialement.

L'utilisation d'aggregate et concatenate me semble inversée. Ce serait pertinent si vous vouliez créer une phrase pour chaque commune avec le nom des points qu'elle contient. Pas pour l'inverse. Là vous voulez une valeur unique. Pas besoin d'aggréger ni de concatener. Facile à faire en SQL, mais la syntaxe de la calculatrice de champs ne m'est pas familière. Je vais tenter de m'y pencher.

Dernière modification par carteq (Thu 16 April 2020 13:48)

Concatenate est essentiel cas d'aggrégation, en effet. Mais là, aucun besoin d'aggréger.

J'ai essayé :

aggregate( layer:='communes', aggregate:='concatenate', expression:="champ_nom_commune_dans_la_couche_de_communes",
               filter:=intersects( $geometry, geometry( @parent ) ) )

Ça fonctionne. J'ai testé 19000 entités. Ça prend un temps ridiculement long comparé au SQL.

Je me demande si votre expression "nom" n'est pas erroné dans votre essai mentionné plus haut ? Est-ce vraiment le nom du champ dans la couche commune ?

Je vais chercher quelque chose de plus direct que aggregate.

Dans la video, il y un lien de un-à-plusieurs. Dans votre cas, c'est un-pour-un : un point n'est que dans une commune. Rien à aggréger. Il faut juste aller chercher la valeur.

@parent :

Il semble qu'on ait pas le choix d'utiliser aggregate si on utilise @parent.

Dernière modification par carteq (Thu 16 April 2020 19:35)

bonjour,

J'aimerais croiser les geometries d'un shp polygone avec celle d'un shp ligne.

Il faudrait que j'attribue dans un nouvelle colonne de mon shp polygone, les linéaires(via : $length) de ligne contenue dans chaque geometrie de polygone.

Cela est t il possible par la calculatrice d'expression?

Bonjour,

Pour les dernières versions de QGIS 3, la fonction est intégrée, exemple :

array_to_string(overlay_within('communes',insee ))
Pour récupérer l'attribut insee de la couche 'commune'.

Pour les versions précédentes, il faut installer l'extension refunction et utiliser ainsi la fonction geomwithin('communes',insee ).

S'il n'y a pas de résultats, plusieurs possibilités :

- Les couches n'ont pas le même système de projection

- L'attribut est mauvais ex : colonne formatée pour recevoir des nombres entiers alors qu'on essaye de mettre du texte.

- Plusieurs polygones se superposent, la fonction ne sait pas lequel prendre

- La formule est mauvaise, il faut remplacer within par touches ou nearest ou contains ou interstects selon l'opération spatiale voulue. (Mais dans votre cas ça devrait être within.

Bonne journée,

Dernière modification par Robin.k (Thu 28 April 2022 15:46)