[GRASS-user] g.region Options Yield Different Values

Glynn Clements glynn at gclements.plus.com
Tue Nov 24 02:19:18 EST 2009


Rich Shepard wrote:

> > Yes. If you create a projected location, you need to enter projected
> > coordinates.
> 
> Glynn,
> 
>    And if they're available only as lat/long geographic coordinates ...?

Create the location anyhow and set the region from a map once you've
imported one.

If for some reason you need suitable bounds before you've imported any
maps you can use use "m.proj -i ..." to convert lat-lon coordinates to
projected coordinates (but projecting the corners of a lat-lon
rectangle may be an underestimate; e.g. for LCC, the southernmost
point will typically be somewhere in the middle of the southern edge,
not one of the corners).

For a more accurate estimate, create a lat-lon location, set the
region, use v.in.region and v.split to create a sampled boundary,
re-project that to the LCC location with v.proj, and use
"g.region vect=..." to set the default region from the projected
boundary.

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


More information about the grass-user mailing list