Integration of Hyperspectral Shortwave and Longwave Infrared Remote-Sensing Data for Mineral Mapping of Makhtesh Ramon in Israel

Hyperspectral remote-sensing in the reflected infrared and thermal infrared regions offers a unique and efficient alternative for mineral mapping, as most minerals exhibit spectral features in these regions, mainly in the shortwave and longwave infrared. Airborne hyperspectral data in both spectral regions, acquired with the AisaFENIX and AisaOWL (Specim) sensors over Makhtesh Ramon in Israel, were analyzed. Calculating the reflectance and emissivity spectra of each pixel in the shortwave infrared and longwave infrared region images, respectively, and determining mineral indices enabled identifying the dominant minerals in this area—kaolinite, calcite, dolomite, quartz, feldspars and gypsum—and mapping their spatial distribution in the surface. The benefit of using hyperspectral data from both reflected infrared and thermal infrared regions to improve mineral identification was demonstrated.


Introduction
Remote-sensing techniques, using spaceborne and airborne sensors in the infrared region, offer an efficient method for mineral mapping (e.g., [1][2][3][4]), as most minerals exhibit spectral features in the reflected infrared (RIR)-mainly shortwave infrared (SWIR, 1.0-2.5 µm), and thermal infrared (TIR)-mainly longwave infrared (LWIR, 8.0-13.0µm) regions.For example, clay minerals, carbonates and gypsum can be detected in both the RIR and TIR regions, whereas feldspars and quartz can be detected only in the TIR, mainly LWIR region.Merging both RIR and TIR regions' data enables to identify a variety of minerals and mapping their distribution in the surface, as shown in several studies (e.g., [5,6]).
In a previous study [7], we developed a procedure based on comparing day and night LWIR region hyperspectral (henceforth LWIR) data acquired with the AisaOWL (Specim) sensor over Makhtesh Ramon in Israel, for mineral mapping.Dominant minerals-quartz, feldspars, clay minerals, gypsum and carbonates-were identified and their distribution was mapped.Pursuant to those results, we extended the study by adding the SWIR region hyperspectral (henceforth SWIR) data acquired with the AisaFENIX (Specim) sensor.Herein we examine the benefit of combining the LWIR and SWIR data for mineral mapping.

The Study Area and Airborne Data
Makhtesh Ramon, an eroded anticline, 40 km long and 7 km wide, with a total area of about 200 km 2 , is located in the Negev desert in Southern Israel, as described by Notesco et al. [7].The surface consists of a variety of minerals and rock types, including clay minerals, calcite, dolomite, gypsum, iron oxides, chert, marl, basalt, tuff, andesite and sandstone.It is an ideal site for examining new remote-sensing technology and ideas, being an arid area with extremely sparse vegetation and cloudless most of the year [8].
The airborne data, two flight lines covering a total area of 15 km 2 (Figure 1), were acquired during two flight campaigns.The LWIR data were acquired with Specim's AisaOWL sensor, covering the 7.7-12 µm spectral range, as described by Notesco et al. [7].The SWIR data were acquired with Specim's AisaFENIX sensor, covering the 0.4-2.5 µm spectral range with 448 spectral bands and a FWHM of 0.0035-0.0055µm [9], on 19 September 2014 at 08:50 UTC.The flight altitude was 1.34 km above ground level, resulting in a spatial resolution of 2 m, as in the LWIR data.

The Study Area and Airborne Data
Makhtesh Ramon, an eroded anticline, 40 km long and 7 km wide, with a total area of about 200 km 2 , is located in the Negev desert in Southern Israel, as described by Notesco et al. [7].The surface consists of a variety of minerals and rock types, including clay minerals, calcite, dolomite, gypsum, iron oxides, chert, marl, basalt, tuff, andesite and sandstone.It is an ideal site for examining new remote-sensing technology and ideas, being an arid area with extremely sparse vegetation and cloudless most of the year [8].
The airborne data, two flight lines covering a total area of 15 km 2 (Figure 1), were acquired during two flight campaigns.The LWIR data were acquired with Specim's AisaOWL sensor, covering the 7.7-12 μm spectral range, as described by Notesco et al. [7].The SWIR data were acquired with Specim's AisaFENIX sensor, covering the 0.4-2.5 μm spectral range with 448 spectral bands and a FWHM of 0.0035-0.0055μm [9], on 19 September 2014 at 08:50 UTC.The flight altitude was 1.34 km above ground level, resulting in a spatial resolution of 2 m, as in the LWIR data.The study area (image source: Esri, DigitalGlobe, GeoEye, Earthstar Geographics, CNES/Airbus DS, USDA, USGS, AEX, Getmapping, Aerogrid, IGN, IGP, swisstopo, and the GIS User Community), the two flight lines (1 and 2) and the location of 18 regions of interest (ROIs).Each ROI, consisting of more than 10 pixels, represents hundreds of square meters of a uniform surface that was sampled and measured with an X-ray diffractometer (XRD) as described by Notesco et al. [7].The XRD analysis results are shown below in Table 1.

Data Analysis
Absorption features in the SWIR reflectance of clay minerals (at 1.4 and 2.2 μm), gypsum (at 1.4-1.5, 1.75, 1.9 and 2.2 μm), and carbonates (at 2.3 μm) are due to vibration modes in the hydroxyl functional groups, water molecules and carbonate ion, respectively.Spectral features in the LWIR emissivity of silicates (quartz, feldspars and clay minerals at 8.3-9.5 μm), gypsum (at 8.6 μm) and carbonates (at 11.3 μm) are due to vibration modes of the silicon-oxygen bond (Si-O), sulfur-oxygen bond (S-O) and carbon-oxygen bond (C-O), respectively [10][11][12][13].Each ROI, consisting of more than 10 pixels, represents hundreds of square meters of a uniform surface that was sampled and measured with an X-ray diffractometer (XRD) as described by Notesco et al. [7].The XRD analysis results are shown below in Table 1.

SWIR Data Analysis
Each radiance image acquired with the AisaFENIX sensor was atmospherically corrected using ACORN-6 software (ImSpec LLC, Advanced Imaging and Spectroscopy), resulting in a reflectance image, followed by continuum removal normalization in specific spectral ranges to emphasize indicative absorption features [14,15], between 2.12-2.24µm and between 2.26-2.37 µm, indicative of clay minerals and carbonates, respectively, as demonstrated in Figure 2. A pixel in the normalized reflectance image between 2.12 and 2.24 µm was classified as clay if NR 2.20µm ď 0.90, where NR is the normalized reflectance value at the indicated wavelength.Furthermore, kaolinite could be identified with a ratio (NR 2.16µm ´NR 2.20µm )/(NR 2.23µm ´NR 2.20µm ) < 1, whereas for other clay minerals this ratio is ě1 [16,17], as demonstrated in Figure 2a for kaolinite, illite and montmorillonite with ratios of 0.28, 1.04 and 1.22, respectively.A pixel in the normalized reflectance image between 2.26 and 2.37 µm was classified as dolomite or calcite based on the absorption at 2.32 or 2.34 µm, respectively, as demonstrated in Figure 2b.

SWIR Data Analysis
Each radiance image acquired with the AisaFENIX sensor was atmospherically corrected using ACORN-6 software (ImSpec LLC, Advanced Imaging and Spectroscopy), resulting in a reflectance image, followed by continuum removal normalization in specific spectral ranges to emphasize indicative absorption features [14,15], between 2.12-2.24μm and between 2.26-2.37 μm, indicative of clay minerals and carbonates, respectively, as demonstrated in Figure 2. A pixel in the normalized reflectance image between 2.12 and 2.24 μm was classified as clay if NR2.20μm ≤ 0.90, where NR is the normalized reflectance value at the indicated wavelength.Furthermore, kaolinite could be identified with a ratio (NR2.16μm− NR2.20μm)/(NR2.23μm− NR2.20μm) < 1, whereas for other clay minerals this ratio is ≥1 [16,17], as demonstrated in Figure 2a for kaolinite, illite and montmorillonite with ratios of 0.28, 1.04 and 1.22, respectively.A pixel in the normalized reflectance image between 2.26 and 2.37 μm was classified as dolomite or calcite based on the absorption at 2.32 or 2.34 μm, respectively, as demonstrated in Figure 2b.(a) Spectra of regions of interest (ROIs) A and B from the normalized reflectance image and of clay minerals from a spectral library (Slib.)[12] for comparison, with the indicative absorption feature at 2.20 μm; (b) Spectra of ROIs C and D from the normalized reflectance image and of carbonates from the spectral library, with indicative absorption features at 2.32 and 2.34 μm of dolomite and calcite, respectively.

LWIR Data Analysis
The LWIR data analysis procedure was described by Notesco et al. [7] as "The Gain Spectrum Approach".Briefly, the at-sensor radiance of each pixel in the image was divided by a tangent blackbody radiation curve, resulting in an Ls (at-sensor radiance)/Lb (fitted blackbody radiation curve) image.The spectrum of each pixel in the Ls/Lb image was divided by a gain-factor spectrum, derived from the Ls/Lb spectrum of an in situ surface, reducing the atmospheric contribution and emphasizing mineral-related features, and resulting in an approximate emissivity image.Dominant minerals-quartz, feldspars, clay minerals, gypsum and carbonates-were identified by using the continuum removal-based indices given by Notesco et al. [7].

Results and Discussion
Mineralogical maps of the surface of Makhtesh Ramon, based on the analyses of the SWIR and LWIR images, are shown in Figure 3. Clay minerals and carbonates-calcite and dolomite-were identified based on the SWIR region data (Figure 3a).More than 95% of the pixels identified as clay minerals resulted in a ratio (NR2.16μm− NR2.20μm)/(NR2.23μm− NR2.20μm) < 0.75, e.g., 0.35 and 0.59 in ROIs A and B (Figure 2a), respectively, and were therefore identified as kaolinite.Identifying gypsum based on the SWIR image was problematic because of its spectral features' proximity to atmospheric

LWIR Data Analysis
The LWIR data analysis procedure was described by Notesco et al. [7] as "The Gain Spectrum Approach".Briefly, the at-sensor radiance of each pixel in the image was divided by a tangent blackbody radiation curve, resulting in an Ls (at-sensor radiance)/Lb (fitted blackbody radiation curve) image.The spectrum of each pixel in the Ls/Lb image was divided by a gain-factor spectrum, derived from the Ls/Lb spectrum of an in situ surface, reducing the atmospheric contribution and emphasizing mineral-related features, and resulting in an approximate emissivity image.Dominant minerals-quartz, feldspars, clay minerals, gypsum and carbonates-were identified by using the continuum removal-based indices given by Notesco et al. [7].

Results and Discussion
Mineralogical maps of the surface of Makhtesh Ramon, based on the analyses of the SWIR and LWIR images, are shown in Figure 3. Clay minerals and carbonates-calcite and dolomite-were identified based on the SWIR region data (Figure 3a).More than 95% of the pixels identified as clay minerals resulted in a ratio (NR 2.16µm ´NR 2.20µm )/(NR 2.23µm ´NR 2.20µm ) < 0.75, e.g., 0.35 and 0.59 in ROIs A and B (Figure 2a), respectively, and were therefore identified as kaolinite.Identifying gypsum based on the SWIR image was problematic because of its spectral features' proximity to atmospheric water vapor (1.40, 1.90 µm) and clay minerals (2.20 µm) features.Quartz, feldspars, clay minerals, gypsum and carbonates were identified based on the LWIR region data (Figure 3b).Classification of the dominant minerals resembled the XRD analysis results of selected rock samples collected from the study area (Table 1).Clay minerals and carbonates were identified in both SWIR and LWIR data (Figure 3a,b).90% of the pixels identified with carbonates were identified in both SWIR (as calcite or dolomite) and LWIR (as carbonates) data sets, strengthening the attribution of the absorption features at 2.32 and 2.34 μm and 11.2 μm (Notesco et al. [7]) to carbonates; however, only 10% of the pixels identified with clay minerals were identified in both SWIR (as kaolinite) and LWIR (as clay minerals) data sets (Figure 4).This is explained by the proximity of the spectral features of kaolinite and quartz in the LWIR emissivity (Figure 5).ROI A, with a relative large amount of kaolinite (according to XRD Classification of the dominant minerals resembled the XRD analysis results of selected rock samples collected from the study area (Table 1).Clay minerals and carbonates were identified in both SWIR and LWIR data (Figure 3a,b).90% of the pixels identified with carbonates were identified in both SWIR (as calcite or dolomite) and LWIR (as carbonates) data sets, strengthening the attribution of the absorption features at 2.32 and 2.34 µm and 11.2 µm (Notesco et al. [7]) to carbonates; however, only 10% of the pixels identified with clay minerals were identified in both SWIR (as kaolinite) and LWIR (as clay minerals) data sets (Figure 4).This is explained by the proximity of the spectral features of kaolinite and quartz in the LWIR emissivity (Figure 5).ROI A, with a relative large amount of kaolinite (according to XRD analysis), was identified in both SWIR and LWIR data whereas ROI B, with a relative large amount of quartz and less kaolinite (according to XRD analysis, Table 1), was identified as kaolinite in the SWIR data and as quartz in the LWIR data.About 10% of the pixels identified with kaolinite in the SWIR data were identified with gypsum in the LWIR data (Figure 4).ROI E, with a relative large amount of gypsum and no kaolinite (according to XRD analysis, Table 1), was identified as kaolinite in the SWIR data and as gypsum in the LWIR data (Figure 5).The proximity of the absorption at 2.20 µm of gypsum and kaolinite may lead to "false kaolinite" identification in the SWIR data.However, gypsum can be identified in the LWIR data with its noticeable feature at 8.63 µm.Thus, when classifying kaolinite pixels in the SWIR image (Figure 3a), we omitted the pixels identified as gypsum in the LWIR image (Figure 3b).Note that the presence of small amounts of kaolinite in such pixels cannot be completely negated.
Remote Sens. 2016, 8, 318 5 of 8 analysis), was identified in both SWIR and LWIR data whereas ROI B, with a relative large amount of quartz and less kaolinite (according to XRD analysis, Table 1), was identified as kaolinite in the SWIR data and as quartz in the LWIR data.About 10% of the pixels identified with kaolinite in the SWIR data were identified with gypsum in the LWIR data (Figure 4).ROI E, with a relative large amount of gypsum and no kaolinite (according to XRD analysis, Table 1), was identified as kaolinite in the SWIR data and as gypsum in the LWIR data (Figure 5).The proximity of the absorption at 2.20 μm of gypsum and kaolinite may lead to "false kaolinite" identification in the SWIR data.However, gypsum can be identified in the LWIR data with its noticeable feature at 8.63 μm.Thus, when classifying kaolinite pixels in the SWIR image (Figure 3a), we omitted the pixels identified as gypsum in the LWIR image (Figure 3b).Note that the presence of small amounts of kaolinite in such pixels cannot be completely negated.The overall mineralogy in each pixel (combining both SWIR and LWIR classification) was determined by giving an arbitrary number to each binary image, representing a specific mineral, and then summing all of these images.
The combined mineral map is shown in Figure 6 and is compared to the geological map of the area.Across flight line 1, the mineralogy changes from gypsum-, carbonates-and kaolinite-dominant rocks (Mohila formation "tr3") in the southwest part, interrupted by intrusion of feldspars-dominant  Remote Sens. 2016, 8, 318 5 of 8 analysis), was identified in both SWIR and LWIR data whereas ROI B, with a relative large amount of quartz and less kaolinite (according to XRD analysis, Table 1), was identified as kaolinite in the SWIR data and as quartz in the LWIR data.About 10% of the pixels identified with kaolinite in the SWIR data were identified with gypsum in the LWIR data (Figure 4).ROI E, with a relative large amount of gypsum and no kaolinite (according to XRD analysis, Table 1), was identified as kaolinite in the SWIR data and as gypsum in the LWIR data (Figure 5).The proximity of the absorption at 2.20 μm of gypsum and kaolinite may lead to "false kaolinite" identification in the SWIR data.However, gypsum can be identified in the LWIR data with its noticeable feature at 8.63 μm.Thus, when classifying kaolinite pixels in the SWIR image (Figure 3a), we omitted the pixels identified as gypsum in the LWIR image (Figure 3b).Note that the presence of small amounts of kaolinite in such pixels cannot be completely negated.The overall mineralogy in each pixel (combining both SWIR and LWIR classification) was determined by giving an arbitrary number to each binary image, representing a specific mineral, and then summing all of these images.
The combined mineral map is shown in Figure 6 and is compared to the geological map of the area.Across flight line 1, the mineralogy changes from gypsum-, carbonates-and kaolinite-dominant rocks (Mohila formation "tr3") in the southwest part, interrupted by intrusion of feldspars-dominant  The overall mineralogy in each pixel (combining both SWIR and LWIR classification) was determined by giving an arbitrary number to each binary image, representing a specific mineral, and then summing all of these images.
The combined mineral map is shown in Figure 6 and is compared to the geological map of the area.Across flight line 1, the mineralogy changes from gypsum-, carbonates-and kaolinite-dominant rocks (Mohila formation "tr3") in the southwest part, interrupted by intrusion of feldspars-dominant basalt rocks ("im"), to carbonates-, quartz-and kaolinite-dominant rocks (Mishor and Ardon formations "jl1"), to quartz-and kaolinite-dominant rocks (Inmar formation "jl2") and quartz-, kaolinite-and carbonates-dominant rocks in the northeastern part (Kurnub group "lck").The addition of carbonates in the northeastern part is from weathering processes occurring in the carbonates-dominant Har Ardon (700 m above sea level) rocks (Hevyon formation "c1", outside the flight line) to the lower (500 m above sea level) quartz-and kaolinite-dominant rocks.Across flight line 2, with a gradual elevation from 420 to 740 m above sea level (southeast to northwest), the mineralogy changes from carbonates-, gypsum-and kaolinite-dominant Mohila formation ("tr3") rocks to kaolinite-, quartz-and carbonates-dominant Alluvium ("q") and Inmar formation ("jl2") rocks.The quartz-dominant "S shape" in the southeast part is related to sedimentation processes along Nahal Ramon (a dry wadi) stream.The carbonates supplement to the quartz and kaolinite in Inmar formation rocks demonstrates the contribution of weathering processes in the elevated carbonates-dominant Hevyon formation ("c1") and En Yorqe'am, Zafit and Avnon formations ("c2") rocks.The northwest part is interrupted by an intrusion of feldspars-dominant basalt rocks ("βlc").The classification of dominant minerals based on the SWIR and LWIR data is compatible with the rock types, demonstrating the ability of combined hyperspectral technology to identify a variety of minerals and mapping their distribution in the surface.
Remote Sens. 2016, 8, 318 6 of 8 basalt rocks ("im"), to carbonates-, quartz-and kaolinite-dominant rocks (Mishor and Ardon formations "jl1"), to quartz-and kaolinite-dominant rocks (Inmar formation "jl2") and quartz-, kaolinite-and carbonates-dominant rocks in the northeastern part (Kurnub group "lck").The addition of carbonates in the northeastern part is from weathering processes occurring in the carbonates-dominant Har Ardon (700 m above sea level) rocks (Hevyon formation "c1", outside the flight line) to the lower (500 m above sea level) quartz-and kaolinite-dominant rocks.Across flight line 2, with a gradual elevation from 420 to 740 m above sea level (southeast to northwest), the mineralogy changes from carbonates-, gypsum-and kaolinite-dominant Mohila formation ("tr3") rocks to kaolinite-, quartz-and carbonates-dominant Alluvium ("q") and Inmar formation ("jl2") rocks.The quartz-dominant "S shape" in the southeast part is related to sedimentation processes along Nahal Ramon (a dry wadi) stream.The carbonates supplement to the quartz and kaolinite in Inmar formation rocks demonstrates the contribution of weathering processes in the elevated carbonates-dominant Hevyon formation ("c1") and En Yorqe'am, Zafit and Avnon formations ("c2") rocks.The northwest part is interrupted by an intrusion of feldspars-dominant basalt rocks ("βlc").
The classification of dominant minerals based on the SWIR and LWIR data is compatible with the rock types, demonstrating the ability of combined hyperspectral technology to identify a variety of minerals and mapping their distribution in the surface.

Summary and Conclusions
This study presents the ability to identify minerals and map their distribution in the surface using hyperspectral remote-sensing data in the SWIR and LWIR regions.Two airborne data sets, acquired with AisaFENIX and AisaOWL sensors over Makhtesh Ramon in Israel, were analyzed.The reflectance spectrum of each pixel in the SWIR images was calculated by applying an atmospheric correction procedure to the radiance images and the 2.0-2.4 µm spectral range was analyzed, enabling identification of kaolinite, calcite and dolomite.The approximate emissivity spectrum of each pixel in the LWIR images was calculated by applying a procedure described by Notesco et al. [7], and the 8.0-12.0µm spectral range was analyzed, enabling identification of quartz and feldspars (minerals that are featureless in the SWIR region), clay minerals, gypsum and carbonates.Combining the results demonstrated the benefit of using both SWIR and LWIR data sets to identify a variety of minerals and improve their classification.For example, pixels with a relative large amount of quartz and less kaolinite (according to XRD analysis) were identified as quartz only in the LWIR data and as kaolinite only in the SWIR data.The ability to identify gypsum in the LWIR data enabled eliminating "kaolinite false-positive" identification in the SWIR data.The resultant mineralogy was in quite good agreement with the mineralogy of selected rock samples, which were collected from the study area as derived from XRD analysis.The resulting combined mineral map, classifying the dominant minerals in each pixel in the images, was in agreement with the geological map of the area.The results of this study emphasize the importance of integrating both RIR and TIR hyperspectral technology for better mineral mapping of the earth's surface.

Figure 1 .
Figure 1.The study area (image source: Esri, DigitalGlobe, GeoEye, Earthstar Geographics, CNES/Airbus DS, USDA, USGS, AEX, Getmapping, Aerogrid, IGN, IGP, swisstopo, and the GIS User Community), the two flight lines (1 and 2) and the location of 18 regions of interest (ROIs).Each ROI, consisting of more than 10 pixels, represents hundreds of square meters of a uniform surface that was sampled and measured with an X-ray diffractometer (XRD) as described by Notesco et al.[7].The XRD analysis results are shown below in Table1.

Figure 1 .
Figure 1.The study area (image source: Esri, DigitalGlobe, GeoEye, Earthstar Geographics, CNES/Airbus DS, USDA, USGS, AEX, Getmapping, Aerogrid, IGN, IGP, swisstopo, and the GIS User Community), the two flight lines (1 and 2) and the location of 18 regions of interest (ROIs).Each ROI, consisting of more than 10 pixels, represents hundreds of square meters of a uniform surface that was sampled and measured with an X-ray diffractometer (XRD) as described by Notesco et al.[7].The XRD analysis results are shown below in Table1.

Figure 2 .
Figure 2.(a) Spectra of regions of interest (ROIs) A and B from the normalized reflectance image and of clay minerals from a spectral library (Slib.)[12]for comparison, with the indicative absorption feature at 2.20 μm; (b) Spectra of ROIs C and D from the normalized reflectance image and of carbonates from the spectral library, with indicative absorption features at 2.32 and 2.34 μm of dolomite and calcite, respectively.

Figure 4 .
Figure 4. Identification of minerals in both SWIR and LWIR data sets or in only one-SWIR or LWIR-data set.

Figure 5 .
Figure 5. SWIR reflectance spectra (a) and LWIR approximate emissivity spectra (b) of ROIs A, B and E. The relevant features are emphasized with thick lines for the identified minerals.Atmospheric water vapor absorptions were omitted from SWIR reflectance spectra.

Figure 4 .
Figure 4. Identification of minerals in both SWIR and LWIR data sets or in only one-SWIR or LWIR-data set.

Figure 4 .
Figure 4. Identification of minerals in both SWIR and LWIR data sets or in only one-SWIR or LWIR-data set.

Figure 5 .
Figure 5. SWIR reflectance spectra (a) and LWIR approximate emissivity spectra (b) of ROIs A, B and E. The relevant features are emphasized with thick lines for the identified minerals.Atmospheric water vapor absorptions were omitted from SWIR reflectance spectra.

Figure 5 .
Figure 5. SWIR reflectance spectra (a) and LWIR approximate emissivity spectra (b) of ROIs A, B and E. The relevant features are emphasized with thick lines for the identified minerals.Atmospheric water vapor absorptions were omitted from SWIR reflectance spectra.