<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="gmail_extra"><br><div class="gmail_quote">On Jun 28, 2018 03:40, "Moritz Lennert" <<a href="mailto:mlennert@club.worldonline.be">mlennert@club.worldonline.be</a>> wrote:<br type="attribution"><blockquote class="quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="quoted-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="quoted-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="elided-text">
> <<a href="mailto:mlennert@club.worldonline.be" target="_blank" rel="noreferrer">mlennert@club.worldonline.be</a> <mailto:<a href="mailto:mlennert@club.worldonline.be" target="_blank" rel="noreferrer">mlennert@club.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="elided-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>