utilisation d'un index ou pas par le planificateur.

Bonjour,
Nous sommes en postgres 12 / postgis 3 : POSTGIS="3.0.3 6660953" [EXTENSION] PGSQL="120" GEOS="3.9.1-CAPI-1.14.2" PROJ="8.0.1" GDAL="GDAL 3.2.3, released 2021/04/27" LIBXML="2.9.7" LIBJSON="0.13.1" LIBPROTOBUF="1.3.0" WAGYU="0.4.3 (Internal)" RASTER
Mon souci viens d'une utilisation de la fonction st_value qui permet de ramener les informations du pixel d'un raster qui correspond à un point (geom). J'aurais pu a poser sur un forum postgis mais le questionnement viens au niveau du planificateur qui refuse d'utiliser un index présent sauf si je mets une clause limit (même avec une très grande valeur) :

Requete 1 (sans clause limit) :
explain
with point_emprise as
(select x,y, ST_SetSRID(st_makepoint(x,y),2154) as geom from extract_complet ),
donnee_source as
(select x,y, st_value(rast,geom)as valeur from geo_raster_france_90m_clc,point_emprise
where st_intersects(rast,geom))
select * from donnee_source

Retour planificateur :
"QUERY PLAN"
"Nested Loop (cost=0.00..161026468635.83 rows=16427665024 width=24)"
" Join Filter: st_intersects(geo_raster_france_90m_clc.rast, st_setsrid(st_makepoint(extract_complet.x, extract_complet.y), 2154), NULL::integer)"
" -> Seq Scan on geo_raster_france_90m_clc (cost=0.00..2350.83 rows=49283 width=264)"
" -> Materialize (cost=0.00..38511.00 rows=1000000 width=16)"
" -> Seq Scan on extract_complet (cost=0.00..28628.00 rows=1000000 width=16)"

Requête 2 (la même avec une clause limit) :
explain
with point_emprise as
(select x,y, ST_SetSRID(st_makepoint(x,y),2154) as geom from extract_complet limit 10000000 ),
donnee_source as
(select x,y, st_value(rast,geom)as valeur from geo_raster_france_90m_clc,point_emprise
where st_intersects(rast,geom))
select * from donnee_source

Retour planificateur :
"QUERY PLAN"
"Nested Loop (cost=0.28..17539324.25 rows=16427665 width=24)"
" -> Limit (cost=0.00..528628.00 rows=1000000 width=48)"
" -> Seq Scan on extract_complet (cost=0.00..528628.00 rows=1000000 width=48)"
" -> Index Scan using convexhull_idx_foretv2 on geo_raster_france_90m_clc (cost=0.28..12.87 rows=2 width=264)"
" Index Cond: ((rast)::geometry && (st_setsrid(st_makepoint(extract_complet.x, extract_complet.y), 2154)))"
" Filter: _st_intersects((st_setsrid(st_makepoint(extract_complet.x, extract_complet.y), 2154)), rast, NULL::integer)"

Le raster est tuilé (près de 50 000 tuiles) et son utilisation est donc la logique même si on ne veut pas parcourir 50000 plus petit rasters au lieu d'un seul (sa non utilisation rend la requête inexploitable).
J'ai bien sûr effectué un vacuum analyze et je ne sais pas quoi faire pour que la requête utilise bien cet index.
Précisions :

· La table extract_tmp est une table temporaire (j'ai aussi essayé d'y passer un vacuum analyze sans succès)

· Le raster geo_raster_france_90m_clc proviens d'une restauration depuis une base postgis 2.5.5
Merci pour vos réponses /conseils.
[INRA]<http://www.inra.fr/>

Alain BENARD

Administrateur Systèmes d'Information

alain(dot)benard(at)inra(dot)fr<mailto:alain(dot)benard(at)inra(dot)fr>

Unité SILVA

Tél. : +33 3 83 39 40 84

Route d'amance

54 280 Champenoux

inra.fr<http://www.inra.fr>

[INRA]<https://youtu.be/xUHmi6JxDI4>