fix(l3): fix missing crs in build_dataset return

j-haacker · Sep 9, 2024 · 8ccf6dc · 8ccf6dc
1 parent 2af9c16
commit 8ccf6dc
Showing 1 changed file with 3 additions and 3 deletions.
diff --git a/cryoswath/l3.py b/cryoswath/l3.py
@@ -309,15 +309,15 @@ def local_closure(roll_iteration):
             del l3_data
             tmp = previously_processed_l3.to_xarray().sortby("x").sortby("y")
             del previously_processed_l3
-            previously_processed_l3 = fill_missing_coords(tmp)
+            previously_processed_l3 = fill_missing_coords(tmp).rio.write_crs(crs)
             tmp.close()
             # save/backup result
             tmp_path = os.path.join(data_path, "tmp", f"{region_id}_l3")
             # try to write new data to file. if anything goes wrong, restore. is this sufficiently safe?
             try:
                 if os.path.isfile(outfilepath):
                     shutil.move(outfilepath, tmp_path)
-                previously_processed_l3.rio.write_crs(crs).to_netcdf(outfilepath)
+                previously_processed_l3.to_netcdf(outfilepath)
             except Exception as err:
                 shutil.move(tmp_path, outfilepath)
                 print("\n")
@@ -326,7 +326,7 @@ def local_closure(roll_iteration):
                     + "message below and continuing. Attempting to save to temporary file.")
                 try:
                     safety_net_tmp_file_path = os.path.join(tmp_path, f"tmp_l3_state__{datetime.datetime.strftime('%dT%H%M%S')}.nc")
-                    previously_processed_l3.rio.write_crs(crs).to_netcdf(safety_net_tmp_file_path)
+                    previously_processed_l3.to_netcdf(safety_net_tmp_file_path)
                 except Exception as err_inner:
                     print("\n")
                     warnings.warn(