<div dir="auto">Thanks moritz.  I will try this. I have about 300 m of elevation to interpolate so i might group by thickness bands of 2 metres etc.  <div dir="auto">Is there a way i can make a batch mode where i prep all the xy p by elevation and then interpolate automatically onw after another all slices?</div></div><br><div class="gmail_quote"><div dir="ltr">On Mon, Jul 9, 2018, 03:01 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">Le Sun, 1 Jul 2018 13:27:03 -0400,<br>
Francois Chartier <<a href="mailto:fra.chartier@gmail.com" target="_blank" rel="noreferrer">fra.chartier@gmail.com</a>> a écrit :<br>
<br>
> Hi,<br>
> <br>
> I am reposting the initial question: <br>
> <br>
> “I am working with a data set that consists of borehole logs with a<br>
> Top of a layer (layer 1) and the top of the underlying layer (layer<br>
> 2) (which is also the bottom of the overlying layer 1).  Everything<br>
> in between the elevation of top of layer 1 and top of layer 2<br>
> correspond to a Layer 1 property.  <br>
> <br>
> The thickness of Layer 1 varies and this layer may not exist<br>
> everywhere (pinches out).  Above the Layer the property is different;<br>
> in other words the property only starts below the Top of layer 1<br>
> until the underlying Top of the next layer. Not sure of the<br>
> capabilities of the interpolation in Grass and working with a very<br>
> large data set (i cannot link every top of layers together), my first<br>
> approach was to create a each Layer property for every elevation<br>
> slice along each Borehole axis, interpolating soil properties at<br>
> every elevation between Boreholes. <br>
> <br>
> The key question is can v.vol.rst (changed this from <a href="http://r3.in.xyz" rel="noreferrer noreferrer" target="_blank">r3.in.xyz</a>)<br>
> interpolate in 3D without a Property at every elevation slice, <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? To provide a bit of background, borehole<br>
> data bases, identify the top of layer as encountered when drilling<br>
> downwards, and provide the elevation of the next layer (pick); in<br>
> between the soil property is the same, however there is no data<br>
> points.  When interpolating, while there is no data point in between<br>
> the two geological picks, the property should still have weight in<br>
> the interpolation process.”<br>
> <br>
<br>
v.vol.rst is a tool to interpolate a continuous value into a 3D space.<br>
So, I do not think that it is the tool you are looking for.<br>
<br>
As already mentioned before:<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>
Have you tried this approach ?<br>
<br>
Moritz<br>
</blockquote></div>