I'm trying to make a coordinate transformation using GDAL/osr python bindings. I have a longitude an a latitude and I want to transform them into the mercator coordinates.
So I define the mercator projection as:
srs = osr.SpatialReference() srs.SetMercator(41.602, 2.403, 1, 0, 0);
and the input point projection as
srs2 = osr.SpatialReference() srs2.SetWellKnownGeogCS("EPSG:4326")
When I transform the point supposed to be 0,0 in the mercator projection, using
transf = CoordinateTransformation(srs2,srs) transf.TransformPoint((2.403,41.602))
The x result is 0, as expected, but the y is a high number: 3799198
What am I doing wrong?