Discussion :

[mcdelay]

Bonjour,

Je cherche à faire une requête spatiale entre 3 couches mais je m'emmêle les pinceaux.
J'ai 3 tables une couche de point (point de collecte), et 2 couches de polygone (zone et batiment). Je cherche a savoir pour un point de collecte, le nombre de bâtiment se trouvant dans la zone la plus proche d'un point de collecte. En gros j'ai 2 requêtes et il faudrait que je les imbrique !

1er : pour connaître le nombre de bâtiment dans une zone

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
 
SELECT COUNT(DISTINCT(bati.ogc_fid)) AS nbreLogement FROM zone left join bati Intersects(bati.the_geom , zone.the_geom)

2ième: pour connaître la zone la plus proche d'un point de collecte

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
 
SELECT zone.ogc_fid, Distance(zone.the_geom ,Pointcollecte.the_geom)   FROM zone, Pointcollecte ORDER BY Distance ASC Limit 1

(Peut être faut-il écrire cette dernière autrement ...)

Si vous avez une idée , merci !

[SQLpro]

Ce que vous faites c'est vraiment du n'importe quoi. Commencez par donner vos tables sous la forme CREATE TABLE ..., ainsi qu'un jeu d'essais sous la forme INSERT ...
Après on pourra commencer à vous aider.

Mais il serait peut être temps de vous forme r à ce qu'est une jointure spatiale !
Mon livre sur SQL possède un chapitre de 60 pages consacré au SIG dans SQL avec des exemples de jointures spatiales...

A +

[mcdelay]

Bonjour,

Je trouve votre réponse un peu arrogante. Je pense avoir été claire, d'ailleurs on m'a donné une réponse sur un autre forum sans que je donne plus d'explications. voici le code :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
with plus_proche as (
   select distinct on(pt.ogc_fid) pt.ogc_fid as pt_id, 
          z.ogc_fid as zone_id,
          z.the_geom as zone_geom
   from zone z, Pointcollecte pt
   where st_dwithin(z.the_geom, pt.the_geom, 500)
   order by pt.ogc_fid, distance(z.the_geom, pt.the_geom)
) select p.pt_id, p.zone_id, count(b.ogc_fid) as nbrelogement
from bati b, plus_proche p
where st_within(b.the_geom, p.zone_geom)
group by p.pt_id, p.zone_id;

... donc un peu d'indulgence pour ceux qui débutent, mais merci d'avoir répondu ;-)