apply_dimension(dim=bands) => do not filter out all nodata tiles, but…

… replace them with 'empty' tile. This avoids issues in subsequent merge_cubes with spatial cube or linear interpolation over time. WorldCereal/worldcereal-classification#141
Open-EO · Sep 19, 2024 · a5a8328 · a5a8328
1 parent a201d12
commit a5a8328
Showing 1 changed file with 6 additions and 2 deletions.
diff --git a/openeo-geotrellis/src/main/scala/org/openeo/geotrellis/OpenEOProcesses.scala b/openeo-geotrellis/src/main/scala/org/openeo/geotrellis/OpenEOProcesses.scala
@@ -561,11 +561,15 @@ class OpenEOProcesses extends Serializable {
     return ContextRDD(new org.apache.spark.rdd.PairRDDFunctions[K,MultibandTile](datacube).mapValues(tile => {
       if (!tile.isInstanceOf[EmptyMultibandTile]) {
         val resultTiles = function(tile.bands)
-        MultibandTile(resultTiles)
+        if(resultTiles.exists(!_.isNoDataTile)){
+          MultibandTile(resultTiles)
+        }else{
+          new EmptyMultibandTile(tile.cols,tile.rows,tile.cellType,tile.bandCount)
+        }
       }else{
         tile
       }
-    }).filter(_._2.bands.exists(!_.isNoDataTile)),datacube.metadata.copy(cellType = scriptBuilder.getOutputCellType()))
+    }),datacube.metadata.copy(cellType = scriptBuilder.getOutputCellType()))
   }
 
   def filterEmptyTile[K:ClassTag](datacube:MultibandTileLayerRDD[K]): RDD[(K, MultibandTile)] with Metadata[TileLayerMetadata[K]]={