changed signs in offset calculation to be consistent with QGIS

EOA-team · Dec 1, 2023 · c9a32a0 · c9a32a0
1 parent 31e4cef
commit c9a32a0
Showing 1 changed file with 4 additions and 4 deletions.
diff --git a/eodal/core/band.py b/eodal/core/band.py
@@ -2230,24 +2230,24 @@ def scale_data(
         scale, offset = self.scale, self.offset
         if self.is_masked_array:
             if pixel_values_to_ignore is None:
-                scaled_array = scale * self.values.data - offset
+                scaled_array = scale * self.values.data + offset
             else:
                 scaled_array = self.values.data.copy().astype(float)
                 scaled_array[~np.isin(scaled_array, pixel_values_to_ignore)] = scale * \
                     scaled_array[~np.isin(scaled_array, pixel_values_to_ignore)] \
-                    - offset
+                    + offset
             # reuse fill value
             fill_value = self.values.fill_value
             scaled_array = np.ma.MaskedArray(
                 data=scaled_array, mask=self.values.mask, fill_value=fill_value)
         elif self.is_ndarray:
             if pixel_values_to_ignore is None:
-                scaled_array = scale * self.values - offset
+                scaled_array = scale * self.values + offset
             else:
                 scaled_array = self.values.copy().astype(float)
                 scaled_array[~np.isin(scaled_array, pixel_values_to_ignore)] = scale * \
                     scaled_array[~np.isin(scaled_array, pixel_values_to_ignore)] \
-                    - offset
+                    + offset
         elif self.is_zarr:
             raise NotImplementedError()