Bonjour,

Dans une base Postgis, j'ai une couche raster et une couche vecteur (ligne).
Je voudrais récupérer les valeurs des pixels du raster qui intersectent ma ligne :
voici ma requête :
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8
9
10
SELECT sampno, 
(ST_Intersection(r.no2, traj.geom)).geom road, 
(ST_Intersection(r.no2, traj.geom)).val elevation,
r.rid
from output_1_t100_27572 as r,
	(SELECT sampno,  ST_Transform(geom, 27572) as geom
	from trajectories.corrected_gps_info_line_t
	where sampno='1108964') as traj
where r.rid=27
and ST_Intersects(traj.geom, r.no2)
r est ma couche raster et traj ma couche vecteur.

Cette requête fonctionne, mais elle met plusieurs minutes à rendre le résultat !
Est-ce qu'il y a un moyen de l'écrire autrement pour optimiser le temps de réponse ? comme c'est pour afficher sur une interface cartographique web, il faudrait que le temps de réponse soit 'raisonnable'

Merci,
Nico