Merge pull request #1484 from OldLipe/feat/reclassify-tile

Add support for multiple tiles in `sits_reclassify()`

Author: gilbertocamara

R/api_gdal.R

@@ -142,6 +142,11 @@
     file <- file[[1L]]
     # Convert to gdal data type
     data_type <- .gdal_data_type[[data_type]]
+    co <- .conf("gdal_presets", "image", "co")
+    if (data_type == "Float64") {
+        co <- setdiff(co, c("COMPRESS=LZW", "PREDICTOR=2"))
+        co <- c("COMPRESS=ZSTD", "PREDICTOR=3", co)
+    }
     # Output file
     file <- .try(
         {
@@ -165,7 +170,7 @@
                         .xmin(bbox), .ymax(bbox), .xmax(bbox), .ymin(bbox)
                     ),
                     "-a_nodata" = miss_value,
-                    "-co" = .conf("gdal_presets", "image", "co")
+                    "-co" = co
                 ),
                 quiet = TRUE
             )

R/api_reclassify.R

@@ -84,18 +84,20 @@
             nlayers = 1L, miss_value = .miss_value(band_conf),
             data_type = .data_type(band_conf)
         )
-        # Copy values from mask cube into mask template
+        # Get mask files
+        paths <-  purrr::map_chr(mask[["file_info"]], .fi_paths)
+        # Copy values into template block
         .gdal_merge_into(
             file = mask_block_file,
-            base_files = .fi_paths(.fi(mask)), multicores = 1L
+            base_files = paths, multicores = 1L
         )
         # Build a new tile for mask based on template
         mask_tile <- .tile_derived_from_file(
             file = mask_block_file,
             band = "class",
             base_tile = .tile(mask),
             derived_class = "class_cube",
-            update_bbox = FALSE
+            update_bbox = TRUE
         )
         # Read and preprocess values
         values <- .tile_read_block(