Thank you very much for the explanations.<br><br>Indeed, with a St_Value, I can assess that the request worked, but without opening the result in ArcGIS (its plant), this notice message make me doubt.<br><br>Samy<br><div class="gmail_quote">
2013/2/20 dustymugs <span dir="ltr"><<a href="mailto:dustymugs@gmail.com" target="_blank">dustymugs@gmail.com</a>></span><br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
The only problem you have in your query is the expression being provided<br>
to ST_MapAlgebraExpr().  That expression should be:<br>
<br>
[rast1.val]*[rast2.val]<br>
<br>
The notice itself can safely be ignored.  In the normal use case, your<br>
source rasters are tiled into smaller chucks.  And between these smaller<br>
chucks, some of these tiles will not intersect.<br>
<br>
+---+---+---+<br>
|   |   |   |<br>
+---+---+---+<br>
|   |   |   |<br>
+---*---*---*---+<br>
|   |   |   |   |<br>
+---*---*---*---+<br>
    |   |   |   |<br>
    +---+---+---+<br>
    |   |   |   |<br>
    +---+---+---+<br>
<br>
If you were to look at the text figure above, the two rasters intersect.<br>
 But some of the smaller chucks when compared would not.<br>
<br>
-bborie<br>
<div><div class="h5"><br>
On 02/19/2013 03:30 PM, Silver Ashen wrote:<br>
> Hello everybody,<br>
><br>
> I'm trying to test the St_MapAlgebraExpr 2 raster version in order to<br>
> multiply the values of a raster by the values of an other one.<br>
><br>
> Here is the request :<br>
><br>
> SELECT St_MapAlgebraExpr(a.rast,b.rast,'(a.rast*b.rast)::double<br>
> precision','32BF','INTERSECTION','NULL','NULL',-1) As rast<br>
> FROM sigma.mnt_5m_2941x6197 as a, sigma.affor_aligned2 as b WHERE<br>
> St_Intersects(a.rast,b.rast);<br>
><br>
> Note that affor_aligned2 was a vector rasterized with St_AsRaster in the<br>
> way that it's fit the alignment of mnt_5m_2941x6197.<br>
><br>
> However, the function above returns the following :<br>
><br>
> "Notice : The two rasters provided have no intersection.  Returning no band<br>
> raster"<br>
><br>
> I don't understand because when I look to the two raster in ArcGIS, there<br>
> is no problem of srid. The two layers are well overlayed...<br>
><br>
>  I don't understand what I did wrong. Any idea ?<br>
><br>
> Thanks,<br>
> Samy<br>
><br>
><br>
><br>
</div></div>> _______________________________________________<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>
><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>