Integration of Absorption Feature Information from Visible to Longwave Infrared Spectral Ranges for Mineral Mapping

Merging hyperspectral data from optical and thermal ranges allows a wider variety of minerals to be mapped and thus allows lithology to be mapped in a more complex way. In contrast, in most of the studies that have taken advantage of the data from the visible (VIS), near-infrared (NIR), shortwave infrared (SWIR) and longwave infrared (LWIR) spectral ranges, these different spectral ranges were analysed and interpreted separately. This limits the complexity of the final interpretation. In this study a presentation is made of how multiple absorption features, which are directly linked to the mineral composition and are present throughout the VIS, NIR, SWIR and LWIR ranges, can be automatically derived and, moreover, how these new datasets can be successfully used for mineral/lithology mapping. The biggest advantage of this approach is that it overcomes the issue of prior definition of endmembers, which is a requested routine employed in all widely used spectral mapping techniques. In this study, two different airborne image datasets were analysed, HyMap (VIS/NIR/SWIR image data) and Airborne Hyperspectral Scanner (AHS, LWIR image data). Both datasets were acquired over the Sokolov lignite open-cast mines in the Czech Republic. It is further demonstrated that even in this case, when the absorption feature information derived from multispectral LWIR data is integrated with the absorption feature information derived from hyperspectral VIS/NIR/SWIR data, an important improvement in terms of more complex mineral mapping is achieved.


•
automatically derived throughout the different spectral ranges • integrated and, moreover, if this new dataset can be successfully used for final mineral/ lithology mapping To map multiple absorption feature parameters automatically, a toolbox was used that was developed using Interactive Data Language (IDL).The biggest advantage of such an approach is that it allows the issue of prior definition of the endmembers to be overcome; this is a requested routine used for all widely-used spectral mapping techniques (e.g., Spectral angle mapping SAM [44], Spectral feature fitting SFF e.g., [45,46] and Spectral unmixing [47,48].Two different airborne image datasets were analysed, HyMap (HyVista Corp., Australian airborne imaging spectrometer, VIS/NIR/SWIR image data) and Airborne Hyperspectral Scanner (AHS, LWIR image data), both datasets were acquired over the Sokolov open-cast lignite mines in the Czech Republic.It is further demonstrated that even in this case, when the absorption feature information derived from multispectral LWIR data is integrated with the absorption feature information derived from hyperspectral VIS/NIR/SWIR data, it is an important contribution and improvement in terms of more complex mineral mapping.

Test Site
The study was performed in the Sokolov basin in the western part of the Czech Republic (Figure 1), in a region affected by long-term extensive lignite mining.The basement of the Sokolov Basin is formed of pre-Variscan and Variscan metamorphic complexes (recorded metamorphism from Devonian to Lower Carboniferous periods) of the Eger, Erzgebirge, Slavkov Forest, Thuring-Vogtland Crystalline Units and granitoids of the Karlovy Vary Pluton.The upper portions of these rocks are frequently weathered to kaolinitic residue.The basal late Eocene Staré Sedlo Formation is formed of well-sorted fluvial sandstones and conglomerates and is overlain by a volcano-sedimentary complex up to 350 m thick, which contains three lignite seams with variable sulphur (S) content.Long-term open cast mining required the removal of up to 180 m of thick overburden (Cypris clays), which was stockpiled and replaced after the lignite was extracted.At the dumps, the material consists mostly of weathered volcanic tuffs and Cypris clays, which can be characterised as well-laminated clays with a dominant kaolinite content; however, different varieties of mineralogical composition are common (e.g., the presence of montmorillonite, illite with admixtures of Ca-Mg-Fe carbonates, sulphates, sulphides, analcite, Mg-micas and bitumen [49].Due to the presence of S in the coal, both active and abandoned lignite mines are affected by acid mine drainage (AMD) [50,51].

Test Site
The study was performed in the Sokolov basin in the western part of the Czech Republic (Figure 1), in a region affected by long-term extensive lignite mining.The basement of the Sokolov Basin is formed of pre-Variscan and Variscan metamorphic complexes (recorded metamorphism from Devonian to Lower Carboniferous periods) of the Eger, Erzgebirge, Slavkov Forest, Thuring-Vogtland Crystalline Units and granitoids of the Karlovy Vary Pluton.The upper portions of these rocks are frequently weathered to kaolinitic residue.The basal late Eocene Staré Sedlo Formation is formed of well-sorted fluvial sandstones and conglomerates and is overlain by a volcanosedimentary complex up to 350 m thick, which contains three lignite seams with variable sulphur (S) content.Long-term open cast mining required the removal of up to 180 m of thick overburden (Cypris clays), which was stockpiled and replaced after the lignite was extracted.At the dumps, the material consists mostly of weathered volcanic tuffs and Cypris clays, which can be characterised as welllaminated clays with a dominant kaolinite content; however, different varieties of mineralogical composition are common (e.g., the presence of montmorillonite, illite with admixtures of Ca-Mg-Fe carbonates, sulphates, sulphides, analcite, Mg-micas and bitumen [49].Due to the presence of S in the coal, both active and abandoned lignite mines are affected by acid mine drainage (AMD) [50,51].Considering the Sokolov site, under various research projects (HypSo, EO-MINERS, DeMinTIR), numerous studies have been published demonstrating how hyperspectral imaging data can be utilised to quantitatively model the substrate pH [6,52], map mineral composition [41], estimate mine water pollution [53] and assess tree health [54][55][56].This study focused on the Lítov dump and the abandoned open pit called Medard (Figure 1), as no human activities were conducted between 2010 and 2011 at these two sites, the years when the two aerial image datasets, HyMap and AHS, were acquired, thus these two sites faced no changes regarding the relief and material transport between 2010 and 2011.

Airborne Imaging Datasets and Their Pre-Processing
Two different airborne image datasets acquired over the Sokolov basin were used in this study.The HyMap image data was acquired in 2010 (August 27) during the HyEUROPE 2010 flight campaign using the HyMap (HyVista Corp., HyVista Corporation Pty Ltd., Baulkham Hills, NSW, Australia) airborne imaging spectrometer.The HyMap sensor records image data in 126 narrow spectral bands covering the entire spectral interval between 0.450 and 2.480 µm of the spectral range with a Full Width Half Maximum (FWHM) of 15 nm and a ground field of view of 4 m.The resulting ground pixel resolution of the image datasets was 5 m.In order to successfully pre-process the hyperspectral data, supportive calibration and validation ground campaigns were organised simultaneously with the HyMap data acquisition.At the selected homogenous targets the ground measurements were acquired by an ASD FieldSpec-3 spectroradiometer to properly calibrate as well as validate the image data and to enable: (i) atmospheric correction of the airborne hyperspectral images and (ii) retrieval of surface reflectance values for further verification.The final atmospheric correction was performed in the ATCOR-4 software package [57] using the MODTRAN 4 physical model of the atmosphere [58].A detailed description of the HyMap data preprocessing can be found in Adar et al. [59].
The second image dataset was acquired by the Airborne Hyperspectral Scanner (AHS) in collaboration with the Spanish Aerospace Institute (INTA) as a set of day (19 July 2011) and night image data (22 July 2011).The AHS is an imaging 80-band line-scanner radiometer with 63 bands in the visible-near infrared (VNIR) and shortwave infrared (SWIR) regions, seven bands in the mid-wave infrared (MWIR) region and 10 bands in the longwave infrared (LWIR) region [60].However, due to cloud cover in the daytime image, it was only possible to use the cloud-free night-time LWIR data.These were acquired after a dry day, on a clear night with no precipitation.The flight lines were acquired at an altitude of 2 km above ground level, resulting in a 5-m pixel size (the same pixel size as the HyMap dataset).The temperature and emissivity were derived from the sensor radiance using the approach described in detail by Notesco et al. [42].
Both datasets were geo-corrected using on-board navigation information.After that both datasets were further georeferenced to the very high spatial resolution aerial orthophotos (pixel size = 0.5 m) achieving sub-pixel positional accuracy.To avoid any spatial misalignments, both datasets were resampled to a 10-m spatial resolution using nearest-neighbour resampling.

Ground Verification Data
Various soil/substrate samples were collected in the field in both years (2010 and 2011) and analysed with a Philips X'Pert X-ray Diffractometer (XRD) at the Czech Geological Survey to resolve their mineralogy.The X-ray powder diffraction patterns were obtained using monochromatic (CuKα) radiation and a graphite secondary monochromator.The whole-sample random patterns were collected in the angular range from 2 • to 70 • (2θ) with steps of 0.05 • (2θ).Oriented clay-fraction specimens (fraction < 2 µm) were prepared by a conventional sedimentation method [61].The oriented clay specimens were analysed after air-drying and after saturation for 10 h with ethylene-glycol vapour at 60 • C. Their diffraction data were acquired in the angular range of 2-50 • (2θ) with steps of 0.05 • (2θ).Mixed-layered minerals were identified by comparing the analysed XRD patterns of the ethylene-glycolated oriented clay fraction with the modelled XRD patterns obtained by NEWMOD code [61].

Absorption Wavelength Mapping
To map absorption features, namely absorption wavelengths and depths, new tools, which are described here for the very first time, were programmed to allow automatic detection of multiple absorption feature parameters.The tools-e.g., called QUANTools-have been created using IDL programming language (ENVI/IDL: version 5.0 and higher, [62]).
The tools can process both spectral libraries and hyperspectral image data.They consist of basic modules (GUIs) (Figure 2) allowing users to: 1.
Define a spectral range within the visible (VIS), near-infrared (NIR), shortwave infrared (SWIR) or thermal (TIR) spectral regions.Different spectral ranges can be defined and analysed consequently, one after another.2.
Employ Continuum Removal (CR)-A standard method to normalise the spectrum, to a departure from the norm [63].

3.
Detect bad spectral bands-A user can use a graphical interface to detect and correct bad (noisy) spectral bands.

4.
Define a number of desired absorption features to be detected within a set spectral range: the user can decide whether to detect an absolute absorption (the most pronounced one) or to define a number of multiple absorption features that can be identified within a set spectral range.

5.
Calculate absorption feature parameters (absorption wavelengths and depths): after correcting noisy bands, the trend of a spectral curve is analysed and saddle points-the local absorption maximum wavelengths (loc_max)-are detected and assigned to an image matrix.The detected absorptions are sorted in ascending order from shorter to longer wavelengths.Additionally, a corresponding absorption depth matrix is also calculated for each absorption feature.
The only decisions made by an operator/expert are to define the spectral region (Figure 2, step 1) and the number of desired absorption features mapped in each spectral region (Figure 2, step 4).All the other steps are done automatically when processing the datasets.A detailed explanation of the processing steps available in the toolbox is given in Section 2.3.2,below.

The Toolbox (QuanTools) Description
QuanTools can be used for spectral absorption band mapping using spectral libraries or hyperspectral image data (Figure 2).A user first defines the spectral range to be analysed within the visible (VIS), near-infrared (NIR), shortwave infrared (SWIR) or longwave infrared (LWIR) regions (Figure 2, step 1).Different spectral ranges can be defined and analysed after one another.Next comes continuum removal (CR, [63]), by which the continuum-that is, a convex hull of straight-line segments, is fitted over a reflectance spectrum and subsequently removed by division or rationing (Figure 2, step 2).The next step is noise detection and removal using local minima of the spectral curve defined as: where ρ CR (λ) and ρ CR (λ − 1) are the values of the spectral curve after Continuum Removal and R is an auxiliary variable.When the local minima are found, it has to be decided whether the value represents noise or not.The decision is based on two sets of values (a defined number of surrounding bands, e.g., called a spectral window)-derived from bands placed before and after the local minimum: where A(λ) represents the set of surrounding bands before the local minimum, B(λ) represents the set after the local minimum, and nn is the number of surrounding bands.Then the minimum values of these sets are detected as: The minimum of the spectral window (A min (λ), B min (λ)) is compared to the local minimum ρ CR (λ).Using this comparison, the noise is detected as follows: Err(λ + 2) . . .
where Err(λ) represents noise ("bad" bands).These error values are recalculated using the values of the surrounding bands (their number is defined by a spectral window size).For sample spectral data, a graphical interphase can be used to set up a spectral window size and check the corrected spectral curves (Figure 3, step 3).This allows the user to compare the results of different settings and finally employ a correction that is optimally tuned for the data under analysis.
After a noise-cleaned image is retained, the operator decides how many absorption features are to be derived (Figure 2, step 4) and the analysis is done automatically (Figure 2, step 5) using the following processing: Then the trend of the spectrum is analysed in a similar way to noise detection.First the saddle points (the local absorption maxima) are detected using these conditions: where R is an auxiliary variable that records increase (1) or decrease (0).During this step, local absorption maximum depths and wavelengths are registered: A set number of desired absorption features is detected as the most pronounced local absorption maxima (a set number of the most pronounced absorptions according to the absorption depths).If the desired number is larger than the number of local maxima detected, the final array is completed by zero values: Finally, the local absorption maximum features are assigned to image matrices.For z absorption features, two matrices with z bands are created; one has the absorption wavelengths assigned, while the second one has corresponding absorption depths assigned, respectively.
Remote Sens. 2017, 9, 1006 7 of 22 where R is an auxiliary variable that records increase (1) or decrease (0).During this step, local absorption maximum depths and wavelengths are registered: A set number of desired absorption features is detected as the most pronounced local absorption maxima (a set number of the most pronounced absorptions according to the absorption depths).If the desired number is larger than the number of local maxima detected, the final array is completed by zero values: Finally, the local absorption maximum features are assigned to image matrices.For z absorption features, two matrices with z bands are created; one has the absorption wavelengths assigned, while the second one has corresponding absorption depths assigned, respectively.

Specific Setting of the QUANTools
The main aim was to test if the major absorption features present though the VIS/NIR/SWIR and LWIR ranges can be (i) derived and integrated into one raster dataset; or (ii) successfully used for a final mineral mapping.The spectral ranges from which the absorption feature parameters were derived were defined as follows: VIS/NIR (HyMap): 0.450-1.200µm, SWIR (HyMap): 2.100-2.400µm and LWIR (AHS): 8.500-12.500µm.Noise detection and correction was employed to the VIS/NIR reflectance (0.450-1.200 µm) and the spectral window was set to 3 (three spectral bands before and after the local minimum were analysed).The number of absorption features to be mapped in each spectral region was set to 2. As a result, two new raster datasets were derived for each spectral range, one having assigned the two major absorption wavelengths and the second raster dataset with corresponding absorption depths.

Specific Setting of the QUANTools
The main aim was to test if the major absorption features present though the VIS/NIR/SWIR and LWIR ranges can be (i) derived and integrated into one raster dataset; or (ii) successfully used for a final mineral mapping.The spectral ranges from which the absorption feature parameters were derived were defined as follows: VIS/NIR (HyMap): 0.450-1.200µm, SWIR (HyMap): 2.100-2.400µm and LWIR (AHS): 8.500-12.500µm.Noise detection and correction was employed to the VIS/NIR reflectance (0.450-1.200 µm) and the spectral window was set to 3 (three spectral bands before and after the local minimum were analysed).The number of absorption features to be mapped in each spectral region was set to 2. As a result, two new raster datasets were derived for each spectral range, one having assigned the two major absorption wavelengths and the second raster dataset with corresponding absorption depths.After the wavelength mapping was employed to different spectral ranges, it was necessary to find a way to further integrate the absorption feature mapping results (Figure 4), more specifically the absorption wavelength and depth matrices.The MNF transformation [64] was employed to compress the data variability of these image matrices (six raster absorption wavelength matrices and six raster absorption depth matrices).
The intention was to test what the advantage will be of adding LWIR data (AHS) to the further mineral classification compared to using only the VIS/NIR/SWIR data (HyMap).Therefore, the MNF transformation was employed for two different scenarios: 1.
the MNF transformation was used so as to be comprised of only eight absorption wavelength/ depth matrices derived on the basis of the HyMap data (VIS/NIR: two absorption wavelength and two absorption depth matrices, SWIR: two absorption wavelength and two absorption depth matrices) 2.
the MNF was employed to comprise of all 12 absorption wavelength/depth matrices derived from both HyMap and AHS datasets (in addition to eight absorption wavelength and depth matrices derived from the HyMap data, two absorption wavelength and two absorption depth matrices derived for the AHS data were added).
For each scenario, the first three MNF bands were visually analysed further to identify training areas (ROIs) representing different material/surfaces.The ROIs were defined as representative pixels of different colour clusters/regions.These were easy to identify when MNF1/MNF2/MNF3 images were displayed as RGB (Figure 5).The ROIs were then used for further supervised classification.In this case the simple non-parametric supervised (parallelepiped) classification was employed [65][66][67], which uses a simple decision rule to classify multispectral data.The parallelepiped classifier uses the threshold of each class signature to determine if a given pixel falls within the class or not.The thresholds specify the dimensions (in standard deviation units) of each side of a parallelepiped surrounding the mean of the class in the feature space.The best results (e.g., class separability, lowest number of unclassified pixels) were achieved when the upper and lower limits of each parallelepiped were set to ±1.5 standard deviations.In the last step, the classifications were overlaid over the original HyMap and AHS image data and the average spectrum from the original HyMap and AHS image data was computed for each class to (i) ensure that each class represents different surface material (is represented by a unique spectral signature) and (ii) be able to describe the mineral composition of each class when interpreting the spectral property.After the wavelength mapping was employed to different spectral ranges, it was necessary to find a way to further integrate the absorption feature mapping results (Figure 4), more specifically the absorption wavelength and depth matrices.The MNF transformation [64] was employed to compress the data variability of these image matrices (six raster absorption wavelength matrices and six raster absorption depth matrices).
The intention was to test what the advantage will be of adding LWIR data (AHS) to the further mineral classification compared to using only the VIS/NIR/SWIR data (HyMap).Therefore, the MNF transformation was employed for two different scenarios: 1. the MNF transformation was used so as to be comprised of only eight absorption wavelength/depth matrices derived on the basis of the HyMap data (VIS/NIR: two absorption wavelength and two absorption depth matrices, SWIR: two absorption wavelength and two absorption depth matrices) 2. the MNF was employed to comprise of all 12 absorption wavelength/depth matrices derived from both HyMap and AHS datasets (in addition to eight absorption wavelength and depth matrices derived from the HyMap data, two absorption wavelength and two absorption depth matrices derived for the AHS data were added).
For each scenario, the first three MNF bands were visually analysed further to identify training areas (ROIs) representing different material/surfaces.The ROIs were defined as representative pixels of different colour clusters/regions.These were easy to identify when MNF1/MNF2/MNF3 images were displayed as RGB (Figure 5).The ROIs were then used for further supervised classification.In this case the simple non-parametric supervised (parallelepiped) classification was employed [65][66][67], which uses a simple decision rule to classify multispectral data.The parallelepiped classifier uses the threshold of each class signature to determine if a given pixel falls within the class or not.The thresholds specify the dimensions (in standard deviation units) of each side of a parallelepiped surrounding the mean of the class in the feature space.The best results (e.g., class separability, lowest number of unclassified pixels) were achieved when the upper and lower limits of each parallelepiped were set to ±1.5 standard deviations.In the last step, the classifications were overlaid over the original HyMap and AHS image data and the average spectrum from the original HyMap and AHS image data was computed for each class to (i) ensure that each class represents different surface material (is represented by a unique spectral signature) and (ii) be able to describe the mineral composition of each class when interpreting the spectral property.

Full-Range (VNIR, SWIR and LWIR) Absorption Wavelength Mapping and Further Classification
The intention was to detect the two most pronounced absorptions within the VIS/NIR, SWIR (both HyMap dataset) and LWIR (AHS dataset) regions.To validate if the wavelengths were detected correctly, the absorption wavelength positions were manually derived from 50 different image pixels of the HyMap and AHS data representing different material/surfaces, and they were compared with the absorption wavelength positions of corresponding pixels detected by QaunTools (Figure 6).Only a small deviation in the NIR range (around 0.900 µm) was detected due to the smoothing employed for the VIS/NIR spectral region prior to wavelength mapping.The wavelength positions of all the other absorption features in the SWIR and LWIR regions were detected without any shifts or errors.

Full-Range (VNIR, SWIR and LWIR) Absorption Wavelength Mapping and Further Classification
The intention was to detect the two most pronounced absorptions within the VIS/NIR, SWIR (both HyMap dataset) and LWIR (AHS dataset) regions.To validate if the wavelengths were detected correctly, the absorption wavelength positions were manually derived from 50 different image pixels of the HyMap and AHS data representing different material/surfaces, and they were compared with the absorption wavelength positions of corresponding pixels detected by QaunTools (Figure 6).Only a small deviation in the NIR range (around 0.900 µm) was detected due to the smoothing employed for the VIS/NIR spectral region prior to wavelength mapping.The wavelength positions of all the other absorption features in the SWIR and LWIR regions were detected without any shifts or errors.

Full-Range (VNIR, SWIR and LWIR) Absorption Wavelength Mapping and Further Classification
The intention was to detect the two most pronounced absorptions within the VIS/NIR, SWIR (both HyMap dataset) and LWIR (AHS dataset) regions.To validate if the wavelengths were detected correctly, the absorption wavelength positions were manually derived from 50 different image pixels of the HyMap and AHS data representing different material/surfaces, and they were compared with the absorption wavelength positions of corresponding pixels detected by QaunTools (Figure 6).Only a small deviation in the NIR range (around 0.900 µm) was detected due to the smoothing employed for the VIS/NIR spectral region prior to wavelength mapping.The wavelength positions of all the other absorption features in the SWIR and LWIR regions were detected without any shifts or errors.As described in the Methods an MNF transformation was employed to the absorption feature mapping results, first when using those derived from the HyMap data only (Scenario 1) and secondly to those derived from both the HyMap and the AHS data (Scenario 2).For each scenario, the ROIs were constructed using RGB colour compositions of the first three MNF bands and these ROIs were consequently used to employ a supervised parallelepiped classification.As a result, eight classes, which were the same for both classifications (Scenarios 1 and 2), were mapped (Figures 7B and 8B).In the case of Scenario 2, where absorption feature parameters derived from both datasets (HyMap and AHS) were used for the consequent mineral mapping, it was possible to map the addition of two classes (Figures 7C and 8C As described in the Methods an MNF transformation was employed to the absorption feature mapping results, first when using those derived from the HyMap data only (Scenario 1) and secondly to those derived from both the HyMap and the AHS data (Scenario 2).For each scenario, the ROIs were constructed using RGB colour compositions of the first three MNF bands and these ROIs were consequently used to employ a supervised parallelepiped classification.As a result, eight classes, which were the same for both classifications (Scenarios 1 and 2), were mapped (Figures 7B and 8B).In the case of Scenario 2, where absorption feature parameters derived from both datasets (HyMap and AHS) were used for the consequent mineral mapping, it was possible to map the addition of two classes (Figures 7C and 8C: *Class 2 and *Class 9).Figures 9 and 10 show the average class spectrum derived from the HyMap data and the AHS data, respectively.The spectral property of these classes is further discussed and linked with the mineralogy in the following text.

Linking the Spectral and Mineral Properties
In the VIS/NIR region, the average class spectra (Figure 9A) reflect the variations in absorption features characterising diverse iron oxy-hydroxides (secondary minerals with Fe 3+ ) and the organic component (organic C and lignite).Secondary minerals with Fe 3+ (hydroxysulfates and oxyhydroxides) exhibit absorption features around 0.500 µm and before 1.00 µm [29,69].The absorption wavelengths and intensities of the absorption features in this region depend on the nature of the crystal field around the Fe atom and on the nature of the bonds around it, because the nature of the magnetic coupling between the Fe 3+ ions (as influenced by the crystal field) facilitates the transition of electrons between energy states [70].In addition, shifts in the wavelength positions can also reflect material mixing, when the Fe 3+ secondary minerals coexist together or with different minerals and an organic component [6].When interpreting the average class spectra (Figure 9A), it can be seen that Classes 1-3, 5 and 7 exhibit the typical absorption features of secondary Fe 3+ -bearing minerals.The wavelength position of the first absorption shifts from 0.480 to 0.540 µm, while the second one varies between 0.900 and 1.080 nm.The shift of the second absorption indicates that Classes 1-3 represent the lithology where jarosite is present at higher amounts or coexist together with other Fe 3+ -bearing minerals, as jarosite exhibits maximum absorption closer to 0.900 µm (shorter wavelengths).On the other hand, the oxyhydroxides (e.g., goethite and hematite) have the second absorption centred around 1.000 µm (longer wavelengths), therefore this mineral is present at high amounts in Classes 4-7 and 10.The absorptions characteristic of organic C and lignite are at 0.550-0.580µm and 0.720-0.770µm, respectively [39], and these two are the most visible in the spectra of Classes 4, 8, 9 and 10.
The overtones and combinations of the fundamental OH and H-O-H vibrations can be mainly observed in the SWIR (Figure 9B).In general, the OH combination bands occur due to the two Al cations in the octahedral sites near 2.210 µm; the spectra of kaolinite shows a distinct absorption doublet at 2.170 and 2.210 µm, whereas muscovite has the main absorption at the shorter wavelength (2.200 µm).The absorption around 2.300 µm characterises carbonates [29]; however, no carbonates were identified at primary, secondary or accessory abundances by the XRD analyses conducted for numerous samples collected in Sokolov.Therefore, in this case, additional smaller absorptions at the longer wavelength (around 2.300 and 2.350 µm) are characteristic of the AlFe-OH and Fe 2 -OH combination bands in phyllosilicates [71].Kaolinite dominates the spectra of Classes 1-2 and 5-6, whereas Classes 3-4 and 7-10 represent mixtures between kaolinite and muscovite as they have less pronounced absorption at 2.170 µm and flat absorption with the same absorption intensities for the wavelength 2.190 and 2.210 µm.
In the LWIR region absorption features, resulting from fundamental molecular vibration modes, show additional information about mineral constituents, such as Si-bearing minerals (mainly quartz and clay minerals).Quartz, due to the molecular vibrations of the Si-O stretching (reststrahlen bands) displays a broad emissivity doublet in the TIR between 8 and 10 µm [72].In addition, it is possible to differentiate among diverse clay minerals (e.g., kaolinite, illite, and montmorillonite) using the LWIR range [34].In Figure 10A the emissivity of some silicates is displayed using the Arizona University Spectral library [68], the original emissivity is displayed together with the equivalent emissivity resampled to the spectral resolution of the AHS data.It shows that the AHS data allows quartz, kaolinite or quartz and kaolinite mixtures to be detected.In the case of the AHS data (Figure 10B), quartz affects the emissivity of the bands centred at 9.25 µm (low emissivity) and 9.75 µm (high emissivity) as well as the slope of the emissivity between these two bands.The distinct absorption features detectable by the AHS bands placed at 9.75 µm and 11.2 µm characterise kaolinite.Muscovite exhibits a wide absorption between 9 and 10 µm, and in general lowers the emissivity between 10-12 µm.It can be concluded that Classes 1 and 9 represent lithologies with a dominant quartz content, Classes 1, 2, 5, 6 and 8 then represent a lithology that has both kaolinite and quartz coexisting together.On the other hand, Classes 4, 7 and 10 predominantly have kaolinite.
The classifications (Figures 7 and 8) were compared to the field documentation and XRD analysis of the samples that fall spatially within each class was carried out.Table 1 shows the mineralogy (XRD analysis) that characterises six out of 10 defined classes (Classes 1-3, 6, 7 and 10)l for Classes 4, 5, 8 and 9 there were no samples (no XRD analysis) that would fall within their spatial extent.When comparing the mineralogy determined by the XRD with the class spectral properties discussed above, there is good agreement.Kaolinite and quartz represent the dominant minerals for Classes 1-3, whereas jarosite is present together with hematite at secondary abundances, as well as muscovite.However, Class 2 has quartz dominating and was described as a quartz-rich crust developed on the tuffs; on the other hand, Class 3 represents less weathered tuffs exposed by erosion.Class 6 represents the fresh clays of the Cypris formation, where kaolinite is the dominant mineral, followed by quartz and muscovite.Class 7 represents the material of backfill overburden, which was described as quartz-rich hard pack material with a clay matrix also containing lignite fragments.Class 10 represents weathered tuffs in which quartz, muscovite and kaolinite are the dominant, minerals whereas lignite and hematite are present at secondary abundances.Although there were no XRD analyses available for Classes 4, 5, 8 and 9, a mineral description (Table 1) was added; however, it should be emphasised that it is based on the interpretation of the class spectral property (Figures 9 and 10) and in this case it was not possible to differentiate between primary and secondary/accessory abundances of different mineral phases.

Discussion
Absorption feature parameters-wavelength position and depth-are the most essential information used in spectroscopy; to put it simply, the maximum absorption wavelength position defines what material it is, whereas absorption depth defines its relative abundance.Different absorption feature mapping/matching techniques have been used by the remote sensing community since hyperspectral data were made available.Among these the spectral feature fitting [73] or its improved version-multi-range spectral feature fitting (MRSFF)-are frequently used [74,75].From the recently developed toolboxes integrating absorption feature matching techniques the Tetracoder [76] or EnGeoMAP 2.0 toolbox [77] should be listed.These approaches are based on the common principle that they compare and statistically assess the fit of the image spectra to the reference spectra.The reference spectra, called endmembers, are scaled to match the image spectra, and can be either laboratory or field spectral measurements or can be extracted directly from the image (image endmembers).Basically, these classifiers require a routine whereby the endmembers (reference spectra) need to be defined prior to spectral mapping.Moreover, if laboratory or field spectra are used, the successful definition of such endmembers usually requires prior knowledge of the material composition and its spatial distribution within the area of interest.If image endmembers are to be used, and an expert with a background in spectroscopy is required to perform image analysis such as pixel purity analysis (PPI, [78,79]).
In contrast, the approach used here does not require a prior definition of the endmembers or any knowledge of the site conditions; this is a clear advantage.Instead, the absorption wavelength positions and depths the major absorption features present in different spectral ranges (VNIR/SWIR/LWIR) are extracted automatically and integrated into new raster datasets (multiple absorption feature wavelength and depth matrices).After compressing the main data variability of these new multi-band datasets, using, for instance, an MNF transformation, it is possible to assess the material spatial variability, define training areas (ROIs) and employ a supervised classification.The approach used here is thus unsupervised at the beginning, as no data or knowledge is required prior to absorption feature mapping.Throughout the automatic processing, there is a point where material variability is visualised in a spatial context; afterwards, the training areas (ROIs) for supervised classification are defined.The selection of the spectral ranges under analysis and the definition of the ROIs are the only expert-dependent parts of the analysis.It is thus recommended to divide spectral ranges into VIS/NIR, SWIR and LWIR regions where the diverse mineral groups exhibit distinct absorptions.The definition of the ROIs is rather intuitive as the expert interprets different colour clusters within the image.
Previously, techniques that utilise different types of interpolations for estimating absorption feature wavelength and depth were proposed.Van der Meer [80] proposed a simple linear interpolation technique in order to derive absorption-band position, depth and asymmetry from hyperspectral image.Such parameters were used to interpret the data in terms of the known alteration phases or to estimate heavy metal contents [81].Rodger et al. [82] proposed a simple quadratic method (SQM) to estimate the wavelengths of absorption features in the shortwave infrared (SWIR) spectral region.The SQM method was tested using spectral data convolved to four different instrument configurations differing in sampling regimes and spectral resolutions.The SQM method was found to estimate feature wavelengths within a reasonable accuracy and to perform well even in noisy environments.Ruitenbeek et al. [83] mapped the wavelength position of the deepest absorption features between 2.100 and 2.400 µm using a second-order polynomial fitting.They concluded that mapping the wavelength position of absorption features between 2.100 and 2.400 µm provided a new method for exploratory analysis of the surface mineralogy and that it would be particularly useful in areas where field validation is sparse and imagery contains shallow spectral absorption features.Moreover, Van der Meer et al. [84] tested two approaches-the 'Wavelength Mapper' [83] and the QUANTools [62] (also used for this study)-and demonstrated for the Rodalquilar epithermal system that deriving absorption feature characteristics, such as the wavelength position and the depth, can be directly linked to mineral type and abundance, and, even more, to subtle changes in mineral chemical composition.
However, as seen from the literature, this is the first demonstration of how the multiple absorption features, and their respective parameters, can be extracted automatically from different spectral ranges (VIS/NIR/SWIR and LWIR), and, furthermore, how, when using this rather simple approach, it is possible to successfully ingrate optical (VNIR and SWIR) and thermal (LWIR) spectral information, gathered in this case by two different sensors.It is also demonstrated how further integration can lead to more complex mineral classification.This is a crucial and highly relevant issue nowadays, when, in addition to the optical sensors, the sensors acquiring the data in the LWIR range are more often available and used (e.g., TASI, AISA OWEL) and even more for the future, when a hyperspectral satellite collecting spectral band through the VIS/NIR/SWIR and LWIR will be operating in orbit (e.g., HySPIRI).
So far the different spectral ranges are more frequently analysed and interpreted separately, without combining the full range of information (VIS/NIR/SWIR/LWIR) into a single integrated data product.Too few approaches have been proposed allowing for real optical and thermal data integration; still, they require an expert decision to be made prior to spectral mapping (e.g., endmember definition).As previously explained, the approach used here does not require any prior definition of the endmembers; moreover, there is no need for any knowledge or ground data from the site under study.
In the approach used here, the level of noise present in the image datasets can be a limitation.However, in an example of HyMap VIS/NIR data, it was demonstrated that this problem can be minimised by employing spectral smoothing and that it can be tailored specifically to the level of noise present in the data.In addition, there might be other variables that may create false alarms, such as varying topography causing brightness differences across the image produced by shadow and slope variations, varying mineral particle size and soil moisture content.QUANTools process the data in such a way that at the beginning they are normalised by employing the continuum removal method, which helps to remove the effect of scattering.Then the absorption feature wavelengths are mapped and, as long as the noisy/false absorptions are eliminated, the absorption wavelength positions do not tend to change spectral locations as regards shadows or slope variations.The other absorption feature parameter-absorption depth-reflects the material quantity but is also sensitive to the sizes of mineral particles.Therefore, varying grain size is also part of a final classification.It is thus important to decide if the grain size is also a criterion for material mapping or not.If only the chemical composition of the targets is requested, then it is possible to use only absorption feature wavelength matrices in a classification.In this study, 2010 HyMap and 2011 AHS image data were used to classify two sites, Lítov and Medard, which faced no changes regarding the relief and material transport between 2010 and 2011; however, some differences regarding changes in microtopography or moisture content may have been present.Both absorption feature parameters-wavelength and depth-were used for classifications and brought results that were in good agreement with the XRD analysis, showing that, at a general level, this approach allows diverse minerals to be mapped, including ones that exhibit multiple absorption features though the VIS/NIR/SWIR and LWIR ranges.There are also other useful absorption feature parameters, such as shoulder positions, symmetry or width, that can be used for gathering the information on minor material components, featureless parameters (e.g., heavy metals) or chemical processes/changes.Mapping these additional absorption feature parameters will be a subject for future QUNATools development.
Using this approach, it was possible to integrate the absorption feature information derived from the VIS/NIR/SWIR regions (HyMap data), together with the absorption feature information derived from the LWIR region (AHS data).This integration led to a mineral classification that differentiated between the presence and abundance of diverse Fe 3+ -bearing minerals and phyllosilicates as well as lignite and quartz contents.In addition to HyMap reflectance, the AHS emissivity data allowed a better discrimination between a quartz-dominating crust and substrates (Classes 2 and 9) from the other classes where quartz did not have such dominating abundances and was present together with other mineral phases at primary or secondary abundance (e.g., Classes 1 and 7).

Conclusions
In this study, it was demonstrated how the multiple absorption features, respectively their parameters such as wavelengths and depths, can be automatically extracted from different spectral ranges and further integrated into one raster dataset.The absorption feature information gathered from the different spectral-ranges of the two different sensors (HyMap and AHS) was integrated and led to a mineral classification that differentiated between the diverse Fe 3+ -bearing minerals and phyllosilicates as well as lignite and quartz contents and overall resulted in a more complex mineral/lithology classification.This is a crucial and highly relevant issue nowadays, when, in addition to optical sensors, sensors acquiring data in the LWIR range are available and more often used (e.g., TASI, AISA OWEL); and even more for the future, when a hyperspectral satellite collecting a spectral band through the VIS/NIR/SWIR and LWIR will be operating in orbit (e.g., HySPIRI).

Figure 1 .
Figure 1.Geographic position of the two sites under study: Lítov dump and Medard Lake, Sokolov basin, Czech Republic.

Figure 1 .
Figure 1.Geographic position of the two sites under study: Lítov dump and Medard Lake, Sokolov basin, Czech Republic.

Figure 3 .
Figure 3. Bad band detection (spectral range 0.450-1.200µm) in QuanTools where the bad bands detected are shown by the green square points and the corrected curve as a dashed line (set spectral window: three neighbouring bands).

Figure 3 .
Figure 3. Bad band detection (spectral range 0.450-1.200µm) in QuanTools where the bad bands detected are shown by the green square points and the corrected curve as a dashed line (set spectral window: three neighbouring bands).

Figure 4 .
Figure 4. Scheme of the processing allowing integration of the absorption feature information detected in VIS/NIR/SWIR and LWIR data and further classification.

Figure 4 .
Figure 4. Scheme of the processing allowing integration of the absorption feature information detected in VIS/NIR/SWIR and LWIR data and further classification.

Figure 5 .
Figure 5. Lítov: MNF1/MNF2/MNF3 images displayed as RGB: (A) Scenario 1 (absorption feature parameters derived only from the HyMap data were used); (B) Scenario 2 (absorption feature parameters derived from both datasets, HyMap and AHS, were used), the MNF bands were further visually analysed to identify training areas (ROIs) representing different material/surfaces.

Figure 6 .
Figure 6.The manually derived absorption wavelength positions from 50 different image pixels of the the HyMap and the AHS image data representing different material/surfaces compared to the absorption wavelength positions derived automatically when using QaunTools.

Figure 5 .
Figure 5. Lítov: MNF1/MNF2/MNF3 images displayed as RGB: (A) Scenario 1 (absorption feature parameters derived only from the HyMap data were used); (B) Scenario 2 (absorption feature parameters derived from both datasets, HyMap and AHS, were used), the MNF bands were further visually analysed to identify training areas (ROIs) representing different material/surfaces.

Figure 5 .
Figure 5. Lítov: MNF1/MNF2/MNF3 images displayed as RGB: (A) Scenario 1 (absorption feature parameters derived only from the HyMap data were used); (B) Scenario 2 (absorption feature parameters derived from both datasets, HyMap and AHS, were used), the MNF bands were further visually analysed to identify training areas (ROIs) representing different material/surfaces.

Figure 6 .
Figure 6.The manually derived absorption wavelength positions from 50 different image pixels of the the HyMap and the AHS image data representing different material/surfaces compared to the absorption wavelength positions derived automatically when using QaunTools.

Figure 6 .
Figure 6.The manually derived absorption wavelength positions from 50 different image pixels of the the HyMap and the AHS image data representing different material/surfaces compared to the absorption wavelength positions derived automatically when using QaunTools.
: *Class 2 and *Class 9).Figures 9 and 10 show the average class spectrum derived from the HyMap data and the AHS data, respectively.The spectral property of these classes is further discussed and linked with the mineralogy in the following text.Remote Sens. 2017, 9, 1006 11 of 22

Figure 7 .
Figure 7. Lítov, mineral mapping: (A) orthophoto showing the Lítov dump; (B) classification using Scenario 1 (absorption feature parameters derived only from the HyMap data were used for the consequent mineral mapping), (C) classification using Scenario 2 (absorption feature parameters derived from both datasets, HyMap and AHS, were used for the consequent mineral mapping), (D) enlargement of the area of interest.The mapped classes correspond to the mineral classes inTable1, *Class 2 and *Class 9: two additional classes mapped when using Scenario 2.

Figure 7 .
Figure 7. Lítov, mineral mapping: (A) orthophoto showing the Lítov dump; (B) classification using Scenario 1 (absorption feature parameters derived only from the HyMap data were used for the consequent mineral mapping), (C) classification using Scenario 2 (absorption feature parameters derived from both datasets, HyMap and AHS, were used for the consequent mineral mapping), (D) enlargement of the area of interest.The mapped classes correspond to the mineral classes inTable1, *Class 2 and *Class 9: two additional classes mapped when using Scenario 2.

Figure 8 .
Figure 8. Medard Lake, mineral mapping: (A) orthophoto showing the Lítov dump, (B) classification using Scenario 1 (absorption feature parameters derived only from the HyMap data were used for the consequent mineral mapping), (C) classification using Scenario 2 (absorption feature parameters derived from both datasets, HyMap and AHS, were used for the consequent mineral mapping), (D) enlargement of the area of interest.The mapped classes correspond to the mineral classes in the Table 1, *Class 2 and *Class 9: two additional classes mapped when using Scenario 2.

Figure 8 . 22 Figure 8 .
Figure 8. Medard Lake, mineral mapping: (A) orthophoto showing the Lítov dump, (B) classification using Scenario 1 (absorption feature parameters derived only from the HyMap data were used for the consequent mineral mapping), (C) classification using Scenario 2 (absorption feature parameters derived from both datasets, HyMap and AHS, were used for the consequent mineral mapping), (D) enlargement of the area of interest.The mapped classes correspond to the mineral classes in the Table 1, *Class 2 and *Class 9: two additional classes mapped when using Scenario 2.

Figure 10 .
Figure 10.(A) Emissivity of some silicates is displayed using the Arizona University Spectral library [68], the original spectra are displayed together with the equivalent spectra resampled to the spectral resolution of the AHS data.(B) The average class spectrum derived from the AHS data.The mapped classes correspond to the mineral classes inTable1, *Class 2 and *Class 9: two additional classes mapped when using Scenario 2.

Figure 10 .
Figure 10.(A) Emissivity of some silicates is displayed using the Arizona University Spectral library [68], the original spectra are displayed together with the equivalent spectra resampled to the spectral resolution of the AHS data.(B) The average class spectrum derived from the AHS data.The mapped classes correspond to the mineral classes inTable1, *Class 2 and *Class 9: two additional classes mapped when using Scenario 2.

Figure 10 .
Figure 10.(A) Emissivity of some silicates is displayed using the Arizona University Spectral library [68], the original spectra are displayed together with the equivalent spectra resampled to the spectral resolution of the AHS data.(B) The average class spectrum derived from the AHS data.The mapped classes correspond to the mineral classes inTable1, *Class 2 and *Class 9: two additional classes mapped when using Scenario 2.
Integration of Absorption Feature Information Detected In VIS/NIR/SWIR and LWIR Data and Further Classification

Table 1 .
Map classes from Figures 7-10 compared to the XRD analysis and field documentation of the samples that fall spatially within each class.For Classes 4, 5, 8 and 9 there were no XRD analyses available; however, the mineral description based on the interpretation of the spectral property was added.