J'essaie de reproduire ce processus ArcGIS dans PostGIS: http://blogs.esri.com/esri/arcgis/2012/11/13/spaghetti_and_meatballs/ . Il décrit comment diviser les points tamponnés en polygones en fonction de leurs intersections, en comptant le nombre de couches et en l'attribuant aux polygones afin de les classer. Je l'utilise pour créer une carte de densité de points approximative avec des vecteurs, et les résultats étaient étonnamment bons pour mon ensemble de données dans ArcGIS. Cependant, j'ai du mal à trouver quelque chose de réalisable dans PostGIS où j'en ai besoin pour produire des couches de densité de points dynamiques pour une carte Web.
Dans ArcGIS, j'ai simplement exécuté l'outil Intersection sur ma couche de points tamponnés pour créer les formes dont j'avais besoin.
Dans PostGIS, j'ai exécuté cette requête:
CREATE TABLE buffer_table AS
SELECT a.gid AS gid, ST_Buffer(a.geo,.003) AS geo
FROM public.pointTable a;
CREATE TABLE intersections AS
SELECT a.gid AS gid_a, b.gid AS gid_b, ST_Intersection(a.geo,b.geo) AS geo
FROM public.pointTable a, public.pointTable b
WHERE ST_Intersects(a.geo, b.geo) AND a.gid < b.gid;
DELETE FROM intersections WHERE id_a = id_b;
La sortie semble à peu près identique à la sortie ArcGIS, sauf qu'elle ne décompose pas les polygones dans la même mesure que celle requise pour une carte de densité significative. Voici des captures d'écran de ce que je veux dire:
ArcGIS est à gauche et PostGIS à droite. C'est un peu difficile à dire, mais l'image ArcGIS montre le polygone «intérieur» créé à l'intersection des 3 tampons. En revanche, la sortie PostGIS ne crée pas ce polygone intérieur et conserve ses composants intacts. Il est donc impossible de fournir une classification uniquement pour cette zone intérieure avec 3 couches les unes sur les autres, contre 1 pour les parties extérieures.
Quelqu'un connaît-il une fonction PostGIS pour décomposer le polygone autant que je le souhaite? Sinon, quelqu'un connaît-il une meilleure façon de produire une carte de densité de points avec des vecteurs dans PostGIS?