Datashading LandSat8 raster satellite imagery#

Datashader is fundamentally a rasterizing library, turning data into rasters (image-like arrays), but it is also useful for already-rasterized data like satellite imagery. For raster data, datashader accepts xarray raster data structures and then re-renders the data into whatever new bounding box and resolution the user requests. The rest of the datashader pipeline can then be used to visualize and analyze the data. This demo shows how to work with a set of raster satellite data, generating images as needed and overlaying them on geographic coordinates using HoloViews (version >1.9), GeoViews, and Bokeh.

import numpy as np
import rioxarray
import holoviews as hv
import geoviews as gv
import datashader as ds
import as ccrs
from holoviews import opts
from holoviews.operation.datashader import rasterize, shade