I am reading
SPS files (see format definition here) from which I gather
Unfortunately, there is no
ESPG reference available nor is there any standard naming scheme for
Ellipsoids so it is easy to get confused and even a well-known
CRS may be better off built manually.
ellipsoidAcronym() method is not capable of recognizing an existing
Ellipsoid if I define my
+rf instead of using the
+ellps. Similarly, the
findMatchingProj() is failing.
In the following example, using
Projection is not recognized (it is actually
EPSG:30732) and I end up creating a useless custom
CRS and cluttering my local database.
def_wo_name = '+proj=utm +zone=32 +a=6378249.145 +rf=293.4663 +towgs84=-209.3622,-87.8162,404.6198,0.0046,3.4784,0.5805,-1.4547 +units=m +no_defs' #def_w_name = '+proj=utm +zone=32 +ellps=clrk80 +towgs84=-209.3622,-87.8162,404.6198,0.0046,3.4784,0.5805,-1.4547 +units=m +no_defs' match = crs.findMatchingProj() crs = QgsCoordinateReferenceSystem.fromProj4(def_wo_name) print(crs.ellipsoidAcronym()) print(match)
The above code output is just the
proj4 string, an empty string (no
Ellipsoid match) and
Is there any method that would allow a better (actual) matching or do I have to check my
+rf against all
Ellipsoids in my database and match manually ?