[Geobd] Soucis d'index spatial

Message posté par : Christophe Vergon (christophe.vergon(a)protonmail.com) ---------------------------------------- Bonjour, La piste de Nicolas est la bonne. Un extrait de code qui permet d'obtenir très rapidement la "vectorisation" d'un raster en carrés et points, plus valeur du raster. Ici le raster est nommé "elevation". Nota il est découpé en tuile de 100*100, conservation de l'index de la tuile (rid). ----------------- Code : with p as (SELECT rid,rast,st_squaregrid(1,st_envelope(rast)) rg FROM elevation LIMIT 100000) SELECT rid,st_value(rast,st_centroid((rg).geom)) val_rast,st_astext((rg).geom)as sgeom, st_astext(st_centroid((rg).geom)) as pgeom FROM p ORDER BY rid DESC ----------------- 100 000 pixels tratés en 2s ... Merci Nicolas. ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=363037#p363037 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