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["EPSG","6326"]],
PRIMEM["Greenwich",0,
AUTHORITY["EPSG","8901"]],
UNIT["degree",0.0174532925199433,
AUTHORITY["EPSG","9122"]],
AUTHORITY["EPSG","4326"]],
PROJECTION["Mercator_1SP"],
PARAMETER["central_meridian",0],
PARAMETER["scale_factor",1],
PARAMETER["false_easting",0],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]],
AXIS["Easting",EAST],
AXIS["Northing",NORTH],
EXTENSION["PROJ4","+proj=merc +a=6378137 +b=6378137 +lat_ts=0 +lon_0=0
+x_0=0 +y_0=0 +k=1 +units=m +nadgrids=@null +wktext +no_defs"]]
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", 0.0174532925199433]],
PROJECTION["Mercator"],
PARAMETER["latitude_of_origin", 0],
PARAMETER["central_meridian", 0],
PARAMETER["false_easting", 0],
PARAMETER["false_northing", 0],
UNIT["meter", 1],
AUTHORITY["EPSG", "3857"]]
and my system parses this and generates:
proj4-library string for coordsys (+over suppressed):
+proj=merc +nadgrids=@null +a=6378137 +b=6378137 +lat_0=0 +lon_0=0 +x_0=0 +y_0=0 +no_defs
Will GDAL accept a "0" for the inverse of flattening to represent a sphere?
Explicitly saying that the datum/geographic coordsys is truly WGS84 is probably technically wrong.
_______________________________________________ gdal-dev mailing list [email protected] https://lists.osgeo.org/mailman/listinfo/gdal-dev
