[GRASS-user] Understanding Regions and Windows

Glynn Clements glynn at gclements.plus.com
Sat Nov 21 17:12:30 EST 2009


Rich Shepard wrote:

>    Starting grass64 from /usr4/grassbase/Oregon I request a listing of the
> region's boundaries and see this:
> 
> GRASS 6.4.0svn (Oregon):/usr4/grassbase/Oregon > g.region -l
> north-west corner: long: 125:18:06.67543W lat: 41:38:42.416382N
> north-east corner: long: 125:18:06.327366W lat: 41:38:42.431669N
> south-east corner: long: 125:18:06.3167W lat: 41:38:42.29503N
> south-west corner: long: 125:18:06.664764W lat: 41:38:42.279743N
> center longitude:  125:18:06.496065W
> center latitude:   41:38:42.355706N
> rows:              1
> cols:              1
> 
>    I've no idea where these values are to be found because in
> /usr4/grassbase/Oregon/rbs/WIND I see this:
> 
> [rshepard at salmo /usr4/grassbase/Oregon]$ less rbs/WIND 
> proj:       99
> zone:       0
> north:      46.15916
> south:      41.933088
> east:       -116.644278
> west:       -124.719106
> cols:       1
> rows:       1
> e-w resol:  8.074828
> n-s resol:  4.226072
> top:        1
> bottom:     0
> cols3:      8
> rows3:      4
> depths:     1
> e-w resol3: 1.0093535
> n-s resol3: 1.056518
> t-b resol:  1
> 
>   Even this is incorrect from the definition of the region and projection
> using the wxPython UI, because the e-w resolution and the n-s resolution are
> both supposed to be in International Feet and were entered as 32.8286760.

The "proj: 99" indicates a cartographic coordinate system, but the
region boundaries seem more appropriate for lat/lon.

>    I'd like to understand where the incorrect corners are found, why the
> region listing does not match the window in the rbs mapset (which matches
> the DEFAULT_WIND and WIND in the PERMANENT mapset, and how to fix this up.

What are the contents of PERMANENT/DEFAULT_WIND and PERMANENT/WIND?

-- 
Glynn Clements <glynn at gclements.plus.com>


More information about the grass-user mailing list