Skip to content
/ hypro Public

Workflow for geometric and radiometric calibration of airborne imaging spectroscopy data, especially using HySpex devices.

License

Notifications You must be signed in to change notification settings

EnSpec/hypro

Repository files navigation

HyPro

© Copyright 2020
Nanfeng Liu
Philip Townsend

Open-source software licensed under GNU GPLv3

Package Overview

Description

HyPro is a Python package that provides tools for processing raw hyperspectral image data from pushbroom sensors.

The main module, HyspexPro, is intended to be run as a script. Its functionality is summarized here.

  1. HyspexPro aims to do geometric and radiometric corrections on HySpex images. Our imaging system consists of two HySpex cameras (VNIR-1800 and SWIR-384) and one GPS-Inertial Measurement Unit (IMU) sensor (iMAR iTraceRT-F400). The VNIR-1800 sensor has 182 spectral bands within 400-1000 nm (spectral sampling = 3.26 nm). The SWIR-384 sensor has 288 bands within 960-2500 nm (spectral sampling = 5.45 nm). For more details about these cameras, please refer to the NEO website. The iTraceRT-F400 sensor records the attitudes (roll, pitch and heading angles) and positions (longitude, latitude and altitude) of the imaging platform.

  2. Geometric correction focuses on removing image geometric distortions mainly caused by platform motions, which can be characterized by roll, pitch and heading angles. We adapted the ray tracing method proposed by Meyer[1] to calculate the map coordinates (easting, northing and elevation) of each raw image pixel (a process called georeferencing). Then, the warp tool provided by GDAL (Geospatial Data Abstraction Library) was used to resample the raw image to regular map grids (a process called georectification).

  3. Radiometric correction focuses on reducing radiometric distortions mainly caused by sensor smile effects (shifts in sensor band center wavelengths and full-widths-at-half-maximum, or FWHMs) and atmospheric effects (absorption and scattering). The spectral matching method proposed by Gao et al.[2] was adapted to detect sensor smile effects. Then, raw radiance spectra were resampled to common center wavelengths using a cubic interpolation. In the atmospheric correction step, two atmospheric parameters (water vapor column and aerosol visibility) were retrieved from image spectra. The water vapor column was estimated via the continuum-interpolated band ratio (CIBR) technique proposed by Kaufman et al.[3]. The estimation of visibility was based on the dense dark vegetation (DDV) method proposed by Kaufman et al.[4]. All radiometric correction steps require an atmospheric lookup table which consists of the at-sensor radiance simulated under different atmospheric conditions and Sun-target-sensor view geometries. An open-source radiative transfer model, libRadtran, was used for this purpose. For more details about libRadtran, please refer to the website.

  4. The input dataset of HyspexPro includes:

    • HySpex raw digital number (DN) images;
    • HySpex IMU & GPS data processed in HySpex NAV;
    • HySpex sensor model files provided by NEO;
    • DEM;
    • HySpex lab radiometric re-calibration data provided by NEO (optional);
    • Ground control points data (optional).
  5. The output dataset of HyspexPro includes:

    • HySpex ground surface reflectance images;
    • View angles (zenith and azimuth);
    • DEM.

Processing Workflow

The standard workflow of HyspexPro is shown diagrammatically below:

DN-to-Reflectance Processing Pipeline

Credits

HyPro is authored by Nanfeng Liu, with contributions from Adam Chlus & Brendan Heberlein.

HyPro is a product of the UW Environmental Spectroscopy Laboratory (EnSpec) and is developed & maintained under Philip Townsend with the support of the University of Wisconsin.

License

HyPro is licensed for public use under GNU General Public License v3. For details of the licensing agreement, refer to ./LICENSE.txt. In brief:

This program is free software: you can redistribute it and/or modify
it under the terms of the GNU General Public License Version 3,
as published by the Free Software Foundation.

This program is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
GNU General Public License for more details.

You should have received a copy of the GNU General Public License Version 3
along with this program. If not, see <https://www.gnu.org/licenses/>.

Ancillary Datasets

Solar Irradiance

HyPro is distributed with solar spectral irradiance data derived from LBLRTM, a radiative transfer model from Verisk Atmospheric & Environmental Research (AER). This is a synthetic spectrum, simulated with a solar source function and database of atomic and molecular lines developed by Dr. Robert Kurucz[5]. The data are derived from LBLRTM 5.21 and are distributed here with permission from AER. This same dataset is used internally by libRadtran.

The spectrum represents top-of-atmosphere spectral irradiance at a Sun-Earth distance of 1 AU, sampled at 0.1-nm intervals and expressed in units of mW/m2/nm.

Note that these are legacy data. Newer, more accurate versions of the LBLRTM solar source function are available. The source code and derived products are subject to licensing terms dictated by AER.

Sun-Earth distance

Also provided are approximate Sun-Earth distances by day of year, expressed in AU. These are public-domain data provided by USGS/NASA Landsat Missions. The eccentricity of Earth's orbit produces variations in total irradiant flux on the order of ±3%. Note that the given distances are only approximate; ephemeris-based calculations may be used for more accurate results.

References

[1] Meyer P (1994). A parametric approach for the geocoding of airborne visible/infrared imaging spectrometer (AVIRIS) data in rugged terrain. Remote Sens Environ 49(2): 118–30. doi:10.1016/0034-4257(94)90048-5

[2] Gao B-C, Montes MJ & Davis CO (2004). Refinement of wavelength calibrations of hyperspectral imaging data using a spectrum-matching technique. Remote Sens Environ 90(4): 424–33. doi:10.1016/j.rse.2003.09.002

[3] Kaufman YJ & Gao B-C (1992). Remote sensing of water vapor in the near IR from EOS/MODIS. IEEE Trans Geosci Remote Sens 30(5): 871–84. doi:10.1109/36.175321

[4] Kaufman YJ, Wald AE, Remer LA, Gao B-C, Li R-R & Flynn L (1997). The MODIS 2.1-/spl mu/m channel-correlation with visible reflectance for use in remote sensing of aerosol. IEEE Trans Geosci Remote Sens 35(5): 1286–98. doi:10.1109/36.628795

[5] Kurucz RL (1992). Synthetic infrared spectra. In Rabin DM, Jefferies JT & Lindsey C (Eds.), Infrared solar physics: proceedings of the 154th symposium of the International Astronomical Union (pp. 523–31). Norwell, MA, USA: Kluwer Academic Publishers. Symposium held in Tucson, AZ, USA, March 2–6, 1992. doi:10.1017/S0074180900124805

About

Workflow for geometric and radiometric calibration of airborne imaging spectroscopy data, especially using HySpex devices.

Resources

License

Stars

Watchers

Forks

Releases

No releases published

Packages

No packages published

Languages