<div dir="auto"><div>So therefore, my approach of generating a soil property at each elevation and interpolating is a correct approach.  <div dir="auto">Is there a way that i can constrain the interpolation horizontally.  I. Cant remember the module name i am using but i am generating the 3d raster from the 3d data points directly. </div><br><br><div class="gmail_quote"><div dir="ltr">On Tue, Jun 26, 2018, 23:38 Moritz Lennert, <<a href="mailto:mlennert@club.worldonline.be">mlennert@club.worldonline.be</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">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 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 correspond to <br>
> a Layer 1 property.<br>
> <br>
> The thickness of Layer 1 varies and this layer may not exist 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 working <br>
> with a very large data set (i cannot link every top of layers together), <br>
> my first approach was to create a each Layer property for every <br>
> elevation slice along each Borehole axis, interpolating soil 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>> 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 layer 3 -<br>
>     can someone confirm this?<br>
> <br>
> To provide a bit of background, borehole data bases, identify the top of <br>
> layer as encountered when drilling downwards, and provide the 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 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>
I don't think <a href="http://r3.in.xyz" rel="noreferrer noreferrer" target="_blank">r3.in.xyz</a> is what you need. This 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>
</blockquote></div></div></div>