> Any suggestions why the results are not the same? Are your nodata values well defined? Might also be, depending on the version you have, the order in which you pass the geometry and the raster to the ST_Intersects() function. Try inverting them. > Any improvements to the query maybe? Not really. Pierre