reflspeckit 0.2.2

Utility package for working with spectral data cubes.

🛰️ reflspeckit

⚙️ A modern toolkit for working with any and all flavors of spectral data with a focus on applications for reflectance/emittance imaging spectroscopy

---

---

🧠 What is reflspeckit?

reflspeckit is a lightweight, modular Python package designed to make analysis of spectral data cubes simple, flexible, and fun. Whether you're exploring planetary hyperspectral data, performing band analysis, or building your own spectral pipelines — this toolkit’s got you covered.

---

⚙️ Package Structure

reflspeckit provides two primary classes for the analysis of spectral data plus a third class specialized for large datasets:

• 📄Spec1D - handles 1-dimensional, single spectrum data
• 📒Spec3D - handles 3-dimensional spectral image cubes
• 🗄️StreamingSpec3D - handles large image cubes using a streaming approach

Each class has equivalent methods, which are listed below:

---

🧰 Available Methods

Method	Description
🚩outlier_removal	Removes anomalous data in the spectral domain
🔊noise_reduction	Provides filtering methods to smooth data in the spectral domain
🚫continuum_removal	Estimates and removes a spectral contrinuum from the data
🖼️make_m3_rgb (3D only)	Creates a standard mafic mineral RGB color-composite image
📈fit_absorption	Performs a least squares polynomial fit over a spectral region

💡 Spectral Utilities

Various spectral utilities are available through the reflspeckit.utils subpackge.

Module	Description
get_nonzero	If you have an empty 3D image array with the first two dimensions being pixels and the third dimension of size N, and each pixel is filled in to a certain depth, M <= N, this function returns a 2D image array that picks out all the pixel values at position M.
rgb_composite	Turns any three 2D arrays into a normalized rgb color composite image. The values of each band will be stretched from 0-255, with values above the 95th percentile being cut off to preserve the color stretch integrity.
wvl_search	Given an estimated wavelength values and a list of real wavelength vales, this module will return the real wavelength values that is closest to the estimate and the index at which this wavelength value is located within the actual wavelength array.

More utilities coming soon! As a work through my Ph.D., I will add all the various utility functions I write for spectral data processing here!

---

🚀 Quick Start

pip install reflspeckit

import reflspeckit as rsk
import matplotlib.pyplot as plt

# Loading in a single spectrum
my_spectrum = rsk.Spec1d(spectrum_array, wavelength_array, unit="nm")
my_spectrum.remove_outliers()
myspectrum.noise_reduction(method="box_filter", filter_width=5)
myspectrum.continuum_removal(method="double_line")
print(myspectrum.filtered)  # Contains filtered spectrum

absorption_feature = myspectrum.fit_absorption(800, 1200, unit="nm")
print(absorption_feature.calculate_ibd())  # Returns integrated band depth.

# Loading in a spectral image cube
my_cube = rsk.Spec1d(cube_array, wavelength_array)
my_cube.remove_outliers()
my_cube.noise_reduction(method="box_filter", filter_width=5)
my_cube.continuum_removal(method="double_line")
print(myspectrum.cube)  # Sequentially replaces myspectrum.cube to save memory.

rgb = my_cube.make_m3_rgb()
plt.imshow(rgb)  # Shows RGB color-composite image.

🔗 Links

• GitHub: https://github.com/z-vig/reflspeckit.git
• Docs: (coming soon!)