Le portail francophone de la géomatique

Mathilde G · Thu 25 September 2014 16:29

Bonjour,

Ca fait maintenant plusieurs jours que je cherche et je ne trouve pas la bonne formule.

Je cherche à créer une vue résultant de l'intersection d'une table a (table de parcelles cadastrales) et deux tables b et c (zones de PLU, et contraintes de PLU).
J'arrive à faire cette requête, mais j'ai réalisé que la table des contraintes (c) superpose quelque fois la table des zones plu (b), je souhaiterai donc sortir les intersections entre les tables b et c de la table c (dans une vue de préférence, pour garder l'info complète) avant d'effectuer ma requête finale avec la table a.

Je crois que passer par une union ne me convient pas car j'ai besoin de garder les attributs (zones U, AU... des PLU) et lorsque il n'y a pas d'intersection avec la table des PLU je souhaite récupéré les attributs de la table c (contraintes).

D'après ce que j'ai compris un ST_difference ne marche pas car il y toujours des géométries de la table B qui sont disjointes de la table C, la requête me retourne donc un nombre démultiplié d'objets.

Je ne vois donc pas bien par quelle requête il faut passer, je vous remercie de m'éclairer un peu et j'espère avoir été claire.

Mathilde

ChristopheV · Fri 26 September 2014 07:54

Bonjour,

j'espère avoir été claire

Là pas trop

la table des contraintes (c) superpose quelque fois la table des zones plu (b), je souhaiterai donc sortir les intersections entre les tables b et c de la table c

superposition et intersection .... ??

Mathilde G · Fri 26 September 2014 09:39

Oui exactement

Mathilde G · Fri 26 September 2014 10:16

Si ça peut aider, j'ai d'abord fait deux vues de mes tables plu et contraintes (prescr) pour sélectionner les attributs qui m'intéressent :

Code:

-- CREATE VIEW selection des zones PLU qui nous intéressent :
CREATE VIEW crozes.v_plu AS
SELECT *
FROM doc_urba
WHERE doc_urba.zone_simpl='AU' OR doc_urba.zone_simpl='U' OR doc_urba.zone_simpl='C' OR doc_urba.zone_simpl='constr';

-- CREATE VIEW selection des zones de prescription du PLU qui nous intéressent :
CREATE VIEW crozes.v_prescr AS
SELECT *
FROM prescription_plu 
WHERE prescription_plu.id_prescription= 1 OR prescription_plu.id_prescription= 2 OR prescription_plu.id_prescription= 6
    ;

Ensuite ma dernière requête ressemblerai à quelque chose comme ça :

Code:

-- Intersection des parcelles perdues avec les documents d'urbanismes et les contraintes du PLU :

CREATE OR REPLACE VIEW crozes.v_intersect_doc_urba AS
(WITH table_temp AS

(SELECT clipped.idu, clipped.type_raison, clipped.sum as supf_plantee, 
ROUND(ST_AREA(clipped_geom)::numeric/10000,4) AS supf_geom, clipped_geom AS geom
FROM (SELECT crozes.v_perdu_plante.gid, crozes.v_perdu_plante.idu, crozes.v_perdu_plante.sum, crozes.v_perdu_plante.type_raison,
    (ST_Dump (ST_Intersection(crozes.v_perdu_plante.geom, v_prescr.geom))).geom AS clipped_geom
    FROM (SELECT *
        FROM crozes.v_prescr) AS v_prescr
        INNER JOIN crozes.v_perdu_plante ON ST_Intersects(crozes.v_perdu_plante.geom, v_prescr.geom)) As clipped


UNION

SELECT clipped.idu, clipped.type_raison, clipped.sum as supf_plantee, 
ROUND(ST_AREA(clipped_geom)::numeric/10000,4) AS supf_geom, clipped_geom AS geom
FROM (SELECT crozes.v_perdu_plante.gid, crozes.v_perdu_plante.idu, crozes.v_perdu_plante.sum, crozes.v_perdu_plante.type_raison,
    (ST_Dump(ST_Intersection(crozes.v_perdu_plante.geom, v_plu.geom))).geom AS clipped_geom
    FROM (SELECT *
        FROM crozes.v_plu) As v_plu
        INNER JOIN crozes.v_perdu_plante ON ST_Intersects(crozes.v_perdu_plante.geom, v_plu.geom)) As clipped
)
SELECT row_number() OVER () AS gid, * FROM table_temp );

Mathilde G · Mon 29 September 2014 09:58

Bonjour,

Je vois que ma question ne vous a pas trop inspirée, alors je vais essayer de reformuler :

J'ai deux tables de polygones (multi), je souhaiterai pouvoir effacer d'une des deux tables toutes les intersections avec la deuxième, dans une vue.
(C'est-à-dire effacer toutes les entités de la première table contenues dans un polygone de la deuxième table et découper les polygones de la première qui intersecteraient la deuxième, de manière à ne garder que les parties extérieures).

J'ai lu qq chose sur ST_relate, je ne sais pas si ça pourrait s'appliquer ici, ou peut-être avez-vous d'autres idées, je vais aussi réessayer un ST_union, mais n'hésitez pas à me donner des pistes
merci de votre coup de pouce,

Mathilde

Nicolas Ribot · Mon 29 September 2014 10:37

Bonjour,

Tel que je le comprends, vous voulez faire la différence entre deux couches, avec trois conditions sur ces couches:
st_intersects(a, b) and not st_contains(b, a) and not st_contains(a, b):

Code:

with tmp as (
    select a.id, b.id as bid, st_difference(a.geom, b.geom)
    from pg1 a, pg2 b
    where st_intersects(a.geom, b.geom) 
      and not st_contains(b.geom, a.geom)
      and not st_contains(a.geom, b.geom)
) select t.* 
from tmp t;

Image: en rouge: couche 1, en bleu-vert: couche2, en noir, le contour du résultat

Nicolas

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

Mathilde G · Mon 29 September 2014 12:47

Bonjour,
et merci pour votre réponse, c'était presque ce que je recherchais ! En tout cas, ça m'a bien aidé. Grâce à votre proposition voici ce que j'ai fait :

J'ai d'abord créé une vue ne contenant un seul polygone car dans ma couche d'intérêt je pouvais avoir un polygone qui intersectait plusieurs polygone de la seconde couche, le ST_difference ne renvoyait donc pas le résultat attendu

Code:

CREATE OR REPLACE VIEW pg2 AS 
SELECT row_number() OVER (ORDER BY 2::integer) AS gid, st_union (geom) as geom
FROM v_plu;

J'ai ensuite fait la requête que vous m'avez proposé en rajoutant une union d'un ST_disjoint qui n'était pas pris en compte, et j'ai ainsi tous les polygones ou parties de polygones en dehors de ceux de la deuxième couche !

Code:

CREATE OR REPLACE VIEW v_diff AS
(with tmp as (
    select a.id, st_difference(a.geom, b.geom) AS geom
    from pg1 a INNER JOIN pg2 b ON
    st_intersects(a.geom, b.geom) 
      and not st_contains(b.geom, a.geom)
      and not st_contains(a.geom, b.geom) 

UNION 


    select a.id, a.geom
    from pg1 a INNER JOIN pg2 b ON
    st_disjoint(a.geom, b.geom)  
) 
SELECT row_number() OVER () AS gid, * FROM tmp );

Encore merci pour votre aide, bonne continuation,

Mathilde

Nicolas Ribot · Mon 29 September 2014 13:47

Ah oui ! J'avais un peu zappé la partie des polygones sans interaction avec l'autre couche...

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 25 September 2014 16:29

POSTGIS: ST_Difference ? Quelle fonction pour l'inverse d'Intersection

#2 Fri 26 September 2014 07:54

Re: POSTGIS: ST_Difference ? Quelle fonction pour l'inverse d'Intersection

#3 Fri 26 September 2014 09:39

Re: POSTGIS: ST_Difference ? Quelle fonction pour l'inverse d'Intersection

#4 Fri 26 September 2014 10:16

Re: POSTGIS: ST_Difference ? Quelle fonction pour l'inverse d'Intersection

Code:

Code:

#5 Mon 29 September 2014 09:58

Re: POSTGIS: ST_Difference ? Quelle fonction pour l'inverse d'Intersection

#6 Mon 29 September 2014 10:37

Re: POSTGIS: ST_Difference ? Quelle fonction pour l'inverse d'Intersection

Code:

#7 Mon 29 September 2014 12:47

Re: POSTGIS: ST_Difference ? Quelle fonction pour l'inverse d'Intersection

Code:

Code:

#8 Mon 29 September 2014 13:47

Re: POSTGIS: ST_Difference ? Quelle fonction pour l'inverse d'Intersection

Pied de page des forums