[GRASS-user] r.proj - Input raster map is outside current region

Kirk Wythers kirk.wythers at gmail.com
Fri Dec 16 13:48:15 EST 2011


I am trying to re-project a raster into a location, and I am getting an error about the region of the raster being outside the region of the new location. 

r.proj input=foresttype location=mnwimi_usfs mapset=PERMANENT output=foresttype
Input Projection Parameters:  +proj=aea +lat_1=29.5 +lat_2=45.5 +lat_0=23 +lon_0=-96 +x_0=0 +y_0=0 +no_defs +a=6378137 +rf=298.257222101 +towgs84=0,0,0,0,0,0,0
Input Unit Factor: 1
Output Projection Parameters:  +proj=aea +lat_1=42.122774 +lat_2=49.01518 +lat_0=45.568977 +lon_0=-84.455955 +x_0=1000000 +y_0=1000000 +no_defs +a=6378137 +rf=298.257222101 +towgs84=0.000,0.000,0.000
Output Unit Factor: 1
ERROR: Input raster map is outside current region

However, running g.region on both the raster to be reprojected and the new locations shows the regions to be the same:

GRASS 6.4.1 (greatlakesalbers):~ > g.region -p
projection: 99 (Albers Equal Area)
zone:       0
datum:      nad83
ellipsoid:  grs80
north:      2960848
south:      2087348
west:       -122110
east:       1141890
nsres:      250
ewres:      250
rows:       3494
cols:       5056
cells:      17665664

GRASS 6.4.1 (mnwimi_usfs):~ > g.region -p rast=foresttype
projection: 99 (Albers Equal Area)
zone:       0
datum:      nad83
ellipsoid:  grs80
north:      2960848
south:      2087348
west:       -122110
east:       1141890
nsres:      250
ewres:      250
rows:       3494
cols:       5056
cells:      17665664

I am most perplexed by the differences in the lat_1 lat_2, etc parameters reported by r.proj, and the region values reported by g.region. Any suggestions would be much appreciated.


More information about the grass-user mailing list