add geotransform and crs

satpy/readers/emit_l1b.py

@@ -22,6 +22,8 @@
 from datetime import datetime
 
 import numpy as np
+from pyproj import CRS
+from affine import Affine
 from pyresample.geometry import SwathDefinition
 from satpy.readers.netcdf_utils import NetCDF4FileHandler
 from satpy.utils import get_legacy_chunk_size
@@ -76,6 +78,16 @@ def sensor(self):
         """Get sensor."""
         return self["/attr/instrument"]
 
+    @property
+    def spatial_ref(self):
+        """Get spatial_ref."""
+        return self["/attr/spatial_ref"]
+
+    @property
+    def geotransform(self):
+        """Get geotransform."""
+        return self["/attr/geotransform"]
+
     @property
     def attrs(self):
         """Return attributes."""
@@ -153,6 +165,10 @@ def get_dataset(self, dataset_id, ds_info):
         self.get_lonlats()
         dataset.attrs['area'] = self.area
 
+        # write crs and geotransform
+        dataset.rio.write_transform(Affine.from_gdal(*self.geotransform), inplace=True)
+        dataset.rio.write_crs(CRS.from_wkt(self.spatial_ref), inplace=True)
+
         return dataset
 
     def get_lonlats(self):