Visualizing AVIRIS data interactively with HyperCoast¶

This notebook demonstrates how to visualize AVIRIS hyperspectral data interactively with HyperCoast. For more information about AVIRIS, please visit the links below:

• https://aviris.jpl.nasa.gov/
• https://aviris.jpl.nasa.gov/dataportal/
• https://popo.jpl.nasa.gov/mmgis-aviris/?s=ujooa
• https://daac.ornl.gov/cgi-bin/dsviewer.pl?ds_id=1988
• https://github.com/ornldaac/deltax_workshop_2022/tree/main
• https://github.com/jjmcnelis/aviris-ng-notebooks/tree/master

In [ ]:

# %pip install "hypercoast[extra]"

# %pip install "hypercoast[extra]"

In [ ]:

import hypercoast

import hypercoast

Download sample dataset from Access the AVIRIS-NG L2 Surface Reflectance product page via ORNL DAAC https://daac.ornl.gov/cgi-bin/dsviewer.pl?ds_id=1988

Download ang20210401t150456_rfl_v2z1.zip and unzip it.

The dataset contains 2 files: ang20210401t150456_rfl_v2z1 and ang20210401t150456_rfl_v2z1.hdr. We will use the ang20210401t150456_rfl_v2z1 file in this notebook.

In [ ]:

filepath = "ang20210401t150456_rfl_v2z1"

filepath = "ang20210401t150456_rfl_v2z1"

Read the AVIRIS data as an xarray.Dataset object.

In [ ]:

ds = hypercoast.read_aviris(filepath)
ds

ds = hypercoast.read_aviris(filepath) ds

Create an interactive map.

In [ ]:

m = hypercoast.Map()
m

m = hypercoast.Map() m

Add the AVIRIS data to the map.

In [ ]:

m.add_aviris(ds, wavelengths=[1000, 700, 400], vmin=0, vmax=0.2)
m.add("spectral")

m.add_aviris(ds, wavelengths=[1000, 700, 400], vmin=0, vmax=0.2) m.add("spectral")