Implement native SQL Alchemy MVT (#2022)

* Implement SQLAlchemy native MVT query

Fix typo

Rename base tile provider

Fix Tile HTML highlighting

Fix missing tile rename

Update CRS conversion

Change fallback for vector tile id

Undo rename tile provider

Reduce diff

Cleanup mvt-postgres provider

Re-add ST_CurveToLine

97ba024

* Add ST_CurveToLine for output geom

* Set id field for vector files

* Add MVT Postgres Provider test

Fix Flake8

* Fix indentation in MVT-Postgres docs

* Fix get_metadata()

- Fix get_metadata function
- Add unit test for get_metadata

* Fix MVT feature properties

- Fix MVT field selection
- Add test for vector tile properties

* Fix flake8

* Add storage_crs information to docs

* Add test_postgresql_mvt_provider to GitHub CI

* Bump GitHub Action Postgis version

Bump GitHub Postgis engine to >3

* Filter SQL connection options

Prevent passing dict to SQL connection which throws an error because it cannot be cached

* Use geoalchemy2 functions

* Amend WorldCRS84Quad Tile bbox

Align with  WorldCRS84Quad spec using #2042

* Support WorldCRS84Quad in default tile html

Add minimal support stretching EPSG3857 TileLayer  to EPSG4326. Really should be a CRS84 based layer like:

L.tileLayer.wms("https://gibs.earthdata.nasa.gov/wms/epsg4326/best/wms.cgi", {
    layers: "BlueMarble_ShadedRelief_Bathymetry",
    attribution: "NASA GIBS",
  })
  .addTo(map);

Author: webb-ben

.github/workflows/main.yml

@@ -24,6 +24,15 @@ jobs:
           - python-version: '3.10'
     env:
       PYGEOAPI_CONFIG: "$(pwd)/pygeoapi-config.yml"
+
+    services:
+      postgres:
+        image: postgis/postgis:14-3.2
+        ports:
+          - 5432:5432
+        env:
+          POSTGRES_DB: test
+          POSTGRES_PASSWORD: ${{ secrets.DatabasePassword || 'postgres' }}
 
     steps:
     - name: Pre-pull Docker Images
@@ -32,7 +41,6 @@ jobs:
         docker pull appropriate/curl:latest &
         docker pull elasticsearch:8.17.0 &
         docker pull opensearchproject/opensearch:2.18.0 &
-        docker pull mdillon/postgis:latest &
         docker pull mongo:8.0.4 &
         docker pull ghcr.io/cgs-earth/sensorthings-action:0.1.0 &
         docker pull postgis/postgis:14-3.2 &
@@ -58,12 +66,6 @@ jobs:
         sudo sysctl -w vm.swappiness=1
         sudo sysctl -w fs.file-max=262144
         sudo sysctl -w vm.max_map_count=262144
-    - name: Install and run PostgreSQL/PostGIS 📦
-      uses: huaxk/postgis-action@v1
-      with:
-        postgresql password: ${{ secrets.DatabasePassword || 'postgres' }}
-        postgresql db: 'test'
-
     - name: "Install and run MySQL 📦"
       uses: mirromutth/mysql-action@v1.1
       with:
@@ -164,6 +166,7 @@ jobs:
         pytest tests/test_oracle_provider.py
         pytest tests/test_parquet_provider.py
         pytest tests/test_postgresql_provider.py
+        pytest tests/test_postgresql_mvt_provider.py
         pytest tests/test_mysql_provider.py
         pytest tests/test_rasterio_provider.py
         pytest tests/test_sensorthings_edr_provider.py 

docs/source/data-publishing/ogcapi-tiles.rst

@@ -140,9 +140,6 @@ MVT-postgresql
 .. note::
    Must have PostGIS installed with protobuf-c support 
 
-.. note::
-   Geometry must be using EPSG:4326
-
 This provider gives support to serving tiles generated using `PostgreSQL <https://www.postgresql.org/>`_ with `PostGIS <https://postgis.net/>`_.
 The tiles are rendered on-the-fly using `ST_AsMVT <https://postgis.net/docs/ST_AsMVT.html>`_ and related methods.
 
@@ -152,24 +149,28 @@ This code block shows how to configure pygeoapi to render Mapbox vector tiles fr
 
    providers:
        - type: tile
-           name: MVT-postgresql
-           data:
-               host: 127.0.0.1
-               port: 3010 # Default 5432 if not provided
-               dbname: test
-               user: postgres
-               password: postgres
-               search_path: [osm, public]
-           id_field: osm_id
-           table: hotosm_bdi_waterways
-           geom_field: foo_geom
-           options:
-               zoom:
-                   min: 0
-                   max: 15
-           format:
-               name: pbf
-               mimetype: application/vnd.mapbox-vector-tile
+         name: MVT-postgresql
+         data:
+             host: 127.0.0.1
+             port: 3010 # Default 5432 if not provided
+             dbname: test
+             user: postgres
+             password: postgres
+             search_path: [osm, public]
+         id_field: osm_id
+         table: hotosm_bdi_waterways
+         geom_field: foo_geom
+         storage_crs: http://www.opengis.net/def/crs/EPSG/0/4326
+         options:
+             zoom:
+                 min: 0
+                 max: 15
+         format:
+             name: pbf
+             mimetype: application/vnd.mapbox-vector-tile
+
+.. tip::
+   Geometry must have correctly defined :ref:`storage_crs<crs>`
 
 PostgreSQL-related connection options can also be added to `options`. Please refer to the :ref:`PostgreSQL OGC Features Provider<PostgreSQL>` documentation for more information.
 

pygeoapi/api/tiles.py

@@ -246,7 +246,7 @@ def get_collection_tiles_data(
         if content is None:
             msg = 'identifier not found'
             return api.get_exception(
-                HTTPStatus.NO_CONTENT, headers, format_, 'NocContent', msg)
+                HTTPStatus.NO_CONTENT, headers, format_, 'NoContent', msg)
         else:
             return headers, HTTPStatus.OK, content