<div dir="ltr">Hi Vaclav,<div><br></div><div>Your script looks very nice (much cleaner than mine!). I definitely need to add some better error handling. I can see how using v.to.rast and r.fillnulls may introduce fewer artifacts than v.surf--looking forward to testing it. Based on the amount of helpful advice we've already received, I'm leaning towards grass-dev for discussing the best methods to use and merging of the scripts, but off-list is also ok. </div>
<div><br></div><div>Eric</div><div><div class="gmail_extra"><br><br><div class="gmail_quote">On Wed, Sep 25, 2013 at 11:26 PM, Vaclav Petras <span dir="ltr"><<a href="mailto:wenzeslaus@gmail.com" target="_blank">wenzeslaus@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi Eric,<div><br></div><div>actually I was working on the same module (named r.local.relief), so we have some duplication now. It is also based on [Hesse2010]. I was about to commit it to GRASS Addons but I need to write documentation first. For now, I'm adding the Python script into a attachment.</div>
<div><br></div><div>I've quickly tested both and they give slightly different results. The other visible difference is the contours-to-really_smooth_elevation step. It seems that usage of r.fillnuls is faster than v.surfs. But this might require more testing as well as the different results do (probably caused by different elevation reconstruction). We can talk about it off-list if you prefer or at least on grass-dev list. And of course, we should discuss some merge of our scripts.</div>
<div><br></div><div>As a side note I want also mention [r.shaded.pca] which I haven't announced yet but which is already in addons. It is based on [Devereux2010].</div><div><br></div><div>Best,</div><div>Vaclav</div>
<div>
<br></div><div><span style="font-family:arial,sans-serif;font-size:13px">[Hesse2010] Hesse, Ralf (2010): LiDAR-derived Local Relief Models - a new tool for archaeological prospection. Archaeological Prospection 17:67-72.</span><br>
</div><div><span style="font-family:arial,sans-serif;font-size:13px">[</span>Devereux2008<span style="font-size:13px;font-family:arial,sans-serif">]</span> Devereux, B. J., Amable, G. S., & Crow, P. P. (2008). Visualisation of LiDAR terrain models for archaeological feature detection. Antiquity, 82(316), 470-479.</div>
<div>[r.shaded.pca] <a href="https://trac.osgeo.org/grass/browser/grass-addons/grass7/raster/r.shaded.pca" target="_blank">https://trac.osgeo.org/grass/browser/grass-addons/grass7/raster/r.shaded.pca</a><br></div></div><div class="gmail_extra">
<br><br><div class="gmail_quote"><div><div class="h5">On Wed, Sep 25, 2013 at 10:50 PM, Eric Goddard <span dir="ltr"><<a href="mailto:egoddard1010@gmail.com" target="_blank">egoddard1010@gmail.com</a>></span> wrote:<br>
</div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5">
<div dir="ltr">Hi all,<div><br></div><div>I've created a grass python script that generates a local relief model (LRM) from a high-resolution (such as from LIDAR) DEM based on Hesse, Ralf (2010): LiDAR-derived Local Relief Models - a new tool for archaeological prospection. Archaeological Prospection 17:67-72. </div>
<div><br></div><div>I've tested it on several data sets of mine and it seems to be working fine, however my data is limited to prehistoric Native American mounds. If anyone else has some data they can test the tool on and report on the results it would be much appreciated. The tool can be downloaded from <a href="https://github.com/egoddard/r.surf.lrm" target="_blank">https://github.com/egoddard/r.surf.lrm</a>.</div>
<div><br></div><div>To install in a linux environment, I followed the instructions at <a href="http://grasswiki.osgeo.org/wiki/GRASS_and_Python#Testing_and_installing_Python_extensions" target="_blank">http://grasswiki.osgeo.org/wiki/GRASS_and_Python#Testing_and_installing_Python_extensions</a>. I haven't been able to find solid information on how users who don't compile themselves can install an addon that isn't available via g.extension. There is <a href="http://grasswiki.osgeo.org/wiki/Compile_and_Install#Scripts" target="_blank">http://grasswiki.osgeo.org/wiki/Compile_and_Install#Scripts</a>, but I'm not sure how that applies to Windows users. Any advice in this area would be much appreciated.</div>
<div><br></div><div>As a final note, I adapted this tool based on Rebecca Bennett's bash script that was published in her dissertation. Thank you, Rebecca, for answering all of my questions and for your efforts in trying to test it out :)</div>
<span><font color="#888888">
<div><br></div><div>Eric</div></font></span></div>
<br></div></div>_______________________________________________<br>
grass-user mailing list<br>
<a href="mailto:grass-user@lists.osgeo.org" target="_blank">grass-user@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/mailman/listinfo/grass-user" target="_blank">http://lists.osgeo.org/mailman/listinfo/grass-user</a><br></blockquote></div><br></div>
</blockquote></div><br></div></div></div>