Le portail francophone de la géomatique

simo lay · Thu 28 February 2013 17:27

Bonjour,
j'ai une table des parcelles sous postgis, mon objectif est de fusionner uniquement les parcelles qui se touchent j'utilise la fonction st_union lorsque les parcelles se touchent (st_touches) le problème c'est que quand j'ai uniquement 2 parcelles qui se touchent je les fusionne et j’obtiens un seul objet sans problème (c'est bon) mais dans le cas où j'ai 3 parcelle (ou plus) qui se touchent la fusion se fait par deux objets (parcelles) et du coup j’obtiens 2 objets au lieu d'avoir un seul.
avec une boucle je pense que c'est possible mais j'ai besoin de quelques pistes

sur l'image jointe le problème est plus claire

merci

Dernière modification par simo lay (Fri 01 March 2013 09:14)

Fichier(s) joint(s) :
Pour accéder aux fichiers vous devez vous inscrire.

ChristopheV · Fri 01 March 2013 20:16

bonjour,

a priori je vois pas trop, pour ma part j'ai solutionné ce problème en passant par la topologie sous postgis 2.0

avec un :

topology.TopoGeo_AddPolygon

pour transformer vos polygone en face, arc, nœud et récupérer la géométrie comme étant la composée des trois faces c'est trivial.

tumasgiu · Sat 02 March 2013 00:51

Bonsoir,
peut être que ça répondra à ta problématique :

Tu peux utiliser la fonction st_Union en prenant un tableau de géométries comme argument.
Cette fonction va te renvoyer un polygone avec un ring pour chaque groupe de parcelle contigües.
Ensuite grâce a St_Dump, tu vas extraire un à un chaque ring de ton polygone.

Code:

SELECT geom into "testUnion" FROM 
    (SELECT (St_dump
                 (St_Union
                      (ARRAY( SELECT geom FROM data)
                      )
                  )
               ).geom AS geom
    ) AS t ;

Dernière modification par tumasgiu (Sat 02 March 2013 00:52)

simo lay · Mon 04 March 2013 09:07

Bonjour,

merci à vous ChristopheV et tumasgiu j'ai testé la solution de tumasgiu et ça répond exactement à ce que je voulais faire.

merci beaucoup

Jkerebel04 · Sat 10 October 2020 11:58

Bonjour
Je reviens vers vous pour une question complémentaire concernant ce sujet. Je souhaiterai ajouter une ou plusieurs conditions indexées sur des colonnes de la table "data". Est-il possible d'ajouter une commande GROUP BY. l'essai ci-dessous le fonctionne pas.

Pourriez-vous m'indiquer des pistes SVP ?

CREATE TABLE test_union AS SELECT geom
FROM
(SELECT (St_dump
   (St_Union
( ARRAY( SELECT geom FROM data )
   )
)
).geom AS geom
   GROUP BY data.num
) AS t ;

JP LLORENS · Mon 12 October 2020 09:35

Bonjour
Votre requêt peut être simplifiée

Code:

SELECT (St_dump(St_Union(geom)).geom AS geom
from  FROM data
GROUP BY data.num

Cordialement
JP

Jkerebel04 · Mon 12 October 2020 19:21

Effectivement ! Cela fonctionne bien mieux. merci.
Est_il possible dans la même requête d'ajouter une condition st_intersect ?

Pour l'instant je réalise la requête en deux temps :

d'abord je met a jour une colonne de la table bâtiment'tab_bat' qui recupère id de la table parcelle 'tab_par' accueillant le batiment
UPDATE tab_bat SET id_par_asso = tab_par.id
FROM tab_par
-- lorsque le centroide sur la surface intersecte la geométrie de tab_bat.
WHERE st_intersects(ST_PointOnSurface(tab_bat.geom),tab_par.geom);

Ensuite je crée la table en unissant tout les batiment qui se touche issus de la même parcelle.
DROP TABLE IF EXISTS tab_union CASCADE;
-- creer la table tab_union avec la géomètrie issue de toutes les partie jointe
CREATE TABLE tab_union AS SELECT
-- Créé la géomètrie unique à partie des différents polygone
(st_dump(st_union(geom))).geom as singlegeom,
id_par_asso
FROM tab_bat
GROUP BY id_par_asso;

A+

JD · Tue 13 October 2020 21:34

Bonsoir,
quelque chose dans le genre :

Code:

CREATE TABLE tab_union AS SELECT
-- Créé la géomètrie unique à partie des différents polygone
    (st_dump(st_union(a.geom))).geom as singlegeom,
    b.id
FROM tab_bat a
LEFT JOIN tab_par b ON  st_intersects(ST_PointOnSurface(a.geom),b.geom);
GROUP BY b.id

Jkerebel04 · Sat 17 October 2020 08:07

Exactement. Et redoutable d'efficacité.
Merci.

Jérome

Communauté	Ressources	L'emploi	Le marché	L'association	Fils RSS
Forum Blogs Wiki Planet Sites amis	Agenda Biblio Fiches AMO-CNIG Paris PCGIS Téléchargements	Offres d'emploi Banque de CV Formations Démarche Métiers	Geo-Entreprises Geo-Communiqué Appels d'offres	L'association Qui sommes-nous ? Nous aider Bienfaiteurs Partenaires Historique Statuts	Tout le site Job Géomatique Services GeoBlogs ...	Faire un don Mentions légales

Le portail francophone de la géomatique

Annonce

#1 Thu 28 February 2013 17:27

[PostGIS] : Fusionner les parcelles qui se touchent

#2 Fri 01 March 2013 20:16

Re: [PostGIS] : Fusionner les parcelles qui se touchent

#3 Sat 02 March 2013 00:51

Re: [PostGIS] : Fusionner les parcelles qui se touchent

Code:

#4 Mon 04 March 2013 09:07

Re: [PostGIS] : Fusionner les parcelles qui se touchent

#5 Sat 10 October 2020 11:58

Re: [PostGIS] : Fusionner les parcelles qui se touchent

#6 Mon 12 October 2020 09:35

Re: [PostGIS] : Fusionner les parcelles qui se touchent

Code:

#7 Mon 12 October 2020 19:21

Re: [PostGIS] : Fusionner les parcelles qui se touchent

#8 Tue 13 October 2020 21:34

Re: [PostGIS] : Fusionner les parcelles qui se touchent

Code:

#9 Sat 17 October 2020 08:07

Re: [PostGIS] : Fusionner les parcelles qui se touchent

Pied de page des forums