<div dir="ltr">Hello,<div><br></div><div>I've been exploring the pmtiles driver a bit, what a great way to read giant blobs of tiles off of cloud storage, really impressed by how fast this is.</div><div><br></div><div>I've been encountering a "Non increasing tile_id" when attempting to read from coarser zoom levels for a large (~300GB) pmtile file:</div><div><br></div><div><a href="https://github.com/OSGeo/gdal/blob/48ddaf28db2826c635243abd469e7c7d2e749bdc/ogr/ogrsf_frmts/pmtiles/ogrpmtilestileiterator.cpp#L329">https://github.com/OSGeo/gdal/blob/48ddaf28db2826c635243abd469e7c7d2e749bdc/ogr/ogrsf_frmts/pmtiles/ogrpmtilestileiterator.cpp#L329</a><br></div><div><br></div><div>For example, this fails for me:</div><div><br></div><div>ds = gdal.OpenEx("</vsis3/...", gdal.OF_VECTOR | gdal.OF_READONLY, open_options=["ZOOM_LEVEL=0"])<br></div><div>lyr = ds.GetLayer()<br></div><div>for feat in lyr:<br>    *** raises RuntimeError: Non increasing tile_id ***</div><div><br>Things seem good when I specify higher zooms (e.g. Z12 or not specified at all), but I can't get features out of the coarser ones.  Using the "pmtiles tiles" command, I can access the 0/0/0 tile just fine though.</div><div><br></div><div>Any thoughts on what could be wrong?  The file was translated from a mbtiles file via the "pmtiles convert" tool.  ogrinfo is happy with it as well:<br><br>INFO: Open of `/vsis3/.../tileset.pmtiles'<br>      using driver `PMTiles' successful.<br>Metadata:<br>  description=mvts for landcover layer<br>  format=pbf<br>  maxzoom=16<br>  minzoom=0<br>  name=mvts<br>  scheme=tms<br>  type=overlay<br>  version=2<br>  ZOOM_LEVEL=16<br>1: landcover<br></div><div><br></div><div>Thanks,</div><div>Patrick</div></div>