Le portail francophone de la géomatique

Bonsoir,

Pour faire ceci en une seule requete, il faut utiliser quelques "tricks" de postgresql, PostGIS < 2.0 ne possedant pas encore d'operateur Nearest Neighbors:

Le principe est de choisir en premier les zones les plus proches des points de collecte, puis de regarder les batiments inclus (ou intersectant, suivant votre cas) dans les zones (non testé dans votre cas):

with plus_proche as (
    select distinct on(z.ogc_fid) ogc_fid, 
    distance(z.the_geom ,pt.the_geom) as d, 
    z.geom as zone_geom
    from zone z, Pointcollecte pt
    where st_dwithin(z.geom, pt.geom, 500)
    order by d
) select pt_ogc_fid, count(b.ogc_fid) as nbrelogement
from bati b
where st_within(b.geom, zone_geom)
group by pt_ogc_fid;

La premiere requete, obtenue avec la clause WITH, tente de lister les zones les plus proches de chaque point et ne garde que le premier candidat => le plus proche.
Pour cela, l'operateur st_dwithin est utilisé avec une distance fixée en dur, c'est la que le bât blesse: il faut connaitre ou estimer la distance max au dela de laquelle aucun point n'est dans une zone.

Ensuite, la deuxieme requete utilise la premiere pour compter le nombre de batiments dans chaque zone la plus proche du point de collecte.

Concernant la difference entre st_ et pas st_, deux choses:
    • L'adoption par Postgis de la norme SQL MM impose de nommer les fonctions st_, pour plus de clarté et les séparer des autres fonctions de la base.
    • Historiquement, dans postgis, l'utilisation des index spatiaux n'etait pas automatique lors de clauses where, mais necessitait d'utiliser l'operateur &&. Lors de la mise en conformité des noms de fonctions, il a ete décidé d'ajouter en plus, dans les fonctions st_*, un appel automatique a l'index, histoire de rendre nos vie plus faciles ( )

Concretement, si vous faites des selects sans preciser de clause where ou de jointure spatiale, aucune différence entre les deux notations.
Dans l'autre cas: jointures spatiales, alors ne pas utiliser les index, c'est prendre le risque de requetes prenant beaucoup de temps (qd je dis beaucoup, un index spatial peut faire tomber le temps de requete de plusieurs heures a 1s).

Dans votre cas, la creation d'index spatiaux sur chaque colonne geo de vos tables est indispensable.

Nicolas

Dernière modification par Nicolas Ribot (Mon 24 October 2011 21:34)

Nicolas,

Merci beaucoup pour ta réponse ! Ca fonctionne très bien et ça m'a permis de faire ma première requête spatiale complexe ;-)

J'ai utilisé la requête de ton deuxième message et je l'ai réadaptée à mes besoins. Petite question pour que je comprenne bien, le order de la première requête permet de classer les zones par distance et quel est l'élément permettant de ne prendre que la zone la plus proche ppur un point de collecte ?

Merci