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

[Glynn Clements]

Rich Shepard wrote:

>    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:

> x_0: 1312336

> and the destination location's PROJ_INFO file contains:

> x_0: 399999.9999999999

Note: 1312336 * 0.3048 = 400000.0128

>    Both have meters as PROJ_UNITS.
> 
>    The original metadata for the 10m DEM shows the following:

>  	  Planar_Distance_Units: User_Defined_Unit

Have your investigated the possibility of a feet/metres mismatch?

I suggest using m.proj in both the source and target locations to
check that both projections are correct. If they aren't, the results
of v.proj and r.proj will be nonsense.

-- 
Glynn Clements <glynn at gclements.plus.com>