Rich Shepard wrote:
The LiDAR source location has x- and y-resolution of 3 feet
(international) and several thousand cells east-west and north-south.
g.region -p for this location displays:
projection: 99 (unnamed)
datum: ** unknown (default: WGS84) **
When in the project location g.region -p produces:
projection: 99 (NAD_1983_HARN_StatePlane_Oregon_North_FIPS_3601_Feet_Intl)
I want to understand why the nsres is now 2561+ feet, the ewres now
feet, and the numbers of rows and columns has decreased so drastically.
Looking forward to learning,
If it is about raster reprojection, there are hints in the r.proj manual
about region settings before applying r.proj.
A simple way to do this is to check the projected bounds of the input map in
the current location's projection using the -p flag. The -g flag reports the
same thing, but in a form which can be directly cut and pasted into a
Some observation: `-g` prints out in shell style, suitable for
`eval`uating. The `eval` utility, for example, constructs a
command, from the given arguments (here, the ouput of `g.region -g`)
which is then executed by the shell. See for example
A ready to copy-n-paste form is obtained via `-f` (this is a recent
addition by Moritz during the code sprint in Bonn).
Maybe the `-f` flag would be a useful addition to `r.proj` as well.
After setting the region in that way you might check the
cell resolution with "g.region -p" then snap it to a regular grid with
g.region's -a flag. E.g. g.region -a res=5 -p. Note that this is just a
A more involved, but more accurate, way to do this is to generate a vector
"box" map of the region in the source location using v.in.region -d. This
"box" map is then reprojected into the target location with v.proj. Next the
region in the target location is set to the extent of the new vector map
with g.region along with the desired raster resolution (g.region -m can be
used in Latitude/Longitude locations to measure the geodetic length of a
pixel). r.proj is then run for the raster map the user wants to reproject.
In this case a little preparation goes a long way.
grass-user mailing list