[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