1. Introduction
WorldView imaging satellites comprise a planned constellation of commercial, orbiting platforms developed by DigitalGlobe Inc. (Longmont, CO, USA) and built by Ball Aerospace & Technologies (Boulder, CO, USA). Financed partially by the National Geospatial-Intelligence Agency (NGA), WorldView satellite systems are launched from Vandenberg Air Force Base on Delta II rockets. WorldView-1 (WV-1) was launched in 2007 with a panchromatic (PAN) imaging system capable of 50-cm spatial resolution [
1]. With no multispectral bands on board, the primary purpose of this single-band PAN system was to rapidly collect high-spatial resolution (or hi-res) imagery, especially suited for generating detailed digital elevation model (DEM) data. The major improvement of WV-1 over other hi-res commercial systems at the time included state-of-the-art geo-locational capability, improved accuracy, and flexibility in programming and tasking that vastly improved the collection of 2-pass data necessary for DEM generation. WV-1 allowed for a 1.7-day revisit time, acquiring 750,000 square kilometers of imagery per day. The next big milestone was encountered with the launch of WorldView-2 (WV-2) in 2009 offering hi-res PAN data at 46-cm pixel size plus visible and near-infrared (VNIR) bands at 1.85-m spatial resolution [
2]. WV-2 was the first system to collect eight hi-res multispectral bands ranging in wavelengths from 0.4 to 1.04 μm, selected primarily for vegetation, coastal, and land-use applications. DigitalGlobe plans to launch WorldView-3 (WV-3) in 2014 offering a similar PAN and VNIR band set, complemented by eight short-wave infrared (SWIR) bands ranging from 1.21 to 2.365 μm (
Table 1) [
3,
4]. It will also have an additional 12 bands included for atmospheric compensation [
3].
Figure 1 shows the proposed WV-3 SWIR band passes. WV-3, with approximately 3.7 m SWIR spatial resolution, will be the first commercial system with both VNIR and SWIR multispectral capability. The SWIR data at 7.5 m spatial resolution will be available for purchase. Commercial pricing has not yet been determined. It is expected that the inclusion of the new SWIR bands will significantly impact surface compositional modeling, and mapping of rock and soil exposures worldwide. Potential applications include improved geologic mapping; environmental/soil/disaster monitoring; exploration for petroleum, minerals, and geothermal resources; as well as other non-renewable resource assessments (especially suited for developing countries). This paper summarizes mineral mapping results achieved by modeling the proposed WV-3 SWIR spectral bands shown in
Figure 1 using the Airborne Visible/Infrared Imaging Spectrometer (AVIRIS), extraction of SWIR spectral signatures and spectral mapping, and comparison to results from NASA’s on-orbit Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) multispectral sensor, which has 6 SWIR bands at 30 m resolution [
5,
6]. This simulation and analysis focuses on the WV-3 SWIR bands, as this is the unique added capability, the VNIR and SWIR measure different physical phenomena, and most spectral features for alteration minerals occur in the SWIR. Combined analysis of the SWIR bands with the WV-3 VNIR bands may also offer additional capabilities, but that evaluation is beyond the scope of this initial investigation.
3. Mineral Mapping
Mineral reflectance spectra of key minerals known to occur at Cuprite at specific locations based on spectral libraries and previous research [
7–
14] were extracted as regions of interest (ROIs) from the AVIRIS, ASTER, and simulated WV-3 data for use in spectral classification and mineral mapping.
Figure 3 shows a comparison of the extracted mean spectra. The AVIRIS spectra are comparable to laboratory-measured library spectra and illustrate the full spectral character of the alteration minerals. The multispectral ASTER data and simulated WV-3 spectra show the discrete sampling of key spectral absorption features by both ASTER (6 bands) and WV-3 (8 bands), generally resulting in good separation for common alteration minerals, however, some ambiguities are apparent. ASTER has demonstrated the viability of detailed multispectral mineral mapping from space [
6,
14,
17–
19]. We expect that WV-3 SWIR will perform similarly, with the added advantage of higher spatial resolution and additional bands improving mineral exploration, mine characterization, and spectral mapping.
The ROIs used to derive the endmember spectra shown in
Figure 3 were used to run a Minimum distance (supervised) classification on the 86 band AVIRIS SWIR data and the simulated WV-3 SWIR data to demonstrate typical information content available using “standard” multispectral classification approaches (
Figure 4). The Minimum Distance classification uses the mean vectors of each class in multi-dimensional space and calculates the Euclidean distance from each image pixel vector to the mean vector for each class [
20]. A standard deviation threshold was used to exclude highly dissimilar pixels. Comparison of the WV-3 results to the AVIRIS data using a confusion matrix approach with the AVIRIS Minimum Distance classification acting as surrogate ground truth shows a excellent match (overall accuracy of 93.34% with a Kappa Coefficient of 0.9081) (
Tables 2 and
3) [
21,
22]. Comparison of the diagonals of each matrix indicates the observed agreement between the two maps for each mineral, with a score of 100% being perfect agreement. The Kappa Coefficient, a statistical measure of agreement taking into account all of the categories, has a value close to zero when the observed agreement is the same as expected by chance and a value approaching one with perfect agreement [
21]. The classifications and the confusion matrix for the Minimum Distance classifications demonstrate that comparable results can be generated from the two datasets using multispectral approaches, but don’t really assess the accuracy of those results with respect to the full mineral information available from the HSI data based on a physical mixing model. Comparison to AVIRIS mineral maps derived using HSI analysis approaches based on spectral matching and unmixing illustrates that there are significant errors in the Minimum Distance classifications (
Figure 4). One of these approaches, the Mixture-Tuned-Matched-Filtering (MTMF) spectral mapping approach [
23] was used to map the spatial distribution of the spectrally predominant mineralogy at the Cuprite, NV site and demonstrates how the WV-3 simulate data should perform for more specific mineral mapping approaches (
Figure 4(Right);
Figure 5). Thus, even though the Minimum Distance classification results are consistent and show high correlation between the AVIRIS and simulated WV-3 data, they do not provide a true measure of WV-3’s viability for mineral mapping using physics-based approaches incorporating physical and mathematical models tied to spectral signatures and mineral mixing.
Hyperspectral data are more typically analyzed using spectral matching or spectral unmixing methods, which take advantage of the full spectral signature for each material and a physical model to account for within pixel mixing of multiple materials. Mixture-Tuned-Matched-Filtering (MTMF) is a method for locating a known signature in the presence of a mixed or unknown background by combining a linear detector and spectral unmixing [
23]. It is ideally suited for use where materials with distinct spectral signatures occur within a single pixel, and allows both identification and abundance estimation. Spectral mixing has been previously observed and reported at the Cuprite site [
13]. The MTMF image-maps produced for this simulation
Figure 4(Right) and
Figure 5(Left) generally correspond to previous mineral mapping results using a variety of multispectral and hyperspectral sensors [
11–
14], however, as no detailed ground truth map actually exists for the site (other than those generated using hyperspectral data), it is difficult to make absolute accuracy comparisons.
For the purposes of this research, the AVIRIS mineral map (
Figure 4(Right) and
Figure 5(Left)) was selected for use as surrogate ground truth, and ASTER data for the Cuprite site and the simulated WV-3 shows spatial patterns for specific minerals that are similar to AVIRIS mineral mapping results; however, in addition to a more noisy appearance likely caused by mineral mixing at the ASTER 30 m pixel, there are some other rather obvious problems with the multispectral mapping results (errors of commission). The simulated WV-3 MTMF results (
Figure 5(Right)) appear to more accurately replicate the AVIRIS mineral patterns for specific minerals, but there are again some obvious errors of omission and commission. In both cases, the multispectral sensors appear to do surprisingly well in defining the mineralogy at the site using the SWIR signatures and the hyperspectral analysis approach, despite having only a few spectral bands compared to the AVIRIS data. Unfortunately, however, it is also clear that the subtle differences between the multispectral signatures, particularly for similar minerals like kaolinite, alunite, muscovite, and silica (all with features near 2.2 μm) result in confusion between similar minerals mapped in the WV-3 simulated data. In addition, the ASTER results are also likely more affected by linear (areal) spectral mixing because of their larger 30 m pixels.
Additional quantitative comparisons of MTMF mineral maps derived from the ASTER data and analysis of the WV-3 simulated data compared to the AVIRIS mineral mapping result were accomplished using the confusion matrix approach (
Table 4). There are significant errors of commission and omission in all cases because of lower spectral resolution and for ASTER, also spectral mixing. It is significant to note that while the statistical measures show problems with the analysis on a per-pixel basis, that visual inspection of the image classes (
Figure 5) indicates that the WV-3 data actually do well at reproducing the main mineral patterns mapped at Cuprite using the AVIRIS data.
The confusion matrices for these data demonstrate that the per-pixel match between the mineral AVIRIS and WV-3 mineral mapping results is only poor to fair using Kappa Coefficient criteria [
21,
22] (
Tables 3 and
4). We note, however, that significant portions of the error are attributable to unclassified areas (
Table 4), meaning that rather than misidentification of materials, that many of the errors are ones where no identification was possible using the multispectral data. This is not an unexpected result, however, as multispectral data are inherently less well suited to mineral mapping than HSI data because they don’t resolve all the spectral features. Improved performance is demonstrated by excluding unclassified pixels from the error analysis (
Table 4). Further examination of the ASTER confusion matrix (
Table 4) and extracted spectra (
Figure 3) shows that performance is higher for calcite and buddingtonite than other minerals because they have spectral features at unique locations (near 2.34 and 2.11 μm respectively) with respect to the other minerals. Confusion between minerals with similar features near 2.2 μm such as alunite, kaolinite, muscovite, and silica is high. This is the result of confusion between similar minerals at lower multispectral spectral resolution, as well as spectral mixing effects. For the simulated WV-3 data, examination of the confusion matrices shows that performance is higher for calcite and buddingtonite; however, there is still confusion between minerals with similar features such as alunite, kaolinite, muscovite, and silica, which at high spectral resolution all have characteristic spectral features near 2.2 μm (
Table 4,
Figure 3).
In both the ASTER and WV-3 cases, ambiguities typically result in “over classification” using the multispectral systems compared to more specific AVIRIS results. That is, minerals are often identified in the MSI data where no specific mineral dominates as indicted by the HSI data analysis. This is mostly likely the result, however, of multiple similar mineral signatures combining to resemble a given mineral due to the reduced MSI spectral resolution. The HSI data are more sensitive to individual minerals in mineral mixtures. The WV-3 multispectral data appear to “fill in” some of the unmapped areas on AVIRIS (these represent additional ambiguities in multispectral images). It is significant to note that both multispectral datasets still produce spatially coherent highly detailed mineral maps that visually match the main spatial patterns in the AVIRIS data, and that the simulated WV-3 data perform significantly better than the ASTER data.
Previous work has demonstrated that that mapping errors caused by similar spectral features can be mitigated somewhat by grouping similar minerals for analysis [
6]. Accordingly, we combined some of the simulated WV-3 results for similar minerals and performed confusion matrix analyses on the combined mineral maps (
Figure 6)
versus the same minerals combined for AVIRIS (not shown). Improved performance is observed in all cases according to the statistics (
Table 5); however, comparing the images with the original AVIRIS mineral map in
Figure 5(Right) demonstrates that significant important information has been lost. The key minerals shown in the AVIRIS map in
Figure 5 are combined for improved per-pixel statistical accuracy, but information about the overall spatial pattern is reduced.
Figure 7 provides a final illustration of the spatial resolution advantage of WV-3
versus ASTER.
Based on these results, and demonstrated real-world ASTER mineral mapping successes, we expect that Worldview-3 will also be successful for mineral mapping using the new SWIR bands and that results will be improved over those determined from ASTER data. Clearly the WV-3 data will provide much sharper images of areas of interest, leading to improved literal analysis and interpretation. In addition, the mineral maps of the simulated WV-3 high spatial resolution data illustrate the advantage of reduced spectral mixing in the higher spatial resolution WV-3 data, leading to more spatially coherent mapping of key minerals. The patterns in the simulated WV-3 data more closely match those of the hyperspectral AVIRIS data compared to the ASTER results, and key mineral patterns are preserved, even where error analysis statistics indicate only fair individual pixel mineral matching performance.
4. Conclusions
This study presents the results of simulated multispectral Worldview-3 (WV-3) data analysis for Short Wave Infrared (SWIR) mineral mapping for hydrothermal ore deposits. WV-3 is a new multispectral remote sensing system to be launched during 2014 by DigitalGlobe, and will provide combined high (7.5 m) spatial resolution and SWIR spectral bands not available from any other source. Because WV-3 is the first of its kind to offer these added SWIR capabilities at high spatial resolution, the interest (particularly from the commercial mineral exploration community) is expected to be extremely high. This is in many ways the sensor that this community has been waiting for—high spatial resolution, global coverage, and increased spectral mapping capabilities. Another multispectral imager, the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER), which has 30 m spatial resolution and had 6 spectral bands in the SWIR (these failed in 2008) was extensively used by the mining industry. The launch of WV-3 is highly anticipated and we expect that the new capabilities of WV-3 will be even more valuable.
Based on these observations and the pending WV-3 launch, simulated SWIR WV-3 imagery was produced utilizing high spatial and spectral resolution Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) (hyperspectral or HSI) data and the proposed WV-3 spatial and spectral responses to validate and quantify the likely SWIR mapping performance, specifically for hydrothermal alteration minerals. AVIRIS data for a well-known altered area (Cuprite, NV, USA) were spectrally and spatially resampled to the WV-3 specifications, and spectral mapping methods were applied to detecting and mapping minerals known to occur for the site. Most multispectral data are analyzed using statistical methods; results shown here for a typical classification method (Minimum Distance supervised classification) illustrate that high consistency can be achieved between the multispectral and hyperspectral classifications using statistical methods. In fact, total agreement was greater than 90% between these two datasets when using this approach. Unfortunately, however, despite this high agreement, spectral signatures and previous mineral mapping demonstrate that the statistical methods do not actually accurately map the known mineral distributions, mainly because they do not take into account the physical mixing model (multiple materials linearly mixing within each pixel).
The idea of utilizing spectral signatures and hyperspectral analysis approaches for multispectral mineral mapping represents a unique approach. Only the 8 well-placed SWIR bands proposed for WV-3 make this worth considering. Confusion matrix analyses comparing the simulated WV-3 mapping results to full-spectral-resolution AVIRIS data indicate that this approach was only partially successful, however, principally because of the similarity of many minerals when using only a few multispectral bands (versus HSI, with up to hundreds of spectral bands). The simulated WV-3 data allowed extraction and mapping of specific mineral spectra; however, spectra for certain similar minerals are only subtly different. These include a broad range of minerals with absorption features near 2.2 μm (such as kaolinite, alunite, muscovite, and silica), which are easily confused in the simulated WV-3 data. Simulated WV-3 performance for specific minerals ranged from approximately 40%–60% match for minerals with similar absorption features to 75%–85% match for minerals with distinctly different absorption features. Even so, standard HSI analysis methods were successfully applied to the WV-3 data to produce detailed mineral maps that separated some of these materials. The AVIRIS data produce a more conservative (and accurate) mineral map taking mineral mixing into consideration, however, the WV-3 data do capture the known alteration patterns at the Cuprite, Nevada site. Additional error analysis demonstrates that improved performance in terms of the classification statistics can be achieved by combining mineral classes with similar spectral signatures. When similar minerals are grouped together, matches between the HSI and WV-3 data improve to better than approximately 75% match. While this does improve the overall accuracy, the images, however, clearly show that there is valuable mineralogic information in the mineral maps produced using the full mineral suite that cannot be replicated in the reduced endmember case. This points towards the need to include additional, narrower multispectral SWIR bands in future sensors.
These simulation results demonstrate that WV-3 will provide significant new capabilities for multispectral mineral mapping from space. The simulated WV-3 data produce similar low spectral resolution spectra to those extracted from reflectance-corrected ASTER data; however WV-3 will significantly outperform ASTER for mineral mapping because of additional SWIR bands and the improvement in spatial resolution to 7.5 m from ASTER’s 30 m pixel size. The limited release (government only) 3.7 m spatial resolution data will be even more capable. The WV-3 simulation does demonstrate reduced capabilities compared to AVIRIS, however, manifested as misclassification and “smoothing” of the spectral mapping because of confusion between similar minerals. This is not unexpected (hyperspectral data are inherently more capable), and despite this fact, the 7.5 m spatial resolution and the simulated spectral resolution appear to generally be well suited to resolving typical mineral occurrences in a hydrothermally altered area. The WV-3 spectral bands should permit identification and mapping of key minerals, while the improved spatial resolution will provide improved discrimination of complex alteration mineral patterns. The WV-3 results are promising and indicate that this sensor will be a significant tool for geologic remote sensing. Additional research is required for additional sites with different mineral distributions. There may also be some benefit to combined analysis of the eight (8) WV-3 visible/near infrared (VNIR) bands with the SWIR bands. This will require further research and case histories.