Water column compensation workflow for hyperspectral imaging data

Our article describes a data processing workflow for hyperspectral imaging data to compensate for the water column in shallow, clear to moderate optical water types. We provide a MATLAB script that can be readily used to implement the described workflow. We break down each code segment of this script so that it is more approachable for use and modification by end users and data providers. The workflow initially implements the method for water column compensation described in Lyzenga (1978) and Lyzenga (1981), generating depth invariant indices from spectral band pairs. Given the high dimensionality of hyperspectral imaging data, an overwhelming number of depth invariant indices are generated in the workflow. As such, a correlation based feature selection methodology is applied to remove redundant depth invariant indices. In a post-processing step, a principal component transformation is applied, extracting features that account for a substantial amount of the variance from the non-redundant depth invariant indices while reducing dimensionality. To fully showcase the developed methodology and its potential for extracting bottom type information, we provide an example output of the water column compensation workflow using hyperspectral imaging data collected over the coast of Philpott's Island in Long Sault Parkway provincial park, Ontario, Canada.•Workflow calculates depth invariant indices for hyperspectral imaging data to compensate for the water column in shallow, clear to moderate optical water types.•The applied principal component transformation generates features that account for a substantial amount of the variance from the depth invariant indices while reducing dimensionality.•The output (both depth invariant index image and principal component image) allows for the analysis of bottom type in shallow, clear to moderate optical water types.