Re: [GRASS-dev] r.stats area question

2013-11-14 Thread Glynn Clements

Eric Goddard wrote:

 My coordinate system is in U.S. Feet (I know, I know...we're still in
 the dark ages). The numbers do appear to be meters, but I just wanted
 to confirm that the output isn't in feet because of the assuming the
 map's coordinate system is in meters part.

Areas are always in square metres, provided that sufficient
information is available to perform conversion.

For lat-lon locations, the conversion is based upon the ellipsoid
parameters and the cell's latitude.

For projected locations, the area per cell is a constant based upon
the location's declared units (i.e. the meters= value in the ouput of
g.proj -g).

This has been the case since at least the oldest version in the
repository (~14 years ago).

AFAICT, the only situation where you could get output in map units
is for an X-Y location or a projected location which lacks a
PROJ_UNITS file.

-- 
Glynn Clements gl...@gclements.plus.com
___
grass-dev mailing list
grass-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-dev


[GRASS-dev] r.stats area question

2013-11-13 Thread Eric Goddard
Hi all,

Quick question regarding the units of area output from r.stats. the -a
flag specifies the output area to be printed in square meters. The
manual, however, says:

Area statistics are given in units of square meters and/or cell counts.

and:

If a single map layer is specified on the command line, a list of
areas in square meters (assuming the map's coordinate system is in
meters) for each category in the raster map layer will be printed...

My coordinate system is in U.S. Feet (I know, I know...we're still in
the dark ages). The numbers do appear to be meters, but I just wanted
to confirm that the output isn't in feet because of the assuming the
map's coordinate system is in meters part.

Thanks,
Eric
___
grass-dev mailing list
grass-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-dev