diff --git a/py_tools_ds/geo/projection.py b/py_tools_ds/geo/projection.py index 4a98a8d94d71403c4f63e95bfe32222ec72a6145..85a3f620d7d0feb6d1d601e2fdda6259fadccaa9 100644 --- a/py_tools_ds/geo/projection.py +++ b/py_tools_ds/geo/projection.py @@ -85,10 +85,16 @@ def prj_equal(prj1, prj2): if prj1 is None and prj2 is None or prj1 == prj2: return True else: - crs1 = CRS.from_user_input(prj1) - crs2 = CRS.from_user_input(prj2) + from pyproj import __version__ as ver + if ver.startswith('2') and int(ver.split('.')[1]) >= 5: + # CRS.equals was added in pyproj 2.5 which does not exist for Python 2.7 in conda-forge channel + crs1 = CRS.from_user_input(prj1) + crs2 = CRS.from_user_input(prj2) - return crs1.equals(crs2) + return crs1.equals(crs2) + + else: + return get_proj4info(proj=prj1) == get_proj4info(proj=prj2) def isProjectedOrGeographic(prj):