Imaging spectrometry in multiple spectral regions for geologic analysis

Integration of multi-wavelength-range imaging spectrometry data improves geologic mapping by combining complementary information from visible, near-IR, shortwave-IR, and longwave-IR spectral regions.
01 October 2015
Fred A. Kruse

Spectroscopy is a well-established method of identifying minerals and other geologic materials in the visible to near-IR (VNIR), shortwave-IR (SWIR), and longwave-IR (LWIR) spectral regions.1–4 This approach forms the basis for remote measurements using imaging spectrometry (also called hyperspectral imagery, or HSI), which comprises imagery in up to hundreds of contiguous spectral bands, with a spectrum at each image pixel.5 Hyperspectral data in the VNIR-SWIR has been available for more than 30 years, and analysis of it for geologic applications is considered mature.6 LWIR multispectral data has also been available for a long time7 and is often used in conjunction with VNIR-SWIR data.8–10 LWIR-HSI data, however, has been more difficult to obtain. It is only now becoming broadly available, and provides new capabilities that are just coming into routine use.11 Although typically analyzed in their separate wavelength ranges, there are some examples of using these individual HSI data modes together.12–14 However, fully integrated analysis is still being explored. The research we describe here is part of ongoing full-range HSI data analysis.15–19

Purchase SPIE Field Guide to Atmospheric OpticsOur objective was to improve accurate identification and mapping of geologic materials using the complementary information available from the ‘ full’ spectral range (0.4–12μm). Using data from NASA's Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) spanning the VNIR-SWIR, and spatially coincident Hyperspectral Thermal Emission Spectrometer (HyTES) data covering the LWIR, we applied a standardized approach20 to analyzing the HSI data. A model-based atmospheric correction converted AVIRIS VNIR and SWIR HSI radiance data to reflectance.21 The HyTES data was corrected using an LWIR atmospheric compensation model22 followed by temperature/emissivity separation.23 We extracted endmembers (key spectra) from each modal dataset and then performed partial spectral unmixing20 (mixture-tuned-matched filtering, or MTMF) on the individual spectral ranges (see Figures 13). AVIRIS directly mapped a variety of specific minerals, while HyTES separated and mapped several igneous rock phases. Some unknown materials require further field verification.

Figure 1. (a) Endmembers (key spectra) from the Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) in the visible to near-IR (VNIR) extracted from the data using the standardized approach. (b) AVIRIS VNIR mineral mapping results overlain on US Geological Survey (USGS) digital orthoquad (DOQ). MTMF: Mixture-tuned matched filtering. UTM: Universal Transverse Mercator.

Figure 2. (a) AVIRIS shortwave-IR (SWIR) endmembers extracted from the data using the standardized approach. (b) AVIRIS SWIR mineral mapping results overlain on USGS DOQ.

Figure 3. (a) Hyperspectral Thermal Emission Spectrometer (HyTES) longwave-IR (LWIR) endmembers extracted from the data using the standardized approach. (b) HyTES LWIR MTMF mapping results overlain on USGS DOQ.

The individual VNIR-SWIR-LWIR results were integrated using class merging, masking, and logical operators. We combined partial unmixing mapping results for individual spectral regions using geologic and spectroscopic directed constraints to reduce variability, reveal additional detail, and clarify associations. For example, close examination of the SWIR mineral spectra and mapping of the carbonate units calcite and dolomite reveals that AVIRIS had some difficulties identifying these materials in shadowed areas. Additionally, HyTES mapping appeared to be less sensitive in low-abundance carbonate areas. Results from both sensors were also negatively affected by spectral mixing between calcite and dolomite. Thus, in some locations materials had the SWIR carbonate spectral features, but no LWIR carbonate features, and vice versa. Combining the two MTMF classifications using a logical OR function (mapping locations where either AVIRIS or HyTES found carbonates) produced a more coherent map of the carbonate units that did not have discontinuities at terrain boundaries. A figure illustrating this concept and additional information are available in a recent paper published elsewhere.19

Taking a similar approach toward mapping individual silicate lithologies and other mineral endmembers, we used geologic context to combine similar materials and to generate mineral assemblages related to hydrothermal alteration (see Figure 4). The variability in similar spectral rock units was combined into individual discrete units using the logical OR operator. We identified alteration assemblages based on the joint occurrence of VNIR-determined minerals (iron oxides) and SWIR-determined muscovite in specific LWIR-defined rock units using a logical AND operator. Muscovite (sericite) alteration was also divided into three types based on joint occurrence of either hematite or goethite in VNIR mapping association. VNIR-mapped goethite areas are also present in SWIR-LWIR mapped carbonates, crosscutting calcite/dolomite boundaries, and possibly indicating additional alteration. We determined certain rock types and areas of skarn (calc-silicates) and silicification solely from their LWIR signatures, and field-verified these. Again, details are available in our recent paper.19

Figure 4. Integrated image map from directed combined analysis of individual VNIR-SWIR-LWIR material maps showing key alteration assemblages for the site overlain on USGS DOQ.

Independent analysis of the VNIR and SWIR AVIRIS and HyTES LWIR hyperspectral data for each spectral range produced complex, information-rich spatial patterns. It was possible to directly relate spectral endmembers and spatial material distributions to physiochemical properties and geologic processes based on specific reflectance absorption features and emissivity minima, and geologic context. However, the integrated mapping approach highlights subtle associations that are not easily observable by mapping the spectrally predominant material in individual ranges. Together, the VNIR-SWIR-LWIR combined modal mapping results lead to improved application-specific characterization and mapping that could not be accomplished using a single spectral range (see Figure 4). The advantage of this approach is that it is possible to tailor the analyses to specific materials of interest. However, this requires the application of subject matter expertise to the individual datasets to combine them in ways that make sense in terms of geologic (or other) context. We are now exploring other integrated analysis approaches for geology and other applications that take advantage of these diverse datasets, while minimizing requirements for additional outside information.15–19

Fred A. Kruse
Naval Postgraduate School
Monterey, CA

Fred A. Kruse received the PhD in geology from the Colorado School of Mines in 1987. He is currently a research professor, exploring the spatial distribution and spectral character of Earth surface materials. He has been involved in remote sensing research, applications, and teaching for more than 30 years.

1. G. R. Hunt, Spectral signatures of particulate minerals in the visible and near infrared, Geophysics 42(3), p. 501-513, 1977.
2. G. R. Hunt, R. P. Ashley, Spectra of altered rocks in the visible and near infrared, Econ. Geol. 74(7), p. 1613-1629, 1979.
3. J. W. Salisbury, L. S. Walter, Thermal infrared (2.5–13.5 micrometers) spectroscopic remote sensing of igneous rock types on particulate planetary surfaces, J. Geophys. Res. 94(B7), p. 9192-9202, 1989.
4. R. N. Clark, T. V. V. King, M. Klejwa, G. Swayze, N. Vergo, High spectral resolution reflectance spectroscopy of minerals, J. Geophys. Res. 95, p. 12653-12680, 1990.
5. A. F. H. Goetz, G. Vane, J. E. Solomon, B. N. Rock, Imaging spectrometry for earth remote sensing, Science 228, p. 1147-1153, 1985.
6. F. D. van der Meer, H. M. A. van der Werff, F. J. A. van Ruitenbeek, C. A. Hecker, W. H. Bakker, M. F. Noomen, M. van der Meijde, E. J. M. Carranza, J. B. de Smeth, T. Woldai, Multi- and hyperspectral geologic remote sensing: a review, Int'l J. Appl. Earth Observ. Geoinfo. 14, p. 112-128, 2012.
7. A. B. Kahle, and A. F. H. Goetz, Mineralogic information from a new airborne thermal infrared multispectral scanner, Science 222(4619), p. 24-27, 1983.
8. S. J. Hook, The combined use of multispectral remotely sensed data from the short wave infrared (SWIR) and thermal infrared (TIR) for lithological mapping and mineral exploration, Proc. Fifth Australas. Rem. Sens. Conf. 1, p. 371-380, 1990.
9. M. Abrams, E. Abbot, A. B. Kahle, Combined use of visible, reflected infrared, and thermal infrared images for mapping Hawaiian lava flows, J. Geophys. Res. 96(B1), p. 475-484, 1991.
10. L. C. Rowan, J. C. Mars, C. J. Simpson, Lithologic mapping of the Mordor, NT, Australia ultramafic complex by using the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER), Rem. Sens. Environ. 99, p. 105-126, 2005.
11. M. J. Abrams, S. Hook, M. C. Abrams, Imaging Spectrometry: Basic Principles and Prospective Applications, ch. Imaging spectrometry in the thermal infrared, p. 283-306, Springer, 2011.
12. T. J. Cudahy, J. Wilson, R. Hewson, P. Linton, P. Harris, M. Sears, K. Okada, J. A. Hackwell, Mapping porphyry-skarn alteration at Yerington, Nevada, using airborne hyperspectral VNIR-SWIR-TIR imaging data, Proc. Int'l Geosci. Rem. Sens. Symp. 2, p. 631-633, 2001.
13. T. Bowers, Analysis of VIS-LWIR hyperspectral image data for detailed geologic mapping, Proc. SPIE 4725, p. 116-127, 2002. doi:10.1117/12.478742
14. R. G. Vaughan, W. M. Calvin, Synthesis of high-spatial resolution hyperspectral VNIR/SWIR and TIR image data for mapping weathering and alteration minerals in Virginia City, Nevada, Proc. IEEE Geosci. Rem. Sens. Symp. 2, p. 1296-1299, 2004.
15. F. A. Kruse, M. L. McDowell, S. R. Cone, Combined analysis of visible-near infrared (VNIR), shortwave infrared (SWIR), and longwave infrared (LWIR) imaging spectrometer data, Proc. Hyperspectr. Imag. Sound. Environ Meet.; Imag. Appl. Opt. Tech. Dig., Opt. Soc. Am., 2015. Paper HW3B.1
16. M. L. McDowell, F. A. Kruse, Integrated visible to near infrared, short wave infrared, and long wave infrared spectral analysis for surface composition mapping near Mountain Pass, California, Proc. SPIE 9472, p. 94721C, 2015. doi:10.1117/12.2176871
17. S. R. Cone, F. A. Kruse, M. L. McDowell, Exploration of integrated visible to near-, shortwave-, and longwave-infrared (full range) hyperspectral data analysis, Proc. SPIE 9472, p. 94721D, 2015. doi:10.1117/12.2086670
18. F. A. Kruse, Comparative analysis of Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) and Hyperspectral Thermal Emission Spectrometer (HyTES) longwave infrared (LWIR) hyperspectral data for geologic mapping, Proc. SPIE 9472, p. 94721F, 2015. doi:10.1117/12.2176646
19. F. A. Kruse, Integrated visible and near infrared, shortwave infrared, and longwave infrared (VNIR-SWIR-LWIR), full-range hyperspectral data analysis for geologic mapping, J. Appl. Rem. Sens. 9, p. 096005, 2015. doi:10.1117/1.JRS.9.096005
20. J. W. Boardman, F. A. Kruse, Analysis of imaging spectrometer data using n-dimensional geometry and a mixture-tuned matched filtering (MTMF) approach, Trans. Geosci. Rem. Sens. 49(11), p. 4138-4152, 2011.
21. F. A. Kruse, Comparison of ATREM, ACORN, and FLAASH atmospheric corrections using low-altitude AVIRIS data of Boulder, Colorado, Proc. 13th JPL Airborne Geosci. Works., Jet Propul. Lab., 2004.
22. S. M. Adler-Golden, P. Conforti, M. Gagnon, P. Tremblay, M. Chamberland, Long-wave infrared surface reflectance spectra retrieved from Telops Hyper-Cam imagery, Proc. SPIE 9088, p. 90880U, 2014. doi:10.1117/12.2050446
23. P. S. Kealy, S. J. Hook, Separating temperature and emissivity in thermal infrared multispectral scanner data: implications for recovering land surface temperatures, IEEE Trans. Geosci. Rem. Sens. 31(6), p. 1155-1164, 1993.
Sign in to read the full article
Create a free SPIE account to get access to
premium articles and original research
Forgot your username?