<div dir="ltr"><div><div><div><div><div><div><div><div>Thanks Daniel,<br><br></div>This really helps. I understand it more, not all ;)<br><br></div>I started creating the slope and aspect files first using:<br></div>r.slope.aspect elevation="w001001@temp" slope="slope" aspect="aspect" format="degrees" prec="float" zfactor=1.0 min_slp_allowed=0.0<br>
</div>But this didn't work. It took 0 seconds and I got empty data files.<br></div>I can view my input file so I assume that it is good.<br></div>Reading the manual doesn't give me any clues what I'm doing wrong.<br>
<br></div>My input file contains several houses so slope and aspect can be calculated.<br><br></div>Any new pointers are much appreciated.<br></div><div class="gmail_extra"><br clear="all"><div><br>
<br>
Paul<br>
<br>
<div style="font-size:13.3px;font-family:Verdana,Arial,Helvetica,sans-serif"><div style="border:1px solid #e5e5e5;background-color:#f8f8f8;padding:0.8em 0.8em 0.3em 0.8em"><div style="margin:0 0 8px 0"><p style="margin:0">
<span><b><img style="float:left;margin-left:5px;margin-right:5px" src="http://www.bontepaarden.nl/bontepaarden/images/newButton.png" alt="" height="120" width="160">Paul Meems </b><br><span>Release manager, configuration manager</span><br>
<span>and forum moderator of MapWindow GIS.</span><br><a href="http://www.mapwindow.org/" target="_blank"><span>www.mapwindow.org</span></a><br><br><span>Owner of MapWindow.nl - Support for</span><br><span>Dutch speaking users.</span><br>
<a href="http://www.mapwindow.nl/" target="_blank"><span>www.mapwindow.nl</span></a><br></span></p>
<p style="margin:0"><b><br></b></p></div></div><img src="http://p1.wisestamp.com/pixel.png?p=mozilla&v=3.9.12&t=1322134056890&u=482737fd0c4ca423" height="1" width="1"></div>
</div>
<br><br><div class="gmail_quote">2012/12/27 Daniel Lee <span dir="ltr"><<a href="mailto:lee@isi-solutions.org" target="_blank">lee@isi-solutions.org</a>></span><br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi Paul,<br><br>Welcome to the list :) Here are a couple of pointers:<br><div class="gmail_extra"><br><br><div class="gmail_quote"><div class="im">2012/12/27 Paul Meems <span dir="ltr"><<a href="mailto:bontepaarden@gmail.com" target="_blank">bontepaarden@gmail.com</a>></span><br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>The 'Input options' tab already raises questions for me. For aspect and slope map in ArcMap FROM_DEM is used so I select my clipped tiff again. My area is located in The Netherlands so I use 90 for aspect. I don't see something in ArcMap that looks like linkein so I leave this blank. The same for albedo. I'm not sure what to use for longin and latin. My tiff is projected (Amersfoort/RD_New). I'm not sure either what to use for coefbh, coefdh, horizon or horizonstep.<br>
</div></div></blockquote></div><div><br>For slope and aspect, you need to use maps of the slope and aspect in your study area. They can both be calculated with r.slope.aspect. Using 90 for aspect because you're in the Netherlands doesn't make sense since the aspect is the orientation of each sloped pixel - a hill in the Netherlands has north, east, south and west sides just like anywhere else.<br>
<br>I've heard that it's possible to just leave these options blank and have r.sun calculate slope and aspect on the fly for you, but I've seen mixed reviews of that and have no personal experience in doing so since my company's programs were written at a time when those options were required and - well, if I can, I don't change a running system unless I've got a good reason to :)<br>
<br>The projection of your tiff itself shouldn't matter because GRASS doesn't use the projection of single maps - you work inside a location and the location has a single location. This forces you to use one single coordinate system for your entire project, removing the temptation to do on-the-fly reprojection, which is a really, really good idea. If your projection doesn't match that of your map, you didn't import it right (you can check the projection using g.proj).<br>
<br>coefbh, coefdh, horizon and horizonstep can be left blank. coefbh and coefdh just help you make more realistic assessments - I'm not sure how realistic you need to be - and horizon and horizonstep are similar. You can read about them in the manual but if you're not sure what to enter, I wouldn't touch it.<br>
</div><div class="im"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>
</div><div>On the 'Output options' tab I'm completely puzzled. I would have expected to set file names but all I can do is select my clipped tif.<br></div></div></blockquote></div><div><br>GRASS doesn't work with external files, it creates maps that it manages in an internal database. You just enter the names of the maps you want. If you need them as files later, you can export them.<br>
</div><div class="im"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div></div><div></div><div>The 'Optional' tab. I want the calculation to use nearby houses, trees and chimneys so I check 's'. 'step' I keep as 0.5, assuming this is the same as ArcMap's 'hour interval'. I don't know what declin or dist should be. For time Mode 1 is mentioned. What are the different modes? <br>
</div></div></blockquote></div><div><br>These are described pretty well in the manual. If you're working with a headless Ubuntu server, just google it. All the GRASS manuals are available online. r.sun's manual is located at <a href="http://grass.osgeo.org/grass64/manuals/r.sun.html" target="_blank">http://grass.osgeo.org/grass64/manuals/r.sun.html</a>.<br>
</div><div class="im"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>
</div><div>For ArcMap I need to pass the Latitude. I can't find that for r.sun.<br></div></div></blockquote></div><div><br>You don't need to in GRASS. ArcMap assumes the same latitude for the entire map (bad). GRASS computes the solar zenith and azimuth dependent on each pixel's latitude (good).<br>
<br>Hope that helps!<span class="HOEnZb"><font color="#888888"><br><br>Daniel<br></font></span></div></div><br></div>
</blockquote></div><br></div>