Jointure spatiale en conservant les doublons

[Georgie]

Hello

Je cherche, sous QGIS, une méthode/plugin me permettant de conserver les doublons lors d'une jointure spatiale.

En effet j'ai remarqué que mes doublons sont systématiquement supprimés, ce qui n'est pas souhaitable dans mon cas.

Y-a-t-il une option que je n'aurais pas vue ?

Merci

[SkåL]

Via SQL (spatialite et Virtual Layers) :

SELECT ma_couche1.*,
ma_couche2.champ1,
ma_couche2.champ2,
ma_couche2.champ3
FROM ma_couche1,
ma_couche2
WHERE ST_INTERSECTS(ma_couche1.geometry, ma_couche2.geometry)

La requête sera encore plus efficace avec des couches stockées dans une base de données spatiale qui possède des index spatiaux.
Pour spatialite, il faudra ajouter quelques petites lignes de plus pour faire fonctionner ces index.

Dans la requête ci-dessus, on récupère la géométrie de la couche1.
On pourrait très bien mettre 'ma_couche2.*' mais il faudrait bien spécifier quelle champ géométrique prendre, ou effectuer une opération de jointure comme ST_UNION unissant les géométries de ma_couche1 et de ma_couche2.
J'ai utilisé une jointure d'intersection, mais il est possible d'utiliser 'st_contains'.

[Georgie]

Hello OK merci effectivement ça marche très bien dans Postgis !

Pour être précis, je cherche à faire une jointure spatiale selon une distance d'accroche (une sorte de zone tampon virtuelle, je fais actuellement avec une zone tampon créée juste avant, mais si je peux m'en passer... C'est mieux !).

Et dans QGIS directement, sur des shapes, pas de fonction simplifiée qui me conserverait les doublons ?

[SkåL]

Dans QGIS, on peut avec des shapefiles, passer directement dans Virtual Layers :
- Menu Couche > Ajouter une couche > Ajouter / Modifier une couche virtuelle ...
(conseil perso) - Menu Base de données > Gestionnaire de base de données > Gestionnaire de base de données > Virtual Layers > QGIS Layers > 2ème Bouton 'Fenêtre SQL'

Et pour la distance d'accroche, rien n'empêche un ST_BUFFER :

SELECT ma_couche1.*,
c2.champ1,
c2.champ2,
c2.champ3
FROM ma_couche1,
(SELECT ma_couche2.champ1,
ma_couche2.champ2,
ma_couche2.champ3,
ST_BUFFER(ma_couche2.geometry, 100) AS geometry)) AS c2
WHERE ST_INTERSECTS(ma_couche1.geometry, c2.geometry)

Sinon, peut-être un plugin qui existe mais pour dupliquer les géométries, je n'utilise que ce moyen.

[Georgie]

Ouch !
Extrêmement puissant ces Virtual Layers ! Je ne connaissais pas du tout, c'est magnifique !
#QGIS

Grand merci !