[GRASS-user] Projection Units and Values Negatively Affect r.topidx Results

[Rich Shepard]

On Thu, 18 Feb 2010, Rich Shepard wrote:

> GRASS 6.4.0svn (oregonM):/usr4/grassbase > r.proj in=demNW loc=Oregon
> mapset=PERMANENT out=demNW10 method=bilinear res=10.0 --o
> Input Projection Parameters:  +proj=lcc +lat_1=43 +lat_2=45.5 +lat_0=41.75
> +lon_0=-120.5 +x_0=400000.00000000 +y_0=0 +no_defs +a=6378137
> +rf=298.257222101 +nadgrids=/usr/local/grass-6.4.0svn/etc/nad/WO
> Input Unit Factor: 0.3048
> Output Projection Parameters:  +proj=lcc +lat_1=43 +lat_2=45.5 +lat_0=41.45
> +lon_0=-120.5 +x_0=400000 +y_0=0 +no_defs +a=6378137 +rf=298.257222101
> +nadgrids=/usr/local/grass-6.4.0svn/etc/nad/WO
> Output Unit Factor: 1
> ERROR: Input raster map is outside current region

   More information:

   From the input location:

GRASS 6.4.0svn (Oregon):/usr4/grassbase > g.region rast=demNW -p
projection: 99 (Lambert Conformal Conic)
zone:       0
datum:      nad83
ellipsoid:  grs80
north:      1735231.43724681
south:      1120657.1324816
west:       143667.18968253
east:       952979.22944945
nsres:      32.82281055
ewres:      32.82281055
rows:       18724
cols:       24657
cells:      461677668

   From the destination location:

GRASS 6.4.0svn (oregonM):/usr4/grassbase > g.region -p
projection: 99 (Lambert Conformal Conic)
zone:       0
datum:      nad83
ellipsoid:  grs80
north:      1735231.43724681
south:      1120657.1324816
west:       143667.18968253
east:       952979.22944945
nsres:      32.82281055
ewres:      32.82281055
rows:       18724
cols:       24657
cells:      461677668

   Since they are identical, what could trigger the error message, "Input
raster map is outside current region" other than the unit factors, and
that's what I'm trying to change (from feet to meters).

Rich