I have following problem:
I want to convert EPSG 3035 coordinates to GPS Latitude & Longitude coordinates (EPSG 4326) via Python.
Therefore I'm using GDAL and following Python Code with EPSG 3035 coordinates: N2628/E4704
pointX = 2628
pointY = 4704
# Spatial Reference System
inputEPSG = 3035
outputEPSG = 4326
# create a geometry from coordinates
point = ogr.Geometry(ogr.wkbPoint)
point.AddPoint(pointX, pointY)
print pointX, pointY
# create coordinate transformation
inSpatialRef = osr.SpatialReference()
inSpatialRef.ImportFromEPSG(inputEPSG)
outSpatialRef = osr.SpatialReference()
outSpatialRef.ImportFromEPSG(outputEPSG)
coordTransform = osr.CoordinateTransformation(inSpatialRef, outSpatialRef)
# transform point
point.Transform(coordTransform)
# print point in EPSG 4326
print point.GetX(), point.GetY()
unfortunately this returns me a a point somewhere in the South Atlantic, but it should be somewhere in Austria.
The Data Vendor of the EPSG 3035 Coordinates mentions following Reference System of the layer:
+proj=laea +lat_0=52 +lon_0=10 +x_0=4321000 +y_0=3210000 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs
What am I doing wrong?
x_0
andy_0
(false easting, false northing), you at least end up in Europe:coordTransform.TransformPoint(4321000+4704,3210000+2628)
>>> (10.068529554343733, 52.02359886909443, 0.0)