ATTENTION! The files in this dataset are designed for streaming, not downloading. For the best experience, please follow the instructions available in the resources. In replacement of the former...

The Medium Resolution Digital Elevation Model (MRDEM) product is a multi-source product that integrates elevation data from the Copernicus DEM acquired during the TanDEM-X Mission, and the High…

As the title suggests, I am trying to generate an mbtiles file with terrarium encoded png or webp tiles from Canada's MRDEM dataset. The data source has provided a VRT file which can be found here: https://canelevation-dem.s3.ca-central-1.amazonaws.com/mrdem-30/mrdem-30-dtm.vrt Based on my readings, I need to re-project the data to web mercator and then from there use rio rgbify to generate the DTM image tiles. For some reason the end result is a tile set that appears to be incorrectly encoded and I'm really not sure where I've gone wrong. The commands I run are first: gdalwarp -te -123.2839 48.9985 -121.8017 49.9344 \ -s_srs EPSG:3979 \ -t_srs EPSG:3857 \ -srcnodata -32767 \ -dstnodata 0 \ -of GTiff mrdem-30-dtm.vrt vancouver-mountains.tif The bounding box is a small mountainous area around Vancouver,BC. I've tried running this command without the dstnodata/srcnodata configuration as well. Then I run: rio rgbify \ --min-z=0 \ --max-z=14 \ --format=webp \ vancouver-mountains.tif \ vancouver-mountains.mbtiles I convert to pmtiles and quickly test: pmtiles convert vancouver-mountains.mbtiles vancouver-mountains.pmtiles pmtiles serve . --port=8900 --cors='*' The end result is a tile that looks like so: By contrast here is the same tile that's generated by the Mapterhorn project: Tile reference is: 10/162/349.webp. When running rgbify I do get some warnings: .pyenv/python3/lib/python3.13/site-packages/rasterio/__init__.py:366: NotGeoreferencedWarning: Dataset has no geotransform, gcps, or rpcs. The identity matrix will be returned. So that could be a clue, but I'm not well versed enough in this area to make sense of it. Finally the gdalinfo output for the VRT file is: Driver: VRT/Virtual Raster Files: mrdem-30-dtm.vrt /vsicurl/https://canelevation-dem.s3.ca-central-1.amazonaws.com/mrdem-30/mrdem-30-dtm.tif Size is 183685, 159653 Coordinate System is: PROJCRS["NAD83(CSRS) / Canada Atlas Lambert", BASEGEOGCRS["NAD83(CSRS)", DATUM["NAD83 Canadian Spatial Reference System", ELLIPSOID["GRS 1980",6378137,298.257222101, LENGTHUNIT["metre",1]]], PRIMEM["Greenwich",0, ANGLEUNIT["degree",0.0174532925199433]], ID["EPSG",4617]], CONVERSION["unnamed", METHOD["Lambert Conic Conformal (2SP)", ID["EPSG",9802]], PARAMETER["Latitude of false origin",49, ANGLEUNIT["degree",0.0174532925199433], ID["EPSG",8821]], PARAMETER["Longitude of false origin",-95, ANGLEUNIT["degree",0.0174532925199433], ID["EPSG",8822]], PARAMETER["Latitude of 1st standard parallel",49, ANGLEUNIT["degree",0.0174532925199433], ID["EPSG",8823]], PARAMETER["Latitude of 2nd standard parallel",77, ANGLEUNIT["degree",0.0174532925199433], ID["EPSG",8824]], PARAMETER["Easting at false origin",0, LENGTHUNIT["metre",1], ID["EPSG",8826]], PARAMETER["Northing at false origin",0, LENGTHUNIT["metre",1], ID["EPSG",8827]]], CS[Cartesian,2], AXIS["easting",east, ORDER[1], LENGTHUNIT["metre",1]], AXIS["northing",north, ORDER[2], LENGTHUNIT["metre",1]], ID["EPSG",3979]] Data axis to CRS axis mapping: 1,2 Origin = (-2453970.000000000000000,3887370.000000000000000) Pixel Size = (30.000000000000000,-30.000000000000000) Corner Coordinates: Upper Left (-2453970.000, 3887370.000) (175d38'54.16"W, 68d 7'32.66"N) Lower Left (-2453970.000, -902220.000) (121d27'10.59"W, 36d35'37.91"N) Upper Right ( 3056580.000, 3887370.000) ( 10d43'19.02"W, 62d45'37.02"N) Lower Right ( 3056580.000, -902220.000) ( 63d 0'40.11"W, 34d21' 7.53"N) Center ( 301305.000, 1492575.000) ( 88d57'23.39"W, 62d31'56.78"N) Band 1 Block=128x128 Type=Float32, ColorInterp=Gray NoData Value=-32767 Overviews: 91842x79826, 45921x39913, 22960x19956, 11480x9978, 5740x4989, 2870x2494, 1435x1247, 717x623