[GRASS-user] [GRASSLIST:1178] how to use cloud cover with r.sun?
Vishal Mehta
vishalm1975 at gmail.com
Mon Aug 28 05:33:38 EDT 2006
so should i change the Proj_units file to:
unit: degree
units: degrees
meters: 111300 * cos(y())
-vishal
On 8/25/06, Glynn Clements <glynn at gclements.plus.com> wrote:
>
>
> Vishal Mehta wrote:
>
> > > This problem can be solved to an extent by using the value from the
> > > PROJ_UNITS file to convert ground distances to metres (although that
> > > requires that the PROJ_UNITS file actually contains the correct
> > > value).
> >
> > I did'nt fully understand this comment Glynn. Could you please elaborate
> > this idea?
> > FYI,
> > GRASS 6.0.1 (cru3):~/CRU > g.proj -w
> > GEOGCS["wgs84",
> > DATUM["WGS_1984",
> > SPHEROID["wgs84",6378137,298.257223563],
> > TOWGS84[0.000,0.000,0.000]],
> > PRIMEM["Greenwich",0],
> > UNIT["degree",0.0174532925199433]]
> >
> > And my Proj_units file in the PERMANENT mapset is:
>
> > unit: degree
> > units: degrees
> > meters: 1.0
>
> The "meters" value in the PROJ_UNITS file is supposed to be a scale
> factor which converts map units to metres (i.e. 1.0 for metres, 0.3048
> for feet). For degrees, it should be roughly 111300 * cos(latitude).
>
> r.slope.aspect uses that value to convert map units to metres; if it's
> set correctly (and the zfactor= option is set correctly if your
> elevations aren't in metres), your slope map should be correct. There
> might still be problems with r.sun itself, though.
>
> > I'm not sure what the last line of g.proj -w means, what does unit there
> > refer to, as opposed to the proj_units 'meters' field?
>
> 0.0174532925199433 is pi/180, i.e. it represents the size of the
> (angular) units in radians.
>
> --
> Glynn Clements <glynn at gclements.plus.com>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.osgeo.org/pipermail/grass-user/attachments/20060828/fcd9d8c1/attachment.html
More information about the grass-user
mailing list