[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