Keep attrs after terrain-correction

bopen · Apr 26, 2022 · 73ae6ba · 73ae6ba
1 parent e1fdbeb
commit 73ae6ba
Showing 1 changed file with 3 additions and 0 deletions.
diff --git a/sarsen/apps.py b/sarsen/apps.py
@@ -187,6 +187,9 @@ def terrain_correction(
 
     logger.info("save output")
 
+    geocoded.attrs.update(beta_nought.attrs)
+    geocoded.x.attrs.update(dem_raster.x.attrs)
+    geocoded.y.attrs.update(dem_raster.y.attrs)
     geocoded.rio.set_crs(dem_raster.rio.crs)
     geocoded.rio.to_raster(
         output_urlpath,