1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
|
"""
This module contains helper methods for rasterio.crs.CRS.
"""
from typing import Any
import rasterio
import rasterio.crs
from packaging import version
from pyproj import CRS
from rasterio.errors import CRSError
def crs_from_user_input(crs_input: Any) -> rasterio.crs.CRS:
"""
Return a rasterio.crs.CRS from user input.
This is to deal with change in rasterio.crs.CRS
as well as to assist in the transition between GDAL 2/3.
Parameters
----------
crs_input: Any
Input to create a CRS.
Returns
-------
rasterio.crs.CRS
"""
if isinstance(crs_input, rasterio.crs.CRS):
return crs_input
try:
# old versions of opendatacube CRS
crs_input = crs_input.wkt
except AttributeError:
pass
try:
return rasterio.crs.CRS.from_user_input(crs_input)
except CRSError:
pass
# use pyproj for edge cases
crs = CRS.from_user_input(crs_input)
if version.parse(rasterio.__gdal_version__) > version.parse("3.0.0"):
return rasterio.crs.CRS.from_wkt(crs.to_wkt())
return rasterio.crs.CRS.from_wkt(crs.to_wkt("WKT1_GDAL"))
|
