Offset bug with "Sphere Azimuthal Equidistant" projection

Description

I saw that the "Sphere Azimuthal Equidistant" was supported from 2.10-beta, but when I combine a WMTS from GeoServer with a tileset from another service they don't match (see screenshot). This works fine when I tested with another projection.

I'm using ESRI:53032: https://epsg.io/53032

I'm adding this projection to epsg.properties:

53032=PROJCS["Sphere_Azimuthal_Equidistant",GEOGCS["GCS_Sphere",DATUM["Not_specified_based_on_Authalic_Sphere",SPHEROID["Sphere",6371000,0]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Azimuthal_Equidistant"],PARAMETER["False_Easting",0],PARAMETER["False_Northing",0],PARAMETER["Central_Meridian",0],PARAMETER["Latitude_Of_Origin",0],UNIT["Meter",1],AUTHORITY["EPSG","53032"]]

I then define a gridset/TileMatrixSet for this projection (see screenshot), and create a WMTS.

I'm adding the same projection to Carto/PostGIS - as I want to combine layers from these two services. The extents and resolutions are the same for both services.

When I add both layers they don't match completely. When I try to do the same with another projection, "North Pole Azimuthal Equidistant" - http://epsg.io/102016, everything lines up perfectly.

I'm not sure if the misplacement origins in GeoServer, GeoTools, PostGIS or Carto...

Environment

OSX

Assignee

Ben Caradoc-Davies

Reporter

Bjørn Sandvik

Triage

None

Fix versions

None

Affects versions

Priority

Medium
Configure