[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