feat: CITM2000 debug tile matrix

[blacha]

Motivation

We are starting to want to view imagery in other projections that do not have a standard Web Quad tile matrixes to work with. It would be helpful to have a way to generate a debug tile matrix that aligns with the WebMercatorQuad standard scales, while also being a quad starting at 1x1 -> 2x2 -> 4x4 grids that increase with powers of two.

Modifications

Added a CITM Debug tile matrix based off the WebMercatorQuad scales.

Added a script to generate tile matrices from a supplied EPSG projection,

The script will download the projection extent from a projjson, then calculate an appropriate WebMercatorQuad scale that will cover the entire projection extent,

This script was only tested with CITM and NZTM it will likely fail if used with other projections.

Verification

Tested basemaps server with QGIS and some CI topo50 mapsheets.

[tawera-manaena]

This script was only tested with CITM and NZTM it will likely fail if used with other projections.

We have a Jira ticket that lists the EPSG codes for each Topo Raster Map Series. I'll run some through the script and see what happens.

[tawera-manaena]

@blacha - I'm hitting the inequality check here. Looks like the xyOrder and axisOrder variables hold values for different coordinate spaces. xyOrder being cartesian and axisOrder being geographic. I don't imagine the inequality check intends for this?

# TargetProjection = 3793;

# xyOrder: 'yx'
# axisOrder: 'n,e'

Error: getXyOrder: yx is not the same as expected ordering from proj json

[blacha]

@blacha - I'm hitting the inequality check here. Looks like the xyOrder and axisOrder variables hold values for different coordinate spaces. xyOrder being cartesian and axisOrder being geographic. I don't imagine the inequality check intends for this?

# TargetProjection = 3793;

# xyOrder: 'yx'
# axisOrder: 'n,e'

Error: getXyOrder: yx is not the same as expected ordering from proj json

This should now be fixed, it only handled when the axis was labelled x or y not n or e

[tawera-manaena]

We have a Jira ticket that lists the EPSG codes for each Topo Raster Map Series. I'll run some through the script and see what happens.

@blacha - I've run the script with the following EPSG codes:

• Antarctic (ANT)

  • McMurdo Sound Lambert Conformal 2000; RSRGD2000

    • 5479

• New Zealand Offshore Islands (NZOI)

  • Antipodes Islands Transverse Mercator 2000; NZGD 2000

    • 3790

  • Campbell Island Transverse Mercator 2000; NZGD 2000

    • 3789

  • Raoul Island Transverse Mercator 2000; NZGD 2000

    • 3791

  • Auckland Islands Transverse Mercator 2000; NZGD 2000

    • 3788

Funnily enough, they are all defined in [y, x]. So, we'll need to update the getXyOrder function to handle this once we support these EPSG codes. Otherwise, I've checked the GeoJSON bounds in QGIS and the generated TileMatrix definitions. Both of which look good. I'm happy to approve this work.