<div dir="auto">Any ideas for this question?</div><br><div class="gmail_quote"><div dir="ltr">On Sun, Jul 1, 2018, 13:27 Francois Chartier, <<a href="mailto:fra.chartier@gmail.com">fra.chartier@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div lang="EN-CA" link="blue" vlink="#954F72"><div class="m_-7514966571087167723WordSection1"><p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Arial",sans-serif;color:#222222;background:white">Hi,<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Arial",sans-serif;color:#222222;background:white"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Arial",sans-serif;color:#222222;background:white">I am reposting the initial question: </span><u></u><u></u></p><div><p class="MsoNormal" style="background:white"><span style="font-size:10.0pt;font-family:"Arial",sans-serif;color:#222222"><u></u> <u></u></span></p><div><p class="MsoNormal" style="background:white"><span style="font-size:10.0pt;font-family:"Arial",sans-serif;color:#222222">“I am working with a data set that consists of borehole logs with a Top of a layer (layer 1) and the top of the underlying layer (layer 2) (which is also the bottom of the overlying layer 1).  Everything in between the elevation of top of layer 1 and top of layer 2 correspond to a Layer 1 property.  <u></u><u></u></span></p></div><div><p class="MsoNormal" style="background:white"><span style="font-size:10.0pt;font-family:"Arial",sans-serif;color:#222222"><u></u> <u></u></span></p></div><div><p class="MsoNormal" style="background:white"><span style="font-size:10.0pt;font-family:"Arial",sans-serif;color:#222222">The thickness of Layer 1 varies and this layer may not exist everywhere (pinches out).  Above the Layer the property is different; in other words the property only starts below the Top of layer 1 until the underlying Top of the next layer.<u></u><u></u></span></p></div><div><p class="MsoNormal" style="background:white"><span style="font-size:10.0pt;font-family:"Arial",sans-serif;color:#222222">Not sure of the capabilities of the interpolation in Grass and working with a very large data set (i cannot link every top of layers together), my first approach was to create a each Layer property for every elevation slice along each Borehole axis, interpolating soil properties at every elevation between Boreholes. <u></u><u></u></span></p></div><div><p class="MsoNormal" style="background:white"><span style="font-size:10.0pt;font-family:"Arial",sans-serif;color:#222222"><br><span style="background:white">The key question is can v.vol.rst </span></span><span style="font-size:10.0pt;font-family:"Arial",sans-serif;color:red;background:white">(changed this from <a href="http://r3.in.xyz" target="_blank" rel="noreferrer">r3.in.xyz</a>)</span><span style="font-size:10.0pt;font-family:"Arial",sans-serif;color:#222222;background:white"> interpolate in 3D without a Property at every elevation slice, </span><span style="font-size:10.0pt;font-family:"Arial",sans-serif;color:#222222"><u></u><u></u></span></p></div><div><p class="MsoNormal" style="margin-left:47.25pt;background:white"><u></u><span style="font-size:10.0pt;font-family:Symbol;color:#222222"><span>·<span style="font:7.0pt "Times New Roman""> </span></span></span><u></u><span style="font-size:10.0pt;font-family:"Arial",sans-serif;color:#222222">while respecting the condition that above the Top of the layer 2, the property corresponds to the Overlying top layer 1, and <u></u><u></u></span></p><p class="MsoNormal" style="margin-left:47.25pt;background:white"><u></u><span style="font-size:10.0pt;font-family:Symbol;color:#222222"><span>·<span style="font:7.0pt "Times New Roman""> </span></span></span><u></u><span style="font-size:10.0pt;font-family:"Arial",sans-serif;color:#222222">that the property is continuous until the next underlying layer 3 - can someone confirm this?  <u></u><u></u></span></p><div><p class="MsoNormal" style="background:white"><span style="font-size:10.0pt;font-family:"Arial",sans-serif;color:#222222">To provide a bit of background, borehole data bases, identify the top of layer as encountered when drilling downwards, and provide the elevation of the next layer (pick); in between the soil property is the same, however there is no data points.  When interpolating, while there is no data point in between the two geological picks, the property should still have weight in the interpolation process.”<u></u><u></u></span></p></div><div><p class="MsoNormal" style="background:white"><span style="font-size:10.0pt;font-family:"Arial",sans-serif;color:#222222"><u></u> <u></u></span></p></div><div><p class="MsoNormal" style="background:white"><span style="font-size:10.0pt;font-family:"Arial",sans-serif;color:#222222">Thank You <u></u><u></u></span></p></div></div></div><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">Sent from <a href="https://go.microsoft.com/fwlink/?LinkId=550986" target="_blank" rel="noreferrer">Mail</a> for Windows 10</p><p class="MsoNormal"><u></u> <u></u></p><div style="border:none;border-top:solid #e1e1e1 1.0pt;padding:3.0pt 0cm 0cm 0cm"><p class="MsoNormal" style="border:none;padding:0cm"><b>From: </b><a href="mailto:wenzeslaus@gmail.com" target="_blank" rel="noreferrer">Vaclav Petras</a><br><b>Sent: </b>June 28, 2018 8:27 PM<br><b>To: </b><a href="mailto:fra.chartier@gmail.com" target="_blank" rel="noreferrer">Francois Chartier</a><br><b>Cc: </b><a href="mailto:mlennert@club.worldonline.be" target="_blank" rel="noreferrer">Moritz Lennert</a>; <a href="mailto:grass-user@lists.osgeo.org" target="_blank" rel="noreferrer">GRASS user list</a><br><b>Subject: </b>Re: [GRASS-user] <a href="http://r3.in.xyz" target="_blank" rel="noreferrer">r3.in.xyz</a></p></div><p class="MsoNormal"><u></u> <u></u></p><div><div><p class="MsoNormal">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)?</p></div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal"><a href="https://grass.osgeo.org/grass74/manuals/r.to.rast3elev.html" target="_blank" rel="noreferrer">https://grass.osgeo.org/grass74/manuals/r.to.rast3elev.html</a></p></div><div><p class="MsoNormal"><a href="https://grass.osgeo.org/grass74/manuals/r3.mapcalc.html" target="_blank" rel="noreferrer">https://grass.osgeo.org/grass74/manuals/r3.mapcalc.html</a></p></div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal">Best,</p></div><div><p class="MsoNormal">Vaclav</p></div></div><div><p class="MsoNormal"><u></u> <u></u></p><div><p class="MsoNormal">On Thu, Jun 28, 2018 at 12:18 PM, Francois Chartier <<a href="mailto:fra.chartier@gmail.com" target="_blank" rel="noreferrer">fra.chartier@gmail.com</a>> wrote:</p><blockquote style="border:none;border-left:solid #cccccc 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-top:5.0pt;margin-right:0cm;margin-bottom:5.0pt"><div><p class="MsoNormal" style="margin-left:4.8pt">Is it possible to define a horizontal constraints so that more weight is given on the interpolation of the property horizontally rather than vertically. </p><div><p class="MsoNormal" style="margin-left:4.8pt">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.  </p></div><div><p class="MsoNormal" style="margin-left:4.8pt">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.</p></div></div><div><div><div><p class="MsoNormal" style="margin-left:4.8pt"><u></u> <u></u></p><div><p class="MsoNormal" style="margin-left:4.8pt">On Jun 28, 2018 03:40, "Moritz Lennert" <<a href="mailto:mlennert@club.worldonline.be" target="_blank" rel="noreferrer">mlennert@club.worldonline.be</a>> wrote:</p><blockquote style="border:none;border-left:solid #cccccc 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-top:5.0pt;margin-right:0cm;margin-bottom:5.0pt"><div><p class="MsoNormal" style="margin-right:0cm;margin-bottom:12.0pt;margin-left:9.6pt">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. </p></div><p class="MsoNormal" style="margin-left:9.6pt">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</p><div><p class="MsoNormal" style="margin-left:9.6pt"><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, </p></div><div><p class="MsoNormal" style="margin-left:9.6pt">> <<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" target="_blank" rel="noreferrer">r3.in.xyz</a> <<a href="http://r3.in.xyz" target="_blank" rel="noreferrer">http://r3.in.xyz</a>><br>>     <<a href="http://r3.in.xyz" target="_blank" rel="noreferrer">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>> </p></div><p class="MsoNormal" style="margin-left:9.6pt">>     I don't think <a href="http://r3.in.xyz" target="_blank" rel="noreferrer">r3.in.xyz</a> <<a href="http://r3.in.xyz" target="_blank" rel="noreferrer">http://r3.in.xyz</a>> is what you need. This</p><div><p class="MsoNormal" style="margin-right:0cm;margin-bottom:12.0pt;margin-left:9.6pt"><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></p></div></blockquote></div><p class="MsoNormal" style="margin-left:4.8pt"><u></u> <u></u></p></div></div></div><p class="MsoNormal" style="margin-left:4.8pt"><br>_______________________________________________<br>grass-user mailing list<br><a href="mailto:grass-user@lists.osgeo.org" target="_blank" rel="noreferrer">grass-user@lists.osgeo.org</a><br><a href="https://lists.osgeo.org/mailman/listinfo/grass-user" target="_blank" rel="noreferrer">https://lists.osgeo.org/mailman/listinfo/grass-user</a></p></blockquote></div></div><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal"><u></u> <u></u></p></div><div id="m_-7514966571087167723DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2">
<br><br>
<hr style="border:none;color:#909090;background-color:#b0b0b0;height:1px;width:99%">
<table style="border-collapse:collapse;border:none">
        <tr>
                <td style="border:none;padding:0px 15px 0px 8px">
                        <a href="https://www.avast.com/antivirus" target="_blank" rel="noreferrer">
                                <img border="0" src="https://static.avast.com/emails/avast-mail-stamp.png" alt="Avast logo">
                        </a>
                </td>
                <td>
                        <p style="color:#3d4d5a;font-family:"Calibri","Verdana","Arial","Helvetica";font-size:12pt">
                                This email has been checked for viruses by Avast antivirus software.
                                <br><a href="https://www.avast.com/antivirus" target="_blank" rel="noreferrer">www.avast.com</a>
                        </p>
                </td>
        </tr>
</table>
<br>
<a href="#m_-7514966571087167723_DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2" width="1" height="1" rel="noreferrer"> </a></div></div></blockquote></div>