<div dir="ltr"><div class="gmail_default" style="font-family:monospace,monospace">Hey,<br></div><div class="gmail_default" style="font-family:monospace,monospace"><br>Map Algebra?<br></div><div class="gmail_default" style="font-family:monospace,monospace">(custom function with a thresholding?)<br></div><div class="gmail_default" style="font-family:monospace,monospace"><br><br>Cheers,<br></div><div class="gmail_default" style="font-family:monospace,monospace">Rémi-C<br></div></div><div class="gmail_extra"><br><div class="gmail_quote">2015-05-27 14:57 GMT+02:00 Gabriel Vatin <span dir="ltr"><<a href="mailto:gabriel.vatin@kinaxia.fr" target="_blank">gabriel.vatin@kinaxia.fr</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hello all,<br>
<br>
I'd like to have some feedbacks on an operation that should be possible with Postgis raster, but I can't find out how to do this.<br>
I have a whole raster of French DEM (elevation data) stored on PostGis, with tiled objects. That makes a table with some 14.000 rows of data: france_mnt (rid, rast, filename)<br>
With ST_Value(france_mnt.rast, geometry), I can retrieve the elevation of a given coordinates in the geometry object.<br>
<br>
I'd like to create a new raster, on the fly, with a maximum value of elevation, such as : all pixels with value < 10<br>
This could be used for mapping water elevation in a area, for instance.<br>
<br>
Anyone has already done this? I'll be happy to have some tips on the function to use.<br>
Thanks!<br>
<br>
Gabriel<br>
<br>
_______________________________________________<br>
postgis-users mailing list<br>
<a href="mailto:postgis-users@lists.osgeo.org" target="_blank">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>