Re: [gdal-dev] "import gdal" broken in 3.2.0: what to do ?

2020-11-06 Thread Kurt Schwehr
+1 for dropping the old way On Thu, Nov 5, 2020, 7:59 AM Kai Mühlbauer wrote: > Hi, > > I'm using GDAL since 2012 and one of the first things I've read is this: > > https://trac.osgeo.org/gdal/wiki/GdalOgrInPython > > The first mention of the deprecation of this import was added 13 years > ago.

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

[gdal-dev] GDAL(Technical support for Extracting Elevations From Large DTED Terrain Datasets)

2020-11-06 Thread Nabil Adechoubou
Greetings, Our team decided to use GDAL in our code to support certain business requirements. Each time we do an analysis we need the elevation profile lines, our analyses run through 10,000s of such cases across a large area. Which means we need to extract these elevation profile lines 10,000s

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=