my .prj file

[anthony farrant]

Hi again,

I've done some translation effort but in vain. I'm using QGis to try to figure
out the correct extent because the one given by gdalinfo is the PIXEL
resolution:

Upper Left  (    0.0,    0.0)
Lower Left  (    0.0, 4251.0)
Upper Right ( 4901.0,    0.0)
Lower Right ( 4901.0, 4251.0)
Center      ( 2450.5, 2125.5)

Of course when projecting a WGS84 point layer on top of this gives a
totally incorrect map with wildly shifted points.

Here is my translation from below .prj values:

+proj=laea +ellps=GRS80 +a=6378137.0 +lat_0=52.0 +lon_0=20.0
+x_0=5071000.0 +y_0=3210000.0 +k_0=0.0174532925199433 +units=m
+no_defs

But with this projection QGis retains the given extent (above), which
is of course wrong because it is the pixel extent while my point layer
has a geographical extent.

Here is the dataset:
http://dataservice.eea.europa.eu/dataservice/metadetails.asp?id=650

I'm really at a loss here and info on the WEB is inconclusive.

Any help appreciated,

Anthony


On 5/22/07, anthony farrant <anthony.k.farrant at gmail.com> wrote:
> Hi people,
>
> I have a tiff with the following weird projection. If anybody could help me out
> to decode this into MapServer-EPSG format, I'd be very grateful. Any link to
> resources are welcome too. I'm a noobie. Thanks!
>
> PROJCS["User_Defined_Lambert_Azimuthal_Equal_Area",
> GEOGCS["GCS_User_Defined",
> DATUM["D_User_Defined",
> SPHEROID["User_Defined_Spheroid",6378137.0,0.0]],
> PRIMEM["Greenwich",0.0],
> UNIT["Degree",0.0174532925199433]],
> PROJECTION["Lambert_Azimuthal_Equal_Area"],
> PARAMETER["False_Easting",5071000.0],
> PARAMETER["False_Northing",3210000.0],
> PARAMETER["Central_Meridian",20.0],
> PARAMETER["Latitude_Of_Origin",52.0],UNIT["Meter",1.0]]
>
> -- A
>