[GRASS-user] trouble transforming corine landcover raster projection Lambert laea to lcc

[Yves Jacolin]

Hi Patrick,

The region for lambert2et seems not correct :
> 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
This not a lambert 2 etendu region. 1 000 000 diffents to the west (I may be 
wrong)

Y.

Le Dimanche 25 Mars 2007 17:17, Patrick Giraudoux a écrit :
> Hi,
>
> I  wonder if this list is the right place to ask for help on such a
> subject since the trouble met can well come from projection definition
> error rather than GRASS itself. Anyway, I have a trouble trying to
> transform projection from one system to another. 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):~ >
>
>
>
>
>
>
>
>
>
>
> _______________________________________________
> grassuser mailing list
> grassuser at grass.itc.it
> http://grass.itc.it/mailman/listinfo/grassuser

-- 
Yves Jacolin
-------------
http://yjacolin.gloobe.org
http://softlibre.gloobe.org