<html><body><div style="color:#000; background-color:#fff; font-family:verdana, helvetica, sans-serif;font-size:10pt"><div><span>Hi Adrian,</span></div><div style="color: rgb(0, 0, 0); font-size: 13.3333px; font-family: verdana,helvetica,sans-serif; background-color: transparent; font-style: normal;"><span><br></span></div><div style="color: rgb(0, 0, 0); font-size: 13.3333px; font-family: verdana,helvetica,sans-serif; background-color: transparent; font-style: normal;"><span>This is the approach I used to calculate SVF in GRASS (as you probably know). You can test the similarity to Kokalj et al's original algorithm by processing another copy from your DTM using their ENVI plugin or standalone tool from here http://iaps.zrc-sazu.si/en/svf#v and running some visual or quantitative assessments <br></span></div><div style="color: rgb(0, 0, 0); font-size: 13.3333px; font-family: verdana,helvetica,sans-serif; background-color: transparent; font-style:
normal;"><br><span></span></div><div style="color: rgb(0, 0, 0); font-size: 13.3333px; font-family: verdana,helvetica,sans-serif; background-color: transparent; font-style: normal;"><span>FWIW I can't remember if I did a quantitative comparison of the horizon images I created against the SFV I made from teh standalone (for me it was sufficient that the process and variables of horizon view creation in GRASS were methodologically the same as the published SVF, and visually they were identical)</span></div><div style="color: rgb(0, 0, 0); font-size: 13.3333px; font-family: verdana,helvetica,sans-serif; background-color: transparent; font-style: normal;"><br><span></span></div><div style="color: rgb(0, 0, 0); font-size: 13.3333px; font-family: verdana,helvetica,sans-serif; background-color: transparent; font-style: normal;"><span>Good luck!</span></div><div style="color: rgb(0, 0, 0); font-size: 13.3333px; font-family: verdana,helvetica,sans-serif;
background-color: transparent; font-style: normal;"><span>Rebecca<br></span></div><div><br><blockquote style="border-left: 2px solid rgb(16, 16, 255); margin-left: 5px; margin-top: 5px; padding-left: 5px;"> <div style="font-family: verdana, helvetica, sans-serif; font-size: 10pt;"> <div style="font-family: times new roman, new york, times, serif; font-size: 12pt;"> <div dir="ltr"> <hr size="1"> <font face="Arial" size="2"> <b><span style="font-weight:bold;">From:</span></b> Adrian Kirilov <kirilov.adrian@gmail.com><br> <b><span style="font-weight: bold;">To:</span></b> Hamish <hamish_b@yahoo.com> <br><b><span style="font-weight: bold;">Cc:</span></b> GRASS user list <grass-user@lists.osgeo.org> <br> <b><span style="font-weight: bold;">Sent:</span></b> Sunday, 19 May 2013, 13:57<br> <b><span style="font-weight: bold;">Subject:</span></b> Re: [GRASS-user] Fwd: Airborne LiDAR data strip adjustment<br> </font> </div> <div
class="y_msg_container"><br><div id="yiv4545953661"><div dir="ltr">I applied the same method actually (after a bit thinking and reevaluation of my first idea). I prefer to use v.in.ascii -t -r -z input= fs= z= out= , because I'm interpolating the points with v.surf.bspline (bilinear, to be more precise) afterwards.<div>
<br></div><div style="">I have some more questions now. Sorry, but I don't know if I should put them in a different e-mail and subject...</div><div style=""><br></div><div style="">I used r.horizon module to build horizon view for 16 directions, consequently I received as a result 8 output raster maps:</div>
<div style="">r.horizon elevin=demb30_5_combined@FlemishArdennes direction=0 horizonstep=45 maxdistance=10 horizon=SVF_demb30_5 <br></div><div style=""><br></div><div style="">After that I combined them witn r.mapcalc but with simple addition only:<br>
</div><div style=""><br></div><div style="">r.mapcalc:</div><div style="">'SVF_demb30_5_0@FlemishArdennes + SVF_demb30_5_1@FlemishArdennes + SVF_demb30_5_2@FlemishArdennes + SVF_demb30_5_4@FlemishArdennes + SVF_demb30_5_5@FlemishArdennes + SVF_demb30_5_6@FlemishArdennes + SVF_demb30_5_7@FlemishArdennes'</div>
<div style=""><br></div><div style="">Is this approach correct?</div><div style=""><br></div><div style="">And how can I change the range value of the output raster after I have set up first color table with r.colors color=grey. I want to change the range from 0.0-1.0 to something else, like 0.7-1.0 for example to get more contrast or brightness. The last is in fact the question from my first e-mail (see up) with this subject.</div>
<div style=""><br></div><div style="">Thanks a lot for the help!</div><div style=""><br></div><div style=""><br></div><div style="">Greetings,</div><div style=""><br></div><div style="">Adrian Kirilov </div></div><div class="yiv4545953661gmail_extra"><br>
<br><div class="yiv4545953661gmail_quote">On Sat, May 18, 2013 at 5:37 AM, Hamish <span dir="ltr"><<a rel="nofollow" ymailto="mailto:hamish_b@yahoo.com" target="_blank" href="mailto:hamish_b@yahoo.com">hamish_b@yahoo.com</a>></span> wrote:<br><blockquote class="yiv4545953661gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div class="yiv4545953661im">Adrian wrote:<br>
> It's a series of points. But how can I export them from all<br>
> imported files (strips) later into a single one ASCII file into<br>
> specific region? I think that v.out.ascii for example is only<br>
> for exporting single file. And I need to export multiple,<br>
> falling into the region. Is this possible in GRASS?<br>
<br>
</div>first set the region and grid size with g.region, check it is<br>
ok with the -p flag.<br>
<br>
then depending on what you want to do:<br>
<br>
cat *.asc | r.in.xyz x= y= z= fs= input=- output=mapname method=<br>
<br>
or<br>
<br>
cat *.asc | v.in.ascii -r -bt x= y= z= fs= output=mapname<br>
<br>
<br>
r.in.xyz discards points outside of the current region, as does<br>
v.in.ascii if you use the -r flag.<br>
<br>
<br>
<br>
if your data is already processed, sorted, and cleaned I'd<br>
go straight for r.in.xyz myself, but that's just me.<br>
<span class="yiv4545953661HOEnZb"><font color="#888888"><br>
<br>
Hamish<br>
</font></span></blockquote></div><br></div></div><br>_______________________________________________<br>grass-user mailing list<br><a ymailto="mailto:grass-user@lists.osgeo.org" href="mailto:grass-user@lists.osgeo.org">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><br><br></div> </div> </div> </blockquote></div> </div></body></html>