Hi Pierre and thanks for your reply.<div><br></div><div>My nodata values are well defined and are excluded in the summarystats operation.</div><div>I also tried inverting the order of (geom, rast) without this changing the results.</div>
<div><br></div><div>Could this be explained by the possibility that Postgis and ArcGIS handels pixels at the border of the polygon differently? i.e. that one includes pixels intersecting, while the other only includes pixels if the centroid of the pixel intersects?</div>
<div><br></div><div>Andreas<br><br><div class="gmail_quote">2012/11/15 Pierre Racine <span dir="ltr"><<a href="mailto:Pierre.Racine@sbf.ulaval.ca" target="_blank">Pierre.Racine@sbf.ulaval.ca</a>></span><br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div class="im">> Any suggestions why the results are not the same?<br>
<br>
</div>Are your nodata values well defined?<br>
<br>
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.<br>
<div class="im"><br>
> Any improvements to the query maybe?<br>
<br>
</div>Not really.<br>
<br>
Pierre<br>
<br>
_______________________________________________<br>
postgis-users mailing list<br>
<a href="mailto:postgis-users@lists.osgeo.org">postgis-users@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/cgi-bin/mailman/listinfo/postgis-users" target="_blank">http://lists.osgeo.org/cgi-bin/mailman/listinfo/postgis-users</a><br>
</blockquote></div><br></div>