Add support for K-means and Spectral Angle Mapper (SAM)

[giswqs]

Compute cyanobacteria band ratios

import hypercoast
from hypercoast.pace import cyano_band_ratios

filepath = "data/PACE_OCI.20240730T181157.L2.OC_AOP.V2_0.NRT.nc"
dataset = hypercoast.read_pace(filepath)
da = cyano_band_ratios(dataset, plot=True)

[github-actions]

🚀 Deployed on https://670b45ceda100bf72605e9cf--opengeos.netlify.app

[giswqs]

K-means

import hypercoast
from hypercoast.pace import apply_kmeans

filepath = "data/PACE_OCI.20240730T181157.L2.OC_AOP.V2_0.NRT.nc"
dataset = hypercoast.read_pace(filepath)

cluster_labels, latitudes, longitudes = apply_kmeans(dataset, n_clusters=6)

[giswqs]

PCA

import hypercoast
from hypercoast.pace import apply_pca

filepath = "data/PACE_OCI.20240730T181157.L2.OC_AOP.V2_0.NRT.nc"
dataset = hypercoast.read_pace(filepath)

pca_data = apply_pca(dataset, n_components=3, x_component=0, y_component=1)

[giswqs]

Apply K-means to selected pixels

import hypercoast
from hypercoast.pace import apply_kmeans

filepath = "data/PACE_OCI.20240730T181157.L2.OC_AOP.V2_0.NRT.nc"
dataset = hypercoast.read_pace(filepath)

da = dataset["Rrs"]

filter_condition = (
    (da.sel(wavelength=650) > da.sel(wavelength=620))
    & (da.sel(wavelength=701) > da.sel(wavelength=681))
    & (da.sel(wavelength=701) > da.sel(wavelength=450))
)
extent = [-95, -85, 27, 33]
colors = ["#e41a1c", "#377eb8", "#4daf4a", "#f781bf", "#a65628", "#984ea3"]

cluster_labels, latitudes, longitudes = apply_kmeans(
    da, 
    n_clusters=6, 
    filter_condition=filter_condition,
    extent=extent,
    colors=colors
    )

[giswqs]

Apply Spectral Angle Mapper

import hypercoast
from hypercoast.pace import apply_sam

filepath = "data/PACE_OCI.20240730T181157.L2.OC_AOP.V2_0.NRT.nc"
dataset = hypercoast.read_pace(filepath)

data, latitudes, longitudes = apply_sam(
    dataset,
    n_components=3, 
    n_clusters=6,
    )

extent = [-95, -85, 27, 33]
colors = ["#377eb8", "#ff7f00", "#4daf4a", "#f781bf", "#a65628", "#984ea3"]
data, latitudes, longitudes = apply_sam(
    dataset,
    n_components=3, 
    n_clusters=6,
    extent=extent,
    colors=colors,
    )

[giswqs]

Add SAM filter condition

da = dataset["Rrs"]

filter_condition = (
    (da.sel(wavelength=650) > da.sel(wavelength=620))
    & (da.sel(wavelength=701) > da.sel(wavelength=681))
    & (da.sel(wavelength=701) > da.sel(wavelength=450))
)
extent = [-95, -85, 27, 33]
colors = ["#e41a1c", "#377eb8", "#4daf4a", "#f781bf", "#a65628", "#984ea3"]

data, latitudes, longitudes = apply_sam(
    dataset,
    n_components=3, 
    n_clusters=6,
    filter_condition=filter_condition,
    extent=extent,
    colors=colors,
    )

[giswqs]

Add support for spectral library

import hypercoast
from hypercoast.pace import apply_sam_spectral

filepath = "data/PACE_OCI.20240730T181157.L2.OC_AOP.V2_0.NRT.nc"
dataset = hypercoast.read_pace(filepath)
url = "https://github.com/opengeos/datasets/releases/download/hypercoast/SAM_spectral_library.zip"
hypercoast.download_file(url)
spectral_library = "./SAM_spectral_library/*.csv"

extent = [-95, -85, 27, 33]
data, latitudes, longitudes = apply_sam_spectral(
    dataset,
    spectral_library=spectral_library,
    extent=extent,
    )

da = dataset["Rrs"]
extent = [-95, -85, 27, 33]
filter_condition = (
    (da.sel(wavelength=650) > da.sel(wavelength=620))
    & (da.sel(wavelength=701) > da.sel(wavelength=681))
    & (da.sel(wavelength=701) > da.sel(wavelength=450))
)
data, latitudes, longitudes = apply_sam_spectral(
    da,
    spectral_library=spectral_library,
    filter_condition=filter_condition,
    extent=extent,
    )