<div dir="ltr">Hi Tim,<div>i am the second mentor in your SoC project, and will try to help you with the 3D raster approach in GRASS.</div><div class="gmail_extra"><br></div><div class="gmail_extra">The 3D raster implementation in GRASS is using a tiled based storage approach.[1]</div>
<div class="gmail_extra">The tiles can be stored in compressed form (zlib in GRASS7) or uncompressed form on the hard disk. Usually you store them in compressed form. Hence in case you have many tiles that are empty or have equal data, the compression is very high. GRASS can handle 3D raster maps up to hundreds of GB in size. I tested 3D raster maps with about 12 GB in size. The main memory is not an issue when processing 3D raster maps, since only a limited number of tiles is read into main memory at runtime. You can specify the tile size of each 3D raster map that you create. Hence you can adjust it to the needs of the soil depths resolution. <br>
In addition you can specify a 3D raster mask that works pretty the same way as masks in the raster approach in GRASS.<br><br>Please don't bother to ask, if you have any questions about the 3D raster implementation.<br>
<br>Best regards</div><div class="gmail_extra">Soeren<br><br>[1] <a href="http://grass.osgeo.org/programming7/raster3dlib.html">http://grass.osgeo.org/programming7/raster3dlib.html</a><br><br><div class="gmail_quote">2013/6/25 Tim Bailey <span dir="ltr"><<a href="mailto:timibly@gmail.com" target="_blank">timibly@gmail.com</a>></span><br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div dir="ltr"><div><div><div><div><br><br></div>Hi Ben, <br>
</div><br></div>All that I meant by mask is, in this case,  an r3 map that defines a subset of space that subsequent operations are constrained to. I guess I was just expressing worry about creating runaway data demands.The region settings act as a three dimensional bounding cube, does not seem adequate to constrain a tiling scheme with relatively sparse data.In the compromise scenario that Dylan suggested with 10 meter xy and 1 cm z resolution and 1 meter depth we would be looking at populating 10000 voxels per hectare for a flat landscape. However if we were using a simple bounding box to define the region and we had 5 meters of relief we would end up with 50000 or 60000 voxels. As the area of coverage gets bigger cost of the range in z values gets worse.  <br>
<span class=""><font color="#888888">
</font></span></div><span class=""><font color="#888888"><div><br></div><div>Tim<br></div></font></span></div>
<br>_______________________________________________<br>
grass-dev mailing list<br>
<a href="mailto:grass-dev@lists.osgeo.org">grass-dev@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/mailman/listinfo/grass-dev" target="_blank">http://lists.osgeo.org/mailman/listinfo/grass-dev</a><br></blockquote></div><br></div></div>