fix rasterize mesh when there is no data

qgis · Jun 10, 2021 · a1ec458 · a1ec458
1 parent 5864629
commit a1ec458
Showing 1 changed file with 17 additions and 13 deletions.
diff --git a/src/analysis/processing/qgsalgorithmexportmesh.cpp b/src/analysis/processing/qgsalgorithmexportmesh.cpp
@@ -834,20 +834,24 @@ QVariantMap QgsMeshRasterizeAlgorithm::processAlgorithm( const QVariantMap &para
     if ( feedback )
       QObject::connect( &rasterBlockFeedBack, &QgsFeedback::canceled, feedback, &QgsFeedback::cancel );
 
-    QgsRasterBlock *block = QgsMeshUtils::exportRasterBlock(
-                              mTriangularMesh,
-                              dataGroup.datasetValues,
-                              dataGroup.activeFaces,
-                              dataGroup.metadata.dataType(),
-                              mTransform,
-                              pixelSize,
-                              extent,
-                              &rasterBlockFeedBack );
-
-
+    if ( dataGroup.datasetValues.isValid() )
+    {
+      QgsRasterBlock *block = QgsMeshUtils::exportRasterBlock(
+                                mTriangularMesh,
+                                dataGroup.datasetValues,
+                                dataGroup.activeFaces,
+                                dataGroup.metadata.dataType(),
+                                mTransform,
+                                pixelSize,
+                                extent,
+                                &rasterBlockFeedBack );
+
+      rasterDataProvider->writeBlock( block, i + 1 );
+      rasterDataProvider->setNoDataValue( i + 1, block->noDataValue() );
+    }
+    else
+      rasterDataProvider->setNoDataValue( i + 1, std::numeric_limits<double>::quiet_NaN() );
 
-    rasterDataProvider->writeBlock( block, i + 1 );
-    rasterDataProvider->setNoDataValue( i + 1, block->noDataValue() );
     if ( feedback )
     {
       if ( feedback->isCanceled() )