[go: up one dir, main page]
More Web Proxy on the site http://driver.im/
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Sign in to use this feature.

Years

Between: -

Subjects

remove_circle_outline
remove_circle_outline
remove_circle_outline

Journals

Article Types

Countries / Regions

Search Results (10)

Search Parameters:
Keywords = AISA Fenix

Order results
Result details
Results per page
Select all
Export citation of selected articles as:
26 pages, 20840 KiB  
Article
Vicarious CAL/VAL Approach for Orbital Hyperspectral Sensors Using Multiple Sites
by Daniela Heller Pearlshtien, Stefano Pignatti and Eyal Ben-Dor
Remote Sens. 2023, 15(3), 771; https://doi.org/10.3390/rs15030771 - 29 Jan 2023
Cited by 5 | Viewed by 2749
Abstract
The hyperspectral (HSR) sensors Earth Surface Mineral Dust Source Investigation (EMIT) of the National Aeronautics and Space Administration (NASA) and Environmental Mapping and Analysis Program (EnMAP) of the German Aerospace Center (DLR) were recently launched. These state-of-the-art sensors have joined the already operational [...] Read more.
The hyperspectral (HSR) sensors Earth Surface Mineral Dust Source Investigation (EMIT) of the National Aeronautics and Space Administration (NASA) and Environmental Mapping and Analysis Program (EnMAP) of the German Aerospace Center (DLR) were recently launched. These state-of-the-art sensors have joined the already operational HSR sensors DESIS (DLR), PRISMA (Italian Space Agency), and HISUI (developed by the Japanese Ministry of Economy, Trade, and Industry METI and Japan Aerospace Exploration Agency JAXA). The launching of more HSR sensors is being planned for the near future (e.g., SBG of NASA, and CHIME of the European Space Agency), and the challenge of monitoring and maintaining their calibration accuracy is becoming more relevant. We proposed two test sites: Amiaz Plain (AP) and Makhtesh Ramon (MR) for spectral, radiometric, and geometric calibration/validation (CAL/VAL). The sites are situated in the arid environment of southern Israel and are in the same overpass coverage. Both test sites have already demonstrated favorable results in assessing an HSR sensor’s performance and were chosen to participate in the EMIT and EnMAP validation stage. We first evaluated the feasibility of using AP and MR as CAL/VAL test sites with extensive datasets and sensors, such as the multispectral sensor Landsat (Landsat5 TM and Landsat8 OLI), the airborne HSR sensor AisaFENIX 1K, and the spaceborne HSR sensors DESIS and PRISMA. Field measurements were taken over time. The suggested methodology integrates reflectance and radiometric CAL/VAL test sites into one operational protocol. The method can highlight degradation in the spectral domain early on, help maintain quantitative applications, adjust the sensor’s radiometric calibration during its mission lifetime, and minimize uncertainties of calibration parameters. A PRISMA sensor case study demonstrates the complete operational protocol, i.e., performance evaluation, quality assessment, and cross-calibration between HSR sensors. These CAL/VAL sites are ready to serve as operational sites for other HSR sensors. Full article
(This article belongs to the Special Issue Accuracy and Quality Control of Remote Sensing Data)
Show Figures

Figure 1

Figure 1
<p>Amiaz Plain (AP) and Makhtesh Ramon (MR) calibration/validation (CAL/VAL) test sites. (<b>a</b>) AP and MR locations in Israel. (<b>b</b>) Work protocol flowchart.</p>
Full article ">Figure 2
<p>MR test sites. Different mineralogies and spectral signatures characterize the 6 test sites. Spectral signatures are from the spectral library acquired using an Analytical Spectral Devices (ASD) spectroradiometer.</p>
Full article ">Figure 3
<p>Location of test sites in (<b>a</b>) MR, and (<b>b</b>) AP. (<b>c</b>) Field spectral measurement protocol; example from the AP test site. (<b>d</b>) Average spectral signal (AVG, solid blue line) and standard deviation (SD, dashed black line) for AP 30 × 30 cube.</p>
Full article ">Figure 4
<p>Spatial change detection result for AP, 2000 vs. 2021. Test site is marked in a red box; detected changed pixels are colored in green (center image).</p>
Full article ">Figure 5
<p>Examination of spectral change in AP. (<b>a</b>) Mean Landsat reflectance signal for each year (2000–2021). (<b>b</b>) SAM, ASDS, and RMSE for each year against the year 2000 (t0).</p>
Full article ">Figure 6
<p>Spatial change detection result for MR, 1996 vs. 2021. Spectral test sites are marked by yellow Xs, detected changed pixels are colored in red.</p>
Full article ">Figure 7
<p>Examination of spectral changes at MR VAL test sites: LANDSAT reflectance signal for each year (1996–2021): (<b>a</b>) brown questa; (<b>b</b>) laccolite; (<b>c</b>) gypsum mine; (<b>d</b>) gypsum fans; (<b>e</b>) kaolinite; (<b>f</b>) calcite. The mean value is marked by black dashed lines.</p>
Full article ">Figure 8
<p>Spectral change examination. SAM, ASDS, and RMSE for each year against the year 1996 (t0).</p>
Full article ">Figure 9
<p>Radiance comparisons. TOA radiance from PRISMA L1 vs. MODTRAN simulation.</p>
Full article ">Figure 10
<p>Evaluation of PRISMA L2D spectra (PRI) vs. ASD at the six test sites (PRISMA—solid line, ASD ground truth—dashed line). BQ, brown questa; L, laccolite; GM, gypsum—mine; GF, gypsum—fans.</p>
Full article ">Figure 11
<p>(<b>a</b>) AP cross-calibration results for PRISMA (PRI); original (before calibration; black line), after SBAF cross-calibration (Cal, blue line) against ASD ground truth (red dotted line), ratio of original (ORI) to calibrated (green triangle). (<b>b</b>) Validation of AP SBAF on MR VNIR test site using brown questa (BQ).</p>
Full article ">Figure 12
<p>The spectral band adjustment factors (SBAFs) for PRISMA in different acquisition months (April—blue; June—orange; November—red; December—green; average SBAF—black).</p>
Full article ">Figure 13
<p>The AP cross-calibration result for PRISMA (PRI) with well-calibrated AisaFENIX; original reflectance (before calibration, black line), after SBAF calibration (blue line) against ASD ground truth (red dotted line).</p>
Full article ">Figure 14
<p>The MR test sites of PRISMA original (before calibration) reflectance (black line) against corrected reflectance (after applying the SBAF, blue line), and the field ASD spectra (dashed black line); (<b>a</b>) brown questa (BQ), (<b>b</b>) laccolite, (<b>c</b>) gypsum mine, (<b>d</b>) gypsum fans, (<b>e</b>) kaolinite, (<b>f</b>) calcite.</p>
Full article ">Figure 15
<p>Comparison of PRISMA’s 2019 and 2022 geometric corrections: (<b>a</b>) 2019 MR image; (<b>b</b>) 2022 MR image; (<b>c</b>) 2022 image overlaid on 2019 image. Several drifts in the geolocation are marked with red arrows (2019—solid; 2022—dashed).</p>
Full article ">Figure 16
<p>Comparison of the same coordinates (red cross) for the years 2019–2022 at three MR test sites: kaolinite (top panels), gypsum soil fans, and laccolite (bottom panels). There is a slight shift in the locations.</p>
Full article ">
15 pages, 4303 KiB  
Article
CO2Flux Model Assessment and Comparison between an Airborne Hyperspectral Sensor and Orbital Multispectral Imagery in Southern Amazonia
by João Lucas Della-Silva, Carlos Antonio da Silva Junior, Mendelson Lima, Paulo Eduardo Teodoro, Marcos Rafael Nanni, Luciano Shozo Shiratsuchi, Larissa Pereira Ribeiro Teodoro, Guilherme Fernando Capristo-Silva, Fabio Henrique Rojo Baio, Gabriel de Oliveira, José Francisco de Oliveira-Júnior and Fernando Saragosa Rossi
Sustainability 2022, 14(9), 5458; https://doi.org/10.3390/su14095458 - 1 May 2022
Cited by 7 | Viewed by 3102
Abstract
In environmental research, remote sensing techniques are mostly based on orbital data, which are characterized by limited acquisition and often poor spectral and spatial resolutions in relation to suborbital sensors. This reflects on carbon patterns, where orbital remote sensing bears devoted sensor systems [...] Read more.
In environmental research, remote sensing techniques are mostly based on orbital data, which are characterized by limited acquisition and often poor spectral and spatial resolutions in relation to suborbital sensors. This reflects on carbon patterns, where orbital remote sensing bears devoted sensor systems for CO2 monitoring, even though carbon observations are performed with natural resources systems, such as Landsat, supported by spectral models such as CO2Flux adapted to multispectral imagery. Based on the considerations above, we have compared the CO2Flux model by using four different imagery systems (Landsat 8, PlanetScope, Sentinel-2, and AisaFenix) in the northern part of the state of Mato Grosso, southern Brazilian Amazonia. The study area covers three different land uses, which are primary tropical forest, bare soil, and pasture. After the atmospheric correction and radiometric calibration, the scenes were resampled to 30 m of spatial resolution, seeking for a parametrized comparison of CO2Flux, as well as NDVI (Normalized Difference Vegetation Index) and PRI (Photochemical Reflectance Index). The results obtained here suggest that PlanetScope, MSI/Sentinel-2, OLI/Landsat-8, and AisaFENIX can be similarly scaled, that is, the data variability along a heterogeneous scene in evergreen tropical forest is similar. We highlight that the spatial-temporal dynamics of rainfall seasonality relation to CO2 emission and uptake should be assessed in future research. Our results provide a better understanding on how the merge and/or combination of different airborne and orbital datasets that can provide reliable estimates of carbon emission and absorption within different terrestrial ecosystems in southern Amazonia. Full article
(This article belongs to the Special Issue Dynamics of Heat Spots and Sustainable Agriculture)
Show Figures

Figure 1

Figure 1
<p>Description of the imaged area in Alta Floresta, Mato Grosso, Brazil.</p>
Full article ">Figure 2
<p>(<b>a</b>) Regions of interest, classified as forest (green), bare soil (red), and pasture (yellow) areas; (<b>b</b>) spectral profile transect over the OLI/Landsat-8 scene.</p>
Full article ">Figure 3
<p>Equipment used in the airborne imaging of the study area.</p>
Full article ">Figure 4
<p>(<b>a</b>) AisaFENIX scene to resampled (<b>b</b>) scene (RGB composition on bands 42, 26, and 12). Scene detail from original (<b>c</b>) and resampled (<b>d</b>) image.</p>
Full article ">Figure 5
<p>Normalized Difference Vegetation Index (NDVI) results for OLI/Landsat-8 (<b>a</b>), MSI/Sentinel-2 (<b>b</b>), PlanetScope (<b>c</b>), and AisaFENIX (<b>d</b>) imagery.</p>
Full article ">Figure 6
<p>Photochemical Reflectance Index (PRI) results for OLI/Landsat-8 (<b>a</b>), MSI/Sentinel-2 (<b>b</b>), PlanetScope (<b>c</b>), and AisaFENIX (<b>d</b>) imagery.</p>
Full article ">Figure 7
<p>CO<sub>2</sub>Flux emission (μmol m<sup>−2</sup> s<sup>−1</sup>) results for OLI/Landsat-8 (<b>a</b>), MSI/Sentinel-2 (<b>b</b>), PlanetScope (<b>c</b>), and AisaFENIX (<b>d</b>) imagery.</p>
Full article ">Figure 8
<p>NDVI, PRI, and CO<sub>2</sub>Flux spectral profiles based on the transect described in <a href="#sustainability-14-05458-f002" class="html-fig">Figure 2</a>b t. Horizontal axis refers to the spatial variation over the transect line (kilometers), while vertical axis represents the index value for NDVI, PRI, and μmol m<sup>−2</sup> s<sup>−</sup><sup>1</sup> for CO<sub>2</sub>Flux.</p>
Full article ">Figure 9
<p>Significant interaction among sensors versus LULC for the vegetation indices NDVI, PRI, sPRI, and CO<sub>2</sub>Flux. Uppercase letters express statistical similarity or difference on by comparing LULC for the same sensor, and lowercase letters express statistical similarity or difference on the same LULC for different sensors.</p>
Full article ">Figure 10
<p>Linear regressions and trend lines based on AisaFENIX versus (<b>a</b>) OLI/Landsat-8, (<b>b</b>) MSI/Sentinel-2, and (<b>c</b>) PlanetScope results of CO<sub>2</sub>Flux.</p>
Full article ">Figure 11
<p>NDVI on bare soil area (blue polygon) based on OLI/Landsat-8 (<b>a</b>), MSI/Sentinel-2 (<b>b</b>), PlanetScope (<b>c</b>), and AisaFENIX (<b>d</b>) images.</p>
Full article ">Figure 12
<p>CO<sub>2</sub>Flux on bare soil area (blue polygon) based on Landsat-8 (<b>a</b>), Sentinel-2 (<b>b</b>), PLANET (<b>c</b>), and AisaFENIX (<b>d</b>) images.</p>
Full article ">Figure 13
<p>CO<sub>2</sub>Flux on tropical forest areas (blue polygons) based on OLI/Landsat-8 (<b>a</b>), Sentinel-2 (<b>b</b>), PLANET (<b>c</b>), and AisaFENIX (<b>d</b>) images.</p>
Full article ">Figure 14
<p>CO<sub>2</sub>Flux on pasture area (blue polygon) based on OLI/Landsat-8 (<b>a</b>), Sentinel-2 (<b>b</b>), PLANET (<b>c</b>), and AisaFENIX (<b>d</b>) images.</p>
Full article ">
19 pages, 5697 KiB  
Article
Mapping Particle Size and Soil Organic Matter in Tropical Soil Based on Hyperspectral Imaging and Non-Imaging Sensors
by Marcos Rafael Nanni, José Alexandre Melo Demattê, Marlon Rodrigues, Glaucio Leboso Alemparte Abrantes dos Santos, Amanda Silveira Reis, Karym Mayara de Oliveira, Everson Cezar, Renato Herrig Furlanetto, Luís Guilherme Teixeira Crusiol and Liang Sun
Remote Sens. 2021, 13(9), 1782; https://doi.org/10.3390/rs13091782 - 3 May 2021
Cited by 21 | Viewed by 4146
Abstract
We evaluated the use of airborne hyperspectral imaging and non-imaging sensors in the Vis—NIR—SWIR spectral region to assess particle size and soil organic matter in the surface layer of tropical soils (Oxisols, Ultisols, Entisols). The study area is near Piracicaba municipality, São Paulo [...] Read more.
We evaluated the use of airborne hyperspectral imaging and non-imaging sensors in the Vis—NIR—SWIR spectral region to assess particle size and soil organic matter in the surface layer of tropical soils (Oxisols, Ultisols, Entisols). The study area is near Piracicaba municipality, São Paulo state, Brazil, in a sugarcane cultivation area of 135 hectares. The study area, with bare soil, was imaged in April 2016 by the AisaFENIX aerotransported hyperspectral sensor, with spectral resolution of 3.5 nm between 380 and 970 nm, and 12 nm between 970 and 2500 nm. We collected 66 surface soil samples. The samples were analyzed for particle size and soil organic matter content. Laboratory spectral measurements were performed using a non-imaging spectroradiometer (ASD FieldSpec 3 Jr). Partial Least Square Regression (PLSR) was used to predict clay, silt, sand and soil organic matter (SOM). The PLSR functions developed were applied to the hyperspectral image of the study area, allowing development of a prediction map of clay, sand, and SOM. The developed PLSR models demonstrated the relationship between the predictor variables at the cross-validation step, both for the non-imaging and imaging sensors, when the highest r and R2 values were obtained for clay, sand, and SOM, with R2 over 0.67. We did not obtain a satisfactory model for silt content. For the non-imaging sensor at the prediction step, R2 values for clay and SOM were over 0.7 and sand was lower than 0.54. The imaging sensor yielded models for clay, sand, and SOM with R2 values of 0.62, 0.66, and 0.67, respectively. Pearson correlation between sensors was greater than 0.849 for the prediction of clay, sand, and SOM. Our study successfully generated, from the imaging sensor, a large-scale and detailed predicted soil maps for particle size and SOM, which are important in the management of tropical soils. Full article
(This article belongs to the Special Issue Remote Sensing of Soil Properties)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Location of the study site and sampling points (+).</p>
Full article ">Figure 2
<p>Topography of the study area and location of the collected samples.</p>
Full article ">Figure 3
<p>Particle-size of the samples plotted in a USDA textural triangle.</p>
Full article ">Figure 4
<p>Spectral curves of soil samples. Curves of the 66 samples obtained by the AisaFENIX Sensor (<b>a</b>) and ASD FieldSpec 3 JR Sensor (<b>b</b>); average of spectral data from sandstone and siltstone obtained by the AisaFENIX Sensor (<b>c</b>) and ASD FieldSpec 3 JR Sensor (<b>d</b>). The gap in the curves of the AisaFENIX sensor is due to the absorption of EMR by the atmosphere.</p>
Full article ">Figure 5
<p>RMSE values of the PLSR model obtained by the imaging and non-imaging sensors in the cross-validation and prediction phases.</p>
Full article ">Figure 6
<p>Regression coefficients of PLSR models for imaging and non-imaging sensors for prediction of clay, sand, SOM and silt.</p>
Full article ">Figure 7
<p>Pearson’s correlation (with an alpha significance level of 5% and <span class="html-italic">p</span> value lower than (0.0001) between the contents predicted of clay (<b>a</b>), sand (<b>b</b>), and SOM (<b>c</b>) by PLSR models for imaging and non-imaging sensors.</p>
Full article ">Figure 8
<p>Estimated map of contents of (<b>top</b>) clay, (<b>middle</b>) sand, and (<b>bottom</b>) SOM of study area soils obtained using an airborne imaging sensor.</p>
Full article ">
23 pages, 7472 KiB  
Article
Mapping Canopy Chlorophyll Content in a Temperate Forest Using Airborne Hyperspectral Data
by J. Malin Hoeppner, Andrew K. Skidmore, Roshanak Darvishzadeh, Marco Heurich, Hsing-Chung Chang and Tawanda W. Gara
Remote Sens. 2020, 12(21), 3573; https://doi.org/10.3390/rs12213573 - 31 Oct 2020
Cited by 22 | Viewed by 4157
Abstract
Chlorophyll content, as the primary pigment driving photosynthesis, is directly affected by many natural and anthropogenic disturbances and stressors. Accurate and timely estimation of canopy chlorophyll content (CCC) is essential for effective ecosystem monitoring to allow for successful management interventions to occur. Hyperspectral [...] Read more.
Chlorophyll content, as the primary pigment driving photosynthesis, is directly affected by many natural and anthropogenic disturbances and stressors. Accurate and timely estimation of canopy chlorophyll content (CCC) is essential for effective ecosystem monitoring to allow for successful management interventions to occur. Hyperspectral remote sensing offers the possibility to accurately estimate and map canopy chlorophyll content. In the past, research has predominantly focused on the use of hyperspectral data on canopy chlorophyll content retrieval of crops and grassland ecosystems. Therefore, in this study, a temperate mixed forest, the Bavarian Forest National Park in Germany, was chosen as the study site. We compared different statistical models (narrowband vegetation indices (VIs), partial least squares regression (PLSR) and random forest (RF)) in their accuracy to predict CCC using airborne hyperspectral data. The airborne hyperspectral imagery was acquired by the AisaFenix sensor (623 bands; 3.5 nm spectral resolution in the visible near-infrared (VNIR) region, and 12 nm spectral resolution in the shortwave infrared (SWIR) region; 3 m spatial resolution) on July 6, 2017. In situ leaf chlorophyll content and leaf area index measurements were sampled from the upper canopy of coniferous, mixed, and deciduous forest stands in July and August 2017. The study yielded the highest retrieval accuracies with PLSR (root mean square error (RMSE) = 0.25 g/m2, R2 = 0.66). It further indicated specific spectral regions within the visible (390–400 nm and 470–540 nm), red edge (680–780 nm), near-infrared (1050–1100 nm) and shortwave infrared regions (2000–2270 nm) that were important for CCC retrieval. The results showed that forest CCC can be mapped with relatively high accuracies using image spectroscopy. Full article
(This article belongs to the Special Issue Remote Sensing for Estimating Leaf Chlorophyll Content in Plants)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Location of Bavarian Forest National Park (source: Park administration office) and hyperspectral transects.</p>
Full article ">Figure 2
<p>Methodological framework of the study. Red rectangles are representing input data, grey analysis steps and blue outputs.</p>
Full article ">Figure 3
<p>3-D correlation matrix of CCC and narrowband MERIS terrestrial chlorophyll index (MTCI) for all combinations of wavelengths that yielded R<sup>2</sup> values of 0.6 or higher (not cross-validated). The arrow in the top right corner points to the location that indicates the index with the optimal wavelength combination. The blue rectangle indicates the band combinations with wavelength 1 in the green region, wavelength 2 in the red edge region and wavelength 3 in one of the listed regions except shortwave infrared. The red circle shows the band combinations with all three wavelengths in the red edge.</p>
Full article ">Figure 4
<p>2-D correlation plot of R<sup>2</sup> of canopy chlorophyll content and narrowband Datt derivative index (not cross-validated). The vertical and horizontal black lines indicate the red edge region between 680–780 nm.</p>
Full article ">Figure 5
<p>Scatterplot of the in situ CCC and predicted values made by (<b>a</b>) nMTCI with optimal wavelength combinations, (<b>b</b>) nDD with optimal wavelength combinations, and (<b>c</b>) Partial least squares regression (PLSR) with five components. The dotted line shows the 1:1 relationship between predicted and measured values. (<b>d</b>) Boxplot of in situ and predicted CCC (g/m<sup>2</sup>) using the three selected models. Mean CCC was 1.64 g/m<sup>2</sup> for all models as well as in situ CCC. The median predicted CCC was 1.67 g/m<sup>2</sup>, 1.69 g/m<sup>2</sup> and 1.65 g/m<sup>2</sup> for nDD, nMTCI and PLSR, respectively. The median in situ CCC was 1.64 g/m<sup>2</sup>.</p>
Full article ">Figure 6
<p>Close up of eastern site of light line 193; (<b>a</b>) Canopy chlorophyll content map derived with PLSR, (<b>b</b>) Land cover classification map, (<b>c</b>) Deadwood classification map, (<b>d</b>) Location plan of flight line 193 and the close-up in relation to the Bavarian Forest National Park. Forest classification map from land cover map. Deadwood classification map from infestation map (obtained from aerial photographs) available from the national park administration [<a href="#B85-remotesensing-12-03573" class="html-bibr">85</a>].</p>
Full article ">Figure A1
<p>(<b>a</b>) Flight line 193: Forest classification map, (<b>b</b>) Flight line 193: Deadwood classification map, (<b>c</b>) Flight line 143: Forest classification map (<b>d</b>) Flight line 143: Deadwood classification map, (<b>e</b>) Location plan of flight line 193 and 143 in relation to the Bavarian Forest National Park. Forest classification map from land cover map available from the national park administration. Deadwood classification map from infestation map (obtained from aerial photographs) available from the national park administration.</p>
Full article ">Figure A2
<p>Flight line 193 mapped CCC (<b>a</b>) using nMTCI with wavelength 1 =1819.52 nm, wavelength 2 = 1996.52 nm and wavelength 3 = 390.07 nm, (<b>b</b>) using nDD with wavelength 1 = 528.83 nm and wavelength 3 = 735.02 nm, and (<b>c</b>) using PLSR.</p>
Full article ">Figure A3
<p>Flight line 149 mapped CCC (<b>a</b>) using nMTCI with wavelength 1=1819.52 nm, wavelength 2 = 1996.52 nm and wavelength 3 = 390.07 nm, (<b>b</b>) using nDD with wavelength 1 = 528.83 nm and wavelength 3 = 735.02 nm, and (<b>c</b>) using PLSR.</p>
Full article ">
33 pages, 12246 KiB  
Article
Monitoring of Canopy Stress Symptoms in New Zealand Kauri Trees Analysed with AISA Hyperspectral Data
by Jane J. Meiforth, Henning Buddenbaum, Joachim Hill and James Shepherd
Remote Sens. 2020, 12(6), 926; https://doi.org/10.3390/rs12060926 - 13 Mar 2020
Cited by 13 | Viewed by 5541
Abstract
The endemic New Zealand kauri trees (Agathis australis) are under threat by the deadly kauri dieback disease (Phytophthora agathidicida (PA)). This study aimed to identify spectral index combinations for characterising visible stress symptoms in the kauri canopy. The analysis is [...] Read more.
The endemic New Zealand kauri trees (Agathis australis) are under threat by the deadly kauri dieback disease (Phytophthora agathidicida (PA)). This study aimed to identify spectral index combinations for characterising visible stress symptoms in the kauri canopy. The analysis is based on an aerial AISA hyperspectral image mosaic and 1258 reference crowns in three study sites in the Waitakere Ranges west of Auckland. A field-based assessment scheme for canopy stress symptoms (classes 1–5) was further optimised for use with RGB aerial images. A combination of four indices with six bands in the spectral range 450–1205 nm resulted in a correlation of 0.93 (mean absolute error 0.27, RMSE 0.48) for all crown sizes. Comparable results were achieved with five indices in the 450–970 nm region. A Random Forest (RF) regression gave the most accurate predictions while a M5P regression tree performed nearly as well and a linear regression resulted in slightly lower correlations. Normalised Difference Vegetation Indices (NDVI) in the near-infrared / red spectral range were the most important index combinations, followed by indices with bands in the near-infrared spectral range from 800 to 1205 nm. A test on different crown sizes revealed that stress symptoms in smaller crowns with denser foliage are best described in combination with pigment-sensitive indices that include bands in the green and blue spectral range. A stratified approach with individual models for pre-segmented low and high forest stands improved the overall performance. The regression models were also tested in a pixel-based analysis. A manual interpretation of the resulting raster map with stress symptom patterns observed in aerial imagery indicated a good match. With bandwidths of 10 nm and a maximum number of six bands, the selected index combinations can be used for large-area monitoring on an airborne multispectral sensor. This study establishes the base for a cost-efficient, objective monitoring method for stress symptoms in kauri canopies, suitable to cover large forest areas with an airborne multispectral sensor. Full article
(This article belongs to the Special Issue Hyperspectral Remote Sensing for Biodiversity Mapping)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Kauri growth classes used in this study, according to the mean crown diameter (cdm) [<a href="#B11-remotesensing-12-00926" class="html-bibr">11</a>]. (Photos [<a href="#B12-remotesensing-12-00926" class="html-bibr">12</a>]).</p>
Full article ">Figure 2
<p>Mature kauri stand in different foliage colour variations in the Waitakere Ranges shown in (<b>a</b>) oblique view and (<b>b</b>) nadir view [<a href="#B13-remotesensing-12-00926" class="html-bibr">13</a>].</p>
Full article ">Figure 3
<p>Study sites in the Waitakere Ranges with the reference crowns marked in orange. The labels give the name of the area and the size in square kilometres. Small map: Location of the Waitakere Ranges on the North Island of New Zealand, west of Auckland City (background maps: [<a href="#B83-remotesensing-12-00926" class="html-bibr">83</a>,<a href="#B84-remotesensing-12-00926" class="html-bibr">84</a>]).</p>
Full article ">Figure 4
<p>Workflow for the preparation of crown-based attributes that were used in the analysis.</p>
Full article ">Figure 5
<p>(<b>a</b>) Crown size classes of the reference crowns (total 1258), used in the analysis per low and high forest stand situation. The crown size classes correspond to the classes used in <a href="#remotesensing-12-00926-t002" class="html-table">Table 2</a>. (<b>b</b>) Low and high forest stands were distinguished with an average height of 21 m on pre-segmented stand polygons based on a LiDAR CHM ([<a href="#B89-remotesensing-12-00926" class="html-bibr">89</a>]).</p>
Full article ">Figure 6
<p>Inter-crown (<b>a</b>,<b>b</b>) and within-crown (<b>c</b>,<b>d</b>) spectral variability for the sunlit part of 189 non-symptomatic small kauri crowns (crown diameter &gt; 3–4.8 m) (<b>a</b>,<b>c</b>) and 337 large kauri crowns with no visible stress symptoms (crown diameter &gt;12.8 m) (<b>b</b>,<b>d</b>). The mean spectra are marked in colour.</p>
Full article ">Figure 7
<p>Mean spectra (bold signatures) and standard deviation (stdev, thin signatures) of kauri in three symptom classes: Non-symptomatic (class 1, green), medium symptoms (class 3, orange) and dead trees (class 5, red). The number of pixels (pix) for the different classes is given in parentheses.</p>
Full article ">Figure 8
<p>Spectra of kauri in three different size classes according to their mean crown diameter (cdm) (S = small 3–4.8 m cdm, M = medium &gt;4.8–12.2 m cdm, L = large &gt;12.2 m cdm) sorted according to three symptom levels: “Non-symptomatic” (green), “medium stress symptoms” (orange) and “dead crowns” (red), which include visible undergrowth and epiphytes. For better readability, the spectra of crowns with medium symptoms and dead trees were offset by 3000 and 6000 units, respectively.</p>
Full article ">Figure 9
<p>Combined box- and jitter-plot diagrams of the predicted values versus the actual reference crown values for crowns in all sizes classes. Figure (<b>a</b>) shows the results of the original scale from 1 = “non-symptomatic” to 5 = “dead”. Figure (<b>b</b>) shows the results for a rescaled range, with the former value 5 for dead trees changed to value 8. The analysis is based on a RF regression with 1000 iterations in a 10-fold cross-validation on the baseline 6-band index combination for the full spectral range.</p>
Full article ">Figure 10
<p>Resulting maps (<b>a</b>,<b>c</b>,<b>e</b>) and corresponding RGB aerial images (<b>b</b>,<b>d</b>,<b>f</b>) (2016) [<a href="#B89-remotesensing-12-00926" class="html-bibr">89</a>] of a pixel-based application of the baseline index combination for two forest stands with marked reference crowns and their reference symptom class values. The analysis was carried out as a Random Forest regression in the EnMAP toolbox [<a href="#B96-remotesensing-12-00926" class="html-bibr">96</a>] on selected indices rasters from the crown based model.</p>
Full article ">Figure A1
<p>Photo illustration of the canopy score scheme for kauri used by Auckland Council in 2015 with five classes [<a href="#B136-remotesensing-12-00926" class="html-bibr">136</a>]: 1 = Healthy crown—no visible signs of dieback, 2 = Foliage/canopy thinning, 3 = Some branch dieback, 4 = Severe dieback, 5 = Dead.</p>
Full article ">
26 pages, 5947 KiB  
Article
Detection of New Zealand Kauri Trees with AISA Aerial Hyperspectral Data for Use in Multispectral Monitoring
by Jane J. Meiforth, Henning Buddenbaum, Joachim Hill, James Shepherd and David A. Norton
Remote Sens. 2019, 11(23), 2865; https://doi.org/10.3390/rs11232865 - 2 Dec 2019
Cited by 4 | Viewed by 5684
Abstract
The endemic New Zealand kauri trees (Agathis australis) are of major importance for the forests in the northern part of New Zealand. The mapping of kauri locations is required for the monitoring of the deadly kauri dieback disease (Phytophthora agathidicida [...] Read more.
The endemic New Zealand kauri trees (Agathis australis) are of major importance for the forests in the northern part of New Zealand. The mapping of kauri locations is required for the monitoring of the deadly kauri dieback disease (Phytophthora agathidicida (PTA)). In this study, we developed a method to identify kauri trees by optical remote sensing that can be applied in an area-wide campaign. Dead and dying trees were separated in one class and the remaining trees with no to medium stress symptoms were defined in the two classes “kauri” and “other”. The reference dataset covers a representative selection of 3165 precisely located crowns of kauri and 21 other canopy species in the Waitakere Ranges west of Auckland. The analysis is based on an airborne hyperspectral AISA Fenix image (437–2337 nm, 1 m2 pixel resolution). The kauri spectra show characteristically steep reflectance and absorption features in the near-infrared (NIR) region with a distinct long descent at 1215 nm, which can be parameterised with a modified Normalised Water Index (mNDWI-Hyp). With a Jeffries–Matusita separability over 1.9, the kauri spectra can be well separated from 21 other canopy vegetation spectra. The Random Forest classifier performed slightly better than Support Vector Machine. A combination of the mNDWI-Hyp index with four additional spectral indices with three red to NIR bands resulted in an overall pixel-based accuracy (OA) of 91.7% for crowns larger 3 m diameter. While the user’s and producer’s accuracies for the class “kauri” with 94.6% and 94.8% are suitable for management purposes, the separation of “dead/dying trees” from “other” canopy vegetation poses the main challenge. The OA can be improved to 93.8% by combining “kauri” and “dead/dying” trees in one class, separate classifications for low and high forest stands and a binning to 10 nm bandwidths. Additional wavelengths and their respective indices only improved the OA up to 0.6%. The method developed in this study allows an accurate location of kauri trees for an area-wide mapping with a five-band multispectral sensor in a representative selection of forest ecosystems. Full article
(This article belongs to the Section Forest Remote Sensing)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Kauri growth classes used in this study, depending on the mean crown diameter (cdm). (Photos: [<a href="#B39-remotesensing-11-02865" class="html-bibr">39</a>]).</p>
Full article ">Figure 2
<p>(<b>a</b>) Location of the Waitakere Ranges on the North Island of New Zealand west of Auckland City. The general area with naturally occurring kauri in New Zealand [<a href="#B2-remotesensing-11-02865" class="html-bibr">2</a>] is marked as hatched. (<b>b</b>) Study sites in the Waitakere Ranges with the reference crowns marked in red (background map: [<a href="#B42-remotesensing-11-02865" class="html-bibr">42</a>]).</p>
Full article ">Figure 3
<p>Reference crowns (total 3165), used in the analysis, per class and diameter.</p>
Full article ">Figure 4
<p>Mean spectra of the target classes “kauri”, “dead/dying” and “other” with standard deviations (stdev).</p>
Full article ">Figure 5
<p>Jeffries–Matusita separability [<a href="#B61-remotesensing-11-02865" class="html-bibr">61</a>] of the three target classes for different spectral ranges. A value larger than 1.9 indicates a high separability. The analysis was based on MNF transformations for all bands in the different spectral ranges.</p>
Full article ">Figure 6
<p>Mean spectra of kauri (thick black line) and six selected other canopy species (grey) that got most easily confused with kauri. The number of pixels (pix) used to generate the mean spectra is given in parentheses. The spectra of these species show the lowest separability from the kauri spectrum in this study (see <a href="#remotesensing-11-02865-t0A2" class="html-table">Table A2</a>).</p>
Full article ">Figure 7
<p>Mean spectra of kauri (black) and five other canopy species (grey) that have the highest separabilities from the kauri spectrum in this study (see <a href="#remotesensing-11-02865-t0A2" class="html-table">Table A2</a>). The number of pixels (pix) used to generate the mean spectra is given in parentheses.</p>
Full article ">Figure 8
<p>Mean spectra of the target classes “kauri”, “dead/dying” and “other” with standard deviations (“stdev”). Below: Band positions of 13 selected spectral indices.</p>
Full article ">Figure 9
<p>Performance of selected indices and index combinations to identify the class “dead/dying” (light grey) and to distinguish between “kauri” and “other vegetation” (dark grey) with an RF classification (five-fold random split, 20 repetitions). Please note that the x-axis starts at 55%.</p>
Full article ">Figure 10
<p>Performance of the final 4–8-band index combinations to distinguish the three target classes “kauri”, “dead/dying” and “other” canopy vegetation. (RF, five-fold random split, 20 repetitions). Please note that the y-axis starts at 89%.</p>
Full article ">Figure 11
<p>RGB images of the three first bands of MNF transformations [<a href="#B49-remotesensing-11-02865" class="html-bibr">49</a>] from: (<b>a</b>) the VIS to NIR1 spectral range (431–970 nm); (<b>b</b>) VIS to NIR2 (431–1327 nm); and (<b>c</b>) the full spectral range from VIS to SWIR (431–2337 nm). The importance of the NIR2 and SWIR spectrum is visible in the higher colour contrast of kauri crowns compared to the VNIR image. The numbers in the kauri polygons indicate the stress symptom class for the crown with 1 = non-symptomatic and 5 = dead.</p>
Full article ">Figure 12
<p>Histograms for selected indices on sunlit pixels for all crown diameters, with the class “kauri” marked in light blue, the class “dead/dying” in red and the class “other” in dark blue. (<b>a</b>) The histogram for the mNDWI-Hyp index, which performed best to separate the class kauri from other vegetation by capturing distinctive features in the NIR2 region, is shown. For the separation of the class “dead/dying”, indices in the RED/NIR1 region are better suited, such as (<b>b</b>) the SR800 and (<b>c</b>) the NDNI index (see <a href="#remotesensing-11-02865-t0A3" class="html-table">Table A3</a> for descriptions of these indices).</p>
Full article ">Figure 13
<p>Overall accuracies for two selected sets of six and eight bands in the visible to NIR1 range. The accuracies are calculated for two and three target classes both with and without an additional CHM layer. The results are based on an RF classification with a three-fold split in 10 repetitions on 94,971 pixel values, including small crowns (&lt;3 m diameter). The standard deviations vary from 0.12 to 0.2.</p>
Full article ">Figure 14
<p>Combined results of 10 RF classifications with a 5-fold stratified random split with different seed values. Overview (left) and detailed maps (right) for the Cascades (<b>a</b>,<b>b</b>), Maungaroa (<b>c</b>,<b>d</b>) and Kauri Grove area (<b>e</b>,<b>f</b>). The numbers indicate the symptom classes in kauri crowns (1 = non-symptomatic, 5 = dead).</p>
Full article ">
19 pages, 7132 KiB  
Article
Sequential PCA-based Classification of Mediterranean Forest Plants using Airborne Hyperspectral Remote Sensing
by Alon Dadon, Moshe Mandelmilch, Eyal Ben-Dor and Efrat Sheffer
Remote Sens. 2019, 11(23), 2800; https://doi.org/10.3390/rs11232800 - 27 Nov 2019
Cited by 19 | Viewed by 4830
Abstract
In recent years, hyperspectral remote sensing (HRS) has become common practice for remote analyses of the physiognomy and composition of forests. Supervised classification is often used for this purpose, but demands intensive sampling and analyses, whereas unsupervised classification often requires information retrieval out [...] Read more.
In recent years, hyperspectral remote sensing (HRS) has become common practice for remote analyses of the physiognomy and composition of forests. Supervised classification is often used for this purpose, but demands intensive sampling and analyses, whereas unsupervised classification often requires information retrieval out of the large HRS datasets, thereby not realizing the full potential of the technology. An improved principal component analysis-based classification (PCABC) scheme is presented and intended to provide accurate and sequential image-based unsupervised classification of Mediterranean forest species. In this study, unsupervised classification and reduction of data size are performed simultaneously by applying binary sequential thresholding to principal components, each time on a spatially reduced subscene that includes the entire spectral range. The methodology was tested on HRS data acquired by the airborne AisaFENIX HRS sensor over a Mediterranean forest in Mount Horshan, Israel. A comprehensive field-validation survey was performed, sampling 257 randomly selected individual plants. The PCABC provided highly improved results compared to the traditional unsupervised classification methodologies, reaching an overall accuracy of 91%. The presented approach may contribute to improved monitoring, management, and conservation of Mediterranean and similar forests. Full article
(This article belongs to the Special Issue Monitoring Forest Change with Remote Sensing)
Show Figures

Figure 1

Figure 1
<p>Mount Horshan, northern Israel, aerial photo (WorldView 2 satellite image) overlaid with the Specim AisaFENIX image. The red dots mark the location of 257 ground-truth validation points.</p>
Full article ">Figure 2
<p>PCABC flowchart presenting all stages of the classification process.</p>
Full article ">Figure 3
<p>Illustration of the process for producing the non-vegetation mask by PCABC method in a subscene of the image (marked in red square). (<b>a</b>) Keren Kayemeth LeIsrael (KKL) orthophoto of the subscene showing that the pixels that were masked out are indeed related to roads and shade (non-vegetation pixels). (<b>b</b>) Specim AisaFENIX RGB. (<b>c</b>) PCA first iteration 1, component 1 (non-vegetation pixels appear in black). (<b>d</b>) Marking of the non-vegetation pixels using DS (non-vegetation pixels appear in cyan, yellow and blue colors; DS values: −51.84 to −5.03). (<b>e</b>) Non-vegetation pixels overlaid on RGB image. (<b>f</b>) Final product of the non-vegetation pixel mask applied to the image.</p>
Full article ">Figure 4
<p>Illustration of the process for classifying plant species using the PCABC method in a subscene of the image (marked with a red square). (<b>a</b>) PCA iteration 1, component 1, with no differences among the vegetation pixels. (<b>b</b>) PCA iteration 2, component 1, with clear differences among the vegetation pixels (higher PC values appear in white). (<b>c</b>) Marking of plant species using the DS tool; vegetation pixels marked in different colors (based on the legend). (<b>d</b>) Example of one of the classes (marked in cyan) identified based on the highest DS values.</p>
Full article ">Figure 5
<p>Illustration of combining the highest DS values into one DS cluster. (<b>a</b>) The final plant cluster that was produced in the second iteration of the PCABC process. (<b>b</b>) The tree crowns marked with different DS values (before the merger). The full DS values of the cluster are shown in the legend.</p>
Full article ">Figure 6
<p>Subscenes of the results of classification using two different classifiers on the PCA components image: (<b>a</b>) K-means and (<b>b</b>) ISODATA.</p>
Full article ">Figure 7
<p>Subscenes of the six classes detected using the PCABC methodology. Different colors represent the locations of the plant species in the image. These were later identified as: (<b>a</b>) <span class="html-italic">P. halepensis.</span> (<b>b</b>) Trees covered by lianas. (<b>c</b>) <span class="html-italic">P. lentiscus</span>. (<b>d</b>) Shrubs. (<b>e</b>) <span class="html-italic">Q. calliprinos</span>. (<b>f</b>) <span class="html-italic">Q. ithaburensis</span>.</p>
Full article ">Figure 8
<p>Thematic map of the PCABC plant species classes.</p>
Full article ">Figure 9
<p>A small <span class="html-italic">Pinus halepensis</span> tree identified by PCABC. (<b>a</b>) <span class="html-italic">P. halepensis</span> class based on PCABC (marked in cyan color); red circle marks a small tree with a height of 1.5 m and canopy diameter of 1 m. (<b>b</b>) The <span class="html-italic">P. halepensis</span> tree in the field.</p>
Full article ">Figure 10
<p>Example of a tree covered by lianas. (<b>a</b>) Lianas class based on PCABC (marked in magenta color); red circle marks the location of the photographed tree in the field. (<b>b</b>) Tree covered by lianas in the field.</p>
Full article ">Figure 11
<p>Illustration of similarity of the spectra of detected classes. (<b>a</b>) Plots of mean spectra for the six plant classes. (<b>b</b>) Mean and standard deviation spectra for <span class="html-italic">Q. calliprinos</span> and <span class="html-italic">Q. ithaburensis</span> classes. (<b>c</b>) Means for <span class="html-italic">Q. calliprinos</span> (blue line) and <span class="html-italic">Q. ithaburensis</span> (orange line). Gray line represents the ratio of the two means. Results indicate high similarity between the two species, reflected by a nearly straight line, with values of above 0.8 (reflectance).</p>
Full article ">
13 pages, 3226 KiB  
Article
Mapping Asphaltic Roads’ Skid Resistance Using Imaging Spectroscopy
by Nimrod Carmon and Eyal Ben-Dor
Remote Sens. 2018, 10(3), 430; https://doi.org/10.3390/rs10030430 - 10 Mar 2018
Cited by 23 | Viewed by 5930
Abstract
The purpose of this study is to evaluate a realistic feasibility of using hyperspectral remote sensing (also termed imaging spectroscopy) airborne data for mapping asphaltic roads’ transportation safety. This is done by quantifying the road-tire friction, an attribute responsible for vehicle control and [...] Read more.
The purpose of this study is to evaluate a realistic feasibility of using hyperspectral remote sensing (also termed imaging spectroscopy) airborne data for mapping asphaltic roads’ transportation safety. This is done by quantifying the road-tire friction, an attribute responsible for vehicle control and emergency stopping. We engaged in a real-life operational scenario, where the roads’ friction was modeled against the reflectance information extracted directly from the image. The asphalt pavement’s dynamic friction coefficient was measured by a standardized technique using a Dynatest 6875H (Dynatest Consulting Inc., Westland, MI, USA) Friction Measuring System, which uses the common test-wheel retardation method. The hyperspectral data was acquired by the SPECIM AisaFenix 1K (Specim, Spectral Imaging Ltd., Oulu, Finland) airborne system, covering the entire optical range (350–2500 nm), over a selected study site, with roads characterized by different aging conditions. The spectral radiance data was processed to provide apparent surface reflectance using ground calibration targets and the ACORN-6 atmospheric correction package. Our final dataset was comprised of 1370 clean asphalt pixels coupled with geo-rectified in situ friction measurement points. We developed a partial least squares regression model using PARACUDA-II spectral data mining engine, which uses an automated outlier detection procedure and dual validation routines—a full cross-validation and an iterative internal validation based on a Latin Hypercube sampling algorithm. Our results show prediction capabilities of R2 = 0.632 for full cross-validation and R2 = 0.702 for the best available model in internal validation, both with significant results (p < 0.0001). Using spectral assignment analysis, we located the spectral bands with the highest weight in the model and discussed their possible physical and chemical assignments. The derived model was applied back on the hyperspectral image to predict and map the friction values of every road pixel in the scene. Combining the standard method with imaging spectroscopy may provide the required expansion of the available data to furnish decision makers with a full picture of the roads’ status. This technique’s limitations originate mainly in compositional variations between different roads, and the requirement for the application of multiple calibrations between scenes. Possible improvements could be achieved by using more spectral regions (e.g., thermal) and additional remote sensing techniques (e.g., LIDAR) as well as new platforms (e.g., UAV). Full article
(This article belongs to the Special Issue Recent Progress and Developments in Imaging Spectroscopy)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Location of the study area in central Israel with a true color image from the AisaFenix 1K hyperspectral sensor, and background greyscale image from sentinel-2 (Copernicus Sentinel Data 2017). Overview satellite image from MODIS (Rapid Response Team, NASA/GSFC).</p>
Full article ">Figure 2
<p>Data exploration plots colored by group: (<b>a</b>) Average spectra by friction group with labeled group ranks and CR range in dotted box; (<b>b</b>) CR average spectra (2150–2400 nm) by friction group rank; (<b>c</b>) Histogram of friction values colored by corresponding group rank (see <a href="#remotesensing-10-00430-t002" class="html-table">Table 2</a> for group definition).</p>
Full article ">Figure 3
<p>Outlier detection module results: (<b>a</b>) PC-1/PC-2 z-score scatter plot of the spectral data with a ±3 confidence circle; (<b>b</b>) histogram of friction value before outlier elimination; (<b>c</b>) histogram of friction value after outlier elimination.</p>
Full article ">Figure 4
<p>Model performance results: (<b>a</b>) Measured vs. predicted friction values for full cross-validation; (<b>b</b>) Measured vs. predicted values for internal validation.</p>
Full article ">Figure 5
<p>Beta coefficient spectrum representing the band weight in the model.</p>
Full article ">Figure 6
<p>Roads predicted friction map of the study site.</p>
Full article ">
5044 KiB  
Article
Atmospheric Correction Performance of Hyperspectral Airborne Imagery over a Small Eutrophic Lake under Changing Cloud Cover
by Lauri Markelin, Stefan G. H. Simis, Peter D. Hunter, Evangelos Spyrakos, Andrew N. Tyler, Daniel Clewley and Steve Groom
Remote Sens. 2017, 9(1), 2; https://doi.org/10.3390/rs9010002 - 23 Dec 2016
Cited by 16 | Viewed by 9213
Abstract
Atmospheric correction of remotely sensed imagery of inland water bodies is essential to interpret water-leaving radiance signals and for the accurate retrieval of water quality variables. Atmospheric correction is particularly challenging over inhomogeneous water bodies surrounded by comparatively bright land surface. We present [...] Read more.
Atmospheric correction of remotely sensed imagery of inland water bodies is essential to interpret water-leaving radiance signals and for the accurate retrieval of water quality variables. Atmospheric correction is particularly challenging over inhomogeneous water bodies surrounded by comparatively bright land surface. We present results of AisaFENIX airborne hyperspectral imagery collected over a small inland water body under changing cloud cover, presenting challenging but common conditions for atmospheric correction. This is the first evaluation of the performance of the FENIX sensor over water bodies. ATCOR4, which is not specifically designed for atmospheric correction over water and does not make any assumptions on water type, was used to obtain atmospherically corrected reflectance values, which were compared to in situ water-leaving reflectance collected at six stations. Three different atmospheric correction strategies in ATCOR4 was tested. The strategy using fully image-derived and spatially varying atmospheric parameters produced a reflectance accuracy of ±0.002, i.e., a difference of less than 15% compared to the in situ reference reflectance. Amplitude and shape of the remotely sensed reflectance spectra were in general accordance with the in situ data. The spectral angle was better than 4.1° for the best cases, in the spectral range of 450–750 nm. The retrieval of chlorophyll-a (Chl-a) concentration using a popular semi-analytical band ratio algorithm for turbid inland waters gave an accuracy of ~16% or 4.4 mg/m3 compared to retrieval of Chl-a from reflectance measured in situ. Using fixed ATCOR4 processing parameters for whole images improved Chl-a retrieval results from ~6 mg/m3 difference to reference to approximately 2 mg/m3. We conclude that the AisaFENIX sensor, in combination with ATCOR4 in image-driven parametrization, can be successfully used for inland water quality observations. This implies that the need for in situ reference measurements is not as strict as has been assumed and a high degree of automation in processing is possible. Full article
(This article belongs to the Special Issue Water Optics and Water Colour Remote Sensing)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>(<b>a</b>) map showing the location of Loch Leven within the UK; (<b>b</b>) Image mosaic from 10 FENIX flight lines collected during the campaign (colours are illustrative only). Locations of the in situ station measurements are marked with red crosses; the location of station ST2 is approximate.</p>
Full article ">Figure 2
<p>(<b>a</b>) In situ E<sub>d</sub>(PAR) and AOT measurements taken during the campaign. In situ observations used for station reference spectra are marked with red diamonds. Orange shading indicates the times of airborne flight lines, and blue shading the times of station measurements; (<b>b</b>) Median and standard deviation of E<sub>d</sub>(PAR) from all 30 measurements recorded at each station; (<b>c</b>) SkyRat based on QC passed measurements for calculation of R<sub>rs</sub> (number of valid observations are shown in horizontal axis label). The SkyRat error bars for stations ST2-ST4 are not shown, as there were &lt;3 valid measurements at these stations.</p>
Full article ">Figure 3
<p>(<b>a</b>) R<sub>w</sub> spectra measured from atmospherically corrected FENIX images. Spectra measured at station locations are shown in colour, grey spectra are an additional 50 measurements representing all lake areas to indicate the homogeneity of the lake; (<b>b</b>) Image-based R<sub>w</sub> spectra of water measured from 21 locations under cloud shadows, the black line indicates station ST6 on image FL6. Wavelength range in the x-axis goes up to 1300 nm to indicate the non-zero reflectance in the NIR-SWIR range.</p>
Full article ">Figure 4
<p>(<b>a</b>–<b>f</b>) In situ station and image-derived R<sub>w</sub> of each matchup. The best matchups in terms of Chi-square and spectral angle are shown in a thicker line in panels (<b>b</b>–<b>d</b>,<b>f</b>); Image R<sub>w</sub> from atmospheric correction strategies AC2 and AC3 are included in panels (<b>b</b>,<b>c</b>); Bracketed values in the panel legends state the number of spectra used to calculate the in situ spectrum and error bars; td = time difference between in situ measurement and image acquisition, sd = distance between in situ and image matchup measurement, if applicable. Standard deviations for the in situ station measurements are plotted as error bars (blue shading) where available.</p>
Full article ">Figure 5
<p>Spectral signal-to-noise ratio (SNR) (blue circles) and at-sensor radiance (green line) of AisaFENIX measured over Loch Leven. Curves are mean of 17 individual measurements from 7 flight lines.</p>
Full article ">Figure 6
<p>(<b>a</b>–<b>f</b>) In situ station and image-derived R<sub>w</sub> spectra standardized between 0 and 1. Image spectra corresponding to the best five comparisons are plotted as thicker lines. Image-derived standardized R<sub>w</sub> from atmospheric correction strategies AC2 and AC3 are included in panels (<b>b</b>,<b>c</b>). In panel legends, td = time difference between in situ measurement and image acquisition, and sd = distance between in situ and image measurement.</p>
Full article ">Figure 6 Cont.
<p>(<b>a</b>–<b>f</b>) In situ station and image-derived R<sub>w</sub> spectra standardized between 0 and 1. Image spectra corresponding to the best five comparisons are plotted as thicker lines. Image-derived standardized R<sub>w</sub> from atmospheric correction strategies AC2 and AC3 are included in panels (<b>b</b>,<b>c</b>). In panel legends, td = time difference between in situ measurement and image acquisition, and sd = distance between in situ and image measurement.</p>
Full article ">Figure 7
<p>(<b>a</b>) Difference between image-derived and in situ R<sub>w</sub>; (<b>b</b>) ratio of in situ over image-derived R<sub>w</sub>; (<b>c</b>) RMS of the best five and all matchups, respectively, in units of reflectance; (<b>d</b>) Relative RMS. Matchup labels for panes (<b>a</b>,<b>b</b>) are given in the legend of panel (<b>b</b>), and the five best matchups are plotted with a solid line.</p>
Full article ">Figure 8
<p>(<b>a</b>) Chl-a values in (mg/m<sup>3</sup>) derived from all measured image R<sub>w</sub>; (<b>b</b>) scatter plot between Chl-a values derived from in situ R<sub>w</sub> and image R<sub>w</sub> at station locations.</p>
Full article ">
6797 KiB  
Technical Note
Integration of Hyperspectral Shortwave and Longwave Infrared Remote-Sensing Data for Mineral Mapping of Makhtesh Ramon in Israel
by Gila Notesco, Yaron Ogen and Eyal Ben-Dor
Remote Sens. 2016, 8(4), 318; https://doi.org/10.3390/rs8040318 - 9 Apr 2016
Cited by 22 | Viewed by 6449
Abstract
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 [...] Read more.
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. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>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 <span class="html-italic">et al.</span> [<a href="#B7-remotesensing-08-00318" class="html-bibr">7</a>]. The XRD analysis results are shown below in <a href="#remotesensing-08-00318-t001" class="html-table">Table 1</a>.</p>
Full article ">Figure 2
<p>(<b>a</b>) Spectra of regions of interest (ROIs) A and B from the normalized reflectance image and of clay minerals from a spectral library (Slib.) [<a href="#B12-remotesensing-08-00318" class="html-bibr">12</a>] for comparison, with the indicative absorption feature at 2.20 µm; (<b>b</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.</p>
Full article ">Figure 3
<p>Mineral map of the surface covered by the two flight lines in Makhtesh Ramon based on the SWIR (<b>a</b>) and LWIR (<b>b</b>) images. K—kaolinite, Ca—calcite, Do—dolomite, Q—quartz, F—feldspars, Cm—clay minerals, G—gypsum, C—carbonates, FQ—feldspars + quartz, GQC—gypsum + quartz (+ carbonates), QC—quartz + carbonates.</p>
Full article ">Figure 4
<p>Identification of minerals in both SWIR and LWIR data sets or in only one—SWIR or LWIR—data set.</p>
Full article ">Figure 5
<p>SWIR reflectance spectra (<b>a</b>) and LWIR approximate emissivity spectra (<b>b</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.</p>
Full article ">Figure 6
<p>Mineral map of the surface of Makhtesh Ramon on the geological map of the area (source: Sneh, A., Bartov, Y., Weissbrod, T. and Rosensaft, M., 1998. Geological Map of Israel, 1:200,000. Isr. Geol. Surv.). The table shows the rock types in the study area and the mineral classification based on the hyperspectral data: Q—quartz, Ca—calcite, Do—dolomite, QC—quartz + carbonates (calcite, dolomite), K—kaolinite, QK—quartz + kaolinite, F—feldspars, FQ—feldspars + quartz, G—gypsum, GQC—gypsum + quartz (+ carbonates).</p>
Full article ">
Back to TopTop