[GRASSLIST:4926] Re: 100KM raduis viewshed
Dan Jacobson
jidanni at dman.ddts.net
Mon Nov 11 03:39:25 EST 2002
>>>>> "G" == Glynn Clements <glynn.clements at virgin.net> writes:
>> I want to concentrate on just 30 degrees of my horizon. That means
>> I must use mapcalc to make a patt_map mask that any cell with the
>> azimuths I want to be =1. A pie slice shaped mask. Hmmm. How to
>> proceed?
G> Probably the simplest method would be to create a vector map
G> consisting of a single polygon, then convert it to a raster map with
G> v.to.rast.
I wouldn't know where to begin there. Alternatively?
G> Alternatively, you could compute the heading relative to the centre
G> of the map (or any other point) using the (undocumented)
G> two-argument form of r.mapcalc's atan() function, e.g.
G> r.mapcalc "eval(heading = atan(x() - $xc, y() - $yc), ...
G> Note: the argument order is the reverse of ANSI C's atan2() function.
$ xc=235577 yc=2675353 ;\
r.mapcalc "eval(heading = atan(x() - $xc, y() - $yc), \
nantou.triang = heading>140 && heading<170)"
you have confused me
I even tried if(heading...)
--
http://jidanni.org/ Taiwan(04)25854780
More information about the grass-user
mailing list