How to rasterize georeferenced images from xarray inference dataset
import xarray as xr
inf = xr.open_dataset("merged_filtered_inference.nc")
dates = inf.time_interpolated_negrar.values
for cantina in ["negrar","laguar"]:
for phase in ["target","predicted"]:
for date in dates[[0,9,18]]:
inf[f"ndvi_{phase}_{cantina}"].sel({f"time_interpolated_{cantina}":date}).rio.set_spatial_dims(x_dim=f"lon_0_{cantina}",y_dim=f"lat_0_{cantina}").rio.to_raster(f"rasters/{cantina}_{str(date)[:10]}_{phase}.tif")