Add .cf.crs attribute

Author: dcherian

cf_xarray/accessor.py

@@ -8,6 +8,7 @@
 from collections.abc import Hashable, Iterable, Mapping, MutableMapping, Sequence
 from datetime import datetime
 from typing import (
+    TYPE_CHECKING,
     Any,
     Callable,
     Literal,
@@ -57,6 +58,9 @@
     parse_cf_standard_name_table,
 )
 
+if TYPE_CHECKING:
+    from pyproj import CRS
+
 FlagParam = namedtuple("FlagParam", ["flag_mask", "flag_value"])
 
 #: Classes wrapped by cf_xarray.
@@ -2342,6 +2346,19 @@ def __getitem__(self, key: Hashable | Iterable[Hashable]) -> DataArray | Dataset
         """
         return _getitem(self, key)
 
+    @property
+    def crs(self) -> CRS:
+        from pyproj import CRS
+
+        try:
+            gm = self["grid_mapping"]
+        except KeyError as e:
+            raise ValueError(
+                "CRS is constructed from the 'grid_mapping' variable, but that could not be detected."
+            ) from e
+
+        return CRS.from_cf(gm.attrs)
+
     @property
     def formula_terms(self) -> dict[Hashable, dict[str, str]]:  # numpydoc ignore=SS06
         """