datacube.utils.geometry.assign_crs#

datacube.utils.geometry.assign_crs(xx, crs=None, crs_coord_name='spatial_ref')[source]#

Assign CRS for a non-georegistered array or dataset.

Returns a new object with CRS information populated.

Can also be called without crs argument on data that already has CRS information but not in the format used by datacube, in this case CRS metadata will be restructured into a shape used by datacube. This format allows for better propagation of CRS information through various computations.

xx = datacube.utils.geometry.assign_crs(xr.open_rasterio("some-file.tif"))
print(xx.geobox)
print(xx.astype('float32').geobox)

Parameters

xx (Union[DataArray, Dataset]) – Dataset or DataArray
crs (Union[str, CRS, CRS, Dict[str, Any], None]) – CRS to assign, if omitted try to guess from attributes
crs_coord_name (str) – how to name crs corodinate (defaults to spatial_ref)

Return type

Union[Dataset, DataArray]