In this section we take a look at adding our computed NDVI metrics to our map, along with methods for publishing our maps as interactive web applications.

from ibis import _
con = ibis.duckdb.connect(extensions=["spatial"])

city = (con
        .read_geo("/vsicurl/https://dsl.richmond.edu/panorama/redlining/static/mappinginequality.gpkg")
        .filter(_.city == "New Haven")
        .execute()
       )
import leafmap.maplibregl as leafmap

m = leafmap.Map(style = "voyager")
m.add_cog_layer("https://github.com/espm-157/static-maps/raw/refs/heads/main/ndvi.tif", colormap_name = "greens", nodata=1)
m.add_gdf(city, "fill",  paint = { "fill-color": ["get", "fill"], "fill-opacity": 0.8})
m