acalcutt / wdb-map-gen

OpenMapTiles basemap generation for WifiDB.net
Other
5 stars 0 forks source link

Globallandcover layer #10

Open webmap2015 opened 2 years ago

webmap2015 commented 2 years ago

Hi, There is a two landcover layer.

landcover Landcover is used to describe the physical material at the surface of the earth. At lower zoom levels this is from Natural Earth data for glaciers and ice shelves and at higher zoom levels the landcover is implied by OSM tags. The most common use case for this layer is to style wood (class=wood) and grass (class=grass) areas.

globallandcover Global landcover is a derivated product from imagery made by ESA as a part of the ESA Climate Change Initiative and in particular its Land Cover project. ESA Landcover v1.1 is processed and vectorized into a data source for global landcover from zoom 0 to zoom 9.

It is vectorized and contains classes as grass, crop, scrub, forest, and so on. Vector tiles are generated up to zoom level 9. Data have six classes: crop, forest, grass, scrub, snow, tree.

https://cloud.maptiler.com/tiles/landcover/ https://data.maptiler.com/downloads/planet/ https://data.maptiler.com/downloads/dataset/landcover/ https://docs.maptiler.com/schema/planet/#globallandcover

https://data.maptiler.com/downloads/tileset/landcover/ maptiler-landcover-2020-03-07-planet.mbtiles 385.6 MiB

How can we produce globallandcover mbtiles? Thanks.

acalcutt commented 2 years ago

I'm not really sure on this. I would start by trying to download the ESA Landcover data from https://www.esa-landcover-cci.org/?q=node/164 ( I tried but i couldn't get their ftp to load )

Then on this page ( http://maps.elie.ucl.ac.be/CCI/viewer/download.php ) they have this example to make a GeoTiff gdalwarp -of Gtiff -co COMPRESS=LZW -co TILED=YES -ot Byte -te -180.0000000 -90.0000000 180.0000000 90.0000000 -tr 0.002777777777778 0.002777777777778 -t_srs EPSG:4326 NETCDF:C3S-LC-L4-LCCS-Map-300m-P1Y-2017-v2.1.1.nc:lccs_class C3S-LC-L4-LCCS-Map-300m-P1Y-2017-v2.1.1.tif

But it looks like you are looking for a mbtiles file, so maybe try something like this

infile=NETCDF:C3S-LC-L4-LCCS-Map-300m-P1Y-2017-v2.1.1.nc:lccs_class
vrtfile=warp.vrt
mbtiles=out.mbtiles

#Convert the source to a EPSG:3857 projection VRT
gdalwarp -co COMPRESS=LZW -ot Byte -t_srs EPSG:3857 ${infile} ${vrtfile}
#Import VRT into MBTiles
gdal_translate ${vrtfile} ${mbtiles} -of MBTILES
#Create MBTiles Overview
gdaladdo ${mbtiles}
#CREATE UNIQUE INDEX tile_index on tiles (zoom_level, tile_column, tile_row);
sqlite3 ${mbtiles} 'CREATE UNIQUE INDEX tile_index on tiles (zoom_level, tile_column, tile_row);'