Message posté par : Nicolas Ribot ---------------------------------------- Merci, Pour le plan, ca serait mieux la sortie texte de pg, plus facile a lire je trouve. Concernant la creation de l'index spatial, je n'ai jamais utilisé les options que tu indiques. Tu as essayé une creation simple, style: ----------------- Code : CREATE INDEX ON public.elevationp USING gist(geom); analyze elevationp; ----------------- Pour voir si ca change le plan ? Nicolas ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=362988#p362988 Pour y répondre : geobd(a)ml.georezo.net ou reply de votre messagerie Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement -- Association GeoRezo - le portail géomatique https://georezo.net

Message posté par : Nicolas Ribot ---------------------------------------- Sinon ce qui marche très vite pour récuperer des valeurs dans un raster à partir de vecteurs, c'est d'utiliser gdallocationinfo sur le raster direct: ca peut choper des centaines de milliers de valeurs par seconde et les stocker dans postgis. J'ai pas mal joué avec ca ces temps-ci et ca permet de faire des croisements massifs entre vecteurs et raster (style des linéaires sur des pays entiers vs des rasters a 30m de réso). Le principe est : • de créer une grille virtuelle correspondant au raster (un pg de la grille = 1 pixel du raster) UNIQUEMENT sous les geometries à traiter (les batiments par ex). Ca va vite en général, sauf pour des super gros volumes, style batiments de france vs raster a 1m par ex. Pour les ~ 110 000 batiments Paris sur une grille 30m, ca prend qq secondes de faire cette grille. • de prendre les centroids de cette grille, de les passer a gdallocationinfo qui renvoie les valeurs dans pg direct (ultra rapide). C'est faisable en une passe avec un peu de shell, style: ----------------- Code : psql -c "select x, y from gridtable order by row_id" | gdallocationinfo -valonly -geoloc dem.tif | psql -c "COPY gridtableval(val) from STDIN" ----------------- gridtableval est la table de résultat à créer avant qui permet de stocker la sortie de gdallocationinfo. En ajoutant une pk generated always as identity, on peut lier cette table de résultat sur la table gridtable pour avoir la valeur raster de chaque pixel de la grille • de faire le SQL qui donne la synthèse par batiment (ici le max d'alti pour pixels de chaque batiment). • S'il faut faire des analyses plus poussées, par ex prendre en compte la partie de pixel qui intersecte la geom, on travaille sur ces pixels uniquement et pas tous les pixels du raster. Avec cette méthode, des process qui prenaient des heures prennent qq minutes maintenant. Nico ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=362993#p362993 Pour y répondre : geobd(a)ml.georezo.net ou reply de votre messagerie Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement -- Association GeoRezo - le portail géomatique https://georezo.net

Message posté par : Christophe Vergon (christophe.vergon(a)protonmail.com) ---------------------------------------- Re, Pour l'index en fait il ne faut pas préciser d'opérateur : ----------------- Code : CREATE INDEX IF NOT EXISTS gix_b ON cadastre.batiment USING gist (the_geom) WITH (buffering=auto) TABLESPACE pg_default; ----------------- Là l'index est pris en compte par le query planner ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=363038#p363038 Pour y répondre : geobd(a)ml.georezo.net ou reply de votre messagerie Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement -- Association GeoRezo - le portail géomatique https://georezo.net

Message posté par : T. Rossini ---------------------------------------- Hello, est ce que tu as fait un vacuum analyze sur les deux tables ? ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=362989#p362989 Pour y répondre : geobd(a)ml.georezo.net ou reply de votre messagerie Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement -- Association GeoRezo - le portail géomatique https://georezo.net

Message posté par : Nicolas Ribot ---------------------------------------- ----------------- tumasgiu écrit : Je me souviens aussi que dans de vieilles versions de postgis, st_intersects ne provoquait pas forcement l'utilisation de l'index et qu'il fallait rajouter la condition d'intersection de bounding box a.geom && b.geom ----------------- Oui au fait, que donne postgis_full_version() ? (je vous recommande pg16 + postgis 3.4 (geos 12): ca marche bien !) ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=362994#p362994 Pour y répondre : geobd(a)ml.georezo.net ou reply de votre messagerie Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement -- Association GeoRezo - le portail géomatique https://georezo.net

Message posté par : Christophe Vergon (christophe.vergon(a)protonmail.com) ---------------------------------------- Bonjour, La solution finale. ----------------- Code : with p as (SELECT idbatiment,dur,rid,the_geom ,rast FROM elevation,cadastre.batiment WHERE st_intersects(the_geom,st_envelope(rast)) AND NOT st_bandisnodata(rast)), p1 as (SELECT idbatiment,dur,rid,rast,the_geom as rid_bat_geom FROM p ORDER BY rid), p2 as (SELECT rid,st_centroid((st_squaregrid(1,st_envelope(rast))).geom) pg FROM elevation) SELECT idbatiment,dur,p2.rid,st_value(rast,pg),pg --INTO elevation_batiment FROM p1,p2 WHERE st_intersects(rid_bat_geom,pg) AND p1.rid=p2.rid ----------------- Ce qui donne un temps d'exécution de 37.35 secondes pour un raster de 4 000 * 6 000 et 4 161 bâtiments dans l'emprise du raster. Et un rendu sympathique cf. PJ NB : Ceci permet aussi de faire du pur SQL sans passer par gdgdallocationinfo() ---------------------------------------- Ce message est accompagné de fichiers, pour les télécharger, suivre le lien ci-dessous. ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=363039#p363039 Pour y répondre : geobd(a)ml.georezo.net ou reply de votre messagerie Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement -- Association GeoRezo - le portail géomatique https://georezo.net