Re: [gdal-dev] Spherical Mercator Projection mismatching major/minor axis

2020-11-06 Thread Even Rouault
Craig, > The way I handle it is (slightly older notation): > > PROJCS["WGS84 / Spherical Mercator", > GEOGCS["WGS84basedSpheric_GCS", > DATUM["WGS84basedSpheric_Datum", > SPHEROID["WGS84based_Sphere", 6378137, 0]], > PRIMEM["Greenwich", 0], > UNIT["degree",

Re: [gdal-dev] Spherical Mercator Projection mismatching major/minor axis

2020-11-06 Thread Craig Bruce
On 11/6/20 7:07 AM, jratike80 wrote: PROJCS["WGS 84 / Pseudo-Mercator", GEOGCS["WGS 84", DATUM["WGS_1984", SPHEROID["WGS 84",6378137,298.257223563, AUTHORITY["EPSG","7030"]], TOWGS84[0,0,0,0,0,0,0], AUTHORITY

Re: [gdal-dev] Spherical Mercator Projection mismatching major/minor axis

2020-11-06 Thread Tobias Wendorff
Am 06.11.2020 um 12:07 schrieb jratike80: This does not really belong to my knowledge area but I'll have a try anyway. To add a bit of my knowledge, EPSG:3857 is the king/queen of unconformness (better unconform-mess): The geographic reference is based on geographic coordinates in WGS84 (el

Re: [gdal-dev] Spherical Mercator Projection mismatching major/minor axis

2020-11-06 Thread jratike80
Hi, This does not really belong to my knowledge area but I'll have a try anyway. Check what you have after reading the proj string instead. Here with Python >>> from osgeo import osr >>> spatialRef = osr.SpatialReference() >>> spatialRef.ImportFromProj4("+proj=merc +a=6378137 +b=6378137 +lat_ts=

[gdal-dev] Spherical Mercator Projection mismatching major/minor axis

2020-11-05 Thread Andreas Roth
Hi, I call OSRImportFromProj4 with "+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +towgs84=0,0,0,0,0,0,0 +wktext +no_defs" and with the returned handle i get the major and minor axis ( OSRGetSemiMajor/OSRGetSemiMinor). I would expect the b