On 21/09/16 09:35, Nikos Alexandris wrote:
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.
AFAIR, the -g flag in r.proj already outputs flat format.
grass-user mailing list