Map Rasters! - Map Something!

license: CC-BY-4.0

Remeber our pretty visualization?

That’s not a map!

IMPORTS¶

from typing import Any, Union
import sys
from IPython.display import Image
from  pathlib import Path
import matplotlib.pyplot as plt
from localtileserver import get_leaflet_tile_layer, TileClient
import ipyleaflet
import localtileserver
import leafmap

CONSTANTS¶

GEOTIFF_FOLDER = '../data/raster'
COG_FOLDER = '../data/raster/cogs'
MAP_CENTER = (39.1780195, -120.17694813)

GEOTIFF_PATHS = [p for p in Path(GEOTIFF_FOLDER).glob(f'*.tif')]
COG_PATHS = [p for p in Path(COG_FOLDER).glob(f'*.tif')]
PATHS = COG_PATHS

for p in PATHS:
    print(p.name)

tahoe-dem.tif
tahoe-mean-s2-20210101_20210301.tif

IPYLEAFLET¶

Try replacing the PATHS = COG_PATHS with the PATHS = GEOTIFF_PATHS ... It slows to a halt

client = TileClient(PATHS[0])
dem_layer = get_leaflet_tile_layer(client)
ipyleaflet_map = ipyleaflet.Map(center=client.center(), zoom=client.default_zoom)
ipyleaflet_map.add(dem_layer)
ipyleaflet_map

LEAFMAP¶

basemaps = list(leafmap.basemaps.keys())
print(len(basemaps))
basemaps[:10]

['OpenStreetMap',
 'Google Maps',
 'Google Satellite',
 'Google Terrain',
 'Google Hybrid',
 'FWS NWI Wetlands',
 'FWS NWI Wetlands Raster',
 'NLCD 2021 CONUS Land Cover',
 'NLCD 2019 CONUS Land Cover',
 'NLCD 2016 CONUS Land Cover']

leafmap_map = leafmap.Map(basemap='NASAGIBS.BlueMarble')
leafmap_map.add_raster(str(PATHS[0]), bands=[1], layer_name='dem', vmin=1000, vmax=3000, cmap='hot')
leafmap_map.add_raster(str(PATHS[1]), layer_name='s2', bands=[4,3,2], vmin=0, vmax=5000)
leafmap_map

MAKE YOUR OWN TILE SERVER¶

using rio-tiler you can generate maptiles from cogs. maptiler/app.py is an example taken from rio-tiler’s docs. i’m using it here

url = 'http://127.0.0.1:8000/8/42/97.png'
from io import BytesIO
from urllib.request import urlopen
im_bytes = BytesIO(urlopen(url).read())
arr = plt.imread(im_bytes)
print(arr.shape, arr.max(), arr.min())
plt.axis('off')
_ = plt.imshow(arr)

(256, 256, 4) 1.0 0.0

leafmap_map_2 = leafmap.Map(zoom=8, height=600, width=1200)
leafmap_map_2.center = MAP_CENTER
leafmap_map_2.add_tile_layer(
    url="http://127.0.0.1:8000/{z}/{x}/{y}.png",
    bands=[4,3,2],
    name="Tile Server",
    attribution="Rio Tiler",
)
leafmap_map_2.add_layer_control(position="topleft")
leafmap_map_2