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 <[EMAIL PROTECTED]> 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
