Visualizing Spectral Data from Field Measurements¶

This notebook demonstrates how to visualize spectral data from field measurements.

In [ ]:

# %pip install "hypercoast[extra]"

# %pip install "hypercoast[extra]"

In [ ]:

import hypercoast
import pandas as pd

import hypercoast import pandas as pd

Download a sample dataset.

In [ ]:

url = "https://github.com/opengeos/datasets/releases/download/hypercoast/pace_sample_points.csv"
data = pd.read_csv(url)
data.head()

url = "https://github.com/opengeos/datasets/releases/download/hypercoast/pace_sample_points.csv" data = pd.read_csv(url) data.head()

Download PACE data.

In [ ]:

url = "https://github.com/opengeos/datasets/releases/download/hypercoast/PACE_OCI.20240730T181157.L2.OC_AOP.V2_0.NRT.nc"
filepath = "data/PACE_OCI.20240730T181157.L2.OC_AOP.V2_0.NRT.nc"
hypercoast.download_file(url, filepath)

url = "https://github.com/opengeos/datasets/releases/download/hypercoast/PACE_OCI.20240730T181157.L2.OC_AOP.V2_0.NRT.nc" filepath = "data/PACE_OCI.20240730T181157.L2.OC_AOP.V2_0.NRT.nc" hypercoast.download_file(url, filepath)

Read the PACE dataset.

In [ ]:

dataset = hypercoast.read_pace(filepath)

dataset = hypercoast.read_pace(filepath)

Run the following cell to show the map. Click on the markers to see the spectral data.

In [ ]:

m = hypercoast.Map(center=[27.235094, -87.791748], zoom=6)

m.add_basemap("Hybrid")
wavelengths = [450, 550, 650]
m.add_pace(
    dataset, wavelengths, indexes=[3, 2, 1], vmin=0, vmax=0.02, layer_name="PACE"
)
m.add("spectral")

m.add_field_data(
    data,
    x_col="wavelength",
    y_col_prefix="(",
    x_label="Wavelength (nm)",
    y_label="Reflectance",
    use_marker_cluster=True,
)
m.set_center(-87.791748, 27.235094, zoom=6)
m

m = hypercoast.Map(center=[27.235094, -87.791748], zoom=6) m.add_basemap("Hybrid") wavelengths = [450, 550, 650] m.add_pace( dataset, wavelengths, indexes=[3, 2, 1], vmin=0, vmax=0.02, layer_name="PACE" ) m.add("spectral") m.add_field_data( data, x_col="wavelength", y_col_prefix="(", x_label="Wavelength (nm)", y_label="Reflectance", use_marker_cluster=True, ) m.set_center(-87.791748, 27.235094, zoom=6) m