<div dir="ltr"><div dir="ltr"><div dir="ltr"><br><br><div class="gmail_quote"><div dir="ltr">On Wed, Oct 10, 2018 at 8:34 PM Francois Chartier <<a href="mailto:fra.chartier@gmail.com">fra.chartier@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Hi Vaclav, <div><br></div><div>Question:</div><div>1-can the value remain constant until the next underlying data point?  </div></div></blockquote><div><br></div><div>v.vol.rst is for interpolation, but a constant value can be achieved with r.to.rast3elev or r.to.rast3, but I think in your case you will likely need some prereprocessing to get the right rasters from your irregular points.<br></div><div><br></div><div><a href="https://grass.osgeo.org/grass74/manuals/r.to.rast3elev.html">https://grass.osgeo.org/grass74/manuals/r.to.rast3elev.html</a></div><div><a href="https://grass.osgeo.org/grass74/manuals/r.to.rast3.html">https://grass.osgeo.org/grass74/manuals/r.to.rast3.html</a></div><div><br></div><div>I think there was some addon or script for working with soil types (which is how your case sounds). Perhaps somebody knows more about it.<br></div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>2-by "the output is a regular grid"  does that mean that the points (vector) can be unstructured, but in the end the raster or mesh is a regular grid?</div></div></blockquote><div><br></div><div>Yes, 3D points (XYZ tuples) and (regular) 3D raster (grid).<br></div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>  if that is the case, how does the value in the raster determined?</div></div></blockquote><div><br></div><div>v.vol.rst is interpolating the value using "regularized spline with 
tension" aka RST, i.e. it is applicable to properties which are a real 
number, but it is not suitable for categorical data.</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div> Is it an interpolated value between the points that fall with in the raster cell footprint?</div></div></blockquote><div><br></div><div>The interpolation happens for the raster cells (or you can think of them as grid locations in the middle of the raster cell). Usually, there would be much more raster cells than vector points, but the RST interpolation works the same way in all cases.<br></div><div><br></div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div><br></div><div>thanks,</div><div>Francois</div></div><br><div class="gmail_quote"><div dir="ltr">Le jeu. 4 oct. 2018 à 22:48, Vaclav Petras <<a href="mailto:wenzeslaus@gmail.com" target="_blank">wenzeslaus@gmail.com</a>> a écrit :<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><br><br><div class="gmail_quote"><div dir="ltr">On Tue, Oct 2, 2018 at 9:38 PM Francois Chartier <<a href="mailto:fra.chartier@gmail.com" target="_blank">fra.chartier@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>For example, in the example on the webpage, there is a 3D raster for soil moisture; if the soil moisture is measured only for the top of a cylinder along your vertical hole, the next moisture measurement is taken at a lower depth, and the intervals (cylinder length) are not the same along that vertical axis?</div></div></blockquote><div><br></div><div><br></div><div>Yes, all are (can be) different. The values are interpolated in between. The input are arbitrary places points. The output is a regular grid. Note, you may need to explore zscale parameter of v.vol.rst because of the how your data looks like.<br></div><div><br></div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><br><div class="gmail_quote"><div dir="ltr">Le lun. 1 oct. 2018 à 22:01, Vaclav Petras <<a href="mailto:wenzeslaus@gmail.com" target="_blank">wenzeslaus@gmail.com</a>> a écrit :<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div class="gmail_quote"><div dir="ltr">On Sat, Sep 29, 2018 at 4:45 PM Francois Chartier <<a href="mailto:fra.chartier@gmail.com" target="_blank">fra.chartier@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Can Grass GIS create a 3D unstructured grid based on a 3D point data set ?</div></blockquote><div><br></div><div>All data and visualizations at the following website were created in GRASS GIS including a 3D regular grid (3D raster in GRASS GIS) interpolated from a set of 3D points (vector which is 3D, i.e. with z, in GRASS GIS). Points were soil properties measurement at different depths, 3D raster is then continuous representation of a selected property.<br></div><div><br></div><div><a href="https://geospatial.ncsu.edu/geoforall/multidim-geovis.html" target="_blank">https://geospatial.ncsu.edu/geoforall/multidim-geovis.html</a></div><div><br></div><div>Module used for interpolation was v.vol.rst:<br></div><div><br></div><div><a href="https://grass.osgeo.org/grass74/manuals/v.vol.rst.html" target="_blank">https://grass.osgeo.org/grass74/manuals/v.vol.rst.html</a><br></div><div><br></div><div>Vaclav<br></div></div></div></div></div>
</blockquote></div>
</blockquote></div></div>
</blockquote></div>
</blockquote></div></div></div></div>