[GRASS-user] r.proj: fixing a failure

Rich Shepard rshepard at appl-ecosys.com
Mon Dec 21 21:42:54 EST 2009


   Trying to bring in the raster DEM to the project location, but r.proj
fails because of a difference in false eastings:

GRASS 6.4.0svn (beaver_lake):/usr4/grassbase > r.proj input=elevation10m
location=dem10m_northwest output=DEM method=cubic
Input Projection Parameters:  +proj=lcc +lat_1=43 +lat_2=45.5 +lat_0=41.75
+lon_0=-120.5 +x_0=1312336 +y_0=0 +no_defs +a=6378137 +rf=298.257222101
+nadgrids=/usr/local/grass-6.4.0svn/etc/nad/WO
Input Unit Factor: 1
Output Projection Parameters:  +proj=lcc +lat_1=43 +lat_2=45.5 +lat_0=41.75
+lon_0=-120.5 +x_0=399999.9999999999 +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

   The input location's PROJ_INFO file contains:

name: Lambert Conformal Conic
proj: lcc
datum: nad83
ellps: grs80
lat_1: 43
lat_2: 45.5
lat_0: 41.75
lon_0: -120.5
x_0: 1312336
y_0: 0
no_defs: defined
nadgrids: WO

and the destination location's PROJ_INFO file contains:

name: Lambert Conformal Conic
proj: lcc
datum: nad83
ellps: grs80
lat_1: 43
lat_2: 45.5
lat_0: 41.75
lon_0: -120.5
x_0: 399999.9999999999
y_0: 0
no_defs: defined
nadgrids: WO

   Both have meters as PROJ_UNITS.

   The original metadata for the 10m DEM shows the following:

Spatial_Reference_Information:
   Horizontal_Coordinate_System_Definition:
     Planar:
       Map_Projection:
         Map_Projection_Name: Lambert Conformal Conic
 	Lambert_Conformal_Conic:
 	  Standard_Parallel: 43.000000
           Standard_Parallel: 45.500000
           Longitude_of_Central_Meridian: -120.500000
 	  Latitude_of_Projection_Origin: 41.750000
 	  False_Easting: 1312336.000000
 	  False_Northing: 0.000000
       Planar_Coordinate_Information:
         Panar_Coordinate_Encoding_Method: row and column
         Coordinate_Representation:
           Abscissa_Resolution: 32.828670
 	  Ordinate_Resolution: 32.828670
 	  Planar_Distance_Units: User_Defined_Unit
 	Geodetic_Model:
           Horizontal_Datum_Name: North American Datum of 1983
 	  Ellipsoid_Name: Geodetic Reference System 80
           Semi-major_Axis: 6378137.000000
   	  Denominator_of_Flattening_Ratio: 298.257222

   What do you folks recommend I do to reconcile the different false easting
for the DEM so it matches that of all the vector maps?

Rich



More information about the grass-user mailing list