<html><body><div style="color:#000; background-color:#fff; font-family:Courier New, courier, monaco, monospace, sans-serif;font-size:10pt"><div><span><div style="font-size: 10pt;"><span>Do you need cumulative viewshed analysis?</span></div><div style="background-color: transparent;"><span>Maybe r.cva is for you:</span></div><div style="font-size: 16px; background-color: transparent;"><span><font size="2">http://www.ucl.ac.uk/~tcrnmar/GIS/r.cva.html</font></span></div></span></div><div><br></div>  <div style="font-family: 'Courier New', courier, monaco, monospace, sans-serif; font-size: 10pt;"> <div style="font-family: 'times new roman', 'new york', times, serif; font-size: 12pt;"> <div dir="ltr"> <font size="2" face="Arial"> <hr size="1">  <b><span style="font-weight:bold;">Απο:</span></b> Brian Sanjeewa Rupasinghe <jinkabs@gmail.com><br> <b><span style="font-weight: bold;">Προς:</span></b> Anna Kratochvílová <kratochanna@gmail.com>;
 grass-user@lists.osgeo.org <br> <b><span style="font-weight: bold;">Στάλθηκε:</span></b> 2:45 μ.μ. Τρίτη, 8 Ιανουαρίου 2013<br> <b><span style="font-weight: bold;">Θέμα:</span></b> Re: [GRASS-user] Using a variable inside GRASS commands with Python code/ grass.mapcalc<br> </font> </div> <br><div id="yiv1327436110">Hi,<div><br></div><div>Many thanks. It worked for me. Then i have a similar problem in GRASS.mapcalc. What i need next is to add viewshed of each observer location cumulatively while looping through the observer locations in order to have final integrated viewshed analysis map. So i used the following command in Python (some code ommitted)</div>
<div><br></div><div>for i in:</div><div>    grass.run_command('r.viewshed', input = rinput, output = 'viewshed', coordinate = [x,y], obs_elev = oelv, tgt_elev = th, memory = 4098, flags = 'b', overwrite = True, quiet = True)<span class="yiv1327436110Apple-tab-span" style="white-space:pre;">     </span></div>
<div><br></div><div>    grass.mapcalc("viewshed_cum = viewshed + viewshed_cum", overwrite = True, quiet = True)</div><div><br></div><div>It seems that expressions like sum = c + sum does not work inside <a rel="nofollow" target="_blank" href="http://grass.mapcalc.is/">grass.mapcalc.Is</a> there any alternative or other way around to get this done? I am using GRASS 7 in Windows 7.</div>
<div><br></div><div>Cheers, Brian<br><br><div class="yiv1327436110gmail_quote">On Tue, Jan 8, 2013 at 1:25 PM, Anna Kratochvílová <span dir="ltr"><<a rel="nofollow" ymailto="mailto:kratochanna@gmail.com" target="_blank" href="mailto:kratochanna@gmail.com">kratochanna@gmail.com</a>></span> wrote:<br>
<blockquote class="yiv1327436110gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">Hi Brian,<br>
<div class="yiv1327436110im"><br>
On Tue, Jan 8, 2013 at 6:51 AM, Brian Sanjeewa Rupasinghe<br>
<<a rel="nofollow" ymailto="mailto:jinkabs@gmail.com" target="_blank" href="mailto:jinkabs@gmail.com">jinkabs@gmail.com</a>> wrote:<br>
> Hi all<br>
><br>
> I am generating viewshed of each observer location by iterating through a<br>
> set of observer coordinate file through<br>
> Python code in Windows 7 . Below is the python command used in grass for<br>
> each iteration. According to this, each time<br>
> loop iterates, output is overwritten with the same name. Now what i need is<br>
> to keep each output of all the observer locations. For this, i need to<br>
> concatanate output name (i.e. 'viewshed' with loop index i which will be<br>
> 'viewshed' + str(i) ). How is that possible within<br>
> this grass.run_command?<br>
><br>
<br>
</div>just set a variable within each loop and use it in the run_command:<br>
for i in ... :<br>
    out =  'viewshed' + str(i)<br>
    grass.run_command('r.viewshed', input = rinput, output = out, ...)<br>
<br>
<br>
Regards,<br>
Anna<br>
<div class="yiv1327436110im"><br>
> grass.run_command('r.viewshed', input = rinput, output = 'viewshed',<br>
> coordinate = [x,y], obs_elev = oelv, tgt_elev = th, memory = 4098, flags =<br>
> 'b', overwrite = True, quiet = True)<br>
><br>
> Cheers, Brian<br>
><br>
</div>> _______________________________________________<br>
> grass-user mailing list<br>
> <a rel="nofollow" ymailto="mailto:grass-user@lists.osgeo.org" target="_blank" href="mailto:grass-user@lists.osgeo.org">grass-user@lists.osgeo.org</a><br>
> <a rel="nofollow" target="_blank" href="http://lists.osgeo.org/mailman/listinfo/grass-user">http://lists.osgeo.org/mailman/listinfo/grass-user</a><br>
><br>
</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></body></html>