Le portail francophone de la géomatique

OlivierRustener · Fri 21 February 2014 17:44

Bonjour,

Novice en requête spatiale et SQL, je cherche à récuprer via Postgis toutes les lignes d'une couche en fonction d'une autre couche de délimitation (ex: extraire toutes les lignes de la routes qui se trouve dans la table Villes).

Sous Postgis avec la fonction St_intersects je parviens presque au résultat que j'espère, cependant le découpage des lignes n'est pas "parfait". En effet, une ligne qui se trouve dans la couche de délimitation et qui se poursuit bien au delà n'est pas découpée.

Voici ma requête : select routes.* from routes, villes where st_intersects(routes.geom, villes.geom)

Pouvez-vous, s'il vous plaît m'aider à construire la bonne requête ?

Nicolas Ribot · Fri 21 February 2014 18:12

Bonsoir,

La fonction st_intersects permet de savoir si deux objets s'intersectent. Ca ne les découpe pas.
Il faut utiliser st_intersection(geom, geom) pour découper des objets:

Code:

select st_intersection(routes.geom, villes.geom) 
from routes, villes 
where st_intersects(routes.geom, villes.geom)

Nicolas

OlivierRustener · Fri 21 February 2014 19:05

Merci beaucoup pour cette réponse.

Après vérification dans Qgis, cela découpe bien les routes exactement en fonction de la limite d'une ville, cependant il n'y a plus aucun attribut.
Que faut-il rajouter dans la requête pour garder les informations de la table routes ?

Merci !!!

Nicolas Ribot · Fri 21 February 2014 20:06

Il faut lister les champs que vous voulez dans la partie select:

select champ1, champ2, st_intersection(...
from ...
where ...

Nicolas

OlivierRustener · Sun 23 February 2014 00:35

Merci beaucoup cela fonctionne bien !!!

Cependant, autre problème du coup... Sur de grosse table (plus de 200mo), la requête met une éternité (plus de 30 minutes facilement...) pour sortir un résultat. Est-ce normal ?

PS 1 : Je réalise mes requêtes sur les données OSM de la région Rhône-Alpes
PS 2 : Je deviens fou :-)

Nicolas Ribot · Sun 23 February 2014 12:24

Bonjour,

La performance dépend de beaucoup de parametres, notamment la présence d'un index spatial sur la table de travail, pour ce genre de requetes.
Le jeu de données à découper est important: ca va forcément prendre du temps.

Je vous invite à lire la doc PostGIS, notamment les premiers chapitres de présentation, pour vous familiariser avec l'outil et ses usages.

Nicolas

bernard lortic · Sun 23 February 2014 18:29

Salut,
ouai, je vais encore me faire censuré; mais c'est vrai, vous avez raison : pourquoi faire simple quand on peut faire compliqué ?

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 Fri 21 February 2014 17:44

Requête d'intersection simple (st_intersetcs) ?

#2 Fri 21 February 2014 18:12

Re: Requête d'intersection simple (st_intersetcs) ?

Code:

#3 Fri 21 February 2014 19:05

Re: Requête d'intersection simple (st_intersetcs) ?

#4 Fri 21 February 2014 20:06

Re: Requête d'intersection simple (st_intersetcs) ?

#5 Sun 23 February 2014 00:35

Re: Requête d'intersection simple (st_intersetcs) ?

#6 Sun 23 February 2014 12:24

Re: Requête d'intersection simple (st_intersetcs) ?

#7 Sun 23 February 2014 18:29

Re: Requête d'intersection simple (st_intersetcs) ?

Pied de page des forums