[GRASS-user] r.proj: Regions Do Not Match

Rich Shepard rshepard at appl-ecosys.com
Wed Mar 2 17:09:37 EST 2011


   I am trying to reproject all data to a single projection: Albers Equal
Area (NAD83) and g.region is:

projection: 99 (Albers Equal Area)
zone:       0
datum:      nad83
ellipsoid:  grs80
north:      43
south:      34
west:       -121
east:       -112
nsres:      1
ewres:      1
rows:       9
cols:       9
cells:      81

   The DEMs were originally in lat/lon (NAD83) and imported to that location.
When I try to reproject them I see that the regions do not match. For
example:

GRASS 6.5.svn (Nevada-aea):~/grassdata > r.proj in=dem_area15
location=Nevada-ll mapset=PERMANENT method=cubic_f
Input Projection Parameters:  +proj=longlat +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=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
Output Unit Factor: 1
ERROR: Input raster map is outside current region

   Running r.proj -p produces:

GRASS 6.5.svn (Nevada-aea):~/grassdata > r.proj -p in=dem_area15
location=Nevada-ll mapset=PERMANENT method=cubic_f
Input Projection Parameters:  +proj=longlat +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=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
Output Unit Factor: 1
Input map <dem_area15 at PERMANENT> in location <Nevada-ll>:
Source cols: 3600
Source rows: 3600
Local north: 2064411.8954923
Local south: 1973182.43184852
Local west: -1787211.15607592
Local east: -1679987.035036

   Is this a result of user error? If so, do I re-import these data? If not,
what do I need to do to get all these data into a single location?

Rich


More information about the grass-user mailing list