Le portail francophone de la géomatique

nawak · Tue 03 November 2009 13:30

salut,

Petite question sur une requête sql sous postgres/postgis :

j'ai une entrée dans une table qui possède la géométrie de la france et de la corse : je souhaite n'avoir plus que la france je dessine donc un polygone contourant la france et je fais :

Code:

CREATE TABLE AS
      SELECT intersection(b.the_geom,a.the_geom) as the_geom 
      FROMfrance a,limite b 
      WHERE a.the_geom && b.the_geom;

mais comment puis-je inverser cette selection et n'avoir que les coordonnées de la Corse sans bien sur recréer un polygone contourant la corse...

Merci

Lionel B · Tue 03 November 2009 14:37

where not (a.the_geom && b.the_geom)

cordialement

Yves · Tue 03 November 2009 14:39

Bonjour,

Quelque chose comme cela :

Code:

CREATE TABLE AS 
    SELECT GeomUnion(a.the_geom) AS the_geom
    FROM france a,limite b
    WHERE !ST_Contains(b.the_geom,a.the_geom);

Par contre je ne sais plus si le symbole "!" fonctionne pour les fonctions SQL-MM dans postgis. Il faudra peut être mettre

Code:

ST_Contains(b.the_geom,a.the_geom)=false

Y.

Lionel B · Tue 03 November 2009 14:42

oups, erreur d'interprétaton de ma part.

tu peux faire une différence entre ton premier résultat et les données de départ.

Code:

CREATE TABLE AS
      SELECT difference(c.the_geom,d.the_geom) as the_geom 
      FROM france c,  (SELECT intersection(b.the_geom,a.the_geom) as the_geom 
      FROM france a,limite b WHERE a.the_geom && b.the_geom) as d

Dernière modification par Lionel B (Tue 03 November 2009 14:45)

nawak · Tue 03 November 2009 15:26

yes merci Lionel ta soluce marche..pour yves merci aussi mais ça ne fonctionne pas

Thomas KOCHERT · Tue 03 November 2009 15:45

Bonjour, pour aller au plus simple, j'aurais fait

SELECT * FROM france INNER JOIN limite ON _st_contains(limite.the_geom,france.the_geom)

et ... ON NOT _st_contains....

pour avoir respectivement la France ou la Corse.

Attention à bien mettre le _ devant l'opérateur st_contains sinon l'analyse se fait sur les rectangles englobants et non sur les géométries.

vincentp · Tue 03 November 2009 17:10

Bonjour,
Si tu n'as que la france et la corse dans ton multipolygone, tu peux simplement utiliser

Code:

st_geometryn(the_geom, 2)

pour récupérer la corse si tant est qu'elle est en deuxième.

Si tu as d'autres iles ou polygones dans ton multipoly, je tenterais plutot qu'une intersection ou une différence complexe de faire un dump du multipolygone et de faire le contains des polygones récupérés avec ton rectangle :

Code:

select * from (
  select
    st_dump(the_geom) as the_geom2
  from 
    france
) as t1
, limite as t2
where 
  not st_contains(t2.the_geom, t1.the_geom2)

@thomas : st_contains fait bien le test sur les géométries, mais dans un deuxième temps. Dans un premier temps il utilise les BBox (et leurs indexes), et pour les positifs une vérification est faite sur les géométries.

Christophe Brun · Fri 06 November 2009 14:38

Bonjour,

Juste une petite précision à propos de :

Thomas KOCHERT a écrit:

Bonjour, pour aller au plus simple, j'aurais fait
Attention à bien mettre le _ devant l'opérateur st_contains sinon l'analyse se fait sur les rectangles englobants et non sur les géométries.

Dans les 2 cas l'analyse se fait bien sur les géométries et non sur les bbox comme le fait l'opérateur &&. La subtilité est que avec ST_Contains l'index spatial est utilisé, donc les bbox sont utilisées pour la recherche; avec _St_Contains, l'analyse se passe de l'index.

Bon week-end,

Christophe

PS: oups j'avais pas vu la fin du message de vincentp ! Désolé pour la redite.

Dernière modification par Christophe Brun (Fri 06 November 2009 14:41)

Thomas KOCHERT · Fri 06 November 2009 16:23

Merci pour cette précision !

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 Tue 03 November 2009 13:30

Postgres/postgis requête sql

Code:

#2 Tue 03 November 2009 14:37

Re: Postgres/postgis requête sql

#3 Tue 03 November 2009 14:39

Re: Postgres/postgis requête sql

Code:

Code:

#4 Tue 03 November 2009 14:42

Re: Postgres/postgis requête sql

Code:

#5 Tue 03 November 2009 15:26

Re: Postgres/postgis requête sql

#6 Tue 03 November 2009 15:45

Re: Postgres/postgis requête sql

#7 Tue 03 November 2009 17:10

Re: Postgres/postgis requête sql

Code:

Code:

#8 Fri 06 November 2009 14:38

Re: Postgres/postgis requête sql

#9 Fri 06 November 2009 16:23

Re: Postgres/postgis requête sql

Pied de page des forums