use projection string from MDAL provider in Mesh Layer

qgis · Sep 24, 2018 · 7ad1b71 · 7ad1b71
1 parent 7a85fb0
commit 7ad1b71
Show file tree

Hide file tree

Showing 2 changed files with 9 additions and 1 deletion.
diff --git a/src/providers/mdal/qgsmdalprovider.cpp b/src/providers/mdal/qgsmdalprovider.cpp
@@ -44,14 +44,20 @@ QString QgsMdalProvider::description() const
 
 QgsCoordinateReferenceSystem QgsMdalProvider::crs() const
 {
-  return QgsCoordinateReferenceSystem();
+  return mCrs;
 }
 
 QgsMdalProvider::QgsMdalProvider( const QString &uri, const ProviderOptions &options )
   : QgsMeshDataProvider( uri, options )
 {
   QByteArray curi = uri.toAscii();
   mMeshH = MDAL_LoadMesh( curi.constData() );
+  if ( mMeshH )
+  {
+    const QString proj = MDAL_M_projection( mMeshH );
+    if ( !proj.isEmpty() )
+      mCrs.createFromString( proj );
+  }
 }
 
 QgsMdalProvider::~QgsMdalProvider()

diff --git a/src/providers/mdal/qgsmdalprovider.h b/src/providers/mdal/qgsmdalprovider.h
@@ -26,6 +26,7 @@
 
 class QMutex;
 class QgsCoordinateTransform;
+class QgsCoordinateReferenceSystem;
 
 /**
   \brief Data provider for MDAL layers.
@@ -69,6 +70,7 @@ class QgsMdalProvider : public QgsMeshDataProvider
   private:
     MeshH mMeshH;
     QStringList mExtraDatasetUris;
+    QgsCoordinateReferenceSystem mCrs;
 };
 
 #endif //QGSMDALPROVIDER_H