[Proj] trouble transforming corine landcover raster projection Lambert laea to lcc

Patrick Giraudoux patrick.giraudoux at univ-fcomte.fr
Sun Mar 25 23:22:29 PDT 2007


Dear listers

I am using GRASS which is grounded, to my knowledge, on proj.4 for 
projection transformation. Within a location called corineurope I have 
extracted a subset of the corine landcover 2000 data set 100m resolution 
(whose projection is given below). This subset corresponds to a window 
including the Franche-Comte region (geographical sense). I have created 
another location in a French coordinate system called "Lambert II 
etendu", the region (GRASS sense) including the whole country.

However,  what I get  is just the following sequence:

GRASS 6.0.2 (franceLIIe):~ > r.proj input=fc location=corineurope 
mapset=giraudoux output=fc

Input Projection Parameters: +proj=laea +lat_0=52.0000000000 
+lon_0=10.0000000000 +a=6378137 +rf=298.257222101 +no_defs
Input Unit Factor: 1

Output Projection Parameters: +proj=lcc +lat_0=46.8000000000 
+lat_1=45.8989188889 +lat_2=47.6960144444 +lon_0=2.3372291667 
+x_0=600000.0000000000 +y_0=2200000.0000000000 +a=6378249.145 
+rf=293.465 +no_defs +towgs84=-168.000,-60.000,320.000
Output Unit Factor: 1
ERROR: Input map is outside current region
GRASS 6.0.2 (franceLIIe):~ >

The error message is clear, but after many checks about projection and 
region definitions, I cannot find what is wrong.

Any suggestion appreciated.

Thanks in advance,

Patrick Giraudoux


Input map:


GRASS 6.0.2 (corineurope):~ > g.proj -p
-PROJ_INFO-------------------------------------------------
name       : Lambert Azimuthal Equal Area
proj       : laea
ellps      : grs80
a          : 6378137.0000000000
es         : 0.0066943800
f          : 298.2572221010
lat_0      : 52.0000000000
lon_0      : 10.0000000000
-PROJ_UNITS------------------------------------------------
unit       : meter
units      : meters
meters     : 1.0
GRASS 6.0.2 (corineurope):~ >


GRASS 6.0.2 (corineurope):~ > g.region -p
projection: 99 (Lambert Azimuthal Equal Area)
zone:       0
datum:      ** unknown (default: WGS84) **
ellipsoid:  grs80
north:      2781245
south:      2571341
west:       3952172
east:       4110683
nsres:      100.00190567
ewres:      100.00694006
rows:       2099
cols:       1585
GRASS 6.0.2 (corineurope):~ >


import location:

GRASS 6.0.2 (franceLIIe):~ > g.proj -p
-PROJ_INFO-------------------------------------------------
name       : Lambert Conformal Conic
datum      : ntf
towgs84    : -168.000,-60.000,320.000
proj       : lcc
ellps      : clark80
a          : 6378249.1449999996
es         : 0.0068035113
f          : 293.4650000000
lat_0      : 46.8000000000
lat_1      : 45.8989188889
lat_2      : 47.6960144444
lon_0      : 2.3372291667
x_0        : 600000.0000000000
y_0        : 2200000.0000000000
-PROJ_UNITS------------------------------------------------
unit       : meter
units      : meters
meters     : 1.0

GRASS 6.0.2 (franceLIIe):~ > g.region -p
projection: 99 (Lambert Conformal Conic)
zone:       0
datum:      ntf
ellipsoid:  clark80
north:      2677442
south:      1620430
west:       47679
east:       1197823
nsres:      100.00113529
ewres:      100.00382575
rows:       10570
cols:       11501
GRASS 6.0.2 (franceLIIe):~ >













More information about the Proj mailing list