<div dir="ltr">Tim Bailey<br>Horizon Based Interpolation module for
GRASS<br>Google Summer of Code Week 1 check in 
<br>Friday, June 21 2013<br><br>Hi everyone, 


<p style="margin-bottom:0in">During this first week of coding I have
been working on input specifications for the horizon descriptions
interpolation module. Several early technical challenges emerged. The
process for importing horzons turns out to be a bit confounding. 
Generating input from horizon databases turns out to be non-trivial
task.  I have done an informal survey of every geologist and soil
scientist that I know concerning how they store horizon datasets.  In
addition I have inspected the storage specifications for several
public soils database formats.  
</p>
<p style="margin-bottom:0in"> 
</p>
<p style="margin-bottom:0in">In parallel to the programming project
I have been building three sandboxes to test the methods developed in
this project. One is a site level wetland delineation.  Another is a
database of geotechnical borehole logs that Bob Moskovitz of the
California Geologic Survey provided.  The third is a random Natural
Resources Conservation Service survey subset. 
</p>

<p style="margin-bottom:0in">My work product next week is going to
focus on the interpolation module.  The abstract workflow for the
module as whole is as follows 
</p>
<p style="margin-bottom:0in"><br>
</p>
<p style="margin-bottom:0in"><b>Process for generation of r3 voxel
grid from a population of xy located 1 dimensional horizon
descriptions</b></p>
<p style="margin-bottom:0in"><br>
</p>
<ol><li><p style="margin-bottom:0in"><b>Import of database horizon
        descriptions</b></p>
        </li><li><p style="margin-bottom:0in"><b>Generation of line vector
        representing the path of the horizon description in xyz</b></p>
        </li><li><p style="margin-bottom:0in"><b>Segmentation of the line vector
        into n number xyz points</b></p>
        </li><li><p style="margin-bottom:0in"><b>query of attribute values to
        points from attribute database</b></p>
        </li><li><p style="margin-bottom:0in"><b>Generation of r3 region </b>
        </p>
        </li><li><p style="margin-bottom:0in"><b>Interpolation of point
        attribute values through geostatistical and/or logical operators
        onto voxel grid locations</b></p>
</li></ol>
<p style="margin-bottom:0in">After extensive consultation with my
programming mentor Ben Ducke, I am now moving on to the interpolation
of vector attributes to voxel cells. One point that Ben articulated
was that the interpolation assignment needs to be capable of
operating on integer data so that it can work on</p>

<p style="margin-bottom:0in"><br></p><p style="margin-bottom:0in">Best wishes on this first day of
summer.  
</p><div class=""><div id=":1on" class="" tabindex="0"><img class="" src="https://mail.google.com/mail/u/0/images/cleardot.gif"><br></div><div id=":1on" class="" tabindex="0">Tim Bailey<br></div></div></div>