In [ ]:
from IPython.display import display, Image
from geonotebook.wrappers import RasterData
import matplotlib.pylab as plt

In [ ]:
EXPECTED="https://data.kitware.com/api/v1/file/586a8e158d777f1e3428d842/download"

M.set_center(-120.32, 47.84, 8)
rd = RasterData('vrts/four_band_weld.vrt')

M.add_layer(rd[1, 2, 3], opacity=1.0, gamma=2.5)

display(Image(EXPECTED, format="png"))

In [ ]:
EXPECTED="https://data.kitware.com/api/v1/item/586a8ec98d777f1e3428d844/download"

M.remove_layer(M.layers[0])

cm = plt.get_cmap('winter', 20)

M.add_layer(rd[4], colormap=cm, interval=(-1.0, 1.0))

display(Image(EXPECTED, format="jpeg"))