+1 for GRASS<div><br></div><div>Hartanto Sanjaya has a nice slide for calculating a basic volume for a lake using GRASS.</div><div><br></div><div>May help if you are unfamiliar with the tools available in GRASS.</div><div>

<br></div><div><a href="http://hartanto.wordpress.com/2007/06/16/analisis-3-dimensi-untuk-volume-air-danau/">http://hartanto.wordpress.com/2007/06/16/analisis-3-dimensi-untuk-volume-air-danau/</a></div><div><br></div><div>

mars<br><div><br><div class="gmail_quote">On Tue, Nov 15, 2011 at 11:01 AM, Saber Razmjooei <span dir="ltr"><<a href="mailto:razmjooeis@faunalia.co.uk">razmjooeis@faunalia.co.uk</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">

I had this problem in the past. Unfortunately there are not many plugins<br>
for raster analysis in QGIS and I usually I end up doing in in GRASS.<br>
The command in GRASS for inspection your raster layers, using polygon<br>
vectors:<br>
<a href="http://grass.osgeo.org/gdp/html_grass64/v.what.rast.html" target="_blank">http://grass.osgeo.org/gdp/html_grass64/v.what.rast.html</a><br>
<br>
You need to run r.null to set the unwanted values to null (i.e. positive<br>
or negative values).<br>
<br>
<br>
Cheers<br>
<font color="#888888">Saber<br>
</font><div><div></div><div class="h5"><br>
On Tue, 2011-11-15 at 18:57 +0100, Bernhard Draeyer wrote:<br>
> I have 2 high resolution raster DEM's before and after a recent flood<br>
> event. The aim would be to get for different areas a) erosion/negative<br>
> volumes b) accumulation/positive volumes and net volume.<br>
><br>
> First of all, i subtracted the actual surface from the 'before'-surface<br>
> to get a difference surface with positive (accumulation) and negative<br>
> (erosion) values.<br>
> Then i created a shape layer with polygons of different areas of<br>
> interest. For each of them i need seperate volumetrics. Ideally, the<br>
> polygons should contain in the end attribute values for accumulation,<br>
> erosion and net volumes.<br>
> Now i was looking without success for a well adapted method to solve my<br>
> problem in QGIS or SAGA.<br>
><br>
> How should i proceed? Is there any method in QGIS to sum up positive and<br>
> negative raster values under polygons and put them in seperate attribute<br>
> fields per polygon? Do i have to convert the raster in a point layer? I<br>
> saw, that in esri's 3d analyst exist a tin polygon volume function, that<br>
> covers this task. Does anything similar exist with an os gis?<br>
><br>
><br>
> Many thanks for your help.<br>
><br>
> Bernie<br>
><br>
> _______________________________________________<br>
> Qgis-user mailing list<br>
> <a href="mailto:Qgis-user@lists.osgeo.org">Qgis-user@lists.osgeo.org</a><br>
> <a href="http://lists.osgeo.org/mailman/listinfo/qgis-user" target="_blank">http://lists.osgeo.org/mailman/listinfo/qgis-user</a><br>
<br>
<br>
_______________________________________________<br>
Qgis-user mailing list<br>
<a href="mailto:Qgis-user@lists.osgeo.org">Qgis-user@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/mailman/listinfo/qgis-user" target="_blank">http://lists.osgeo.org/mailman/listinfo/qgis-user</a><br>
</div></div></blockquote></div><br></div></div>