From 62e46b9b0b0166b690a02d5f4237290cbbd20616 Mon Sep 17 00:00:00 2001
From: Aurore Dupuis <aurore.dupuis@gmail.com>
Date: Tue, 21 Jul 2020 08:54:24 +0000
Subject: [PATCH] Merge branch '69-add-compression-to-syntheses' into 'develop'

Resolve "add compression to syntheses"

Closes #69

See merge request remote_modules/let-it-snow!61

(cherry picked from commit e34015cc05eb867f9088cb7e2ecf420683db7b3d)

3c3b599a add compression to snow annual map
bd079a61 add compression to SMOD, SOD
02d3832b Update compute_NOBS.py
---
 python/s2snow/compute_NOBS.py     | 1 +
 python/s2snow/compute_SOD_SMOD.py | 1 +
 python/s2snow/snow_annual_map.py  | 2 +-
 3 files changed, 3 insertions(+), 1 deletion(-)

diff --git a/python/s2snow/compute_NOBS.py b/python/s2snow/compute_NOBS.py
index d8140c87..0a1245aa 100644
--- a/python/s2snow/compute_NOBS.py
+++ b/python/s2snow/compute_NOBS.py
@@ -57,6 +57,7 @@ def compute_NOBS(input_file, output_file=None):
         profile.update(
             dtype=rasterio.uint16,
             driver='GTiff',
+            compress='deflate',
             count=1)
 
         with rasterio.open(output_file, 'w', **profile) as dst:
diff --git a/python/s2snow/compute_SOD_SMOD.py b/python/s2snow/compute_SOD_SMOD.py
index 42bb4e14..0671bc28 100644
--- a/python/s2snow/compute_SOD_SMOD.py
+++ b/python/s2snow/compute_SOD_SMOD.py
@@ -71,6 +71,7 @@ def compute_SOD_SMOD(input_file, sod_file=None, smod_file=None):
         profile = src.profile
         profile.update(
             dtype=rasterio.uint16,
+            compress='deflate',
             count=1)
 
         with rasterio.open(smod_file, 'w', **profile) as dst:
diff --git a/python/s2snow/snow_annual_map.py b/python/s2snow/snow_annual_map.py
index d868dbde..bd93fcb6 100644
--- a/python/s2snow/snow_annual_map.py
+++ b/python/s2snow/snow_annual_map.py
@@ -102,7 +102,7 @@ def compute_SCD(binary_snow_mask_list, multitemp_cloud_vrt, input_dates_filename
     expression = "+".join(["im1b" + str(i) for i in band_index])
     logging.debug("expression: {}".format(expression))
     bandMathApp = band_math([gapfilled_timeserie],
-                            snow_coverage_duration,
+                            snow_coverage_duration  + "?&gdal:co:COMPRESS=DEFLATE",
                             expression,
                             ram,
                             otb.ImagePixelType_uint16)
-- 
GitLab