I am in the process of rendering a series of map tiles based on the OS OpenData products using the script (and an updated version that uses all cores on the machine to speed things up).  The different raster products are rendered at different scales and then displayed using LeafletJS and OpenLayers applications as simple demonstrations.

The following command generates the tiles for the zoom levels I need:

python -z '7-9' -e -p raster -r average osvmd.vrt osvmdtiles/

This creates zoom levels 7, 8 and 9 using the average resampling method from the osvmd.vrt virtual raster and outputs the tiles in the osvmdtiles directory.

Zoom level 7
Zoom level 7
Zoom Level 8
Zoom Level 8
Zoom Level 9
Zoom level 9

The issues with zoom level 7 appear to be random.  It could be anything to do with the source raster, the projection.  It is just odd that the same virtual raster works for the other zoom levels and not that one.

If you have any answers, please post them here:

Leave a Reply

Your email address will not be published. Required fields are marked *

This site uses Akismet to reduce spam. Learn how your comment data is processed.