def open_dataset(file_path): ds = xr.open_dataset(file_path, mask_and_scale=False, drop_variables='crs') ds.attrs['crs'] = datacube.utils.geometry.CRS('EPSG:32754') return ds