I have tried it and it works. Thank you Pierre.<br><br><div class="gmail_quote">On 22 November 2011 15:33, Pierre Racine <span dir="ltr"><<a href="mailto:Pierre.Racine@sbf.ulaval.ca">Pierre.Racine@sbf.ulaval.ca</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><div class="im">> I have used ST_MapAlgebraExpr but I get the error below.<br>
><br>
><br>
> select * from (select ST_MapAlgebraExpr(kb.rast, kb1.rast, 'kb1.rast+kb.rast') as<br>
> rast from kb_dsm_clip kb, kb_dsm_clip100 kb1) as r1 limit 5;<br>
> ERROR: missing FROM-clause entry for table "kb1"<br>
> LINE 1: SELECT (KB1.RAST+KB.RAST)::double precision ^<br>
> QUERY: SELECT (KB1.RAST+KB.RAST)::double precision<br>
<br>
</div>Try this instead:<br>
<br>
SELECT ST_MapAlgebraExpr(kb.rast, kb1.rast, 'rast1 + rast2') as rast<br>
FROM kb_dsm_clip kb, kb_dsm_clip100 kb1<br>
WHERE st_intersects(kb.rast, kb1.rast)<br>
<br>
1) There is no need for a subquery<br>
<br>
2) You must refer to the raster using 'rast1' and 'rast2', not their real names<br>
<br>
3) I don't know if your rasters are tiled but it is always a good idea to restrict the process to intersecting footprints.<br>
<span class="HOEnZb"><font color="#888888"><br>
Pierre<br>
</font></span><div class="HOEnZb"><div class="h5">_______________________________________________<br>
postgis-users mailing list<br>
<a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a><br>
<a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_blank">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>
</div></div></blockquote></div><br>