Failing rendering with non-square pixels

[cbdj]

Description

When rendering an EPSG:4326 COG in a Web Mercator viewport, tiles are reprojected wrongly and miss the lower part of the COG:

When dezoomed, the missing lower part are now displayed but there is still a reprojection issue:

Using a custom epsgResolver based on the @developmentseed/epsg/all.csv.gz instead of the default remote epsg.io yields the same result.

Reproduction

Here is a minimal reproducing example app:

reprojection_issue.zip

[kylebarron]

What version of packages are you using? Your package.json only says latest and you don't have a lockfile.

[kylebarron]

Also you need to ensure you're pinned to the latest version of proj4.js. proj4js/proj4js#545

[cbdj]

Here are the updated package.json and pnpm-lock.yaml.
I'm using @developmentseed/deck.gl-geotiff@0.4.0 and proj4@2.20.4
pnpm-lock.yaml
package.json

[kylebarron]

I spent a bit of time debugging this and it seems down to the fact that your pixels are non-square.

If you check gdalinfo, you see

Origin = (-5.499844444444444,49.035866666666664)
Pixel Size = (0.000366210937500,-0.000224784207161)

So each pixel is 63% wider than it is tall.

The problem is that currently we use OGC TileMatrixSet as the intermediate abstraction for traversing tiles, and that only supports a single pixel size across height and width.

This is planned to be refactored a bit anyways as we add support for (Geo)Zarr in the next month.