Le portail francophone de la géomatique

Merci de ta réponse cquest.
J'ai dû mal m'exprimer et je pense que ce que je n'ai pas pu faire comprendre ce que je cherche à faire.
Je vais essayer de le mettre en "image".
Soit une table contenant des points relevés dans des villes de France. Certaines villes contiennent beaucoup de points, d'autres n'en contiennent qu'un.
Admettons, pour l'image, qu'il y ait 1000 points à Marseille, 10.000 points à Paris, 5000 points à Renne, 1 point à Brest, 5000 points à Lyon, ...
Je ne tiens pas compte des autres données que les coordonnées.
Ce que je cherche à faire est dans PostGIS, une vue qui calculerait le point moyen d'un nuage d'un rayon défini (disons 50km de rayon pour l'image des points de France).
À ce que j'ai vu, ST_ClusterWithin part d'un point connu et retournerait tous les points qui sont dans le rayon donné. Pour chaque ville, il me semble (à moins que j'ai mal compris. Ce qui est tout à fait possible) qu'il retournerait un nuage autour de chaque point. Par exemple, pour Marseille j'aurais 1000 nuages avec sensiblement les mêmes points.
Or c'est l'inverse que je cherche. Je cherche à déterminer les nuages et à en calculer le point moyen.
Le ST_Centroid me retournerai le point moyen de chaque nuage, mais je n'arrive justement pas à déterminer les nuages.
En somme, je cherche à faire une requête requête pourrait déterminer qu'à Lyon il y a un nuage de points distants de moins de 50km, de même qu'à Paris, Marseille, etc ...
La requête retournait des ensembles de points (nuages) distants de moins de la tolérance (le 50km de l'exemple) à partir de la table des points de France. En entrée de requête j'aurais la table des points de France, et en sortie, chaque ligne représenterait un nuage de points. Par exemple une géométrie multipoint (par exemple) pour Marseille, une pour Paris, une pour Renne, une pour Brest et une pour Lyon.
Ainsi je pourrais applique un ST_Centroid sur chaque ligne et trouver le point moyen de chaque nuage.

Merci de votre aide.

Bonsoir,

Je pense que cquest a raison: regardez du coté de st_clusterWithin:

Ca construit un array de geometries qui sont a une distance données les unes des autres: le centroid de cette geometrie vous en donne le centre:

select st_centroid(unnest(st_clusterWithin(geometry, 2000)))
from testpt;

Sur l'image: 35000 points sur des communes, "résumés" en 651 points, centres des clusters initiaux à 2km.

Nicolas

Bonjour,

Une autre solution consiste à calculer le barycentre d'un ensemble de points.
X = somme de 1 à n des (xi) divisé par n (idem Y).

Il suffit de définir l'ensemble en fonction d'une géométrie, d'un rayon, d'un tampon bref ce que l'on veut.

L'avantage de st_clusterwithin est de proposer une agrégation automatique à partir d'un rayon, mais comme le souligne Tumasgiu il est possible et même certain (en fonction du rayon) qu'un point appartienne à plusieurs grappes. Il faut noter aussi qu'il n'est pas possible a priori de récupérer les id des points "grappés", donc pas possible de remonter aux informations liées aux points.

Petite précision la formule donnée plus haut est en fait
G(X,Y) = (somme1àn (AiXi)/somme(Ai) , somme1àn (AiYi)/somme(Ai))

Où Ai sont les poids attribués aux points. (autre avantage de cette méthode).

(sorry tomasgiu: j'ai envoyé ma réponse hier sans rafraîchir la page...)