[Proj] Re: trouble transforming corine landcover raster projection

Melita Kennedy mkennedy2 at earthlink.net
Mon Mar 26 16:47:11 PDT 2007


Patrick,

For LAEA, I have two more parameters: 

false easting: 4321000.0
false northing: 3210000.0

When I use that definition in ArcMap, the extents you report are in France, at least!

Melita

-- 
Melita Kennedy
Product Specialist
ESRI, Inc.

>----------------------------------------------------------------------
>
>Date: Mon, 26 Mar 2007 08:22:29 +0200
>From: Patrick Giraudoux <patrick.giraudoux at univ-fcomte.fr>
>Subject: [Proj] trouble transforming corine landcover raster
>	projection Lambert laea to lcc
>To: proj at lists.maptools.org
>
>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