Le portail francophone de la géomatique

Merci pour le retour,

Effectivement il s'agit d'une requête sous PostgreSQL/PostGIS et il manquait bien une parenthèse..

Par contre le problème de temps de réponse provient bien à mon avis de la fonction du ST_PointOnSurface, car le ST_Centroid, ou plus simplement comme tu propose le ST_Intersects sont quasi instantanés.

"[ST_PointOnSurface] is substantially more computationally expensive than the centroid operation". cf : https://gis.stackexchange.com/questions … calculated

Et je dispose bien dans mes tables d'index géographiques du type.

CREATE INDEX geo_parcelle_geom_idx
ON temp_majic_cca_2018.geo_parcelle
USING gist (geom);

Or c'est bien la fonction ST_PointOnSurface qui m'interesse car ST_Centroid n'est pas assez précis dans le cas de parcelles dont le centroide ne tombe pas dans la parcelle

Excellente idée! merci!

13 secondes pour créer la nouvelle géométrie des parcelles

ALTER TABLE temp_majic_cca_2018.geo_parcelle ADD COLUMN geom_pos geometry(Geometry,2154);
UPDATE temp_majic_cca_2018.geo_parcelle SET geom_pos = ST_PointOnSurface (geo_parcelle.geom);
CREATE INDEX geo_parcelle_geom_pos_idx ON temp_majic_cca_2018.geo_parcelle USING gist(geom_pos);

9 secondes pour créer l'intersection entre cette nouvelle géométrie et mes unités foncières

SELECT 
  min(geo_parcelle.geo_parcelle) as id_parc,
  min(geo_unite_fonciere.id) as id_uf
FROM 
  temp_majic_cca_2018.geo_parcelle, 
  temp_majic_cca_2018.geo_unite_fonciere                     
WHERE 
  ST_Intersects(geo_parcelle.geom_pos, geo_unite_fonciere.geom)
GROUP BY geo_parcelle

Merci encore!

Dernière modification par boblebob (Mon 17 September 2018 10:46)