morecantile 5.2.3__py3-none-any.whl → 5.3.1__py3-none-any.whl

morecantile/__init__.py

@@ -8,7 +8,7 @@ Refs:
 
 """
 
-__version__ = "5.2.3"
+__version__ = "5.3.1"
 
 from .commons import BoundingBox, Coords, Tile  # noqa
 from .defaults import TileMatrixSets, tms  # noqa

morecantile/models.py

@@ -31,6 +31,7 @@ from morecantile.utils import (
     _parse_tile_arg,
     bbox_to_feature,
     check_quadkey_support,
+    lons_contain_antimeridian,
     meters_per_unit,
     point_in_bbox,
     to_rasterio_crs,
@@ -660,6 +661,7 @@ class TileMatrixSet(BaseModel, arbitrary_types_allowed=True):
         ordered_axes: Optional[List[str]] = None,
         geographic_crs: pyproj.CRS = WGS84_CRS,
         screen_pixel_size: float = 0.28e-3,
+        decimation_base: int = 2,
         **kwargs: Any,
     ):
         """
@@ -692,8 +694,12 @@ class TileMatrixSet(BaseModel, arbitrary_types_allowed=True):
             Tile Matrix Set identifier
         geographic_crs: pyproj.CRS
             Geographic (lat,lon) coordinate reference system (default is EPSG:4326)
+        ordered_axes: list of str, optional
+            Override Axis order (e.g `["N", "S"]`) else default to CRS's metadata
         screen_pixel_size: float, optional
             Rendering pixel size. 0.28 mm was the actual pixel size of a common display from 2005 and considered as standard by OGC.
+        decimation_base: int, optional
+            How tiles are divided at each zoom level (default is 2). Must be greater than 1.
         kwargs: Any
             Attributes to forward to the TileMatrixSet
 
@@ -713,6 +719,11 @@ class TileMatrixSet(BaseModel, arbitrary_types_allowed=True):
             transform = pyproj.Transformer.from_crs(extent_crs, crs, always_xy=True)
             extent = transform.transform_bounds(*extent, densify_pts=21)
 
+        if decimation_base <= 1:
+            raise ValueError(
+                "Custom TileMatrixSet requires a decimation base that is greater than 1."
+            )
+
         bbox = BoundingBox(*extent)
         x_origin = bbox.left if not is_inverted else bbox.top
         y_origin = bbox.top if not is_inverted else bbox.left
@@ -723,8 +734,10 @@ class TileMatrixSet(BaseModel, arbitrary_types_allowed=True):
         tile_matrices: List[TileMatrix] = []
         for zoom in range(minzoom, maxzoom + 1):
             res = max(
-                width / (tile_width * matrix_scale[0]) / 2.0**zoom,
-                height / (tile_height * matrix_scale[1]) / 2.0**zoom,
+                width / (tile_width * matrix_scale[0]) / float(decimation_base) ** zoom,
+                height
+                / (tile_height * matrix_scale[1])
+                / float(decimation_base) ** zoom,
             )
             tile_matrices.append(
                 TileMatrix(
@@ -735,8 +748,8 @@ class TileMatrixSet(BaseModel, arbitrary_types_allowed=True):
                         "pointOfOrigin": [x_origin, y_origin],
                         "tileWidth": tile_width,
                         "tileHeight": tile_height,
-                        "matrixWidth": matrix_scale[0] * 2**zoom,
-                        "matrixHeight": matrix_scale[1] * 2**zoom,
+                        "matrixWidth": matrix_scale[0] * decimation_base**zoom,
+                        "matrixHeight": matrix_scale[1] * decimation_base**zoom,
                     }
                 )
             )
@@ -1240,9 +1253,11 @@ class TileMatrixSet(BaseModel, arbitrary_types_allowed=True):
 
         for w, s, e, n in bboxes:
             # Clamp bounding values.
-            w = max(self.bbox.left, w)
+            ws_contain_180th = lons_contain_antimeridian(w, self.bbox.left)
+            es_contain_180th = lons_contain_antimeridian(e, self.bbox.right)
+            w = min(self.bbox.left, w) if ws_contain_180th else max(self.bbox.left, w)
             s = max(self.bbox.bottom, s)
-            e = min(self.bbox.right, e)
+            e = max(self.bbox.right, e) if es_contain_180th else min(self.bbox.right, e)
             n = min(self.bbox.top, n)
 
             for z in zooms:

morecantile/utils.py

@@ -40,6 +40,30 @@ def _parse_tile_arg(*args) -> Tile:
         )
 
 
+def lons_contain_antimeridian(lon1: float, lon2: float) -> bool:
+    """
+    Check if the antimeridian (180th meridian) is between two longitude points
+
+    Parameters
+    ----------
+    lon1: float
+        The first longitude.
+    lon2: float
+        The second longitude
+
+    Returns
+    -------
+    A bool representing whether two longs contain the 180th meridian.
+    """
+    lon1_clipped = max(-180.0, min(lon1, 180))
+    lon2_clipped = max(-180.0, min(lon2, 180))
+    lon1_converted = (lon1_clipped + 360) % 360
+    lon2_converted = (lon2_clipped + 360) % 360
+    ws = [lon1_converted, lon2_converted]
+    sorted(ws)
+    return ws[0] < 180 < ws[1]
+
+
 def meters_per_unit(crs: CRS) -> float:
     """
     Coefficient to convert the coordinate reference system (CRS)

{morecantile-5.2.3.dist-info → morecantile-5.3.1.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: morecantile
-Version: 5.2.3
+Version: 5.3.1
 Summary: Construct and use map tile grids (a.k.a TileMatrixSet / TMS).
 Keywords: GIS,TMS,TileMatrixSet,Map Tile
 Author-email: Vincent Sarago <vincent@developmentseed.com>
@@ -18,6 +18,7 @@ Requires-Dist: attrs
 Requires-Dist: pyproj~=3.1
 Requires-Dist: pydantic~=2.0
 Requires-Dist: pre-commit ; extra == "dev"
+Requires-Dist: bump-my-version ; extra == "dev"
 Requires-Dist: mkdocs ; extra == "docs"
 Requires-Dist: mkdocs-material ; extra == "docs"
 Requires-Dist: pygments ; extra == "docs"

{morecantile-5.2.3.dist-info → morecantile-5.3.1.dist-info}/RECORD

@@ -1,10 +1,10 @@
-morecantile/__init__.py,sha256=wcwgjah4uYI6ix_t1YBpqqfE2dkDhfI18E2cuy4qYMg,436
+morecantile/__init__.py,sha256=xp0X6IQ8GGTQkJFJ9m2nrxtfkWW4njt-9MjFRNVKgUA,436
 morecantile/commons.py,sha256=nOj5f5K5EbDZiC0QTK2f14Mpkz6YIZyCaaUI0S9aWIw,1015
 morecantile/defaults.py,sha256=9lvHmZ9WM_J5uNCS9VWEy8O5_x2dA2Spvk4H5aHv6UQ,1804
 morecantile/errors.py,sha256=rhtdpNglfEz5nC8I-BJKUr_gkOwAPwVi1vhLFJ2StYw,907
-morecantile/models.py,sha256=LssuG3SB5sNSn5JkT_FlRo0hi6U3PCWuQ63jOxo9zms,50806
+morecantile/models.py,sha256=mRCIDQ2VjMjNigF9mUiPy2oLehGmSSwD3VdVO_OH0Wc,51649
 morecantile/py.typed,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-morecantile/utils.py,sha256=J37ixKq_lJC23bqJSAf01mfS-ZEHizvt2uqVE47dPRE,3364
+morecantile/utils.py,sha256=-JUTcrwhT_TEWAPWV8YQBQJoYU2HfVkqQEJaggtnKJg,4020
 morecantile/data/CDB1GlobalGrid.json,sha256=VDc2ukAWtTQeCdOG8YMqsGO-D7eSzSCN1TEaBn2qLZI,36115
 morecantile/data/CanadianNAD83_LCC.json,sha256=RuR9z5MjxKzU4KE-CutrzjpWBz9c-XXykKAxcD83veU,7267
 morecantile/data/EuropeanETRS89_LAEAQuad.json,sha256=pR2Q7eAlLRG74SRdzYhXj19bf5qBvHrE6kzJ5cvko1Q,4679
@@ -21,8 +21,8 @@ morecantile/data/WorldCRS84Quad.json,sha256=l0Wf2faYwzRwO0mG-Tea1Vydpi7V1Oel5L9l
 morecantile/data/WorldMercatorWGS84Quad.json,sha256=JaqoBSu5qVJmL-J7oSVU1etP2OPWx_aVaTc68dGX0Ec,7001
 morecantile/scripts/__init__.py,sha256=-CJncfgWDnSZ8au-SJtgX-OFgCddlf7___d91qQcqQM,23
 morecantile/scripts/cli.py,sha256=zSE15udt2Mb8IR9XB01gm8k_L65ON4_WAkpcJmQazlI,16859
-morecantile-5.2.3.dist-info/entry_points.txt,sha256=keoXuYgnX-mdrGWwXIHqEkSZjQK_iAmtmdrtZ3sCT24,59
-morecantile-5.2.3.dist-info/LICENSE,sha256=18IxFIta7rF_RcVSIgLUUd_alyfQ9bGoZKQm5vOarGU,1073
-morecantile-5.2.3.dist-info/WHEEL,sha256=EZbGkh7Ie4PoZfRQ8I0ZuP9VklN_TvcZ6DSE5Uar4z4,81
-morecantile-5.2.3.dist-info/METADATA,sha256=XlcFVU9Hd8h14s_eho6zY5QrECpXHKgk9eQm3czGfK0,6890
-morecantile-5.2.3.dist-info/RECORD,,
+morecantile-5.3.1.dist-info/entry_points.txt,sha256=keoXuYgnX-mdrGWwXIHqEkSZjQK_iAmtmdrtZ3sCT24,59
+morecantile-5.3.1.dist-info/LICENSE,sha256=18IxFIta7rF_RcVSIgLUUd_alyfQ9bGoZKQm5vOarGU,1073
+morecantile-5.3.1.dist-info/WHEEL,sha256=EZbGkh7Ie4PoZfRQ8I0ZuP9VklN_TvcZ6DSE5Uar4z4,81
+morecantile-5.3.1.dist-info/METADATA,sha256=1HUu9x4xTg7y2cmoqGghUW127eupnemilcY7gK2UYlE,6938
+morecantile-5.3.1.dist-info/RECORD,,