diff --git a/python/s2snow/compute_NOBS.py b/python/s2snow/compute_NOBS.py
index d8140c87233cf47989adac06fc756fc28f8ea633..0a1245aaaa768ee18f361677e8e8d96309687ed5 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: