<div dir="ltr"><div>Hi, I'm not following this discussion fully, but don't you need just r.to.rast3elev followed by some raster algebra with ifs (r3.mapcalc)?<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><br></div><div><a href="https://grass.osgeo.org/grass74/manuals/r3.mapcalc.html">https://grass.osgeo.org/grass74/manuals/r3.mapcalc.html</a></div><div><br></div><div>Best,</div><div>Vaclav<br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Jun 28, 2018 at 12:18 PM, Francois Chartier <span dir="ltr"><<a href="mailto:fra.chartier@gmail.com" target="_blank">fra.chartier@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="auto">Is it possible to define a horizontal constraints so that more weight is given on the interpolation of the property horizontally rather than vertically. <div dir="auto">Similarly in kriging we can define a semi variogram with an orientation in 2d, and i would like to know if it can be done in 3d. </div><div dir="auto">Also what about the question regarding the elevations between the top and bottom of the soil formation. Is there a possibilty for v vol rst to consider data points in between top and bottom.</div></div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><div class="gmail_quote">On Jun 28, 2018 03:40, "Moritz Lennert" <<a href="mailto:mlennert@club.worldonline.be" target="_blank">mlennert@club.worldonline.be</a>> wrote:<br type="attribution"><blockquote class="m_-1109762564626777105quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="m_-1109762564626777105quoted-text">On 27/06/18 16:32, Francois Chartier wrote:<br>
> So therefore, my approach of generating a soil property at each <br>
> elevation and interpolating is a correct approach.<br>
> Is there a way that i can constrain the interpolation horizontally. <br>
<br></div>
Not sure what you mean by constraining the interpolation horizontally.<br>
<br>
You can create a 2D elevation map by interpolating all your height <br>
information of layer 1 (for example using v.surf.rst). Then you create a <br>
second elevation map interpolating the layer 2 height information, etc.<br>
<br>
Moritz<div class="m_-1109762564626777105quoted-text"><br>
<br>
<br>
> I. <br>
> Cant remember the module name i am using but i am generating the 3d <br>
> raster from the 3d data points directly.<br>
> <br>
> <br>
> On Tue, Jun 26, 2018, 23:38 Moritz Lennert, <br></div><div class="m_-1109762564626777105elided-text">
> <<a href="mailto:mlennert@club.worldonline.be" rel="noreferrer" target="_blank">mlennert@club.worldonline.be</a> <mailto:<a href="mailto:mlennert@club.worldonline.be" rel="noreferrer" target="_blank">mlennert@club.<wbr>worldonline.be</a>>> wrote:<br>
> <br>
> On 27/06/18 03:59, Francois Chartier wrote:<br>
> > Hi,<br>
> ><br>
> > I am working with a data set that consists of borehole logs with<br>
> a Top<br>
> > of a layer (layer 1) and the top of the underlying layer (layer 2)<br>
> > (which is also the bottom of the overlying layer 1). Everything in<br>
> > between the elevation of top of layer 1 and top of layer 2<br>
> correspond to<br>
> > a Layer 1 property.<br>
> ><br>
> > The thickness of Layer 1 varies and this layer may not exist<br>
> everywhere<br>
> > (pinches out). Above the Layer the property is different; in other<br>
> > words the property only starts below the Top of layer 1 until the<br>
> > underlying Top of the next layer.<br>
> > Not sure of the capabilities of the interpolation in Grass and<br>
> working<br>
> > with a very large data set (i cannot link every top of layers<br>
> together),<br>
> > my first approach was to create a each Layer property for every<br>
> > elevation slice along each Borehole axis, interpolating soil<br>
> properties<br>
> > at every elevation between Boreholes.<br>
> ><br>
> > I read on another forum that <a href="http://r3.in.xyz" rel="noreferrer noreferrer" target="_blank">r3.in.xyz</a> <<a href="http://r3.in.xyz" rel="noreferrer noreferrer" target="_blank">http://r3.in.xyz</a>><br>
> <<a href="http://r3.in.xyz" rel="noreferrer noreferrer" target="_blank">http://r3.in.xyz</a>> can<br>
> > interpolate in 3D without a Property at every elevation slice,<br>
> ><br>
> > * while respecting the condition that above the Top of the layer 2,<br>
> > the property corresponds to the Overlying top layer 1, and<br>
> > * that the property is continuous until the next underlying<br>
> layer 3 -<br>
> > can someone confirm this?<br>
> ><br>
> > To provide a bit of background, borehole data bases, identify the<br>
> top of<br>
> > layer as encountered when drilling downwards, and provide the<br>
> elevation<br>
> > of the next layer (pick); in between the soil property is the same,<br>
> > however there is no data points. When interpolating, while there<br>
> is no<br>
> > data point in between the two geological picks, the property should<br>
> > still have weight in the interpolation process.<br>
> ><br>
> <br></div>
> I don't think <a href="http://r3.in.xyz" rel="noreferrer noreferrer" target="_blank">r3.in.xyz</a> <<a href="http://r3.in.xyz" rel="noreferrer noreferrer" target="_blank">http://r3.in.xyz</a>> is what you need. This<div class="m_-1109762564626777105elided-text"><br>
> module aggregates 3D<br>
> point data into voxels.<br>
> <br>
> You probably want to use something like r.to.rast3. This would mean<br>
> interpolating each layer separately into 2D elevation maps and then<br>
> assemble them into 3D using r.to.rast3.<br>
> <br>
> Moritz<br>
> <br>
<br>
<br>
</div></blockquote></div><br></div>
</div></div><br>______________________________<wbr>_________________<br>
grass-user mailing list<br>
<a href="mailto:grass-user@lists.osgeo.org">grass-user@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/grass-user" rel="noreferrer" target="_blank">https://lists.osgeo.org/<wbr>mailman/listinfo/grass-user</a><br></blockquote></div><br></div>