<html xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40"><head><meta http-equiv=Content-Type content="text/html; charset=utf-8"><meta name=Generator content="Microsoft Word 15 (filtered medium)"><style><!--
/* Font Definitions */
@font-face
{font-family:Wingdings;
panose-1:5 0 0 0 0 0 0 0 0 0;}
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri",sans-serif;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:#954F72;
text-decoration:underline;}
.MsoChpDefault
{mso-style-type:export-only;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:72.0pt 72.0pt 72.0pt 72.0pt;}
div.WordSection1
{page:WordSection1;}
/* List Definitions */
@list l0
{mso-list-id:757559197;
mso-list-template-ids:-1;}
@list l0:level1
{mso-level-number-format:bullet;
mso-level-text:\F0B7;
mso-level-tab-stop:36.0pt;
mso-level-number-position:left;
text-indent:-18.0pt;
mso-ansi-font-size:10.0pt;
font-family:Symbol;}
@list l0:level2
{mso-level-number-format:bullet;
mso-level-text:o;
mso-level-tab-stop:72.0pt;
mso-level-number-position:left;
text-indent:-18.0pt;
mso-ansi-font-size:10.0pt;
font-family:"Courier New";
mso-bidi-font-family:"Times New Roman";}
@list l0:level3
{mso-level-number-format:bullet;
mso-level-text:\F0A7;
mso-level-tab-stop:108.0pt;
mso-level-number-position:left;
text-indent:-18.0pt;
mso-ansi-font-size:10.0pt;
font-family:Wingdings;}
@list l0:level4
{mso-level-number-format:bullet;
mso-level-text:\F0A7;
mso-level-tab-stop:144.0pt;
mso-level-number-position:left;
text-indent:-18.0pt;
mso-ansi-font-size:10.0pt;
font-family:Wingdings;}
@list l0:level5
{mso-level-number-format:bullet;
mso-level-text:\F0A7;
mso-level-tab-stop:180.0pt;
mso-level-number-position:left;
text-indent:-18.0pt;
mso-ansi-font-size:10.0pt;
font-family:Wingdings;}
@list l0:level6
{mso-level-number-format:bullet;
mso-level-text:\F0A7;
mso-level-tab-stop:216.0pt;
mso-level-number-position:left;
text-indent:-18.0pt;
mso-ansi-font-size:10.0pt;
font-family:Wingdings;}
@list l0:level7
{mso-level-number-format:bullet;
mso-level-text:\F0A7;
mso-level-tab-stop:252.0pt;
mso-level-number-position:left;
text-indent:-18.0pt;
mso-ansi-font-size:10.0pt;
font-family:Wingdings;}
@list l0:level8
{mso-level-number-format:bullet;
mso-level-text:\F0A7;
mso-level-tab-stop:288.0pt;
mso-level-number-position:left;
text-indent:-18.0pt;
mso-ansi-font-size:10.0pt;
font-family:Wingdings;}
@list l0:level9
{mso-level-number-format:bullet;
mso-level-text:\F0A7;
mso-level-tab-stop:324.0pt;
mso-level-number-position:left;
text-indent:-18.0pt;
mso-ansi-font-size:10.0pt;
font-family:Wingdings;}
ol
{margin-bottom:0cm;}
ul
{margin-bottom:0cm;}
--></style></head><body lang=EN-CA link=blue vlink="#954F72"><div class=WordSection1><p class=MsoNormal><span style='font-size:10.0pt;font-family:"Arial",sans-serif;color:#222222;background:white'>Hi,<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:10.0pt;font-family:"Arial",sans-serif;color:#222222;background:white'><o:p> </o:p></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><o:p></o:p></p><div><p class=MsoNormal style='background:white'><span style='font-size:10.0pt;font-family:"Arial",sans-serif;color:#222222'><o:p> </o:p></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. <o:p></o:p></span></p></div><div><p class=MsoNormal style='background:white'><span style='font-size:10.0pt;font-family:"Arial",sans-serif;color:#222222'><o:p> </o:p></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.<o:p></o:p></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. <o:p></o:p></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 r3.in.xyz)</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'><o:p></o:p></span></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:47.25pt;text-indent:-18.0pt;mso-list:l0 level1 lfo1;background:white'><![if !supportLists]><span style='font-size:10.0pt;font-family:Symbol;color:#222222'><span style='mso-list:Ignore'>·<span style='font:7.0pt "Times New Roman"'> </span></span></span><![endif]><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 <o:p></o:p></span></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:47.25pt;text-indent:-18.0pt;mso-list:l0 level1 lfo1;background:white'><![if !supportLists]><span style='font-size:10.0pt;font-family:Symbol;color:#222222'><span style='mso-list:Ignore'>·<span style='font:7.0pt "Times New Roman"'> </span></span></span><![endif]><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? <o:p></o:p></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.”<o:p></o:p></span></p></div><div><p class=MsoNormal style='background:white'><span style='font-size:10.0pt;font-family:"Arial",sans-serif;color:#222222'><o:p> </o:p></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 <o:p></o:p></span></p></div></div></div><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Sent from <a href="https://go.microsoft.com/fwlink/?LinkId=550986">Mail</a> for Windows 10</p><p class=MsoNormal><o:p> </o:p></p><div style='mso-element:para-border-div;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">Vaclav Petras</a><br><b>Sent: </b>June 28, 2018 8:27 PM<br><b>To: </b><a href="mailto:fra.chartier@gmail.com">Francois Chartier</a><br><b>Cc: </b><a href="mailto:mlennert@club.worldonline.be">Moritz Lennert</a>; <a href="mailto:grass-user@lists.osgeo.org">GRASS user list</a><br><b>Subject: </b>Re: [GRASS-user] r3.in.xyz</p></div><p class=MsoNormal><o:p> </o:p></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><o:p> </o:p></p></div><div><p class=MsoNormal><a href="https://grass.osgeo.org/grass74/manuals/r.to.rast3elev.html">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">https://grass.osgeo.org/grass74/manuals/r3.mapcalc.html</a></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal>Best,</p></div><div><p class=MsoNormal>Vaclav</p></div></div><div><p class=MsoNormal><o:p> </o:p></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">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'><o:p> </o:p></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">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='mso-margin-top-alt:0cm;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">mlennert@club.worldonline.be</a> <mailto:<a href="mailto:mlennert@club.worldonline.be" target="_blank">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">r3.in.xyz</a> <<a href="http://r3.in.xyz" target="_blank">http://r3.in.xyz</a>><br>> <<a href="http://r3.in.xyz" 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>> </p></div><p class=MsoNormal style='margin-left:9.6pt'>> I don't think <a href="http://r3.in.xyz" target="_blank">r3.in.xyz</a> <<a href="http://r3.in.xyz" target="_blank">http://r3.in.xyz</a>> is what you need. This</p><div><p class=MsoNormal style='mso-margin-top-alt:0cm;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'><o:p> </o:p></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">grass-user@lists.osgeo.org</a><br><a href="https://lists.osgeo.org/mailman/listinfo/grass-user" target="_blank">https://lists.osgeo.org/mailman/listinfo/grass-user</a></p></blockquote></div></div><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><o:p> </o:p></p></div><div id="DAB4FAD8-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">
<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">www.avast.com</a>
</p>
</td>
</tr>
</table>
<br />
<a href="#DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2" width="1" height="1"> </a></div></body></html>