qgscoordinatereferencesystem.cpp

Go to the documentation of this file.

/***************************************************************************
                          qgscoordinatereferencesystem.cpp
 
                             -------------------
    begin                : 2007
    copyright            : (C) 2007 by Gary E. Sherman
    email                : [email protected]
***************************************************************************/
 
/***************************************************************************
 *                                                                         *
 *   This program is free software; you can redistribute it and/or modify  *
 *   it under the terms of the GNU General Public License as published by  *
 *   the Free Software Foundation; either version 2 of the License, or     *
 *   (at your option) any later version.                                   *
 *                                                                         *
 ***************************************************************************/
#include "qgscoordinatereferencesystem.h"
 
#include <cmath>
#include <proj.h>
#include <proj_constants.h>
#include <proj_experimental.h>
#include <sqlite3.h>
 
#include "qgis.h"
#include "qgsapplication.h"
#include "qgscoordinatereferencesystem_legacy_p.h"
#include "qgscoordinatereferencesystem_p.h"
#include "qgscoordinatereferencesystemregistry.h"
#include "qgscoordinatereferencesystemutils.h"
#include "qgsdatums.h"
#include "qgsellipsoidutils.h"
#include "qgslocalec.h"
#include "qgslogger.h"
#include "qgsmessagelog.h"
#include "qgsogcutils.h"
#include "qgsogrutils.h"
#include "qgsprojectionfactors.h"
#include "qgsprojoperation.h"
#include "qgsprojutils.h"
#include "qgsreadwritelocker.h"
#include "qgssettings.h"
 
#include <QDir>
#include <QDomElement>
#include <QDomNode>
#include <QFile>
#include <QFileInfo>
#include <QRegularExpression>
#include <QString>
#include <QTextStream>
 
#include "moc_qgscoordinatereferencesystem.cpp"
 
using namespace Qt::StringLiterals;
 
//gdal and ogr includes (needed for == operator)
#include <ogr_srs_api.h>
#include <cpl_error.h>
#include <cpl_conv.h>
#include <cpl_csv.h>
 
CUSTOM_CRS_VALIDATION QgsCoordinateReferenceSystem::sCustomSrsValidation = nullptr;
 
typedef QHash< long, QgsCoordinateReferenceSystem > SrIdCrsCacheHash;
typedef QHash< QString, QgsCoordinateReferenceSystem > StringCrsCacheHash;
 
Q_GLOBAL_STATIC( QReadWriteLock, sSrIdCacheLock )
Q_GLOBAL_STATIC( SrIdCrsCacheHash, sSrIdCache )
bool QgsCoordinateReferenceSystem::sDisableSrIdCache = false;
 
Q_GLOBAL_STATIC( QReadWriteLock, sOgcLock )
Q_GLOBAL_STATIC( StringCrsCacheHash, sOgcCache )
bool QgsCoordinateReferenceSystem::sDisableOgcCache = false;
 
Q_GLOBAL_STATIC( QReadWriteLock, sProj4CacheLock )
Q_GLOBAL_STATIC( StringCrsCacheHash, sProj4Cache )
bool QgsCoordinateReferenceSystem::sDisableProjCache = false;
 
Q_GLOBAL_STATIC( QReadWriteLock, sCRSWktLock )
Q_GLOBAL_STATIC( StringCrsCacheHash, sWktCache )
bool QgsCoordinateReferenceSystem::sDisableWktCache = false;
 
Q_GLOBAL_STATIC( QReadWriteLock, sCRSSrsIdLock )
Q_GLOBAL_STATIC( SrIdCrsCacheHash, sSrsIdCache )
bool QgsCoordinateReferenceSystem::sDisableSrsIdCache = false;
 
Q_GLOBAL_STATIC( QReadWriteLock, sCrsStringLock )
Q_GLOBAL_STATIC( StringCrsCacheHash, sStringCache )
bool QgsCoordinateReferenceSystem::sDisableStringCache = false;
 
QString getFullProjString( PJ *obj )
{
  // see https://lists.osgeo.org/pipermail/proj/2019-May/008565.html, it's not sufficient to just
  // use proj_as_proj_string
  QgsProjUtils::proj_pj_unique_ptr boundCrs( proj_crs_create_bound_crs_to_WGS84( QgsProjContext::get(), obj, nullptr ) );
  if ( boundCrs )
  {
    if ( const char *proj4src = proj_as_proj_string( QgsProjContext::get(), boundCrs.get(), PJ_PROJ_4, nullptr ) )
    {
      return QString( proj4src );
    }
  }
 
  return QString( proj_as_proj_string( QgsProjContext::get(), obj, PJ_PROJ_4, nullptr ) );
}
//--------------------------
 
QgsCoordinateReferenceSystem::QgsCoordinateReferenceSystem()
{
  static QgsCoordinateReferenceSystem nullCrs = QgsCoordinateReferenceSystem( QString() );
 
  d = nullCrs.d;
}
 
QgsCoordinateReferenceSystem::QgsCoordinateReferenceSystem( const QString &definition )
{
  d = new QgsCoordinateReferenceSystemPrivate();
  createFromString( definition );
}
 
QgsCoordinateReferenceSystem::QgsCoordinateReferenceSystem( const long id, CrsType type )
{
  d = new QgsCoordinateReferenceSystemPrivate();
  Q_NOWARN_DEPRECATED_PUSH
  createFromId( id, type );
  Q_NOWARN_DEPRECATED_POP
}
 
QgsCoordinateReferenceSystem::QgsCoordinateReferenceSystem( const QgsCoordinateReferenceSystem &srs ) //NOLINT
  : d( srs.d )
  , mValidationHint( srs.mValidationHint )
  , mNativeFormat( srs.mNativeFormat )
{}
 
QgsCoordinateReferenceSystem &QgsCoordinateReferenceSystem::operator=( const QgsCoordinateReferenceSystem &srs ) //NOLINT
{
  if ( &srs == this )
    return *this;
 
  d = srs.d;
  mValidationHint = srs.mValidationHint;
  mNativeFormat = srs.mNativeFormat;
  return *this;
}
 
QList<long> QgsCoordinateReferenceSystem::validSrsIds()
{
  QList<long> results;
  // check both standard & user defined projection databases
  QStringList dbs = QStringList() << QgsApplication::srsDatabaseFilePath() << QgsApplication::qgisUserDatabaseFilePath();
 
  const auto constDbs = dbs;
  for ( const QString &db : constDbs )
  {
    QFileInfo myInfo( db );
    if ( !myInfo.exists() )
    {
      QgsDebugError( "failed : " + db + " does not exist!" );
      continue;
    }
 
    sqlite3_database_unique_ptr database;
    sqlite3_statement_unique_ptr statement;
 
    //check the db is available
    int result = openDatabase( db, database );
    if ( result != SQLITE_OK )
    {
      QgsDebugError( "failed : " + db + " could not be opened!" );
      continue;
    }
 
    QString sql = u"select srs_id from tbl_srs"_s;
    int rc;
    statement = database.prepare( sql, rc );
    while ( true )
    {
      // this one is an infinitive loop, intended to fetch any row
      int ret = statement.step();
 
      if ( ret == SQLITE_DONE )
      {
        // there are no more rows to fetch - we can stop looping
        break;
      }
 
      if ( ret == SQLITE_ROW )
      {
        results.append( statement.columnAsInt64( 0 ) );
      }
      else
      {
        QgsMessageLog::logMessage( QObject::tr( "SQLite error: %2\nSQL: %1" ).arg( sql, sqlite3_errmsg( database.get() ) ), QObject::tr( "SpatiaLite" ) );
        break;
      }
    }
  }
  std::sort( results.begin(), results.end() );
  return results;
}
 
QgsCoordinateReferenceSystem QgsCoordinateReferenceSystem::fromOgcWmsCrs( const QString &ogcCrs )
{
  QgsCoordinateReferenceSystem crs;
  crs.createFromOgcWmsCrs( ogcCrs );
  return crs;
}
 
QgsCoordinateReferenceSystem QgsCoordinateReferenceSystem::fromEpsgId( long epsg )
{
  QgsCoordinateReferenceSystem res = fromOgcWmsCrs( "EPSG:" + QString::number( epsg ) );
  if ( res.isValid() )
    return res;
 
  // pre proj6 builds allowed use of ESRI: codes here (e.g. 54030), so we need to keep compatibility
  res = fromOgcWmsCrs( "ESRI:" + QString::number( epsg ) );
  if ( res.isValid() )
    return res;
 
  return QgsCoordinateReferenceSystem();
}
 
QgsCoordinateReferenceSystem QgsCoordinateReferenceSystem::fromProj4( const QString &proj4 )
{
  return fromProj( proj4 );
}
 
QgsCoordinateReferenceSystem QgsCoordinateReferenceSystem::fromProj( const QString &proj )
{
  QgsCoordinateReferenceSystem crs;
  crs.createFromProj( proj );
  return crs;
}
 
QgsCoordinateReferenceSystem QgsCoordinateReferenceSystem::fromWkt( const QString &wkt )
{
  QgsCoordinateReferenceSystem crs;
  crs.createFromWkt( wkt );
  return crs;
}
 
QgsCoordinateReferenceSystem QgsCoordinateReferenceSystem::fromSrsId( long srsId )
{
  QgsCoordinateReferenceSystem crs;
  crs.createFromSrsId( srsId );
  return crs;
}
 
QgsCoordinateReferenceSystem QgsCoordinateReferenceSystem::createCompoundCrs( const QgsCoordinateReferenceSystem &horizontalCrs, const QgsCoordinateReferenceSystem &verticalCrs, QString &error )
{
  error.clear();
  PJ *horizontalObj = horizontalCrs.projObject();
  PJ *verticalObj = verticalCrs.projObject();
  if ( horizontalObj && verticalObj )
  {
    QStringList errors;
    QgsProjUtils::proj_pj_unique_ptr compoundCrs = QgsProjUtils::createCompoundCrs( horizontalObj, verticalObj, &errors );
    if ( compoundCrs )
      return QgsCoordinateReferenceSystem::fromProjObject( compoundCrs.get() );
 
    QStringList formattedErrorList;
    for ( const QString &rawError : std::as_const( errors ) )
    {
      QString formattedError = rawError;
      formattedError.replace( "proj_create_compound_crs: "_L1, QString() );
      formattedErrorList.append( formattedError );
    }
    error = formattedErrorList.join( '\n' );
  }
  return QgsCoordinateReferenceSystem();
}
 
QgsCoordinateReferenceSystem QgsCoordinateReferenceSystem::createGeocentricCrs( const QString &ellipsoid )
{
  const QgsEllipsoidUtils::EllipsoidParameters ellipsoidParams = QgsEllipsoidUtils::ellipsoidParameters( ellipsoid );
  if ( !ellipsoidParams.valid )
    return QgsCoordinateReferenceSystem();
 
  QgsProjUtils::proj_pj_unique_ptr crs( proj_create_geocentric_crs(
    QgsProjContext::get(),
    /*crs_name*/ nullptr,
    /*datum_name*/ nullptr,
    /*ellps_name*/ nullptr,
    /*semi_major_metre*/ ellipsoidParams.semiMajor,
    /*inv_flattening*/ ellipsoidParams.inverseFlattening,
    /*prime_meridian_name*/ nullptr,
    /*prime_meridian_offset*/ 0,
    /*angular_units*/ nullptr, // "NULL for degrees"
    /*angular_units_conv*/ 0,  // "0 for degrees if angular_units == NULL"
    /*linear_units*/ nullptr,  // "NULL for meter"
    /*linear_units_conv*/ 0    // "0 for Metre if linear_units == NULL"
  ) );
  return QgsCoordinateReferenceSystem::fromProjObject( crs.get() );
}
 
QgsCoordinateReferenceSystem::~QgsCoordinateReferenceSystem() //NOLINT
{}
 
bool QgsCoordinateReferenceSystem::createFromId( const long id, CrsType type )
{
  bool result = false;
  switch ( type )
  {
    case InternalCrsId:
      result = createFromSrsId( id );
      break;
    case PostgisCrsId:
      Q_NOWARN_DEPRECATED_PUSH
      result = createFromSrid( id );
      Q_NOWARN_DEPRECATED_POP
      break;
    case EpsgCrsId:
      result = createFromOgcWmsCrs( u"EPSG:%1"_s.arg( id ) );
      break;
    default:
      //THIS IS BAD...THIS PART OF CODE SHOULD NEVER BE REACHED...
      QgsDebugError( u"Unexpected case reached!"_s );
  };
  return result;
}
 
bool QgsCoordinateReferenceSystem::createFromString( const QString &definition )
{
  if ( definition.isEmpty() )
    return false;
 
  QgsReadWriteLocker locker( *sCrsStringLock(), QgsReadWriteLocker::Read );
  if ( !sDisableStringCache )
  {
    QHash< QString, QgsCoordinateReferenceSystem >::const_iterator crsIt = sStringCache()->constFind( definition );
    if ( crsIt != sStringCache()->constEnd() )
    {
      // found a match in the cache
      *this = crsIt.value();
      return d->mIsValid;
    }
  }
  locker.unlock();
 
  bool result = false;
  const thread_local QRegularExpression reCrsId( u"^(epsg|esri|osgeo|ignf|ogc|nkg|zangi|iau_2015|iau2000|postgis|internal|user)\\:(\\w+)$"_s, QRegularExpression::CaseInsensitiveOption );
  QRegularExpressionMatch match = reCrsId.match( definition );
  if ( match.capturedStart() == 0 )
  {
    QString authName = match.captured( 1 ).toLower();
    if ( authName == "epsg"_L1 )
    {
      result = createFromOgcWmsCrs( definition );
    }
    else if ( authName == "postgis"_L1 )
    {
      const long id = match.captured( 2 ).toLong();
      Q_NOWARN_DEPRECATED_PUSH
      result = createFromSrid( id );
      Q_NOWARN_DEPRECATED_POP
    }
    else if ( authName == "esri"_L1
              || authName == "osgeo"_L1
              || authName == "ignf"_L1
              || authName == "zangi"_L1
              || authName == "iau2000"_L1
              || authName == "ogc"_L1
              || authName == "nkg"_L1
              || authName == "iau_2015"_L1 )
    {
      result = createFromOgcWmsCrs( definition );
    }
    else
    {
      const long id = match.captured( 2 ).toLong();
      Q_NOWARN_DEPRECATED_PUSH
      result = createFromId( id, InternalCrsId );
      Q_NOWARN_DEPRECATED_POP
    }
  }
  else
  {
    const thread_local QRegularExpression reCrsStr( u"^(?:(wkt|proj4|proj)\\:)?(.+)$"_s, QRegularExpression::CaseInsensitiveOption );
    match = reCrsStr.match( definition );
    if ( match.capturedStart() == 0 )
    {
      if ( match.captured( 1 ).startsWith( "proj"_L1, Qt::CaseInsensitive ) )
      {
        result = createFromProj( match.captured( 2 ) );
      }
      else
      {
        result = createFromWkt( match.captured( 2 ) );
      }
    }
  }
 
  locker.changeMode( QgsReadWriteLocker::Write );
  if ( !sDisableStringCache )
    sStringCache()->insert( definition, *this );
  return result;
}
 
bool QgsCoordinateReferenceSystem::createFromUserInput( const QString &definition )
{
  if ( definition.isEmpty() )
    return false;
 
  QString userWkt;
  OGRSpatialReferenceH crs = OSRNewSpatialReference( nullptr );
 
  if ( OSRSetFromUserInput( crs, definition.toLocal8Bit().constData() ) == OGRERR_NONE )
  {
    userWkt = QgsOgrUtils::OGRSpatialReferenceToWkt( crs );
    OSRDestroySpatialReference( crs );
  }
  //QgsDebugMsgLevel( "definition: " + definition + " wkt = " + wkt, 2 );
  return createFromWkt( userWkt );
}
 
void QgsCoordinateReferenceSystem::setupESRIWktFix()
{
  // make sure towgs84 parameter is loaded if gdal >= 1.9
  // this requires setting GDAL_FIX_ESRI_WKT=GEOGCS (see qgis bug #5598 and gdal bug #4673)
  const char *configOld = CPLGetConfigOption( "GDAL_FIX_ESRI_WKT", "" );
  const char *configNew = "GEOGCS";
  // only set if it was not set, to let user change the value if needed
  if ( strcmp( configOld, "" ) == 0 )
  {
    CPLSetConfigOption( "GDAL_FIX_ESRI_WKT", configNew );
    if ( strcmp( configNew, CPLGetConfigOption( "GDAL_FIX_ESRI_WKT", "" ) ) != 0 )
      QgsLogger::warning( u"GDAL_FIX_ESRI_WKT could not be set to %1 : %2"_s.arg( configNew, CPLGetConfigOption( "GDAL_FIX_ESRI_WKT", "" ) ) );
    QgsDebugMsgLevel( u"set GDAL_FIX_ESRI_WKT : %1"_s.arg( configNew ), 4 );
  }
  else
  {
    QgsDebugMsgLevel( u"GDAL_FIX_ESRI_WKT was already set : %1"_s.arg( configNew ), 4 );
  }
}
 
bool QgsCoordinateReferenceSystem::createFromOgcWmsCrs( const QString &crs )
{
  if ( crs.isEmpty() )
    return false;
 
  QgsReadWriteLocker locker( *sOgcLock(), QgsReadWriteLocker::Read );
  if ( !sDisableOgcCache )
  {
    QHash< QString, QgsCoordinateReferenceSystem >::const_iterator crsIt = sOgcCache()->constFind( crs );
    if ( crsIt != sOgcCache()->constEnd() )
    {
      // found a match in the cache
      *this = crsIt.value();
      return d->mIsValid;
    }
  }
  locker.unlock();
 
  QString wmsCrs = crs;
 
  QString authority;
  QString code;
  const QgsOgcCrsUtils::CRSFlavor crsFlavor = QgsOgcCrsUtils::parseCrsName( crs, authority, code );
  const QString authorityLower = authority.toLower();
  if ( crsFlavor == QgsOgcCrsUtils::CRSFlavor::AUTH_CODE && ( authorityLower == "user"_L1 || authorityLower == "custom"_L1 || authorityLower == "qgis"_L1 ) )
  {
    if ( createFromSrsId( code.toInt() ) )
    {
      locker.changeMode( QgsReadWriteLocker::Write );
      if ( !sDisableOgcCache )
        sOgcCache()->insert( crs, *this );
      return d->mIsValid;
    }
  }
  else if ( crsFlavor != QgsOgcCrsUtils::CRSFlavor::UNKNOWN )
  {
    wmsCrs = authority + ':' + code;
  }
 
  // first chance for proj 6 - scan through legacy systems and try to use authid directly
  const QString legacyKey = wmsCrs.toLower();
  for ( auto it = sAuthIdToQgisSrsIdMap.constBegin(); it != sAuthIdToQgisSrsIdMap.constEnd(); ++it )
  {
    if ( it.key().compare( legacyKey, Qt::CaseInsensitive ) == 0 )
    {
      const QStringList parts = it.key().split( ':' );
      const QString auth = parts.at( 0 );
      const QString code = parts.at( 1 );
      if ( loadFromAuthCode( auth, code ) )
      {
        locker.changeMode( QgsReadWriteLocker::Write );
        if ( !sDisableOgcCache )
          sOgcCache()->insert( crs, *this );
        return d->mIsValid;
      }
    }
  }
 
  if ( loadFromDatabase( QgsApplication::srsDatabaseFilePath(), u"lower(auth_name||':'||auth_id)"_s, wmsCrs.toLower() ) )
  {
    locker.changeMode( QgsReadWriteLocker::Write );
    if ( !sDisableOgcCache )
      sOgcCache()->insert( crs, *this );
    return d->mIsValid;
  }
 
  // NAD27
  if ( wmsCrs.compare( "CRS:27"_L1, Qt::CaseInsensitive ) == 0 || wmsCrs.compare( "OGC:CRS27"_L1, Qt::CaseInsensitive ) == 0 )
  {
    // TODO: verify same axis orientation
    return createFromOgcWmsCrs( u"EPSG:4267"_s );
  }
 
  // NAD83
  if ( wmsCrs.compare( "CRS:83"_L1, Qt::CaseInsensitive ) == 0 || wmsCrs.compare( "OGC:CRS83"_L1, Qt::CaseInsensitive ) == 0 )
  {
    // TODO: verify same axis orientation
    return createFromOgcWmsCrs( u"EPSG:4269"_s );
  }
 
  // WGS84
  if ( wmsCrs.compare( "CRS:84"_L1, Qt::CaseInsensitive ) == 0 || wmsCrs.compare( "OGC:CRS84"_L1, Qt::CaseInsensitive ) == 0 )
  {
    if ( loadFromDatabase( QgsApplication::srsDatabaseFilePath(), u"lower(auth_name||':'||auth_id)"_s, u"epsg:4326"_s ) )
    {
      d->mAxisInverted = false;
      d->mAxisInvertedDirty = false;
    }
 
    locker.changeMode( QgsReadWriteLocker::Write );
    if ( !sDisableOgcCache )
      sOgcCache()->insert( crs, *this );
 
    return d->mIsValid;
  }
 
  // Try loading from Proj's db using authority and code
  // While this CRS wasn't found in QGIS' srs db, it may be present in proj's
  if ( !authority.isEmpty() && !code.isEmpty() && loadFromAuthCode( authority, code ) )
  {
    locker.changeMode( QgsReadWriteLocker::Write );
    if ( !sDisableOgcCache )
      sOgcCache()->insert( crs, *this );
    return d->mIsValid;
  }
 
  locker.changeMode( QgsReadWriteLocker::Write );
  if ( !sDisableOgcCache )
    sOgcCache()->insert( crs, QgsCoordinateReferenceSystem() );
  return d->mIsValid;
}
 
// Misc helper functions -----------------------
 
 
void QgsCoordinateReferenceSystem::validate()
{
  if ( d->mIsValid || !sCustomSrsValidation )
    return;
 
  // try to validate using custom validation routines
  if ( sCustomSrsValidation )
    sCustomSrsValidation( *this );
}
 
bool QgsCoordinateReferenceSystem::createFromSrid( const long id )
{
  QgsReadWriteLocker locker( *sSrIdCacheLock(), QgsReadWriteLocker::Read );
  if ( !sDisableSrIdCache )
  {
    QHash< long, QgsCoordinateReferenceSystem >::const_iterator crsIt = sSrIdCache()->constFind( id );
    if ( crsIt != sSrIdCache()->constEnd() )
    {
      // found a match in the cache
      *this = crsIt.value();
      return d->mIsValid;
    }
  }
  locker.unlock();
 
  // first chance for proj 6 - scan through legacy systems and try to use authid directly
  for ( auto it = sAuthIdToQgisSrsIdMap.constBegin(); it != sAuthIdToQgisSrsIdMap.constEnd(); ++it )
  {
    if ( it.value().endsWith( u",%1"_s.arg( id ) ) )
    {
      const QStringList parts = it.key().split( ':' );
      const QString auth = parts.at( 0 );
      const QString code = parts.at( 1 );
      if ( loadFromAuthCode( auth, code ) )
      {
        locker.changeMode( QgsReadWriteLocker::Write );
        if ( !sDisableSrIdCache )
          sSrIdCache()->insert( id, *this );
 
        return d->mIsValid;
      }
    }
  }
 
  bool result = loadFromDatabase( QgsApplication::srsDatabaseFilePath(), u"srid"_s, QString::number( id ) );
 
  locker.changeMode( QgsReadWriteLocker::Write );
  if ( !sDisableSrIdCache )
    sSrIdCache()->insert( id, *this );
 
  return result;
}
 
bool QgsCoordinateReferenceSystem::createFromSrsId( const long id )
{
  QgsReadWriteLocker locker( *sCRSSrsIdLock(), QgsReadWriteLocker::Read );
  if ( !sDisableSrsIdCache )
  {
    QHash< long, QgsCoordinateReferenceSystem >::const_iterator crsIt = sSrsIdCache()->constFind( id );
    if ( crsIt != sSrsIdCache()->constEnd() )
    {
      // found a match in the cache
      *this = crsIt.value();
      return d->mIsValid;
    }
  }
  locker.unlock();
 
  // first chance for proj 6 - scan through legacy systems and try to use authid directly
  for ( auto it = sAuthIdToQgisSrsIdMap.constBegin(); it != sAuthIdToQgisSrsIdMap.constEnd(); ++it )
  {
    if ( it.value().startsWith( QString::number( id ) + ',' ) )
    {
      const QStringList parts = it.key().split( ':' );
      const QString auth = parts.at( 0 );
      const QString code = parts.at( 1 );
      if ( loadFromAuthCode( auth, code ) )
      {
        locker.changeMode( QgsReadWriteLocker::Write );
        if ( !sDisableSrsIdCache )
          sSrsIdCache()->insert( id, *this );
        return d->mIsValid;
      }
    }
  }
 
  bool result = loadFromDatabase( id < Qgis::USER_CRS_START_ID ? QgsApplication::srsDatabaseFilePath() : QgsApplication::qgisUserDatabaseFilePath(), u"srs_id"_s, QString::number( id ) );
 
  locker.changeMode( QgsReadWriteLocker::Write );
  if ( !sDisableSrsIdCache )
    sSrsIdCache()->insert( id, *this );
  return result;
}
 
bool QgsCoordinateReferenceSystem::loadFromDatabase( const QString &db, const QString &expression, const QString &value )
{
  d.detach();
 
  QgsDebugMsgLevel( "load CRS from " + db + " where " + expression + " is " + value, 3 );
  d->mIsValid = false;
  d->mWktPreferred.clear();
 
  QFileInfo myInfo( db );
  if ( !myInfo.exists() )
  {
    QgsDebugError( "failed : " + db + " does not exist!" );
    return d->mIsValid;
  }
 
  sqlite3_database_unique_ptr database;
  sqlite3_statement_unique_ptr statement;
  int myResult;
  //check the db is available
  myResult = openDatabase( db, database );
  if ( myResult != SQLITE_OK )
  {
    return d->mIsValid;
  }
 
  /*
    srs_id INTEGER PRIMARY KEY,
    description text NOT NULL,
    projection_acronym text NOT NULL,
    ellipsoid_acronym NOT NULL,
    parameters text NOT NULL,
    srid integer NOT NULL,
    auth_name varchar NOT NULL,
    auth_id integer NOT NULL,
    is_geo integer NOT NULL);
  */
 
  QString mySql = "select srs_id,description,projection_acronym,"
                  "ellipsoid_acronym,parameters,srid,auth_name||':'||auth_id,is_geo,wkt "
                  "from tbl_srs where "
                  + expression
                  + '='
                  + QgsSqliteUtils::quotedString( value )
                  + " order by deprecated";
  statement = database.prepare( mySql, myResult );
  QString wkt;
  // XXX Need to free memory from the error msg if one is set
  if ( myResult == SQLITE_OK && statement.step() == SQLITE_ROW )
  {
    d->mSrsId = statement.columnAsText( 0 ).toLong();
    d->mDescription = statement.columnAsText( 1 );
    d->mProjectionAcronym = statement.columnAsText( 2 );
    d->mEllipsoidAcronym.clear();
    d->mProj4 = statement.columnAsText( 4 );
    d->mWktPreferred.clear();
    d->mSRID = statement.columnAsText( 5 ).toLong();
    d->mAuthId = statement.columnAsText( 6 );
    d->mIsGeographic = statement.columnAsText( 7 ).toInt() != 0;
    wkt = statement.columnAsText( 8 );
    d->mAxisInvertedDirty = true;
 
    if ( d->mSrsId >= Qgis::USER_CRS_START_ID && ( d->mAuthId.isEmpty() || d->mAuthId == QChar( ':' ) ) )
    {
      d->mAuthId = u"USER:%1"_s.arg( d->mSrsId );
    }
    else if ( !d->mAuthId.startsWith( "USER:"_L1, Qt::CaseInsensitive ) )
    {
      QStringList parts = d->mAuthId.split( ':' );
      QString auth = parts.at( 0 );
      QString code = parts.at( 1 );
 
      {
        QgsProjUtils::proj_pj_unique_ptr crs( proj_create_from_database( QgsProjContext::get(), auth.toLatin1(), code.toLatin1(), PJ_CATEGORY_CRS, false, nullptr ) );
        d->setPj( QgsProjUtils::unboundCrs( crs.get() ) );
      }
 
      d->mIsValid = d->hasPj();
      setMapUnits();
    }
 
    if ( !d->mIsValid )
    {
      if ( !wkt.isEmpty() )
      {
        setWktString( wkt );
        // set WKT string resets the description to that description embedded in the WKT, so manually overwrite this back to the
        // value from the user DB
        d->mDescription = statement.columnAsText( 1 );
      }
      else
        setProjString( d->mProj4 );
    }
  }
  else
  {
    QgsDebugMsgLevel( "failed : " + mySql, 4 );
  }
  return d->mIsValid;
}
 
void QgsCoordinateReferenceSystem::removeFromCacheObjectsBelongingToCurrentThread( PJ_CONTEXT *pj_context )
{
  // Not completely sure about object order destruction after main() has
  // exited. So it is safer to check sDisableCache before using sCacheLock
  // in case sCacheLock would have been destroyed before the current TLS
  // QgsProjContext object that has called us...
 
  if ( !sDisableSrIdCache )
  {
    QgsReadWriteLocker locker( *sSrIdCacheLock(), QgsReadWriteLocker::Write );
    if ( !sDisableSrIdCache )
    {
      for ( auto it = sSrIdCache()->begin(); it != sSrIdCache()->end(); )
      {
        auto &v = it.value();
        if ( v.d->removeObjectsBelongingToCurrentThread( pj_context ) )
          it = sSrIdCache()->erase( it );
        else
          ++it;
      }
    }
  }
  if ( !sDisableOgcCache )
  {
    QgsReadWriteLocker locker( *sOgcLock(), QgsReadWriteLocker::Write );
    if ( !sDisableOgcCache )
    {
      for ( auto it = sOgcCache()->begin(); it != sOgcCache()->end(); )
      {
        auto &v = it.value();
        if ( v.d->removeObjectsBelongingToCurrentThread( pj_context ) )
          it = sOgcCache()->erase( it );
        else
          ++it;
      }
    }
  }
  if ( !sDisableProjCache )
  {
    QgsReadWriteLocker locker( *sProj4CacheLock(), QgsReadWriteLocker::Write );
    if ( !sDisableProjCache )
    {
      for ( auto it = sProj4Cache()->begin(); it != sProj4Cache()->end(); )
      {
        auto &v = it.value();
        if ( v.d->removeObjectsBelongingToCurrentThread( pj_context ) )
          it = sProj4Cache()->erase( it );
        else
          ++it;
      }
    }
  }
  if ( !sDisableWktCache )
  {
    QgsReadWriteLocker locker( *sCRSWktLock(), QgsReadWriteLocker::Write );
    if ( !sDisableWktCache )
    {
      for ( auto it = sWktCache()->begin(); it != sWktCache()->end(); )
      {
        auto &v = it.value();
        if ( v.d->removeObjectsBelongingToCurrentThread( pj_context ) )
          it = sWktCache()->erase( it );
        else
          ++it;
      }
    }
  }
  if ( !sDisableSrsIdCache )
  {
    QgsReadWriteLocker locker( *sCRSSrsIdLock(), QgsReadWriteLocker::Write );
    if ( !sDisableSrsIdCache )
    {
      for ( auto it = sSrsIdCache()->begin(); it != sSrsIdCache()->end(); )
      {
        auto &v = it.value();
        if ( v.d->removeObjectsBelongingToCurrentThread( pj_context ) )
          it = sSrsIdCache()->erase( it );
        else
          ++it;
      }
    }
  }
  if ( !sDisableStringCache )
  {
    QgsReadWriteLocker locker( *sCrsStringLock(), QgsReadWriteLocker::Write );
    if ( !sDisableStringCache )
    {
      for ( auto it = sStringCache()->begin(); it != sStringCache()->end(); )
      {
        auto &v = it.value();
        if ( v.d->removeObjectsBelongingToCurrentThread( pj_context ) )
          it = sStringCache()->erase( it );
        else
          ++it;
      }
    }
  }
}
 
bool QgsCoordinateReferenceSystem::hasAxisInverted() const
{
  if ( d->mAxisInvertedDirty )
  {
    d->mAxisInverted = QgsProjUtils::axisOrderIsSwapped( d->threadLocalProjObject() );
    d->mAxisInvertedDirty = false;
  }
 
  return d->mAxisInverted;
}
 
QList<Qgis::CrsAxisDirection> QgsCoordinateReferenceSystem::axisOrdering() const
{
  if ( type() == Qgis::CrsType::Compound )
    return horizontalCrs().axisOrdering() + verticalCrs().axisOrdering();
 
  const PJ *projObject = d->threadLocalProjObject();
  if ( !projObject )
    return {};
 
  PJ_CONTEXT *context = QgsProjContext::get();
  QgsProjUtils::proj_pj_unique_ptr pjCs( proj_crs_get_coordinate_system( context, projObject ) );
  if ( !pjCs )
    return {};
 
  const thread_local QMap< Qgis::CrsAxisDirection, QString > mapping = {
    { Qgis::CrsAxisDirection::North, u"north"_s },
    { Qgis::CrsAxisDirection::NorthNorthEast, u"northNorthEast"_s },
    { Qgis::CrsAxisDirection::NorthEast, u"northEast"_s },
    { Qgis::CrsAxisDirection::EastNorthEast, u"eastNorthEast"_s },
    { Qgis::CrsAxisDirection::East, u"east"_s },
    { Qgis::CrsAxisDirection::EastSouthEast, u"eastSouthEast"_s },
    { Qgis::CrsAxisDirection::SouthEast, u"southEast"_s },
    { Qgis::CrsAxisDirection::SouthSouthEast, u"southSouthEast"_s },
    { Qgis::CrsAxisDirection::South, u"south"_s },
    { Qgis::CrsAxisDirection::SouthSouthWest, u"southSouthWest"_s },
    { Qgis::CrsAxisDirection::SouthWest, u"southWest"_s },
    { Qgis::CrsAxisDirection::WestSouthWest, u"westSouthWest"_s },
    { Qgis::CrsAxisDirection::West, u"west"_s },
    { Qgis::CrsAxisDirection::WestNorthWest, u"westNorthWest"_s },
    { Qgis::CrsAxisDirection::NorthWest, u"northWest"_s },
    { Qgis::CrsAxisDirection::NorthNorthWest, u"northNorthWest"_s },
    { Qgis::CrsAxisDirection::GeocentricX, u"geocentricX"_s },
    { Qgis::CrsAxisDirection::GeocentricY, u"geocentricY"_s },
    { Qgis::CrsAxisDirection::GeocentricZ, u"geocentricZ"_s },
    { Qgis::CrsAxisDirection::Up, u"up"_s },
    { Qgis::CrsAxisDirection::Down, u"down"_s },
    { Qgis::CrsAxisDirection::Forward, u"forward"_s },
    { Qgis::CrsAxisDirection::Aft, u"aft"_s },
    { Qgis::CrsAxisDirection::Port, u"port"_s },
    { Qgis::CrsAxisDirection::Starboard, u"starboard"_s },
    { Qgis::CrsAxisDirection::Clockwise, u"clockwise"_s },
    { Qgis::CrsAxisDirection::CounterClockwise, u"counterClockwise"_s },
    { Qgis::CrsAxisDirection::ColumnPositive, u"columnPositive"_s },
    { Qgis::CrsAxisDirection::ColumnNegative, u"columnNegative"_s },
    { Qgis::CrsAxisDirection::RowPositive, u"rowPositive"_s },
    { Qgis::CrsAxisDirection::RowNegative, u"rowNegative"_s },
    { Qgis::CrsAxisDirection::DisplayRight, u"displayRight"_s },
    { Qgis::CrsAxisDirection::DisplayLeft, u"displayLeft"_s },
    { Qgis::CrsAxisDirection::DisplayUp, u"displayUp"_s },
    { Qgis::CrsAxisDirection::DisplayDown, u"displayDown"_s },
    { Qgis::CrsAxisDirection::Future, u"future"_s },
    { Qgis::CrsAxisDirection::Past, u"past"_s },
    { Qgis::CrsAxisDirection::Towards, u"towards"_s },
    { Qgis::CrsAxisDirection::AwayFrom, u"awayFrom"_s },
  };
 
  QList< Qgis::CrsAxisDirection > res;
  const int axisCount = proj_cs_get_axis_count( context, pjCs.get() );
  if ( axisCount > 0 )
  {
    res.reserve( axisCount );
 
    for ( int i = 0; i < axisCount; ++i )
    {
      const char *outDirection = nullptr;
      proj_cs_get_axis_info( context, pjCs.get(), i, nullptr, nullptr, &outDirection, nullptr, nullptr, nullptr, nullptr );
      // get first word of direction only
      const thread_local QRegularExpression rx( u"([^\\s]+).*"_s );
      const QRegularExpressionMatch match = rx.match( QString( outDirection ) );
      if ( !match.hasMatch() )
        continue;
 
      const QString direction = match.captured( 1 );
      Qgis::CrsAxisDirection dir = Qgis::CrsAxisDirection::Unspecified;
      for ( auto it = mapping.constBegin(); it != mapping.constEnd(); ++it )
      {
        if ( it.value().compare( direction, Qt::CaseInsensitive ) == 0 )
        {
          dir = it.key();
          break;
        }
      }
 
      res.append( dir );
    }
  }
  return res;
}
 
bool QgsCoordinateReferenceSystem::createFromWkt( const QString &wkt )
{
  return createFromWktInternal( wkt, QString() );
}
 
bool QgsCoordinateReferenceSystem::createFromWktInternal( const QString &wkt, const QString &description )
{
  if ( wkt.isEmpty() )
    return false;
 
  d.detach();
 
  QgsReadWriteLocker locker( *sCRSWktLock(), QgsReadWriteLocker::Read );
  if ( !sDisableWktCache )
  {
    QHash< QString, QgsCoordinateReferenceSystem >::const_iterator crsIt = sWktCache()->constFind( wkt );
    if ( crsIt != sWktCache()->constEnd() )
    {
      // found a match in the cache
      *this = crsIt.value();
 
      if ( !description.isEmpty() && d->mDescription.isEmpty() )
      {
        // now we have a name for a previously unknown CRS! Update the cached CRS accordingly, so that we use the name from now on...
        d->mDescription = description;
        locker.changeMode( QgsReadWriteLocker::Write );
        sWktCache()->insert( wkt, *this );
      }
      return d->mIsValid;
    }
  }
  locker.unlock();
 
  d->mIsValid = false;
  d->mProj4.clear();
  d->mWktPreferred.clear();
  if ( wkt.isEmpty() )
  {
    QgsDebugMsgLevel( u"theWkt is uninitialized, operation failed"_s, 4 );
    return d->mIsValid;
  }
 
  // try to match against user crs
  QgsCoordinateReferenceSystem::RecordMap record = getRecord( "select * from tbl_srs where wkt=" + QgsSqliteUtils::quotedString( wkt ) + " order by deprecated" );
  if ( !record.empty() )
  {
    long srsId = record[u"srs_id"_s].toLong();
    if ( srsId > 0 )
    {
      createFromSrsId( srsId );
    }
  }
  else
  {
    setWktString( wkt );
    if ( !description.isEmpty() )
    {
      d->mDescription = description;
    }
    if ( d->mSrsId == 0 )
    {
      // lastly, try a tolerant match of the created proj object against all user CRSes (allowing differences in parameter order during the comparison)
      long id = matchToUserCrs();
      if ( id >= Qgis::USER_CRS_START_ID )
      {
        createFromSrsId( id );
      }
    }
  }
 
  locker.changeMode( QgsReadWriteLocker::Write );
  if ( !sDisableWktCache )
    sWktCache()->insert( wkt, *this );
 
  return d->mIsValid;
  //setMapunits will be called by createfromproj above
}
 
bool QgsCoordinateReferenceSystem::isValid() const
{
  return d->mIsValid;
}
 
bool QgsCoordinateReferenceSystem::createFromProj4( const QString &proj4String )
{
  return createFromProj( proj4String );
}
 
bool QgsCoordinateReferenceSystem::createFromProj( const QString &projString, const bool identify )
{
  if ( projString.isEmpty() )
    return false;
 
  d.detach();
 
  if ( projString.trimmed().isEmpty() )
  {
    d->mIsValid = false;
    d->mProj4.clear();
    d->mWktPreferred.clear();
    return false;
  }
 
  QgsReadWriteLocker locker( *sProj4CacheLock(), QgsReadWriteLocker::Read );
  if ( !sDisableProjCache )
  {
    QHash< QString, QgsCoordinateReferenceSystem >::const_iterator crsIt = sProj4Cache()->constFind( projString );
    if ( crsIt != sProj4Cache()->constEnd() )
    {
      // found a match in the cache
      *this = crsIt.value();
      return d->mIsValid;
    }
  }
  locker.unlock();
 
  //
  // Examples:
  // +proj=tmerc +lat_0=0 +lon_0=-62 +k=0.999500 +x_0=400000 +y_0=0
  // +ellps=clrk80 +towgs84=-255,-15,71,0,0,0,0 +units=m +no_defs
  //
  // +proj=lcc +lat_1=46.8 +lat_0=46.8 +lon_0=2.337229166666664 +k_0=0.99987742
  // +x_0=600000 +y_0=2200000 +a=6378249.2 +b=6356515.000000472 +units=m +no_defs
  //
  QString myProj4String = projString.trimmed();
  myProj4String.remove( u"+type=crs"_s );
  myProj4String = myProj4String.trimmed();
 
  d->mIsValid = false;
  d->mWktPreferred.clear();
 
  if ( identify )
  {
    // first, try to use proj to do this for us...
    const QString projCrsString = myProj4String + ( myProj4String.contains( u"+type=crs"_s ) ? QString() : u" +type=crs"_s );
    QgsProjUtils::proj_pj_unique_ptr crs( proj_create( QgsProjContext::get(), projCrsString.toLatin1().constData() ) );
    if ( crs )
    {
      QString authName;
      QString authCode;
      if ( QgsProjUtils::identifyCrs( crs.get(), authName, authCode, QgsProjUtils::FlagMatchBoundCrsToUnderlyingSourceCrs ) )
      {
        const QString authid = u"%1:%2"_s.arg( authName, authCode );
        if ( createFromOgcWmsCrs( authid ) )
        {
          locker.changeMode( QgsReadWriteLocker::Write );
          if ( !sDisableProjCache )
            sProj4Cache()->insert( projString, *this );
          return d->mIsValid;
        }
      }
    }
 
    // try a direct match against user crses
    QgsCoordinateReferenceSystem::RecordMap myRecord = getRecord( "select * from tbl_srs where parameters=" + QgsSqliteUtils::quotedString( myProj4String ) + " order by deprecated" );
    long id = 0;
    if ( !myRecord.empty() )
    {
      id = myRecord[u"srs_id"_s].toLong();
      if ( id >= Qgis::USER_CRS_START_ID )
      {
        createFromSrsId( id );
      }
    }
    if ( id < Qgis::USER_CRS_START_ID )
    {
      // no direct matches, so go ahead and create a new proj object based on the proj string alone.
      setProjString( myProj4String );
 
      // lastly, try a tolerant match of the created proj object against all user CRSes (allowing differences in parameter order during the comparison)
      id = matchToUserCrs();
      if ( id >= Qgis::USER_CRS_START_ID )
      {
        createFromSrsId( id );
      }
    }
  }
  else
  {
    setProjString( myProj4String );
  }
 
  locker.changeMode( QgsReadWriteLocker::Write );
  if ( !sDisableProjCache )
    sProj4Cache()->insert( projString, *this );
 
  return d->mIsValid;
}
 
//private method meant for internal use by this class only
QgsCoordinateReferenceSystem::RecordMap QgsCoordinateReferenceSystem::getRecord( const QString &sql )
{
  QString myDatabaseFileName;
  QgsCoordinateReferenceSystem::RecordMap myMap;
  QString myFieldName;
  QString myFieldValue;
  sqlite3_database_unique_ptr database;
  sqlite3_statement_unique_ptr statement;
  int myResult;
 
  // Get the full path name to the sqlite3 spatial reference database.
  myDatabaseFileName = QgsApplication::srsDatabaseFilePath();
  QFileInfo myInfo( myDatabaseFileName );
  if ( !myInfo.exists() )
  {
    QgsDebugError( "failed : " + myDatabaseFileName + " does not exist!" );
    return myMap;
  }
 
  //check the db is available
  myResult = openDatabase( myDatabaseFileName, database );
  if ( myResult != SQLITE_OK )
  {
    return myMap;
  }
 
  statement = database.prepare( sql, myResult );
  // XXX Need to free memory from the error msg if one is set
  if ( myResult == SQLITE_OK && statement.step() == SQLITE_ROW )
  {
    int myColumnCount = statement.columnCount();
    //loop through each column in the record adding its expression name and value to the map
    for ( int myColNo = 0; myColNo < myColumnCount; myColNo++ )
    {
      myFieldName = statement.columnName( myColNo );
      myFieldValue = statement.columnAsText( myColNo );
      myMap[myFieldName] = myFieldValue;
    }
    if ( statement.step() != SQLITE_DONE )
    {
      QgsDebugMsgLevel( u"Multiple records found in srs.db"_s, 4 );
      //be less fussy on proj 6 -- the db has MANY more entries!
    }
  }
  else
  {
    QgsDebugMsgLevel( "failed :  " + sql, 4 );
  }
 
  if ( myMap.empty() )
  {
    myDatabaseFileName = QgsApplication::qgisUserDatabaseFilePath();
    QFileInfo myFileInfo;
    myFileInfo.setFile( myDatabaseFileName );
    if ( !myFileInfo.exists() )
    {
      QgsDebugError( u"user qgis.db not found"_s );
      return myMap;
    }
 
    //check the db is available
    myResult = openDatabase( myDatabaseFileName, database );
    if ( myResult != SQLITE_OK )
    {
      return myMap;
    }
 
    statement = database.prepare( sql, myResult );
    // XXX Need to free memory from the error msg if one is set
    if ( myResult == SQLITE_OK && statement.step() == SQLITE_ROW )
    {
      int myColumnCount = statement.columnCount();
      //loop through each column in the record adding its field name and value to the map
      for ( int myColNo = 0; myColNo < myColumnCount; myColNo++ )
      {
        myFieldName = statement.columnName( myColNo );
        myFieldValue = statement.columnAsText( myColNo );
        myMap[myFieldName] = myFieldValue;
      }
 
      if ( statement.step() != SQLITE_DONE )
      {
        QgsDebugMsgLevel( u"Multiple records found in srs.db"_s, 4 );
        myMap.clear();
      }
    }
    else
    {
      QgsDebugMsgLevel( "failed :  " + sql, 4 );
    }
  }
  return myMap;
}
 
// Accessors -----------------------------------
 
long QgsCoordinateReferenceSystem::srsid() const
{
  return d->mSrsId;
}
 
long QgsCoordinateReferenceSystem::postgisSrid() const
{
  return d->mSRID;
}
 
QString QgsCoordinateReferenceSystem::authid() const
{
  return d->mAuthId;
}
 
QString QgsCoordinateReferenceSystem::description() const
{
  if ( d->mDescription.isNull() )
  {
    return QString();
  }
  else
  {
    return d->mDescription;
  }
}
 
QString QgsCoordinateReferenceSystem::userFriendlyIdentifier( Qgis::CrsIdentifierType type ) const
{
  QString id;
  if ( !authid().isEmpty() )
  {
    if ( type != Qgis::CrsIdentifierType::ShortString && !description().isEmpty() )
      id = u"%1 - %2"_s.arg( authid(), description() );
    else
      id = authid();
  }
  else if ( !description().isEmpty() )
    id = description();
  else if ( type == Qgis::CrsIdentifierType::ShortString )
    id = isValid() ? QObject::tr( "Custom CRS" ) : QObject::tr( "Unknown CRS" );
  else if ( !toWkt( Qgis::CrsWktVariant::Preferred ).isEmpty() )
    id = QObject::tr( "Custom CRS: %1" )
           .arg( type == Qgis::CrsIdentifierType::MediumString ? ( toWkt( Qgis::CrsWktVariant::Preferred ).left( 50 ) + QString( QChar( 0x2026 ) ) ) : toWkt( Qgis::CrsWktVariant::Preferred ) );
  else if ( !toProj().isEmpty() )
    id = QObject::tr( "Custom CRS: %1" ).arg( type == Qgis::CrsIdentifierType::MediumString ? ( toProj().left( 50 ) + QString( QChar( 0x2026 ) ) ) : toProj() );
  if ( !id.isEmpty() && !std::isnan( d->mCoordinateEpoch ) )
    id += u" @ %1"_s.arg( qgsDoubleToString( d->mCoordinateEpoch, 3 ) );
 
  return id;
}
 
QString QgsCoordinateReferenceSystem::projectionAcronym() const
{
  if ( d->mProjectionAcronym.isNull() )
  {
    return QString();
  }
  else
  {
    return d->mProjectionAcronym;
  }
}
 
QString QgsCoordinateReferenceSystem::ellipsoidAcronym() const
{
  if ( d->mEllipsoidAcronym.isNull() )
  {
    if ( PJ *obj = d->threadLocalProjObject() )
    {
      QgsProjUtils::proj_pj_unique_ptr ellipsoid( proj_get_ellipsoid( QgsProjContext::get(), obj ) );
      if ( ellipsoid )
      {
        const QString ellipsoidAuthName( proj_get_id_auth_name( ellipsoid.get(), 0 ) );
        const QString ellipsoidAuthCode( proj_get_id_code( ellipsoid.get(), 0 ) );
        if ( !ellipsoidAuthName.isEmpty() && !ellipsoidAuthCode.isEmpty() )
          d->mEllipsoidAcronym = u"%1:%2"_s.arg( ellipsoidAuthName, ellipsoidAuthCode );
        else
        {
          double semiMajor, semiMinor, invFlattening;
          int semiMinorComputed = 0;
          if ( proj_ellipsoid_get_parameters( QgsProjContext::get(), ellipsoid.get(), &semiMajor, &semiMinor, &semiMinorComputed, &invFlattening ) )
          {
            d->mEllipsoidAcronym = u"PARAMETER:%1:%2"_s.arg( qgsDoubleToString( semiMajor ), qgsDoubleToString( semiMinor ) );
          }
          else
          {
            d->mEllipsoidAcronym.clear();
          }
        }
      }
    }
    return d->mEllipsoidAcronym;
  }
  else
  {
    return d->mEllipsoidAcronym;
  }
}
 
QString QgsCoordinateReferenceSystem::toProj4() const
{
  return toProj();
}
 
QString QgsCoordinateReferenceSystem::toProj() const
{
  if ( !d->mIsValid )
    return QString();
 
  if ( d->mProj4.isEmpty() )
  {
    if ( PJ *obj = d->threadLocalProjObject() )
    {
      d->mProj4 = getFullProjString( obj );
    }
  }
  // Stray spaces at the end?
  return d->mProj4.trimmed();
}
 
std::string QgsCoordinateReferenceSystem::toJsonString( bool multiline, int indentationWidth, const QString &schema ) const
{
  if ( !d->mIsValid )
    return {};
 
  if ( PJ *obj = d->threadLocalProjObject() )
  {
    const QByteArray multiLineOption = u"MULTILINE=%1"_s.arg( multiline ? u"YES"_s : u"NO"_s ).toLocal8Bit();
    const QByteArray indentatationWidthOption = u"INDENTATION_WIDTH=%1"_s.arg( multiline ? QString::number( indentationWidth ) : u"0"_s ).toLocal8Bit();
    const QByteArray schemaOption = u"SCHEMA=%1"_s.arg( schema ).toLocal8Bit();
    const char *const options[] = { multiLineOption.constData(), indentatationWidthOption.constData(), schemaOption.constData(), nullptr };
 
    const char *json = proj_as_projjson( QgsProjContext::get(), obj, options );
    return json ? std::string { json } : std::string {};
  }
  else
  {
    return {};
  }
}
 
Qgis::CrsType QgsCoordinateReferenceSystem::type() const
{
  // NOLINTBEGIN(bugprone-branch-clone)
  switch ( d->mProjType )
  {
    case PJ_TYPE_UNKNOWN:
      return Qgis::CrsType::Unknown;
 
    case PJ_TYPE_ELLIPSOID:
    case PJ_TYPE_PRIME_MERIDIAN:
    case PJ_TYPE_GEODETIC_REFERENCE_FRAME:
    case PJ_TYPE_DYNAMIC_GEODETIC_REFERENCE_FRAME:
    case PJ_TYPE_VERTICAL_REFERENCE_FRAME:
    case PJ_TYPE_DYNAMIC_VERTICAL_REFERENCE_FRAME:
    case PJ_TYPE_DATUM_ENSEMBLE:
    case PJ_TYPE_CONVERSION:
    case PJ_TYPE_TRANSFORMATION:
    case PJ_TYPE_CONCATENATED_OPERATION:
    case PJ_TYPE_OTHER_COORDINATE_OPERATION:
    case PJ_TYPE_TEMPORAL_DATUM:
    case PJ_TYPE_ENGINEERING_DATUM:
    case PJ_TYPE_PARAMETRIC_DATUM:
      return Qgis::CrsType::Other;
 
    case PJ_TYPE_CRS:
    case PJ_TYPE_GEOGRAPHIC_CRS:
      //not possible
      return Qgis::CrsType::Other;
 
    case PJ_TYPE_GEODETIC_CRS:
      return Qgis::CrsType::Geodetic;
    case PJ_TYPE_GEOCENTRIC_CRS:
      return Qgis::CrsType::Geocentric;
    case PJ_TYPE_GEOGRAPHIC_2D_CRS:
      return Qgis::CrsType::Geographic2d;
    case PJ_TYPE_GEOGRAPHIC_3D_CRS:
      return Qgis::CrsType::Geographic3d;
    case PJ_TYPE_VERTICAL_CRS:
      return Qgis::CrsType::Vertical;
    case PJ_TYPE_PROJECTED_CRS:
      return Qgis::CrsType::Projected;
    case PJ_TYPE_COMPOUND_CRS:
      return Qgis::CrsType::Compound;
    case PJ_TYPE_TEMPORAL_CRS:
      return Qgis::CrsType::Temporal;
    case PJ_TYPE_ENGINEERING_CRS:
      return Qgis::CrsType::Engineering;
    case PJ_TYPE_BOUND_CRS:
      return Qgis::CrsType::Bound;
    case PJ_TYPE_OTHER_CRS:
      return Qgis::CrsType::Other;
#if PROJ_VERSION_MAJOR > 9 || ( PROJ_VERSION_MAJOR == 9 && PROJ_VERSION_MINOR >= 2 )
    case PJ_TYPE_DERIVED_PROJECTED_CRS:
      return Qgis::CrsType::DerivedProjected;
    case PJ_TYPE_COORDINATE_METADATA:
      return Qgis::CrsType::Other;
#endif
  }
  return Qgis::CrsType::Unknown;
  // NOLINTEND(bugprone-branch-clone)
}
 
bool QgsCoordinateReferenceSystem::isDeprecated() const
{
  const PJ *pj = projObject();
  if ( !pj )
    return false;
 
  return proj_is_deprecated( pj );
}
 
bool QgsCoordinateReferenceSystem::isGeographic() const
{
  return d->mIsGeographic;
}
 
bool QgsCoordinateReferenceSystem::isDynamic() const
{
  const PJ *pj = projObject();
  if ( !pj )
    return false;
 
  return QgsProjUtils::isDynamic( pj );
}
 
QString QgsCoordinateReferenceSystem::celestialBodyName() const
{
  const PJ *pj = projObject();
  if ( !pj )
    return QString();
 
  PJ_CONTEXT *context = QgsProjContext::get();
 
  return QString( proj_get_celestial_body_name( context, pj ) );
}
 
void QgsCoordinateReferenceSystem::setCoordinateEpoch( double epoch )
{
  if ( d->mCoordinateEpoch == epoch )
    return;
 
  // detaching clears the proj object, so we need to clone the existing one first
  QgsProjUtils::proj_pj_unique_ptr clone( proj_clone( QgsProjContext::get(), projObject() ) );
  d.detach();
  d->mCoordinateEpoch = epoch;
  d->setPj( std::move( clone ) );
}
 
double QgsCoordinateReferenceSystem::coordinateEpoch() const
{
  return d->mCoordinateEpoch;
}
 
QgsDatumEnsemble QgsCoordinateReferenceSystem::datumEnsemble() const
{
  QgsDatumEnsemble res;
  res.mValid = false;
 
  const PJ *pj = projObject();
  if ( !pj )
    return res;
 
  PJ_CONTEXT *context = QgsProjContext::get();
 
  QgsProjUtils::proj_pj_unique_ptr ensemble = QgsProjUtils::crsToDatumEnsemble( pj );
  if ( !ensemble )
    return res;
 
  res.mValid = true;
  res.mName = QString( proj_get_name( ensemble.get() ) );
  res.mAuthority = QString( proj_get_id_auth_name( ensemble.get(), 0 ) );
  res.mCode = QString( proj_get_id_code( ensemble.get(), 0 ) );
  res.mRemarks = QString( proj_get_remarks( ensemble.get() ) );
  res.mScope = QString( proj_get_scope( ensemble.get() ) );
  res.mAccuracy = proj_datum_ensemble_get_accuracy( context, ensemble.get() );
 
  const int memberCount = proj_datum_ensemble_get_member_count( context, ensemble.get() );
  for ( int i = 0; i < memberCount; ++i )
  {
    QgsProjUtils::proj_pj_unique_ptr member( proj_datum_ensemble_get_member( context, ensemble.get(), i ) );
    if ( !member )
      continue;
 
    QgsDatumEnsembleMember details;
    details.mName = QString( proj_get_name( member.get() ) );
    details.mAuthority = QString( proj_get_id_auth_name( member.get(), 0 ) );
    details.mCode = QString( proj_get_id_code( member.get(), 0 ) );
    details.mRemarks = QString( proj_get_remarks( member.get() ) );
    details.mScope = QString( proj_get_scope( member.get() ) );
 
    res.mMembers << details;
  }
  return res;
}
 
QgsProjectionFactors QgsCoordinateReferenceSystem::factors( const QgsPoint &point ) const
{
  QgsProjectionFactors res;
 
  // we have to make a transformation object corresponding to the crs
  QString projString = toProj();
  projString.replace( "+type=crs"_L1, QString() );
 
  QgsProjUtils::proj_pj_unique_ptr transformation( proj_create( QgsProjContext::get(), projString.toUtf8().constData() ) );
  if ( !transformation )
    return res;
 
  PJ_COORD coord = proj_coord( 0, 0, 0, HUGE_VAL );
  coord.uv.u = point.x() * M_PI / 180.0;
  coord.uv.v = point.y() * M_PI / 180.0;
 
  proj_errno_reset( transformation.get() );
  const PJ_FACTORS pjFactors = proj_factors( transformation.get(), coord );
  if ( proj_errno( transformation.get() ) )
  {
    return res;
  }
 
  res.mIsValid = true;
  res.mMeridionalScale = pjFactors.meridional_scale;
  res.mParallelScale = pjFactors.parallel_scale;
  res.mArealScale = pjFactors.areal_scale;
  res.mAngularDistortion = pjFactors.angular_distortion;
  res.mMeridianParallelAngle = pjFactors.meridian_parallel_angle * 180 / M_PI;
  res.mMeridianConvergence = pjFactors.meridian_convergence * 180 / M_PI;
  res.mTissotSemimajor = pjFactors.tissot_semimajor;
  res.mTissotSemiminor = pjFactors.tissot_semiminor;
  res.mDxDlam = pjFactors.dx_dlam;
  res.mDxDphi = pjFactors.dx_dphi;
  res.mDyDlam = pjFactors.dy_dlam;
  res.mDyDphi = pjFactors.dy_dphi;
  return res;
}
 
QgsProjOperation QgsCoordinateReferenceSystem::operation() const
{
  if ( !d->mIsValid )
    return QgsProjOperation();
 
  QgsProjOperation res;
 
  // we have to make a transformation object corresponding to the crs
  QString projString = toProj();
  projString.replace( "+type=crs"_L1, QString() );
  if ( projString.isEmpty() )
    return QgsProjOperation();
 
  QgsProjUtils::proj_pj_unique_ptr transformation( proj_create( QgsProjContext::get(), projString.toUtf8().constData() ) );
  if ( !transformation )
    return res;
 
  PJ_PROJ_INFO info = proj_pj_info( transformation.get() );
 
  if ( info.id )
  {
    return QgsApplication::coordinateReferenceSystemRegistry()->projOperations().value( QString( info.id ) );
  }
 
  return res;
}
 
Qgis::DistanceUnit QgsCoordinateReferenceSystem::mapUnits() const
{
  if ( !d->mIsValid )
    return Qgis::DistanceUnit::Unknown;
 
  return d->mMapUnits;
}
 
QgsRectangle QgsCoordinateReferenceSystem::bounds() const
{
  if ( !d->mIsValid )
    return QgsRectangle();
 
  PJ *obj = d->threadLocalProjObject();
  if ( !obj )
    return QgsRectangle();
 
  double westLon = 0;
  double southLat = 0;
  double eastLon = 0;
  double northLat = 0;
 
  if ( !proj_get_area_of_use( QgsProjContext::get(), obj, &westLon, &southLat, &eastLon, &northLat, nullptr ) )
    return QgsRectangle();
 
 
  // don't use the constructor which normalizes!
  QgsRectangle rect;
  rect.setXMinimum( westLon );
  rect.setYMinimum( southLat );
  rect.setXMaximum( eastLon );
  rect.setYMaximum( northLat );
  return rect;
}
 
QString QgsCoordinateReferenceSystem::toOgcUri() const
{
  const auto parts { authid().split( ':' ) };
  if ( parts.length() == 2 )
  {
    if ( parts[0] == "EPSG"_L1 )
      return u"http://www.opengis.net/def/crs/EPSG/0/%1"_s.arg( parts[1] );
    else if ( parts[0] == "OGC"_L1 )
    {
      return u"http://www.opengis.net/def/crs/OGC/1.3/%1"_s.arg( parts[1] );
    }
    else
    {
      QgsMessageLog::logMessage( u"Error converting published CRS to URI %1: (not OGC or EPSG)"_s.arg( authid() ), u"CRS"_s, Qgis::MessageLevel::Critical );
    }
  }
  else
  {
    QgsMessageLog::logMessage( u"Error converting published CRS to URI: %1"_s.arg( authid() ), u"CRS"_s, Qgis::MessageLevel::Critical );
  }
  return QString();
}
 
QString QgsCoordinateReferenceSystem::toOgcUrn() const
{
  const auto parts { authid().split( ':' ) };
  if ( parts.length() == 2 )
  {
    if ( parts[0] == "EPSG"_L1 )
      return u"urn:ogc:def:crs:EPSG::%1"_s.arg( parts[1] );
    else if ( parts[0] == "OGC"_L1 )
    {
      return u"urn:ogc:def:crs:OGC:1.3:%1"_s.arg( parts[1] );
    }
    else if ( parts[0].startsWith( "IAU"_L1 ) )
    {
      if ( parts[0].contains( "_"_L1 ) )
      {
        const auto subParts = parts[0].split( '_' );
        if ( subParts.length() == 2 )
        {
          return u"urn:ogc:def:crs:%1:%2:%3"_s.arg( subParts[0], subParts[1], parts[1] );
        }
      }
      else
      {
        return u"urn:ogc:def:crs:%1::%2"_s.arg( parts[0], parts[1] );
      }
    }
    else
    {
      QgsMessageLog::logMessage( u"Error converting published CRS to URN %1: (not OGC or EPSG)"_s.arg( authid() ), u"CRS"_s, Qgis::MessageLevel::Critical );
    }
  }
  else
  {
    QgsMessageLog::logMessage( u"Error converting published CRS to URN: %1"_s.arg( authid() ), u"CRS"_s, Qgis::MessageLevel::Critical );
  }
  return QString();
}
 
 
void QgsCoordinateReferenceSystem::updateDefinition()
{
  if ( !d->mIsValid )
    return;
 
  if ( d->mSrsId >= Qgis::USER_CRS_START_ID )
  {
    // user CRS, so update to new definition
    createFromSrsId( d->mSrsId );
  }
  else
  {
    // nothing to do -- only user CRS definitions can be changed
  }
}
 
void QgsCoordinateReferenceSystem::setProjString( const QString &proj4String )
{
  d.detach();
  d->mProj4 = proj4String;
  d->mWktPreferred.clear();
 
  QgsLocaleNumC l;
  QString trimmed = proj4String.trimmed();
 
  trimmed += " +type=crs"_L1;
  PJ_CONTEXT *ctx = QgsProjContext::get();
 
  {
    d->setPj( QgsProjUtils::proj_pj_unique_ptr( proj_create( ctx, trimmed.toLatin1().constData() ) ) );
  }
 
  if ( !d->hasPj() )
  {
#ifdef QGISDEBUG
    const int errNo = proj_context_errno( ctx );
    QgsDebugError( u"proj string rejected: %1"_s.arg( proj_context_errno_string( ctx, errNo ) ) );
#endif
    d->mIsValid = false;
  }
  else
  {
    d->mEllipsoidAcronym.clear();
    d->mIsValid = true;
  }
 
  setMapUnits();
}
 
bool QgsCoordinateReferenceSystem::setWktString( const QString &wkt )
{
  bool res = false;
  d->mIsValid = false;
  d->mWktPreferred.clear();
 
  PROJ_STRING_LIST warnings = nullptr;
  PROJ_STRING_LIST grammarErrors = nullptr;
  {
    d->setPj( QgsProjUtils::proj_pj_unique_ptr( proj_create_from_wkt( QgsProjContext::get(), wkt.toLatin1().constData(), nullptr, &warnings, &grammarErrors ) ) );
  }
 
  res = d->hasPj();
  if ( !res )
  {
    QgsDebugMsgLevel( u"\n---------------------------------------------------------------"_s, 2 );
    QgsDebugMsgLevel( u"This CRS could *** NOT *** be set from the supplied Wkt "_s, 2 );
    QgsDebugMsgLevel( "INPUT: " + wkt, 2 );
    for ( auto iter = warnings; iter && *iter; ++iter )
    {
      QgsDebugMsgLevel( *iter, 2 );
    }
    for ( auto iter = grammarErrors; iter && *iter; ++iter )
    {
      QgsDebugMsgLevel( *iter, 2 );
    }
    QgsDebugMsgLevel( u"---------------------------------------------------------------\n"_s, 2 );
  }
  proj_string_list_destroy( warnings );
  proj_string_list_destroy( grammarErrors );
 
  QgsReadWriteLocker locker( *sProj4CacheLock(), QgsReadWriteLocker::Unlocked );
  if ( !res )
  {
    locker.changeMode( QgsReadWriteLocker::Write );
    if ( !sDisableWktCache )
      sWktCache()->insert( wkt, *this );
    return d->mIsValid;
  }
 
  // try to match to a known authority
  if ( d->hasPj() )
  {
    // try 1 - maybe we can directly grab the auth name and code from the crs already?
    QString authName( proj_get_id_auth_name( d->threadLocalProjObject(), 0 ) );
    QString authCode( proj_get_id_code( d->threadLocalProjObject(), 0 ) );
 
    if ( authName.isEmpty() || authCode.isEmpty() )
    {
      // try 2, use proj's identify method and see if there's a nice candidate we can use
      QgsProjUtils::identifyCrs( d->threadLocalProjObject(), authName, authCode );
    }
 
    if ( !authName.isEmpty() && !authCode.isEmpty() )
    {
      QgsCoordinateReferenceSystem fromAuthCode;
      if ( fromAuthCode.loadFromAuthCode( authName, authCode ) )
      {
        *this = fromAuthCode;
        locker.changeMode( QgsReadWriteLocker::Write );
        if ( !sDisableWktCache )
          sWktCache()->insert( wkt, *this );
        return d->mIsValid;
      }
    }
 
    // Still a valid CRS, just not a known one
    d->mIsValid = true;
    d->mDescription = QString( proj_get_name( d->threadLocalProjObject() ) );
    setMapUnits();
  }
 
  return d->mIsValid;
}
 
void QgsCoordinateReferenceSystem::setMapUnits()
{
  if ( !d->mIsValid )
  {
    d->mMapUnits = Qgis::DistanceUnit::Unknown;
    return;
  }
 
  if ( !d->hasPj() )
  {
    d->mMapUnits = Qgis::DistanceUnit::Unknown;
    return;
  }
 
  PJ_CONTEXT *context = QgsProjContext::get();
  // prefer horizontal CRS units, if present
  QgsProjUtils::proj_pj_unique_ptr crs( QgsProjUtils::crsToHorizontalCrs( d->threadLocalProjObject() ) );
  if ( !crs )
    crs = QgsProjUtils::unboundCrs( d->threadLocalProjObject() );
 
  if ( !crs )
  {
    d->mMapUnits = Qgis::DistanceUnit::Unknown;
    return;
  }
 
  QgsProjUtils::proj_pj_unique_ptr coordinateSystem( proj_crs_get_coordinate_system( context, crs.get() ) );
  if ( !coordinateSystem )
  {
    d->mMapUnits = Qgis::DistanceUnit::Unknown;
    return;
  }
 
  const int axisCount = proj_cs_get_axis_count( context, coordinateSystem.get() );
  if ( axisCount > 0 )
  {
    const char *outUnitName = nullptr;
    // Read only first axis
    proj_cs_get_axis_info( context, coordinateSystem.get(), 0, nullptr, nullptr, nullptr, nullptr, &outUnitName, nullptr, nullptr );
 
    const QString unitName( outUnitName );
 
    // proj unit names are freeform -- they differ from authority to authority :(
    // see https://lists.osgeo.org/pipermail/proj/2019-April/008444.html
    if ( unitName.compare( "degree"_L1, Qt::CaseInsensitive ) == 0
         || unitName.compare( "degree minute second"_L1, Qt::CaseInsensitive ) == 0
         || unitName.compare( "degree minute second hemisphere"_L1, Qt::CaseInsensitive ) == 0
         || unitName.compare( "degree minute"_L1, Qt::CaseInsensitive ) == 0
         || unitName.compare( "degree hemisphere"_L1, Qt::CaseInsensitive ) == 0
         || unitName.compare( "degree minute hemisphere"_L1, Qt::CaseInsensitive ) == 0
         || unitName.compare( "hemisphere degree"_L1, Qt::CaseInsensitive ) == 0
         || unitName.compare( "hemisphere degree minute"_L1, Qt::CaseInsensitive ) == 0
         || unitName.compare( "hemisphere degree minute second"_L1, Qt::CaseInsensitive ) == 0
         || unitName.compare( "degree (supplier to define representation)"_L1, Qt::CaseInsensitive ) == 0 )
      d->mMapUnits = Qgis::DistanceUnit::Degrees;
    else if ( unitName.compare( "metre"_L1, Qt::CaseInsensitive ) == 0 || unitName.compare( 'm'_L1, Qt::CaseInsensitive ) == 0 || unitName.compare( "meter"_L1, Qt::CaseInsensitive ) == 0 )
      d->mMapUnits = Qgis::DistanceUnit::Meters;
    else if ( unitName.compare( "US survey foot"_L1, Qt::CaseInsensitive ) == 0 )
      d->mMapUnits = Qgis::DistanceUnit::FeetUSSurvey;
    else if ( unitName.compare( "foot"_L1, Qt::CaseInsensitive ) == 0 )
      d->mMapUnits = Qgis::DistanceUnit::Feet;
    else if ( unitName.compare( "British yard (Sears 1922)"_L1, Qt::CaseInsensitive ) == 0 )
      d->mMapUnits = Qgis::DistanceUnit::YardsBritishSears1922;
    else if ( unitName.compare( "British yard (Sears 1922 truncated)"_L1, Qt::CaseInsensitive ) == 0 )
      d->mMapUnits = Qgis::DistanceUnit::YardsBritishSears1922Truncated;
    else if ( unitName.compare( "British foot (Sears 1922)"_L1, Qt::CaseInsensitive ) == 0 )
      d->mMapUnits = Qgis::DistanceUnit::FeetBritishSears1922;
    else if ( unitName.compare( "British foot (Sears 1922 truncated)"_L1, Qt::CaseInsensitive ) == 0 )
      d->mMapUnits = Qgis::DistanceUnit::FeetBritishSears1922Truncated;
    else if ( unitName.compare( "British chain (Sears 1922)"_L1, Qt::CaseInsensitive ) == 0 )
      d->mMapUnits = Qgis::DistanceUnit::ChainsBritishSears1922;
    else if ( unitName.compare( "British chain (Sears 1922 truncated)"_L1, Qt::CaseInsensitive ) == 0 )
      d->mMapUnits = Qgis::DistanceUnit::ChainsBritishSears1922Truncated;
    else if ( unitName.compare( "British link (Sears 1922)"_L1, Qt::CaseInsensitive ) == 0 )
      d->mMapUnits = Qgis::DistanceUnit::LinksBritishSears1922;
    else if ( unitName.compare( "British link (Sears 1922 truncated)"_L1, Qt::CaseInsensitive ) == 0 )
      d->mMapUnits = Qgis::DistanceUnit::LinksBritishSears1922Truncated;
    else if ( unitName.compare( "British yard (Benoit 1895 A)"_L1, Qt::CaseInsensitive ) == 0 )
      d->mMapUnits = Qgis::DistanceUnit::YardsBritishBenoit1895A;
    else if ( unitName.compare( "British foot (Benoit 1895 A)"_L1, Qt::CaseInsensitive ) == 0 )
      d->mMapUnits = Qgis::DistanceUnit::FeetBritishBenoit1895A;
    else if ( unitName.compare( "British chain (Benoit 1895 A)"_L1, Qt::CaseInsensitive ) == 0 )
      d->mMapUnits = Qgis::DistanceUnit::ChainsBritishBenoit1895A;
    else if ( unitName.compare( "British link (Benoit 1895 A)"_L1, Qt::CaseInsensitive ) == 0 )
      d->mMapUnits = Qgis::DistanceUnit::LinksBritishBenoit1895A;
    else if ( unitName.compare( "British yard (Benoit 1895 B)"_L1, Qt::CaseInsensitive ) == 0 )
      d->mMapUnits = Qgis::DistanceUnit::YardsBritishBenoit1895B;
    else if ( unitName.compare( "British foot (Benoit 1895 B)"_L1, Qt::CaseInsensitive ) == 0 )
      d->mMapUnits = Qgis::DistanceUnit::FeetBritishBenoit1895B;
    else if ( unitName.compare( "British chain (Benoit 1895 B)"_L1, Qt::CaseInsensitive ) == 0 )
      d->mMapUnits = Qgis::DistanceUnit::ChainsBritishBenoit1895B;
    else if ( unitName.compare( "British link (Benoit 1895 B)"_L1, Qt::CaseInsensitive ) == 0 )
      d->mMapUnits = Qgis::DistanceUnit::LinksBritishBenoit1895B;
    else if ( unitName.compare( "British foot (1865)"_L1, Qt::CaseInsensitive ) == 0 )
      d->mMapUnits = Qgis::DistanceUnit::FeetBritish1865;
    else if ( unitName.compare( "British foot (1936)"_L1, Qt::CaseInsensitive ) == 0 )
      d->mMapUnits = Qgis::DistanceUnit::FeetBritish1936;
    else if ( unitName.compare( "Indian foot"_L1, Qt::CaseInsensitive ) == 0 )
      d->mMapUnits = Qgis::DistanceUnit::FeetIndian;
    else if ( unitName.compare( "Indian foot (1937)"_L1, Qt::CaseInsensitive ) == 0 )
      d->mMapUnits = Qgis::DistanceUnit::FeetIndian1937;
    else if ( unitName.compare( "Indian foot (1962)"_L1, Qt::CaseInsensitive ) == 0 )
      d->mMapUnits = Qgis::DistanceUnit::FeetIndian1962;
    else if ( unitName.compare( "Indian foot (1975)"_L1, Qt::CaseInsensitive ) == 0 )
      d->mMapUnits = Qgis::DistanceUnit::FeetIndian1975;
    else if ( unitName.compare( "Indian yard"_L1, Qt::CaseInsensitive ) == 0 )
      d->mMapUnits = Qgis::DistanceUnit::YardsIndian;
    else if ( unitName.compare( "Indian yard (1937)"_L1, Qt::CaseInsensitive ) == 0 )
      d->mMapUnits = Qgis::DistanceUnit::YardsIndian1937;
    else if ( unitName.compare( "Indian yard (1962)"_L1, Qt::CaseInsensitive ) == 0 )
      d->mMapUnits = Qgis::DistanceUnit::YardsIndian1962;
    else if ( unitName.compare( "Indian yard (1975)"_L1, Qt::CaseInsensitive ) == 0 )
      d->mMapUnits = Qgis::DistanceUnit::YardsIndian1975;
    else if ( unitName.compare( "Gold Coast foot"_L1, Qt::CaseInsensitive ) == 0 )
      d->mMapUnits = Qgis::DistanceUnit::FeetGoldCoast;
    else if ( unitName.compare( "Clarke's foot"_L1, Qt::CaseInsensitive ) == 0 )
      d->mMapUnits = Qgis::DistanceUnit::FeetClarkes;
    else if ( unitName.compare( "Clarke's yard"_L1, Qt::CaseInsensitive ) == 0 )
      d->mMapUnits = Qgis::DistanceUnit::YardsClarkes;
    else if ( unitName.compare( "Clarke's chain"_L1, Qt::CaseInsensitive ) == 0 )
      d->mMapUnits = Qgis::DistanceUnit::ChainsClarkes;
    else if ( unitName.compare( "Clarke's link"_L1, Qt::CaseInsensitive ) == 0 )
      d->mMapUnits = Qgis::DistanceUnit::LinksClarkes;
    else if ( unitName.compare( "kilometre"_L1, Qt::CaseInsensitive ) == 0 ) //#spellok
      d->mMapUnits = Qgis::DistanceUnit::Kilometers;
    else if ( unitName.compare( "centimetre"_L1, Qt::CaseInsensitive ) == 0 ) //#spellok
      d->mMapUnits = Qgis::DistanceUnit::Centimeters;
    else if ( unitName.compare( "millimetre"_L1, Qt::CaseInsensitive ) == 0 ) //#spellok
      d->mMapUnits = Qgis::DistanceUnit::Millimeters;
    else if ( unitName.compare( "Statute mile"_L1, Qt::CaseInsensitive ) == 0 )
      d->mMapUnits = Qgis::DistanceUnit::Miles;
    else if ( unitName.compare( "nautical mile"_L1, Qt::CaseInsensitive ) == 0 )
      d->mMapUnits = Qgis::DistanceUnit::NauticalMiles;
    else if ( unitName.compare( "yard"_L1, Qt::CaseInsensitive ) == 0 )
      d->mMapUnits = Qgis::DistanceUnit::Yards;
    else if ( unitName.compare( "fathom"_L1, Qt::CaseInsensitive ) == 0 )
      d->mMapUnits = Qgis::DistanceUnit::Fathoms;
    else if ( unitName.compare( "US survey chain"_L1, Qt::CaseInsensitive ) == 0 )
      d->mMapUnits = Qgis::DistanceUnit::ChainsUSSurvey;
    else if ( unitName.compare( "chain"_L1, Qt::CaseInsensitive ) == 0 )
      d->mMapUnits = Qgis::DistanceUnit::ChainsInternational;
    else if ( unitName.compare( "link"_L1, Qt::CaseInsensitive ) == 0 )
      d->mMapUnits = Qgis::DistanceUnit::LinksInternational;
    else if ( unitName.compare( "US survey link"_L1, Qt::CaseInsensitive ) == 0 )
      d->mMapUnits = Qgis::DistanceUnit::LinksUSSurvey;
    else if ( unitName.compare( "US survey mile"_L1, Qt::CaseInsensitive ) == 0 )
      d->mMapUnits = Qgis::DistanceUnit::MilesUSSurvey;
    else if ( unitName.compare( "German legal metre"_L1, Qt::CaseInsensitive ) == 0 )
      d->mMapUnits = Qgis::DistanceUnit::MetersGermanLegal;
    // TODO - maybe more values to handle here?
    else
      d->mMapUnits = Qgis::DistanceUnit::Unknown;
    return;
  }
  else
  {
    d->mMapUnits = Qgis::DistanceUnit::Unknown;
    return;
  }
}
 
 
long QgsCoordinateReferenceSystem::findMatchingProj()
{
  if ( d->mEllipsoidAcronym.isNull() || d->mProjectionAcronym.isNull() || !d->mIsValid )
  {
    QgsDebugMsgLevel(
      "QgsCoordinateReferenceSystem::findMatchingProj will only "
      "work if prj acr ellipsoid acr and proj4string are set"
      " and the current projection is valid!",
      4
    );
    return 0;
  }
 
  sqlite3_database_unique_ptr database;
  sqlite3_statement_unique_ptr statement;
  int myResult;
 
  // Set up the query to retrieve the projection information
  // needed to populate the list
  QString mySql = QString(
                    "select srs_id,parameters from tbl_srs where "
                    "projection_acronym=%1 and ellipsoid_acronym=%2 order by deprecated"
  )
                    .arg( QgsSqliteUtils::quotedString( d->mProjectionAcronym ), QgsSqliteUtils::quotedString( d->mEllipsoidAcronym ) );
  // Get the full path name to the sqlite3 spatial reference database.
  QString myDatabaseFileName = QgsApplication::srsDatabaseFilePath();
 
  //check the db is available
  myResult = openDatabase( myDatabaseFileName, database );
  if ( myResult != SQLITE_OK )
  {
    return 0;
  }
 
  statement = database.prepare( mySql, myResult );
  if ( myResult == SQLITE_OK )
  {
    while ( statement.step() == SQLITE_ROW )
    {
      QString mySrsId = statement.columnAsText( 0 );
      QString myProj4String = statement.columnAsText( 1 );
      if ( toProj() == myProj4String.trimmed() )
      {
        return mySrsId.toLong();
      }
    }
  }
 
  //
  // Try the users db now
  //
 
  myDatabaseFileName = QgsApplication::qgisUserDatabaseFilePath();
  //check the db is available
  myResult = openDatabase( myDatabaseFileName, database );
  if ( myResult != SQLITE_OK )
  {
    return 0;
  }
 
  statement = database.prepare( mySql, myResult );
 
  if ( myResult == SQLITE_OK )
  {
    while ( statement.step() == SQLITE_ROW )
    {
      QString mySrsId = statement.columnAsText( 0 );
      QString myProj4String = statement.columnAsText( 1 );
      if ( toProj() == myProj4String.trimmed() )
      {
        return mySrsId.toLong();
      }
    }
  }
 
  return 0;
}
 
bool QgsCoordinateReferenceSystem::operator==( const QgsCoordinateReferenceSystem &srs ) const
{
  // shortcut
  if ( d == srs.d )
    return true;
 
  if ( !d->mIsValid && !srs.d->mIsValid )
    return true;
 
  if ( !d->mIsValid || !srs.d->mIsValid )
    return false;
 
  if ( !qgsNanCompatibleEquals( d->mCoordinateEpoch, srs.d->mCoordinateEpoch ) )
    return false;
 
  const bool isUser = d->mSrsId >= Qgis::USER_CRS_START_ID;
  const bool otherIsUser = srs.d->mSrsId >= Qgis::USER_CRS_START_ID;
  if ( isUser != otherIsUser )
    return false;
 
  // we can't directly compare authid for user crses -- the actual definition of these may have changed
  if ( !isUser && ( !d->mAuthId.isEmpty() || !srs.d->mAuthId.isEmpty() ) )
    return d->mAuthId == srs.d->mAuthId;
 
  return toWkt( Qgis::CrsWktVariant::Preferred ) == srs.toWkt( Qgis::CrsWktVariant::Preferred );
}
 
bool QgsCoordinateReferenceSystem::operator!=( const QgsCoordinateReferenceSystem &srs ) const
{
  return !( *this == srs );
}
 
QString QgsCoordinateReferenceSystem::toWkt( Qgis::CrsWktVariant variant, bool multiline, int indentationWidth ) const
{
  if ( PJ *obj = d->threadLocalProjObject() )
  {
    const bool isDefaultPreferredFormat = variant == Qgis::CrsWktVariant::Preferred && !multiline;
    if ( isDefaultPreferredFormat && !d->mWktPreferred.isEmpty() )
    {
      // can use cached value
      return d->mWktPreferred;
    }
 
    PJ_WKT_TYPE type = PJ_WKT1_GDAL;
    switch ( variant )
    {
      case Qgis::CrsWktVariant::Wkt1Gdal:
        type = PJ_WKT1_GDAL;
        break;
      case Qgis::CrsWktVariant::Wkt1Esri:
        type = PJ_WKT1_ESRI;
        break;
      case Qgis::CrsWktVariant::Wkt2_2015:
        type = PJ_WKT2_2015;
        break;
      case Qgis::CrsWktVariant::Wkt2_2015Simplified:
        type = PJ_WKT2_2015_SIMPLIFIED;
        break;
      case Qgis::CrsWktVariant::Wkt2_2019:
        type = PJ_WKT2_2019;
        break;
      case Qgis::CrsWktVariant::Wkt2_2019Simplified:
        type = PJ_WKT2_2019_SIMPLIFIED;
        break;
    }
 
    const QByteArray multiLineOption = u"MULTILINE=%1"_s.arg( multiline ? u"YES"_s : u"NO"_s ).toLocal8Bit();
    const QByteArray indentatationWidthOption = u"INDENTATION_WIDTH=%1"_s.arg( multiline ? QString::number( indentationWidth ) : u"0"_s ).toLocal8Bit();
    const char *const options[] = { multiLineOption.constData(), indentatationWidthOption.constData(), nullptr };
    QString res = proj_as_wkt( QgsProjContext::get(), obj, type, options );
 
    if ( isDefaultPreferredFormat )
    {
      // cache result for later use
      d->mWktPreferred = res;
    }
 
    return res;
  }
  return QString();
}
 
bool QgsCoordinateReferenceSystem::readXml( const QDomNode &node )
{
  d.detach();
  bool result = true;
  QDomNode srsNode = node.namedItem( u"spatialrefsys"_s );
 
  if ( !srsNode.isNull() )
  {
    bool initialized = false;
 
    bool ok = false;
    long srsid = srsNode.namedItem( u"srsid"_s ).toElement().text().toLong( &ok );
 
    QDomNode node;
 
    if ( ok && srsid > 0 && srsid < Qgis::USER_CRS_START_ID )
    {
      node = srsNode.namedItem( u"authid"_s );
      if ( !node.isNull() )
      {
        createFromOgcWmsCrs( node.toElement().text() );
        if ( isValid() )
        {
          initialized = true;
        }
      }
 
      if ( !initialized )
      {
        node = srsNode.namedItem( u"epsg"_s );
        if ( !node.isNull() )
        {
          operator=( QgsCoordinateReferenceSystem::fromEpsgId( node.toElement().text().toLong() ) );
          if ( isValid() )
          {
            initialized = true;
          }
        }
      }
    }
 
    // if wkt is present, prefer that since it's lossless (unlike proj4 strings)
    if ( !initialized )
    {
      // before doing anything, we grab and set the stored CRS name (description).
      // this way if the stored CRS doesn't match anything available locally (i.e. from Proj's db
      // or the user's custom CRS list), then we will correctly show the CRS with its original
      // name (instead of just "custom crs")
      const QString description = srsNode.namedItem( u"description"_s ).toElement().text();
 
      const QString wkt = srsNode.namedItem( u"wkt"_s ).toElement().text();
      initialized = createFromWktInternal( wkt, description );
    }
 
    if ( !initialized )
    {
      node = srsNode.namedItem( u"proj4"_s );
      const QString proj4 = node.toElement().text();
      initialized = createFromProj( proj4 );
    }
 
    if ( !initialized )
    {
      // Setting from elements one by one
      node = srsNode.namedItem( u"proj4"_s );
      const QString proj4 = node.toElement().text();
      if ( !proj4.trimmed().isEmpty() )
        setProjString( node.toElement().text() );
 
      node = srsNode.namedItem( u"srsid"_s );
      d->mSrsId = node.toElement().text().toLong();
 
      node = srsNode.namedItem( u"srid"_s );
      d->mSRID = node.toElement().text().toLong();
 
      node = srsNode.namedItem( u"authid"_s );
      d->mAuthId = node.toElement().text();
 
      node = srsNode.namedItem( u"description"_s );
      d->mDescription = node.toElement().text();
 
      node = srsNode.namedItem( u"projectionacronym"_s );
      d->mProjectionAcronym = node.toElement().text();
 
      node = srsNode.namedItem( u"ellipsoidacronym"_s );
      d->mEllipsoidAcronym = node.toElement().text();
 
      node = srsNode.namedItem( u"geographicflag"_s );
      d->mIsGeographic = node.toElement().text() == "true"_L1;
 
      d->mWktPreferred.clear();
 
      //make sure the map units have been set
      setMapUnits();
    }
 
    const QString epoch = srsNode.toElement().attribute( u"coordinateEpoch"_s );
    if ( !epoch.isEmpty() )
    {
      bool epochOk = false;
      d->mCoordinateEpoch = epoch.toDouble( &epochOk );
      if ( !epochOk )
        d->mCoordinateEpoch = std::numeric_limits< double >::quiet_NaN();
    }
    else
    {
      d->mCoordinateEpoch = std::numeric_limits< double >::quiet_NaN();
    }
 
    mNativeFormat = qgsEnumKeyToValue<Qgis::CrsDefinitionFormat>( srsNode.toElement().attribute( u"nativeFormat"_s ), Qgis::CrsDefinitionFormat::Wkt );
 
    const QDomNode topoBaseCrsNode = srsNode.namedItem( u"topocentricBaseCrs"_s );
    if ( !topoBaseCrsNode.isNull() )
    {
      d->mTopocentricBaseCrs = std::make_unique<QgsCoordinateReferenceSystem>();
      d->mTopocentricBaseCrs->readXml( topoBaseCrsNode );
    }
  }
  else
  {
    // Return empty CRS if none was found in the XML.
    d = new QgsCoordinateReferenceSystemPrivate();
    result = false;
  }
  return result;
}
 
bool QgsCoordinateReferenceSystem::writeXml( QDomNode &node, QDomDocument &doc ) const
{
  QDomElement layerNode = node.toElement();
  QDomElement srsElement = doc.createElement( u"spatialrefsys"_s );
 
  srsElement.setAttribute( u"nativeFormat"_s, qgsEnumValueToKey<Qgis::CrsDefinitionFormat>( mNativeFormat ) );
 
  if ( std::isfinite( d->mCoordinateEpoch ) )
  {
    srsElement.setAttribute( u"coordinateEpoch"_s, d->mCoordinateEpoch );
  }
 
  QDomElement wktElement = doc.createElement( u"wkt"_s );
  wktElement.appendChild( doc.createTextNode( toWkt( Qgis::CrsWktVariant::Preferred ) ) );
  srsElement.appendChild( wktElement );
 
  QDomElement proj4Element = doc.createElement( u"proj4"_s );
  proj4Element.appendChild( doc.createTextNode( toProj() ) );
  srsElement.appendChild( proj4Element );
 
  QDomElement srsIdElement = doc.createElement( u"srsid"_s );
  srsIdElement.appendChild( doc.createTextNode( QString::number( srsid() ) ) );
  srsElement.appendChild( srsIdElement );
 
  QDomElement sridElement = doc.createElement( u"srid"_s );
  sridElement.appendChild( doc.createTextNode( QString::number( postgisSrid() ) ) );
  srsElement.appendChild( sridElement );
 
  QDomElement authidElement = doc.createElement( u"authid"_s );
  authidElement.appendChild( doc.createTextNode( authid() ) );
  srsElement.appendChild( authidElement );
 
  QDomElement descriptionElement = doc.createElement( u"description"_s );
  descriptionElement.appendChild( doc.createTextNode( description() ) );
  srsElement.appendChild( descriptionElement );
 
  QDomElement projectionAcronymElement = doc.createElement( u"projectionacronym"_s );
  projectionAcronymElement.appendChild( doc.createTextNode( projectionAcronym() ) );
  srsElement.appendChild( projectionAcronymElement );
 
  QDomElement ellipsoidAcronymElement = doc.createElement( u"ellipsoidacronym"_s );
  ellipsoidAcronymElement.appendChild( doc.createTextNode( ellipsoidAcronym() ) );
  srsElement.appendChild( ellipsoidAcronymElement );
 
  QDomElement geographicFlagElement = doc.createElement( u"geographicflag"_s );
  QString geoFlagText = u"false"_s;
  if ( isGeographic() )
  {
    geoFlagText = u"true"_s;
  }
 
  geographicFlagElement.appendChild( doc.createTextNode( geoFlagText ) );
  srsElement.appendChild( geographicFlagElement );
 
  if ( d->mTopocentricBaseCrs && d->mTopocentricBaseCrs->isValid() )
  {
    QDomElement topoBaseCrsElement = doc.createElement( u"topocentricBaseCrs"_s );
    d->mTopocentricBaseCrs->writeXml( topoBaseCrsElement, doc );
    srsElement.appendChild( topoBaseCrsElement );
  }
 
  layerNode.appendChild( srsElement );
 
  return true;
}
 
//
// Static helper methods below this point only please!
//
 
 
// Returns the whole proj4 string for the selected srsid
//this is a static method! NOTE I've made it private for now to reduce API clutter TS
QString QgsCoordinateReferenceSystem::projFromSrsId( const int srsId )
{
  QString myDatabaseFileName;
  QString myProjString;
  QString mySql = u"select parameters from tbl_srs where srs_id = %1 order by deprecated"_s.arg( srsId );
 
  //
  // Determine if this is a user projection or a system on
  // user projection defs all have srs_id >= 100000
  //
  if ( srsId >= Qgis::USER_CRS_START_ID )
  {
    myDatabaseFileName = QgsApplication::qgisUserDatabaseFilePath();
    QFileInfo myFileInfo;
    myFileInfo.setFile( myDatabaseFileName );
    if ( !myFileInfo.exists() ) //its unlikely that this condition will ever be reached
    {
      QgsDebugError( u"users qgis.db not found"_s );
      return QString();
    }
  }
  else //must be  a system projection then
  {
    myDatabaseFileName = QgsApplication::srsDatabaseFilePath();
  }
 
  sqlite3_database_unique_ptr database;
  sqlite3_statement_unique_ptr statement;
 
  int rc;
  rc = openDatabase( myDatabaseFileName, database );
  if ( rc )
  {
    return QString();
  }
 
  statement = database.prepare( mySql, rc );
 
  if ( rc == SQLITE_OK )
  {
    if ( statement.step() == SQLITE_ROW )
    {
      myProjString = statement.columnAsText( 0 );
    }
  }
 
  return myProjString;
}
 
int QgsCoordinateReferenceSystem::openDatabase( const QString &path, sqlite3_database_unique_ptr &database, bool readonly )
{
  int myResult;
  if ( readonly )
    myResult = database.open_v2( path, SQLITE_OPEN_READONLY, nullptr );
  else
    myResult = database.open( path );
 
  if ( myResult != SQLITE_OK )
  {
    QgsDebugError( "Can't open database: " + database.errorMessage() );
    // XXX This will likely never happen since on open, sqlite creates the
    //     database if it does not exist.
    // ... unfortunately it happens on Windows
    QgsMessageLog::logMessage( QObject::tr( "Could not open CRS database %1\nError(%2): %3" ).arg( path ).arg( myResult ).arg( database.errorMessage() ), QObject::tr( "CRS" ) );
  }
  return myResult;
}
 
void QgsCoordinateReferenceSystem::setCustomCrsValidation( CUSTOM_CRS_VALIDATION f )
{
  sCustomSrsValidation = f;
}
 
CUSTOM_CRS_VALIDATION QgsCoordinateReferenceSystem::customCrsValidation()
{
  return sCustomSrsValidation;
}
 
void QgsCoordinateReferenceSystem::debugPrint()
{
  QgsDebugMsgLevel( u"***SpatialRefSystem***"_s, 1 );
  QgsDebugMsgLevel( "* Valid : " + ( d->mIsValid ? QString( "true" ) : QString( "false" ) ), 1 );
  QgsDebugMsgLevel( "* SrsId : " + QString::number( d->mSrsId ), 1 );
  QgsDebugMsgLevel( "* Proj4 : " + toProj(), 1 );
  QgsDebugMsgLevel( "* WKT   : " + toWkt( Qgis::CrsWktVariant::Preferred ), 1 );
  QgsDebugMsgLevel( "* Desc. : " + d->mDescription, 1 );
  if ( mapUnits() == Qgis::DistanceUnit::Meters )
  {
    QgsDebugMsgLevel( u"* Units : meters"_s, 1 );
  }
  else if ( mapUnits() == Qgis::DistanceUnit::Feet )
  {
    QgsDebugMsgLevel( u"* Units : feet"_s, 1 );
  }
  else if ( mapUnits() == Qgis::DistanceUnit::Degrees )
  {
    QgsDebugMsgLevel( u"* Units : degrees"_s, 1 );
  }
}
 
void QgsCoordinateReferenceSystem::setValidationHint( const QString &html )
{
  mValidationHint = html;
}
 
QString QgsCoordinateReferenceSystem::validationHint() const
{
  return mValidationHint;
}
 
long QgsCoordinateReferenceSystem::saveAsUserCrs( const QString &name, Qgis::CrsDefinitionFormat nativeFormat )
{
  return QgsApplication::coordinateReferenceSystemRegistry()->addUserCrs( *this, name, nativeFormat );
}
 
void QgsCoordinateReferenceSystem::setNativeFormat( Qgis::CrsDefinitionFormat format )
{
  mNativeFormat = format;
}
 
Qgis::CrsDefinitionFormat QgsCoordinateReferenceSystem::nativeFormat() const
{
  return mNativeFormat;
}
 
long QgsCoordinateReferenceSystem::getRecordCount()
{
  sqlite3_database_unique_ptr database;
  sqlite3_statement_unique_ptr statement;
  int myResult;
  long myRecordCount = 0;
  //check the db is available
  myResult = database.open_v2( QgsApplication::qgisUserDatabaseFilePath(), SQLITE_OPEN_READONLY, nullptr );
  if ( myResult != SQLITE_OK )
  {
    QgsDebugError( u"Can't open database: %1"_s.arg( database.errorMessage() ) );
    return 0;
  }
  // Set up the query to retrieve the projection information needed to populate the ELLIPSOID list
  QString mySql = u"select count(*) from tbl_srs"_s;
  statement = database.prepare( mySql, myResult );
  if ( myResult == SQLITE_OK )
  {
    if ( statement.step() == SQLITE_ROW )
    {
      QString myRecordCountString = statement.columnAsText( 0 );
      myRecordCount = myRecordCountString.toLong();
    }
  }
  return myRecordCount;
}
 
bool testIsGeographic( PJ *crs )
{
  PJ_CONTEXT *pjContext = QgsProjContext::get();
  bool isGeographic = false;
 
  // check horizontal CRS units
  QgsProjUtils::proj_pj_unique_ptr horizontalCrs( QgsProjUtils::crsToHorizontalCrs( crs ) );
  if ( !horizontalCrs )
    return false;
 
  QgsProjUtils::proj_pj_unique_ptr coordinateSystem( proj_crs_get_coordinate_system( pjContext, horizontalCrs.get() ) );
  if ( coordinateSystem )
  {
    const int axisCount = proj_cs_get_axis_count( pjContext, coordinateSystem.get() );
    if ( axisCount > 0 )
    {
      const char *outUnitAuthName = nullptr;
      const char *outUnitAuthCode = nullptr;
      // Read only first axis
      proj_cs_get_axis_info( pjContext, coordinateSystem.get(), 0, nullptr, nullptr, nullptr, nullptr, nullptr, &outUnitAuthName, &outUnitAuthCode );
 
      if ( outUnitAuthName && outUnitAuthCode )
      {
        const char *unitCategory = nullptr;
        if ( proj_uom_get_info_from_database( pjContext, outUnitAuthName, outUnitAuthCode, nullptr, nullptr, &unitCategory ) )
        {
          isGeographic = QString( unitCategory ).compare( "angular"_L1, Qt::CaseInsensitive ) == 0;
        }
      }
    }
  }
  return isGeographic;
}
 
void getOperationAndEllipsoidFromProjString( const QString &proj, QString &operation, QString &ellipsoid )
{
  thread_local const QRegularExpression projRegExp( u"\\+proj=(\\S+)"_s );
  const QRegularExpressionMatch projMatch = projRegExp.match( proj );
  if ( !projMatch.hasMatch() )
  {
    QgsDebugMsgLevel( u"no +proj argument found [%2]"_s.arg( proj ), 2 );
    return;
  }
  operation = projMatch.captured( 1 );
 
  const QRegularExpressionMatch ellipseMatch = projRegExp.match( proj );
  if ( ellipseMatch.hasMatch() )
  {
    ellipsoid = ellipseMatch.captured( 1 );
  }
  else
  {
    // satisfy not null constraint on ellipsoid_acronym field
    // possibly we should drop the constraint, yet the crses with missing ellipsoid_acronym are malformed
    // and will result in oddities within other areas of QGIS (e.g. project ellipsoid won't be correctly
    // set for these CRSes). Better just hack around and make the constraint happy for now,
    // and hope that the definitions get corrected in future.
    ellipsoid = "";
  }
}
 
 
bool QgsCoordinateReferenceSystem::loadFromAuthCode( const QString &auth, const QString &code )
{
  if ( !QgsApplication::coordinateReferenceSystemRegistry()->authorities().contains( auth.toLower() ) )
    return false;
 
  d.detach();
  d->mIsValid = false;
  d->mWktPreferred.clear();
 
  PJ_CONTEXT *pjContext = QgsProjContext::get();
  QgsProjUtils::proj_pj_unique_ptr crs( proj_create_from_database( pjContext, auth.toUtf8().constData(), code.toUtf8().constData(), PJ_CATEGORY_CRS, false, nullptr ) );
  if ( !crs )
  {
    return false;
  }
 
  crs = QgsProjUtils::unboundCrs( crs.get() );
 
  QString proj4 = getFullProjString( crs.get() );
  proj4.replace( "+type=crs"_L1, QString() );
  proj4 = proj4.trimmed();
 
  d->mIsValid = true;
  d->mProj4 = proj4;
  d->mWktPreferred.clear();
  d->mDescription = QString( proj_get_name( crs.get() ) );
  d->mAuthId = u"%1:%2"_s.arg( auth, code );
  d->mIsGeographic = testIsGeographic( crs.get() );
  d->mAxisInvertedDirty = true;
  QString operation;
  QString ellipsoid;
  getOperationAndEllipsoidFromProjString( proj4, operation, ellipsoid );
  d->mProjectionAcronym = operation;
  d->mEllipsoidAcronym.clear();
  d->setPj( std::move( crs ) );
 
  const QString dbVals = sAuthIdToQgisSrsIdMap.value( u"%1:%2"_s.arg( auth, code ).toUpper() );
  if ( !dbVals.isEmpty() )
  {
    const QStringList parts = dbVals.split( ',' );
    d->mSrsId = parts.at( 0 ).toInt();
    d->mSRID = parts.at( 1 ).toInt();
  }
 
  setMapUnits();
 
  return true;
}
 
QList<long> QgsCoordinateReferenceSystem::userSrsIds()
{
  QList<long> results;
  // check user defined projection database
  const QString db = QgsApplication::qgisUserDatabaseFilePath();
 
  QFileInfo myInfo( db );
  if ( !myInfo.exists() )
  {
    QgsDebugError( "failed : " + db + " does not exist!" );
    return results;
  }
 
  sqlite3_database_unique_ptr database;
  sqlite3_statement_unique_ptr statement;
 
  //check the db is available
  int result = openDatabase( db, database );
  if ( result != SQLITE_OK )
  {
    QgsDebugError( "failed : " + db + " could not be opened!" );
    return results;
  }
 
  QString sql = u"select srs_id from tbl_srs where srs_id >= %1"_s.arg( Qgis::USER_CRS_START_ID );
  int rc;
  statement = database.prepare( sql, rc );
  while ( true )
  {
    int ret = statement.step();
 
    if ( ret == SQLITE_DONE )
    {
      // there are no more rows to fetch - we can stop looping
      break;
    }
 
    if ( ret == SQLITE_ROW )
    {
      results.append( statement.columnAsInt64( 0 ) );
    }
    else
    {
      QgsMessageLog::logMessage( QObject::tr( "SQLite error: %2\nSQL: %1" ).arg( sql, sqlite3_errmsg( database.get() ) ), QObject::tr( "SpatiaLite" ) );
      break;
    }
  }
 
  return results;
}
 
long QgsCoordinateReferenceSystem::matchToUserCrs() const
{
  PJ *obj = d->threadLocalProjObject();
  if ( !obj )
    return 0;
 
  const QList< long > ids = userSrsIds();
  for ( long id : ids )
  {
    QgsCoordinateReferenceSystem candidate = QgsCoordinateReferenceSystem::fromSrsId( id );
    if ( candidate.projObject() && proj_is_equivalent_to( obj, candidate.projObject(), PJ_COMP_EQUIVALENT ) )
    {
      return id;
    }
  }
  return 0;
}
 
static void sync_db_proj_logger( void * /* user_data */, int level, const char *message )
{
#ifndef QGISDEBUG
  Q_UNUSED( message )
#endif
  if ( level == PJ_LOG_ERROR )
  {
    QgsDebugMsgLevel( u"PROJ: %1"_s.arg( message ), 2 );
  }
  else if ( level == PJ_LOG_DEBUG )
  {
    QgsDebugMsgLevel( u"PROJ: %1"_s.arg( message ), 3 );
  }
}
 
int QgsCoordinateReferenceSystem::syncDatabase()
{
  setlocale( LC_ALL, "C" );
  QString dbFilePath = QgsApplication::srsDatabaseFilePath();
 
  int inserted = 0, updated = 0, deleted = 0, errors = 0;
 
  QgsDebugMsgLevel( u"Load srs db from: %1"_s.arg( QgsApplication::srsDatabaseFilePath().toLocal8Bit().constData() ), 4 );
 
  sqlite3_database_unique_ptr database;
  if ( database.open( dbFilePath ) != SQLITE_OK )
  {
    QgsDebugError( u"Could not open database: %1 (%2)\n"_s.arg( QgsApplication::srsDatabaseFilePath(), database.errorMessage() ) );
    return -1;
  }
 
  if ( sqlite3_exec( database.get(), "BEGIN TRANSACTION", nullptr, nullptr, nullptr ) != SQLITE_OK )
  {
    QgsDebugError( u"Could not begin transaction: %1 (%2)\n"_s.arg( QgsApplication::srsDatabaseFilePath(), database.errorMessage() ) );
    return -1;
  }
 
  sqlite3_statement_unique_ptr statement;
  int result;
  char *errMsg = nullptr;
 
  bool createdTypeColumn = false;
  if ( sqlite3_exec( database.get(), "ALTER TABLE tbl_srs ADD COLUMN srs_type text", nullptr, nullptr, nullptr ) == SQLITE_OK )
  {
    createdTypeColumn = true;
    if ( sqlite3_exec( database.get(), "CREATE INDEX srs_type ON tbl_srs(srs_type)", nullptr, nullptr, nullptr ) != SQLITE_OK )
    {
      QgsDebugError( u"Could not create index for srs_type"_s );
      return -1;
    }
  }
 
  if ( sqlite3_exec( database.get(), "create table tbl_info (proj_major INT, proj_minor INT, proj_patch INT)", nullptr, nullptr, nullptr ) == SQLITE_OK )
  {
    QString sql = u"INSERT INTO tbl_info(proj_major, proj_minor, proj_patch) VALUES (%1, %2,%3)"_s
                    .arg( QString::number( PROJ_VERSION_MAJOR ), QString::number( PROJ_VERSION_MINOR ), QString::number( PROJ_VERSION_PATCH ) );
    if ( sqlite3_exec( database.get(), sql.toUtf8(), nullptr, nullptr, &errMsg ) != SQLITE_OK )
    {
      QgsDebugError( u"Could not execute: %1 [%2/%3]\n"_s.arg( sql, database.errorMessage(), errMsg ? errMsg : "(unknown error)" ) );
      if ( errMsg )
        sqlite3_free( errMsg );
      return -1;
    }
  }
  else
  {
    // retrieve last update details
    QString sql = u"SELECT proj_major, proj_minor, proj_patch FROM tbl_info"_s;
    statement = database.prepare( sql, result );
    if ( result != SQLITE_OK )
    {
      QgsDebugError( u"Could not prepare: %1 [%2]\n"_s.arg( sql, database.errorMessage() ) );
      return -1;
    }
    if ( statement.step() == SQLITE_ROW )
    {
      int major = statement.columnAsInt64( 0 );
      int minor = statement.columnAsInt64( 1 );
      int patch = statement.columnAsInt64( 2 );
      if ( !createdTypeColumn && major == PROJ_VERSION_MAJOR && minor == PROJ_VERSION_MINOR && patch == PROJ_VERSION_PATCH )
        // yay, nothing to do!
        return 0;
    }
    else
    {
      QgsDebugError( u"Could not retrieve previous CRS sync PROJ version number"_s );
      return -1;
    }
  }
 
  PJ_CONTEXT *pjContext = QgsProjContext::get();
  // silence proj warnings
  proj_log_func( pjContext, nullptr, sync_db_proj_logger );
 
  PROJ_STRING_LIST authorities = proj_get_authorities_from_database( pjContext );
 
  int nextSrsId = 67218;
  int nextSrId = 520007218;
  for ( auto authIter = authorities; authIter && *authIter; ++authIter )
  {
    const QString authority( *authIter );
    QgsDebugMsgLevel( u"Loading authority '%1'"_s.arg( authority ), 2 );
    PROJ_STRING_LIST codes = proj_get_codes_from_database( pjContext, *authIter, PJ_TYPE_CRS, true );
 
    QStringList allCodes;
 
    for ( auto codesIter = codes; codesIter && *codesIter; ++codesIter )
    {
      const QString code( *codesIter );
      allCodes << QgsSqliteUtils::quotedString( code );
      QgsDebugMsgLevel( u"Loading code '%1'"_s.arg( code ), 4 );
      QgsProjUtils::proj_pj_unique_ptr crs( proj_create_from_database( pjContext, *authIter, *codesIter, PJ_CATEGORY_CRS, false, nullptr ) );
      if ( !crs )
      {
        QgsDebugError( u"Could not load '%1:%2'"_s.arg( authority, code ) );
        continue;
      }
 
      const PJ_TYPE pjType = proj_get_type( crs.get() );
 
      QString srsTypeString;
      // NOLINTBEGIN(bugprone-branch-clone)
      switch ( pjType )
      {
        // don't need these in the CRS db
        case PJ_TYPE_ELLIPSOID:
        case PJ_TYPE_PRIME_MERIDIAN:
        case PJ_TYPE_GEODETIC_REFERENCE_FRAME:
        case PJ_TYPE_DYNAMIC_GEODETIC_REFERENCE_FRAME:
        case PJ_TYPE_VERTICAL_REFERENCE_FRAME:
        case PJ_TYPE_DYNAMIC_VERTICAL_REFERENCE_FRAME:
        case PJ_TYPE_DATUM_ENSEMBLE:
        case PJ_TYPE_CONVERSION:
        case PJ_TYPE_TRANSFORMATION:
        case PJ_TYPE_CONCATENATED_OPERATION:
        case PJ_TYPE_OTHER_COORDINATE_OPERATION:
        case PJ_TYPE_TEMPORAL_DATUM:
        case PJ_TYPE_ENGINEERING_DATUM:
        case PJ_TYPE_PARAMETRIC_DATUM:
        case PJ_TYPE_UNKNOWN:
          continue;
 
        case PJ_TYPE_CRS:
        case PJ_TYPE_GEOGRAPHIC_CRS:
          continue; // not possible
 
        case PJ_TYPE_GEODETIC_CRS:
          srsTypeString = qgsEnumValueToKey( Qgis::CrsType::Geodetic );
          break;
 
        case PJ_TYPE_GEOCENTRIC_CRS:
          srsTypeString = qgsEnumValueToKey( Qgis::CrsType::Geocentric );
          break;
 
        case PJ_TYPE_GEOGRAPHIC_2D_CRS:
          srsTypeString = qgsEnumValueToKey( Qgis::CrsType::Geographic2d );
          break;
 
        case PJ_TYPE_GEOGRAPHIC_3D_CRS:
          srsTypeString = qgsEnumValueToKey( Qgis::CrsType::Geographic3d );
          break;
 
        case PJ_TYPE_PROJECTED_CRS:
          srsTypeString = qgsEnumValueToKey( Qgis::CrsType::Projected );
          break;
 
        case PJ_TYPE_COMPOUND_CRS:
          srsTypeString = qgsEnumValueToKey( Qgis::CrsType::Compound );
          break;
 
        case PJ_TYPE_TEMPORAL_CRS:
          srsTypeString = qgsEnumValueToKey( Qgis::CrsType::Temporal );
          break;
 
        case PJ_TYPE_ENGINEERING_CRS:
          srsTypeString = qgsEnumValueToKey( Qgis::CrsType::Engineering );
          break;
 
        case PJ_TYPE_BOUND_CRS:
          srsTypeString = qgsEnumValueToKey( Qgis::CrsType::Bound );
          break;
 
        case PJ_TYPE_VERTICAL_CRS:
          srsTypeString = qgsEnumValueToKey( Qgis::CrsType::Vertical );
          break;
 
#if PROJ_VERSION_MAJOR > 9 || ( PROJ_VERSION_MAJOR == 9 && PROJ_VERSION_MINOR >= 2 )
        case PJ_TYPE_DERIVED_PROJECTED_CRS:
          srsTypeString = qgsEnumValueToKey( Qgis::CrsType::DerivedProjected );
          break;
        case PJ_TYPE_COORDINATE_METADATA:
          continue;
#endif
        case PJ_TYPE_OTHER_CRS:
          srsTypeString = qgsEnumValueToKey( Qgis::CrsType::Other );
          break;
      }
      // NOLINTEND(bugprone-branch-clone)
 
      crs = QgsProjUtils::unboundCrs( crs.get() );
 
      QString proj4 = getFullProjString( crs.get() );
      proj4.replace( "+type=crs"_L1, QString() );
      proj4 = proj4.trimmed();
 
      if ( proj4.isEmpty() )
      {
        QgsDebugMsgLevel( u"No proj4 for '%1:%2'"_s.arg( authority, code ), 2 );
        // satisfy not null constraint
        proj4 = "";
      }
 
      // there's a not-null constraint on these columns, so we must use empty strings instead
      QString operation = "";
      QString ellps = "";
      getOperationAndEllipsoidFromProjString( proj4, operation, ellps );
 
      const QString translatedOperation = QgsCoordinateReferenceSystemUtils::translateProjection( operation );
      if ( translatedOperation.isEmpty() && !operation.isEmpty() )
      {
        std::cout << u"Operation needs translation in QgsCoordinateReferenceSystemUtils::translateProjection: %1"_s.arg( operation ).toLocal8Bit().constData() << std::endl;
        qFatal( "aborted" );
      }
 
      const bool deprecated = proj_is_deprecated( crs.get() );
      const QString name( proj_get_name( crs.get() ) );
 
      QString sql = u"SELECT parameters,description,deprecated,srs_type,projection_acronym FROM tbl_srs WHERE auth_name='%1' AND auth_id='%2'"_s.arg( authority, code );
      statement = database.prepare( sql, result );
      if ( result != SQLITE_OK )
      {
        QgsDebugError( u"Could not prepare: %1 [%2]\n"_s.arg( sql, database.errorMessage() ) );
        continue;
      }
 
      QString dbSrsProj4;
      QString dbSrsDesc;
      QString dbSrsType;
      QString dbOperation;
      bool dbSrsDeprecated = deprecated;
      if ( statement.step() == SQLITE_ROW )
      {
        dbSrsProj4 = statement.columnAsText( 0 );
        dbSrsDesc = statement.columnAsText( 1 );
        dbSrsDeprecated = statement.columnAsText( 2 ).toInt() != 0;
        dbSrsType = statement.columnAsText( 3 );
        dbOperation = statement.columnAsText( 4 );
      }
 
      if ( !dbSrsProj4.isEmpty() || !dbSrsDesc.isEmpty() )
      {
        if ( proj4 != dbSrsProj4 || name != dbSrsDesc || deprecated != dbSrsDeprecated || dbSrsType != srsTypeString || dbOperation != operation )
        {
          errMsg = nullptr;
          sql = u"UPDATE tbl_srs SET parameters=%1,description=%2,deprecated=%3, srs_type=%4,projection_acronym=%5 WHERE auth_name=%6 AND auth_id=%7"_s.arg( QgsSqliteUtils::quotedString( proj4 ) )
                  .arg( QgsSqliteUtils::quotedString( name ) )
                  .arg( deprecated ? 1 : 0 )
                  .arg( QgsSqliteUtils::quotedString( srsTypeString ), QgsSqliteUtils::quotedString( operation ), QgsSqliteUtils::quotedString( authority ), QgsSqliteUtils::quotedString( code ) );
 
          if ( sqlite3_exec( database.get(), sql.toUtf8(), nullptr, nullptr, &errMsg ) != SQLITE_OK )
          {
            QgsDebugError( u"Could not execute: %1 [%2/%3]\n"_s.arg( sql, database.errorMessage(), errMsg ? errMsg : "(unknown error)" ) );
            if ( errMsg )
              sqlite3_free( errMsg );
            errors++;
          }
          else
          {
            updated++;
          }
        }
      }
      else
      {
        const bool isGeographic = testIsGeographic( crs.get() );
 
        // work out srid and srsid
        const QString dbVals = sAuthIdToQgisSrsIdMap.value( u"%1:%2"_s.arg( authority, code ) );
        QString srsId;
        QString srId;
        if ( !dbVals.isEmpty() )
        {
          const QStringList parts = dbVals.split( ',' );
          srsId = parts.at( 0 );
          srId = parts.at( 1 );
        }
        if ( srId.isEmpty() )
        {
          srId = QString::number( nextSrId );
          nextSrId++;
        }
        if ( srsId.isEmpty() )
        {
          srsId = QString::number( nextSrsId );
          nextSrsId++;
        }
 
        if ( !srsId.isEmpty() )
        {
          sql = u"INSERT INTO tbl_srs(srs_id, description,projection_acronym,ellipsoid_acronym,parameters,srid,auth_name,auth_id,is_geo,deprecated,srs_type) VALUES (%1, %2,%3,%4,%5,%6,%7,%8,%9,%10,%11)"_s
                  .arg( srsId )
                  .arg( QgsSqliteUtils::quotedString( name ), QgsSqliteUtils::quotedString( operation ), QgsSqliteUtils::quotedString( ellps ), QgsSqliteUtils::quotedString( proj4 ) )
                  .arg( srId )
                  .arg( QgsSqliteUtils::quotedString( authority ) )
                  .arg( QgsSqliteUtils::quotedString( code ) )
                  .arg( isGeographic ? 1 : 0 )
                  .arg( deprecated ? 1 : 0 )
                  .arg( QgsSqliteUtils::quotedString( srsTypeString ) );
        }
        else
        {
          sql = u"INSERT INTO tbl_srs(description,projection_acronym,ellipsoid_acronym,parameters,srid,auth_name,auth_id,is_geo,deprecated,srs_type) VALUES (%1,%2,%3,%4,%5,%6,%7,%8,%9,%10)"_s
                  .arg( QgsSqliteUtils::quotedString( name ), QgsSqliteUtils::quotedString( operation ), QgsSqliteUtils::quotedString( ellps ), QgsSqliteUtils::quotedString( proj4 ) )
                  .arg( srId )
                  .arg( QgsSqliteUtils::quotedString( authority ) )
                  .arg( QgsSqliteUtils::quotedString( code ) )
                  .arg( isGeographic ? 1 : 0 )
                  .arg( deprecated ? 1 : 0 )
                  .arg( QgsSqliteUtils::quotedString( srsTypeString ) );
        }
 
        errMsg = nullptr;
        if ( sqlite3_exec( database.get(), sql.toUtf8(), nullptr, nullptr, &errMsg ) == SQLITE_OK )
        {
          inserted++;
        }
        else
        {
          qCritical( "Could not execute: %s [%s/%s]\n", sql.toLocal8Bit().constData(), sqlite3_errmsg( database.get() ), errMsg ? errMsg : "(unknown error)" );
          errors++;
 
          if ( errMsg )
            sqlite3_free( errMsg );
        }
      }
    }
 
    proj_string_list_destroy( codes );
 
    const QString sql = u"DELETE FROM tbl_srs WHERE auth_name='%1' AND NOT auth_id IN (%2)"_s.arg( authority, allCodes.join( ',' ) );
    if ( sqlite3_exec( database.get(), sql.toUtf8(), nullptr, nullptr, nullptr ) == SQLITE_OK )
    {
      deleted = sqlite3_changes( database.get() );
    }
    else
    {
      errors++;
      qCritical( "Could not execute: %s [%s]\n", sql.toLocal8Bit().constData(), sqlite3_errmsg( database.get() ) );
    }
  }
  proj_string_list_destroy( authorities );
 
  QString sql = u"UPDATE tbl_info set proj_major=%1,proj_minor=%2,proj_patch=%3"_s.arg( QString::number( PROJ_VERSION_MAJOR ), QString::number( PROJ_VERSION_MINOR ), QString::number( PROJ_VERSION_PATCH ) );
  if ( sqlite3_exec( database.get(), sql.toUtf8(), nullptr, nullptr, &errMsg ) != SQLITE_OK )
  {
    QgsDebugError( u"Could not execute: %1 [%2/%3]\n"_s.arg( sql, database.errorMessage(), errMsg ? errMsg : "(unknown error)" ) );
    if ( errMsg )
      sqlite3_free( errMsg );
    return -1;
  }
 
  if ( sqlite3_exec( database.get(), "COMMIT", nullptr, nullptr, nullptr ) != SQLITE_OK )
  {
    QgsDebugError( u"Could not commit transaction: %1 [%2]\n"_s.arg( QgsApplication::srsDatabaseFilePath(), sqlite3_errmsg( database.get() ) ) );
    return -1;
  }
 
#ifdef QGISDEBUG
  QgsDebugMsgLevel( u"CRS update (inserted:%1 updated:%2 deleted:%3 errors:%4)"_s.arg( inserted ).arg( updated ).arg( deleted ).arg( errors ), 4 );
#else
  Q_UNUSED( deleted )
#endif
 
  if ( errors > 0 )
    return -errors;
  else
    return updated + inserted;
}
 
const QHash<QString, QgsCoordinateReferenceSystem> &QgsCoordinateReferenceSystem::stringCache()
{
  return *sStringCache();
}
 
const QHash<QString, QgsCoordinateReferenceSystem> &QgsCoordinateReferenceSystem::projCache()
{
  return *sProj4Cache();
}
 
const QHash<QString, QgsCoordinateReferenceSystem> &QgsCoordinateReferenceSystem::ogcCache()
{
  return *sOgcCache();
}
 
const QHash<QString, QgsCoordinateReferenceSystem> &QgsCoordinateReferenceSystem::wktCache()
{
  return *sWktCache();
}
 
const QHash<long, QgsCoordinateReferenceSystem> &QgsCoordinateReferenceSystem::srIdCache()
{
  return *sSrIdCache();
}
 
const QHash<long, QgsCoordinateReferenceSystem> &QgsCoordinateReferenceSystem::srsIdCache()
{
  return *sSrsIdCache();
}
 
QgsCoordinateReferenceSystem QgsCoordinateReferenceSystem::toGeographicCrs() const
{
  if ( isGeographic() )
  {
    return *this;
  }
 
  if ( PJ *obj = d->threadLocalProjObject() )
  {
    PJ_CONTEXT *pjContext = QgsProjContext::get();
    QgsProjUtils::proj_pj_unique_ptr geoCrs( proj_crs_get_geodetic_crs( pjContext, obj ) );
    if ( !geoCrs )
      return QgsCoordinateReferenceSystem();
 
    const PJ_TYPE pjType = proj_get_type( geoCrs.get() );
    if ( pjType == PJ_TYPE_GEOCENTRIC_CRS )
    {
      // special case: while a geocentric crs IS a geodetic CRS, this particular QGIS method advertises
      // that it will always return a geographic latitude/longitude based CRS. So we build a geographic
      // CRS using the same datum as the original CRS
      QgsProjUtils::proj_pj_unique_ptr cs( proj_create_ellipsoidal_2D_cs( pjContext, PJ_ELLPS2D_LONGITUDE_LATITUDE, "Degree", 1.0 ) );
      QgsProjUtils::proj_pj_unique_ptr datum( proj_crs_get_datum( pjContext, geoCrs.get() ) );
      QgsProjUtils::proj_pj_unique_ptr datumEnsemble( proj_crs_get_datum_ensemble( pjContext, geoCrs.get() ) );
      QgsProjUtils::proj_pj_unique_ptr geoGraphicCrs( proj_create_geographic_crs_from_datum( pjContext, nullptr, datum ? datum.get() : datumEnsemble.get(), cs.get() ) );
      if ( !geoGraphicCrs )
        return QgsCoordinateReferenceSystem();
      return QgsCoordinateReferenceSystem::fromProjObject( geoGraphicCrs.get() );
    }
 
    if ( !testIsGeographic( geoCrs.get() ) )
      return QgsCoordinateReferenceSystem();
 
    QgsProjUtils::proj_pj_unique_ptr normalized( proj_normalize_for_visualization( pjContext, geoCrs.get() ) );
    if ( !normalized )
      return QgsCoordinateReferenceSystem();
 
    return QgsCoordinateReferenceSystem::fromProjObject( normalized.get() );
  }
  else
  {
    return QgsCoordinateReferenceSystem();
  }
}
 
QgsCoordinateReferenceSystem QgsCoordinateReferenceSystem::toGeocentricCrs() const
{
  if ( type() == Qgis::CrsType::Geocentric )
  {
    return *this;
  }
 
  if ( PJ *obj = d->threadLocalProjObject() )
  {
    PJ_CONTEXT *pjContext = QgsProjContext::get();
 
    // we need the horizontal, unbound crs in order to extract the datum:
    QgsProjUtils::proj_pj_unique_ptr horizontalCrs = QgsProjUtils::crsToHorizontalCrs( obj );
    if ( !horizontalCrs )
    {
      return QgsCoordinateReferenceSystem();
    }
 
    QgsProjUtils::proj_pj_unique_ptr datum( proj_crs_get_datum( pjContext, horizontalCrs.get() ) );
    QgsProjUtils::proj_pj_unique_ptr datumEnsemble( proj_crs_get_datum_ensemble( pjContext, horizontalCrs.get() ) );
    if ( !datum && !datumEnsemble )
      return QgsCoordinateReferenceSystem();
 
    QgsProjUtils::proj_pj_unique_ptr crs( proj_create_geocentric_crs_from_datum(
      pjContext,
      /*crs_name*/ nullptr,
      /*datum_or_datum_ensemble*/ datumEnsemble ? datumEnsemble.get() : datum.get(),
      /*linear_units*/ nullptr, // "NULL for meter"
      /*linear_units_conv*/ 0   // "0 for meter if linear_units == NULL"
    ) );
    if ( crs )
      return QgsCoordinateReferenceSystem::fromProjObject( crs.get() );
  }
 
  return QgsCoordinateReferenceSystem();
}
 
QgsCoordinateReferenceSystem QgsCoordinateReferenceSystem::horizontalCrs() const
{
  switch ( type() )
  {
    case Qgis::CrsType::Unknown:
    case Qgis::CrsType::Geodetic:
    case Qgis::CrsType::Geocentric:
    case Qgis::CrsType::Geographic2d:
    case Qgis::CrsType::Projected:
    case Qgis::CrsType::Temporal:
    case Qgis::CrsType::Engineering:
    case Qgis::CrsType::Bound:
    case Qgis::CrsType::Other:
    case Qgis::CrsType::DerivedProjected:
    case Qgis::CrsType::Geographic3d:
      return *this;
 
    case Qgis::CrsType::Vertical:
      return QgsCoordinateReferenceSystem();
 
    case Qgis::CrsType::Compound:
      break;
  }
 
  if ( PJ *obj = d->threadLocalProjObject() )
  {
    QgsProjUtils::proj_pj_unique_ptr hozCrs = QgsProjUtils::crsToHorizontalCrs( obj );
    if ( hozCrs )
      return QgsCoordinateReferenceSystem::fromProjObject( hozCrs.get() );
  }
  return QgsCoordinateReferenceSystem();
}
 
QgsCoordinateReferenceSystem QgsCoordinateReferenceSystem::verticalCrs() const
{
  switch ( type() )
  {
    case Qgis::CrsType::Unknown:
    case Qgis::CrsType::Geodetic:
    case Qgis::CrsType::Geocentric:
    case Qgis::CrsType::Geographic2d:
    case Qgis::CrsType::Projected:
    case Qgis::CrsType::Temporal:
    case Qgis::CrsType::Engineering:
    case Qgis::CrsType::Bound:
    case Qgis::CrsType::Other:
    case Qgis::CrsType::DerivedProjected:
    case Qgis::CrsType::Geographic3d:
      return QgsCoordinateReferenceSystem();
 
    case Qgis::CrsType::Vertical:
      return *this;
 
    case Qgis::CrsType::Compound:
      break;
  }
 
  if ( PJ *obj = d->threadLocalProjObject() )
  {
    QgsProjUtils::proj_pj_unique_ptr vertCrs = QgsProjUtils::crsToVerticalCrs( obj );
    if ( vertCrs )
      return QgsCoordinateReferenceSystem::fromProjObject( vertCrs.get() );
  }
  return QgsCoordinateReferenceSystem();
}
 
bool QgsCoordinateReferenceSystem::hasVerticalAxis() const
{
  if ( PJ *obj = d->threadLocalProjObject() )
  {
    return QgsProjUtils::hasVerticalAxis( obj );
  }
  return false;
}
 
QString QgsCoordinateReferenceSystem::geographicCrsAuthId() const
{
  if ( isGeographic() )
  {
    return d->mAuthId;
  }
  else if ( PJ *obj = d->threadLocalProjObject() )
  {
    QgsProjUtils::proj_pj_unique_ptr geoCrs( proj_crs_get_geodetic_crs( QgsProjContext::get(), obj ) );
    return geoCrs ? u"%1:%2"_s.arg( proj_get_id_auth_name( geoCrs.get(), 0 ), proj_get_id_code( geoCrs.get(), 0 ) ) : QString();
  }
  else
  {
    return QString();
  }
}
 
bool QgsCoordinateReferenceSystem::topocentricOrigin( double &latDeg, double &lonDeg ) const
{
  if ( !isValid() )
    return false;
 
  PJ_CONTEXT *ctx = QgsProjContext::get();
  const PJ *pj = projObject();
  if ( !pj )
    return false;
 
  QgsProjUtils::proj_pj_unique_ptr conversion( proj_crs_get_coordoperation( ctx, pj ) );
  if ( !conversion )
    return false;
 
  const int paramCount = proj_coordoperation_get_param_count( ctx, conversion.get() );
  bool hasLat = false, hasLon = false;
 
  for ( int i = 0; i < paramCount; i++ )
  {
    const char *paramName = nullptr, *authName = nullptr, *code = nullptr, *valueString = nullptr, *unitName = nullptr, *unitAuthName = nullptr, *unitCode = nullptr, *unitCategory = nullptr;
    double value = 0.0;
    double unitConvFactor = 1.0;
    if ( !proj_coordoperation_get_param( ctx, conversion.get(), i, &paramName, &authName, &code, &value, &valueString, &unitConvFactor, &unitName, &unitAuthName, &unitCode, &unitCategory ) )
      continue;
 
    const int paramCode = QString( code ).toInt();
 
    if ( paramCode == EPSG_CODE_PARAMETER_LATITUDE_TOPOGRAPHIC_ORIGIN )
    {
      latDeg = value * unitConvFactor * 180.0 / M_PI;
      hasLat = true;
    }
    else if ( paramCode == EPSG_CODE_PARAMETER_LONGITUDE_TOPOGRAPHIC_ORIGIN )
    {
      lonDeg = value * unitConvFactor * 180.0 / M_PI;
      hasLon = true;
    }
  }
 
  return hasLat && hasLon;
}
 
 
QgsCoordinateReferenceSystem QgsCoordinateReferenceSystem::toTopocentricCrs( double latitude, double longitude ) const
{
  if ( !isValid() )
    return QgsCoordinateReferenceSystem();
 
 
  PJ_CONTEXT *ctx = QgsProjContext::get();
  const PJ *pj = projObject();
  if ( !pj )
    return QgsCoordinateReferenceSystem();
 
  QgsProjUtils::proj_pj_unique_ptr baseCrs;
  double lat = 0.0, lon = 0.0;
  const bool isTopocentric = topocentricOrigin( lat, lon );
  if ( isTopocentric )
  {
    baseCrs.reset( proj_get_source_crs( ctx, pj ) );
    if ( baseCrs )
      pj = baseCrs.get();
  }
 
  QgsProjUtils::proj_pj_unique_ptr geodeticCrs( proj_crs_get_geodetic_crs( ctx, pj ) );
  if ( !geodeticCrs )
    return QgsCoordinateReferenceSystem();
 
  if ( proj_get_type( geodeticCrs.get() ) == PJ_TYPE_GEOGRAPHIC_2D_CRS )
  {
    QgsProjUtils::proj_pj_unique_ptr geodetic3DCrs( proj_crs_promote_to_3D( ctx, nullptr, geodeticCrs.get() ) );
    if ( geodetic3DCrs )
      geodeticCrs = std::move( geodetic3DCrs );
  }
 
  if ( proj_get_type( geodeticCrs.get() ) == PJ_TYPE_GEOCENTRIC_CRS )
  {
    QgsProjUtils::proj_pj_unique_ptr cs( proj_create_ellipsoidal_2D_cs( ctx, PJ_ELLPS2D_LONGITUDE_LATITUDE, "Degree", 0 ) );
    QgsProjUtils::proj_pj_unique_ptr datum( proj_crs_get_datum( ctx, geodeticCrs.get() ) );
    QgsProjUtils::proj_pj_unique_ptr datumEnsemble( proj_crs_get_datum_ensemble( ctx, geodeticCrs.get() ) );
    if ( !datum && !datumEnsemble )
      return QgsCoordinateReferenceSystem();
    geodeticCrs.reset( proj_create_geographic_crs_from_datum( ctx, nullptr, datum ? datum.get() : datumEnsemble.get(), cs.get() ) );
    if ( !geodeticCrs )
      return QgsCoordinateReferenceSystem();
  }
 
  const char *projWkt = proj_as_wkt( ctx, geodeticCrs.get(), PJ_WKT2_2019, nullptr );
  if ( !projWkt )
    return QgsCoordinateReferenceSystem();
 
  QString baseWkt = QString::fromUtf8( projWkt );
  if ( baseWkt.startsWith( "GEOGCRS["_L1 ) )
    baseWkt.replace( 0, 8, "BASEGEOGCRS["_L1 );
  else
    return QgsCoordinateReferenceSystem();
 
  const QString wkt = u"GEODCRS[\"Topocentric\","
                      "%1,"
                      "DERIVINGCONVERSION[\"Geographic/topocentric\","
                      "  METHOD[\"Geographic/topocentric conversions\",ID[\"EPSG\",9837]],"
                      "  PARAMETER[\"Latitude of topocentric origin\",%2,"
                      "    ANGLEUNIT[\"degree\",0.0174532925199433],ID[\"EPSG\",8834]],"
                      "  PARAMETER[\"Longitude of topocentric origin\",%3,"
                      "    ANGLEUNIT[\"degree\",0.0174532925199433],ID[\"EPSG\",8835]],"
                      "  PARAMETER[\"Ellipsoidal height of topocentric origin\",0,"
                      "    LENGTHUNIT[\"metre\",1],ID[\"EPSG\",8836]]],"
                      "CS[Cartesian,3],"
                      "AXIS[\"(X)\",east,ORDER[1],LENGTHUNIT[\"metre\",1]],"
                      "AXIS[\"(Y)\",north,ORDER[2],LENGTHUNIT[\"metre\",1]],"
                      "AXIS[\"(Z)\",up,ORDER[3],LENGTHUNIT[\"metre\",1]]]"_s.arg( baseWkt )
                        .arg( qgsDoubleToString( latitude ) )
                        .arg( qgsDoubleToString( longitude ) );
 
  QgsProjUtils::proj_pj_unique_ptr topocentric( proj_create_from_wkt( ctx, wkt.toUtf8().constData(), nullptr, nullptr, nullptr ) );
  if ( !topocentric )
    return QgsCoordinateReferenceSystem();
 
  QgsCoordinateReferenceSystem crs = QgsCoordinateReferenceSystem::fromProjObject( topocentric.get() );
  if ( isTopocentric )
    crs.d->mTopocentricBaseCrs = std::make_unique<QgsCoordinateReferenceSystem>( *d->mTopocentricBaseCrs );
  else
    crs.d->mTopocentricBaseCrs = std::make_unique<QgsCoordinateReferenceSystem>( *this );
 
  return crs;
}
 
QgsCoordinateReferenceSystem QgsCoordinateReferenceSystem::topocentricBaseCrs() const
{
  return d->mTopocentricBaseCrs ? *d->mTopocentricBaseCrs : QgsCoordinateReferenceSystem();
}
 
PJ *QgsCoordinateReferenceSystem::projObject() const
{
  return d->threadLocalProjObject();
}
 
QgsCoordinateReferenceSystem QgsCoordinateReferenceSystem::fromProjObject( PJ *object )
{
  QgsCoordinateReferenceSystem crs;
  crs.createFromProjObject( object );
  return crs;
}
 
bool QgsCoordinateReferenceSystem::createFromProjObject( PJ *object )
{
  d.detach();
  d->mIsValid = false;
  d->mProj4.clear();
  d->mWktPreferred.clear();
 
  if ( !object )
  {
    return false;
  }
 
  switch ( proj_get_type( object ) )
  {
    case PJ_TYPE_GEODETIC_CRS:
    case PJ_TYPE_GEOCENTRIC_CRS:
    case PJ_TYPE_GEOGRAPHIC_CRS:
    case PJ_TYPE_GEOGRAPHIC_2D_CRS:
    case PJ_TYPE_GEOGRAPHIC_3D_CRS:
    case PJ_TYPE_VERTICAL_CRS:
    case PJ_TYPE_PROJECTED_CRS:
    case PJ_TYPE_COMPOUND_CRS:
    case PJ_TYPE_TEMPORAL_CRS:
    case PJ_TYPE_ENGINEERING_CRS:
    case PJ_TYPE_BOUND_CRS:
    case PJ_TYPE_OTHER_CRS:
      break;
 
    default:
      return false;
  }
 
  d->setPj( QgsProjUtils::unboundCrs( object ) );
 
  if ( !d->hasPj() )
  {
    return d->mIsValid;
  }
  else
  {
    // maybe we can directly grab the auth name and code from the crs
    const QString authName( proj_get_id_auth_name( d->threadLocalProjObject(), 0 ) );
    const QString authCode( proj_get_id_code( d->threadLocalProjObject(), 0 ) );
    if ( !authName.isEmpty() && !authCode.isEmpty() && createFromOgcWmsCrs( u"%1:%2"_s.arg( authName, authCode ) ) )
    {
      return d->mIsValid;
    }
    else
    {
      // Still a valid CRS, just not a known one
      d->mIsValid = true;
      d->mDescription = QString( proj_get_name( d->threadLocalProjObject() ) );
      setMapUnits();
      d->mIsGeographic = testIsGeographic( d->threadLocalProjObject() );
    }
  }
 
  return d->mIsValid;
}
 
QStringList QgsCoordinateReferenceSystem::recentProjections()
{
  QStringList projections;
  const QList<QgsCoordinateReferenceSystem> res = QgsApplication::coordinateReferenceSystemRegistry()->recentCrs();
  projections.reserve( res.size() );
  for ( const QgsCoordinateReferenceSystem &crs : res )
  {
    projections << QString::number( crs.srsid() );
  }
  return projections;
}
 
QList<QgsCoordinateReferenceSystem> QgsCoordinateReferenceSystem::recentCoordinateReferenceSystems()
{
  return QgsApplication::coordinateReferenceSystemRegistry()->recentCrs();
}
 
void QgsCoordinateReferenceSystem::pushRecentCoordinateReferenceSystem( const QgsCoordinateReferenceSystem &crs )
{
  QgsApplication::coordinateReferenceSystemRegistry()->pushRecent( crs );
}
 
void QgsCoordinateReferenceSystem::removeRecentCoordinateReferenceSystem( const QgsCoordinateReferenceSystem &crs )
{
  QgsApplication::coordinateReferenceSystemRegistry()->removeRecent( crs );
}
 
void QgsCoordinateReferenceSystem::clearRecentCoordinateReferenceSystems()
{
  QgsApplication::coordinateReferenceSystemRegistry()->clearRecent();
}
 
void QgsCoordinateReferenceSystem::invalidateCache( bool disableCache )
{
  sSrIdCacheLock()->lockForWrite();
  if ( !sDisableSrIdCache )
  {
    if ( disableCache )
      sDisableSrIdCache = true;
    sSrIdCache()->clear();
  }
  sSrIdCacheLock()->unlock();
 
  sOgcLock()->lockForWrite();
  if ( !sDisableOgcCache )
  {
    if ( disableCache )
      sDisableOgcCache = true;
    sOgcCache()->clear();
  }
  sOgcLock()->unlock();
 
  sProj4CacheLock()->lockForWrite();
  if ( !sDisableProjCache )
  {
    if ( disableCache )
      sDisableProjCache = true;
    sProj4Cache()->clear();
  }
  sProj4CacheLock()->unlock();
 
  sCRSWktLock()->lockForWrite();
  if ( !sDisableWktCache )
  {
    if ( disableCache )
      sDisableWktCache = true;
    sWktCache()->clear();
  }
  sCRSWktLock()->unlock();
 
  sCRSSrsIdLock()->lockForWrite();
  if ( !sDisableSrsIdCache )
  {
    if ( disableCache )
      sDisableSrsIdCache = true;
    sSrsIdCache()->clear();
  }
  sCRSSrsIdLock()->unlock();
 
  sCrsStringLock()->lockForWrite();
  if ( !sDisableStringCache )
  {
    if ( disableCache )
      sDisableStringCache = true;
    sStringCache()->clear();
  }
  sCrsStringLock()->unlock();
}
 
bool QgsCoordinateReferenceSystem::isEarthCrs() const
{
  return celestialBodyName() == "Earth"_L1;
}
 
bool QgsCoordinateReferenceSystem::isSameCelestialBody( const QgsCoordinateReferenceSystem &other ) const
{
  return celestialBodyName() == other.celestialBodyName();
}
 
// invalid < regular < user
bool operator>( const QgsCoordinateReferenceSystem &c1, const QgsCoordinateReferenceSystem &c2 )
{
  if ( c1.d == c2.d )
    return false;
 
  if ( !c1.d->mIsValid && !c2.d->mIsValid )
    return false;
 
  if ( !c1.d->mIsValid && c2.d->mIsValid )
    return false;
 
  if ( c1.d->mIsValid && !c2.d->mIsValid )
    return true;
 
  const bool c1IsUser = c1.d->mSrsId >= Qgis::USER_CRS_START_ID;
  const bool c2IsUser = c2.d->mSrsId >= Qgis::USER_CRS_START_ID;
 
  if ( c1IsUser && !c2IsUser )
    return true;
 
  if ( !c1IsUser && c2IsUser )
    return false;
 
  if ( !c1IsUser && !c2IsUser && !c1.d->mAuthId.isEmpty() && !c2.d->mAuthId.isEmpty() )
  {
    if ( c1.d->mAuthId != c2.d->mAuthId )
      return c1.d->mAuthId > c2.d->mAuthId;
  }
 
  const QString wkt1 = c1.toWkt( Qgis::CrsWktVariant::Preferred );
  const QString wkt2 = c2.toWkt( Qgis::CrsWktVariant::Preferred );
  if ( wkt1 != wkt2 )
    return wkt1 > wkt2;
 
  if ( c1.d->mCoordinateEpoch == c2.d->mCoordinateEpoch )
    return false;
 
  if ( std::isnan( c1.d->mCoordinateEpoch ) && std::isnan( c2.d->mCoordinateEpoch ) )
    return false;
 
  if ( std::isnan( c1.d->mCoordinateEpoch ) && !std::isnan( c2.d->mCoordinateEpoch ) )
    return false;
 
  if ( !std::isnan( c1.d->mCoordinateEpoch ) && std::isnan( c2.d->mCoordinateEpoch ) )
    return true;
 
  return c1.d->mCoordinateEpoch > c2.d->mCoordinateEpoch;
}
 
bool operator<( const QgsCoordinateReferenceSystem &c1, const QgsCoordinateReferenceSystem &c2 )
{
  if ( c1.d == c2.d )
    return false;
 
  if ( !c1.d->mIsValid && !c2.d->mIsValid )
    return false;
 
  if ( c1.d->mIsValid && !c2.d->mIsValid )
    return false;
 
  if ( !c1.d->mIsValid && c2.d->mIsValid )
    return true;
 
  const bool c1IsUser = c1.d->mSrsId >= Qgis::USER_CRS_START_ID;
  const bool c2IsUser = c2.d->mSrsId >= Qgis::USER_CRS_START_ID;
 
  if ( !c1IsUser && c2IsUser )
    return true;
 
  if ( c1IsUser && !c2IsUser )
    return false;
 
  if ( !c1IsUser && !c2IsUser && !c1.d->mAuthId.isEmpty() && !c2.d->mAuthId.isEmpty() )
  {
    if ( c1.d->mAuthId != c2.d->mAuthId )
      return c1.d->mAuthId < c2.d->mAuthId;
  }
 
  const QString wkt1 = c1.toWkt( Qgis::CrsWktVariant::Preferred );
  const QString wkt2 = c2.toWkt( Qgis::CrsWktVariant::Preferred );
  if ( wkt1 != wkt2 )
    return wkt1 < wkt2;
 
  if ( c1.d->mCoordinateEpoch == c2.d->mCoordinateEpoch )
    return false;
 
  if ( std::isnan( c1.d->mCoordinateEpoch ) && std::isnan( c2.d->mCoordinateEpoch ) )
    return false;
 
  if ( !std::isnan( c1.d->mCoordinateEpoch ) && std::isnan( c2.d->mCoordinateEpoch ) )
    return false;
 
  if ( std::isnan( c1.d->mCoordinateEpoch ) && !std::isnan( c2.d->mCoordinateEpoch ) )
    return true;
 
  return c1.d->mCoordinateEpoch < c2.d->mCoordinateEpoch;
}
 
bool operator>=( const QgsCoordinateReferenceSystem &c1, const QgsCoordinateReferenceSystem &c2 )
{
  return !( c1 < c2 );
}
bool operator<=( const QgsCoordinateReferenceSystem &c1, const QgsCoordinateReferenceSystem &c2 )
{
  return !( c1 > c2 );
}