WMTS BoundingBox fix and WMTS template adjustment (#1052)

* fixes a bug where the WMTS BoundingBox coordinates flipped incorrectly
if boundingbox element is used, as the coords must match the CRS order

updates WMTS XML template to avoid placing extra new lines in various
places

* fixed typo in import

* updated changelog, removed experimental change in wmts template that is needed for another pending pr

* reformatting

---------

Co-authored-by: vincentsarago <vincent.sarago@gmail.com>

Author: AndrewAnnex

CHANGES.md

@@ -14,6 +14,11 @@
 * update `/healthz` endpoint to return dependencies versions (titiler, rasterio, gdal, ...) (author @scottyhq, https://github.com/developmentseed/titiler/pull/1056)
 * migrate `templates/index.html` to bootstrap5, remove unused css, reuse bs classes (author @PratapVardhan, https://github.com/developmentseed/titiler/pull/1048)
 
+### Misc
+
+* Updated WMTS Capabilities template to avoid inserting extra new lines (author @AndrewAnnex, https://github.com/developmentseed/titiler/pull/1052).
+* Updated WMTS endpoint in titiler.mosaic and titiler.core to return layer bounds in coordinate ordering matching CRS order if WGS84 is not used (author @AndrewAnnex, https://github.com/developmentseed/titiler/pull/1052).
+
 ## 0.19.2 (2024-11-28)
 
 ### Misc

src/titiler/core/titiler/core/factory.py

@@ -26,6 +26,7 @@
 from morecantile import TileMatrixSet
 from morecantile import tms as morecantile_tms
 from morecantile.defaults import TileMatrixSets
+from morecantile.models import crs_axis_inverted
 from pydantic import Field
 from rio_tiler.colormap import ColorMaps
 from rio_tiler.colormap import cmap as default_cmap
@@ -1120,6 +1121,10 @@ def wmts(
             if tms.rasterio_geographic_crs != WGS84_CRS:
                 bbox_crs_type = "BoundingBox"
                 bbox_crs_uri = CRS_to_uri(tms.rasterio_geographic_crs)
+                # WGS88BoundingBox is always xy ordered, but BoundingBox must match the CRS order
+                if crs_axis_inverted(tms.geographic_crs):
+                    # match the bounding box coordinate order to the CRS
+                    bounds = [bounds[1], bounds[0], bounds[3], bounds[2]]
 
             return self.templates.TemplateResponse(
                 request,

src/titiler/core/titiler/core/templates/wmts.xml

@@ -33,7 +33,7 @@
         </ows:Operation>
     </ows:OperationsMetadata>
     <Contents>
-    {% for layer in layers %}
+    {% for layer in layers -%}
         <Layer>
             <ows:Title>{{ layer.title }}</ows:Title>
             <ows:Identifier>{{ layer.name }}</ows:Identifier>
@@ -51,13 +51,13 @@
             </TileMatrixSetLink>
             <ResourceURL format="{{ media_type }}" resourceType="tile" template="{{ layer.tiles_url }}?{{ layer.query_string | escape }}" />
         </Layer>
-    {% endfor %}
+    {%- endfor %}
         <TileMatrixSet>
             <ows:Identifier>{{ tileMatrixSetId }}</ows:Identifier>
             <ows:SupportedCRS>{{ supported_crs }}</ows:SupportedCRS>
-            {% for item in tileMatrix %}
+            {% for item in tileMatrix -%}
             {{ item | safe }}
-            {% endfor %}
+            {%- endfor %}
         </TileMatrixSet>
     </Contents>
     <ServiceMetadataURL xlink:href="{{ request.url | escape }}" />

src/titiler/mosaic/titiler/mosaic/factory.py

@@ -14,6 +14,7 @@
 from geojson_pydantic.geometries import Polygon
 from morecantile import tms as morecantile_tms
 from morecantile.defaults import TileMatrixSets
+from morecantile.models import crs_axis_inverted
 from pydantic import Field
 from rio_tiler.constants import MAX_THREADS, WGS84_CRS
 from rio_tiler.io import BaseReader, MultiBandReader, MultiBaseReader, Reader
@@ -907,9 +908,9 @@ def wmts(
 
             layers = [
                 {
-                    "title": src_path
-                    if isinstance(src_path, str)
-                    else "TiTiler Mosaic",
+                    "title": (
+                        src_path if isinstance(src_path, str) else "TiTiler Mosaic"
+                    ),
                     "name": "default",
                     "tiles_url": tiles_url,
                     "query_string": urlencode(qs, doseq=True) if qs else None,
@@ -922,6 +923,10 @@ def wmts(
             if tms.rasterio_geographic_crs != WGS84_CRS:
                 bbox_crs_type = "BoundingBox"
                 bbox_crs_uri = CRS_to_uri(tms.rasterio_geographic_crs)
+                # WGS88BoundingBox is always xy ordered, but BoundingBox must match the CRS order
+                if crs_axis_inverted(tms.geographic_crs):
+                    # match the bounding box coordinate order to the CRS
+                    bounds = [bounds[1], bounds[0], bounds[3], bounds[2]]
 
             return self.templates.TemplateResponse(
                 request,