<div dir="ltr">Thanks Even!<div><br></div><div>I got a much smaller reproduce case with data ready at <a href="https://github.com/OSGeo/gdal/issues/9288">https://github.com/OSGeo/gdal/issues/9288</a></div><div><br></div><div>I created the test case with `pmtiles extract --maxzoom 7 ...`, interestingly enough using a coarser zoom didn't surface the issue!</div><div><br></div><div>Patrick</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Feb 23, 2024 at 1:37 AM Even Rouault <<a href="mailto:even.rouault@spatialys.com">even.rouault@spatialys.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hi,<br>
<br>
this is likely a bug. I would need access to the file to investigate, <br>
ideally through /vsicurl, /vsis3, etc. . I don't have storage to host <br>
300 GB. And possibly try to reduce the file size by removing upper zoom <br>
levels<br>
<br>
Even<br>
<br>
Le 23/02/2024 à 06:54, Patrick Young via gdal-dev a écrit :<br>
> Hello,<br>
><br>
> I've been exploring the pmtiles driver a bit, what a great way to read <br>
> giant blobs of tiles off of cloud storage, really impressed by how <br>
> fast this is.<br>
><br>
> I've been encountering a "Non increasing tile_id" when attempting to <br>
> read from coarser zoom levels for a large (~300GB) pmtile file:<br>
><br>
> <a href="https://github.com/OSGeo/gdal/blob/48ddaf28db2826c635243abd469e7c7d2e749bdc/ogr/ogrsf_frmts/pmtiles/ogrpmtilestileiterator.cpp#L329" rel="noreferrer" target="_blank">https://github.com/OSGeo/gdal/blob/48ddaf28db2826c635243abd469e7c7d2e749bdc/ogr/ogrsf_frmts/pmtiles/ogrpmtilestileiterator.cpp#L329</a><br>
><br>
> For example, this fails for me:<br>
><br>
> ds = gdal.OpenEx("</vsis3/...", gdal.OF_VECTOR | gdal.OF_READONLY, <br>
> open_options=["ZOOM_LEVEL=0"])<br>
> lyr = ds.GetLayer()<br>
> for feat in lyr:<br>
>     *** raises RuntimeError: Non increasing tile_id ***<br>
><br>
> Things seem good when I specify higher zooms (e.g. Z12 or not <br>
> specified at all), but I can't get features out of the coarser ones.  <br>
> Using the "pmtiles tiles" command, I can access the 0/0/0 tile just <br>
> fine though.<br>
><br>
> Any thoughts on what could be wrong?  The file was translated from a <br>
> mbtiles file via the "pmtiles convert" tool.  ogrinfo is happy with it <br>
> 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>
><br>
> Thanks,<br>
> Patrick<br>
><br>
> _______________________________________________<br>
> gdal-dev mailing list<br>
> <a href="mailto:gdal-dev@lists.osgeo.org" target="_blank">gdal-dev@lists.osgeo.org</a><br>
> <a href="https://lists.osgeo.org/mailman/listinfo/gdal-dev" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/gdal-dev</a><br>
<br>
-- <br>
<a href="http://www.spatialys.com" rel="noreferrer" target="_blank">http://www.spatialys.com</a><br>
My software is free, but my time generally not.<br>
<br>
</blockquote></div>