Fix NIR reflectance computation to use XArray

pytroll · Jun 5, 2018 · d7bc3e6 · d7bc3e6
1 parent 14114a5
commit d7bc3e6
Showing 1 changed file with 4 additions and 2 deletions.
diff --git a/satpy/composites/__init__.py b/satpy/composites/__init__.py
@@ -469,9 +469,11 @@ def __call__(self, projectables, optional_datasets=None, **info):
         """
         self._init_refl3x(projectables)
         _nir, _ = projectables
-        proj = Dataset(self._get_reflectance(projectables, optional_datasets) * 100, **_nir.info)
+        refl = self._get_reflectance(projectables, optional_datasets) * 100
+        proj = xr.DataArray(refl.filled(np.nan), dims=_nir.dims,
+                            coords=_nir.coords, attrs=_nir.attrs)
 
-        proj.info['units'] = '%'
+        proj.attrs['units'] = '%'
         self.apply_modifier_info(_nir, proj)
 
         return proj