Fixed georeferencing bug.

preprocessing/calculate_ndvi.py

@@ -199,6 +199,7 @@ def download_ndvi_imagery(config_file: str, interp_chunk: dict = {'x': 400, 'y':
         else:
             file_chunksize = (1, ndvi.dims['y'], ndvi.dims['x'])
         
+        ndvi.attrs['crs'] = ndvi.rio.crs.to_string()
         write_job = ndvi.to_netcdf(ndvi_cube_path, encoding = {"NDVI": {"dtype": "u1", "_FillValue": 0, "chunksizes": file_chunksize}}, compute = False)
         write_job = write_job.persist()
         

preprocessing/lib_era5_land_pixel.py

@@ -856,6 +856,7 @@ def era5Land_daily_to_yearly_pixel(list_era5land_files: List[str], output_file:
         final_weather_ds.rio.write_crs(target_crs, inplace = True)
         final_weather_ds.rio.write_transform(inplace = True)
         final_weather_ds['spatial_ref'] = spatial_ref
+        final_weather_ds.attrs['crs'] = final_weather_ds.rio.crs.to_string()
         final_weather_ds = final_weather_ds.set_coords('spatial_ref')
 
         # Save empty output