switch from rasterio/gdal to pyproj

morecantile/models.py

@@ -7,6 +7,7 @@
 from pydantic import AnyHttpUrl, BaseModel, Field, PrivateAttr, validator
 from pyproj import CRS, Transformer
 from pyproj.enums import WktVersion
+from pyproj.exceptions import ProjError
 
 from .commons import BoundingBox, Coords, Tile
 from .errors import (
@@ -28,7 +29,6 @@
     from rasterio.crs import CRS as rasterioCRS
     from rasterio.env import GDALVersion
 except ModuleNotFoundError:
-    rasterio = None
     rasterioCRS = None
     GDALVersion = None
 
@@ -149,12 +149,21 @@ def __init__(self, **data):
         """Create PyProj transforms and check if TileMatrixSet supports quadkeys."""
         super().__init__(**data)
         self._is_quadtree = check_quadkey_support(self.tileMatrix)
-        self._to_wgs84 = Transformer.from_crs(
-            self.supportedCRS, WGS84_CRS, always_xy=True
-        )
-        self._from_wgs84 = Transformer.from_crs(
-            WGS84_CRS, self.supportedCRS, always_xy=True
-        )
+        try:
+            self._to_wgs84 = Transformer.from_crs(
+                self.supportedCRS, WGS84_CRS, always_xy=True
+            )
+            self._from_wgs84 = Transformer.from_crs(
+                WGS84_CRS, self.supportedCRS, always_xy=True
+            )
+        except ProjError:
+            warnings.warn(
+                "Could not create coordinate Transformer from input CRS to WGS84"
+                "tms methods might not be available.",
+                UserWarning,
+            )
+            self._to_wgs84 = None
+            self._from_wgs84 = None
 
     @validator("tileMatrix")
     def sort_tile_matrices(cls, v):

tests/test_models.py

@@ -217,8 +217,6 @@ def test_nztm_quad_is_quad():
 def test_nztm_quad_scales():
     nztm_tms = morecantile.tms.get("NZTM2000Quad")
     google_tms = morecantile.tms.get("WebMercatorQuad")
-    print(dir(google_tms))
-
     for z in range(2, nztm_tms.maxzoom + 2):
         assert (
             round(
@@ -278,7 +276,10 @@ def test_schema():
         "+proj=stere +lat_0=90 +lon_0=0 +k=2 +x_0=0 +y_0=0 +R=3396190 +units=m +no_defs"
     )
     extent = [-13584760.000, -13585240.000, 13585240.000, 13584760.000]
-    tms = morecantile.TileMatrixSet.custom(extent, crs, identifier="MarsNPolek2MOLA5k")
+    with pytest.warns(UserWarning):
+        tms = morecantile.TileMatrixSet.custom(
+            extent, crs, identifier="MarsNPolek2MOLA5k"
+        )
     assert tms.schema()
     assert tms.schema_json()
     assert tms.dict(exclude_none=True)