[go: up one dir, main page]
More Web Proxy on the site http://driver.im/
Next Issue
Volume 14, July-1
Previous Issue
Volume 14, June-1
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
remotesensing-logo

Journal Browser

Journal Browser

Remote Sens., Volume 14, Issue 12 (June-2 2022) – 229 articles

Cover Story (view full-size image): Carpinteria Salt Marsh Reserve, like many coastal salt marshes, is a dynamic and productive ecosystem that provides a wide variety of ecosystem services for coastal environments and communities. However, disturbances such as debris flows and sea level rise have the potential to degrade those services. In the absence of field data, we employed Sentinel-2 imagery and random forest classification to quantify landcover change associated with the Montecito Debris Flows of 2018. While total vegetated area remained constant after debris flow, the proportion of the high marsh vegetation community decreased, accompanied by a potential loss of species diversity. Such plant community shifts, identifiable through post-classification change detection, may negatively impact marsh function and resilience, especially in deposition-prone wetlands. View this paper
  • Issues are regarded as officially published after their release is announced to the table of contents alert mailing list.
  • You may sign up for e-mail alerts to receive table of contents of newly released issues.
  • PDF is the official format for papers published in both, html and pdf forms. To view the papers in pdf format, click on the "PDF Full-text" link, and use the free Adobe Reader to open them.
Order results
Result details
Section
Select all
Export citation of selected articles as:
14 pages, 7560 KiB  
Article
Evaluation of Tidal Effect in Long-Strip DInSAR Measurements Based on GPS Network and Tidal Models
by Wei Peng, Qijie Wang, Yunmeng Cao, Xuemin Xing and Wenjie Hu
Remote Sens. 2022, 14(12), 2954; https://doi.org/10.3390/rs14122954 - 20 Jun 2022
Cited by 6 | Viewed by 2160
Abstract
A long-strip differential interferometric synthetic aperture radar (DInSAR) measurement based on multi-frame image mosaicking is currently the realizable approach to measure large-scale ground deformation. As the spatial range of the mosaicked images increases, the nonlinear variation of ground ocean tidal loading (OTL) displacements [...] Read more.
A long-strip differential interferometric synthetic aperture radar (DInSAR) measurement based on multi-frame image mosaicking is currently the realizable approach to measure large-scale ground deformation. As the spatial range of the mosaicked images increases, the nonlinear variation of ground ocean tidal loading (OTL) displacements is more significant, and using plane fitting to remove the large-scale errors will produce large tidal displacement residuals in a region with a complex coastline. To conveniently evaluate the ground tidal effect on mosaic DInSAR interferograms along the west coast of the U.S., a three-dimensional ground OTL displacements grid is generated by integrating tidal constituents’ estimation of the GPS reference station network and global/regional ocean tidal models. Meanwhile, a solid earth tide (SET) model based on IERS conventions is used to estimate the high-precision SET displacements. Experimental results show that the OTL and SET in a long-strip interferogram can reach 77.5 mm, which corresponds to a 19.3% displacement component. Furthermore, the traditional bilinear ramp fitting methods will cause 7.2~20.3 mm residual tidal displacement in the mosaicked interferograms, and the integrated tidal constituents displacements calculation method can accurately eliminate the tendency of tidal displacement in the long-strip interferograms. Full article
Show Figures

Figure 1

Figure 1
<p>The standard deviation map of global ground OTL displacement with an interval of 300 s in the U, N, and E directions based on the FES2014b model.</p>
Full article ">Figure 2
<p>Distribution of 1038 continuous GPS sites and multiple Sentinel-1 SLC image ranges on the west coast of the U.S.</p>
Full article ">Figure 3
<p>Comparative analysis of the phasor of tidal constituents estimated by the kinematic PPP technique and FES2014b+osu.usawest models.</p>
Full article ">Figure 4
<p>StdDev maps of the difference of the OTL displacements estimated by the proposed tide calculation method and FES2014b+osu.usawest model in U, N, and E directions, and the areas with lager difference (red box).</p>
Full article ">Figure 5
<p>The OTL displacement StdDev along the LOS direction (Incident angle is 39°) of the ascending (Heading direction is −13°) and descending (Heading direction is 193°) Sentinel-1 measurements with an interval of 12 days based on the proposed OTL calculation method.</p>
Full article ">Figure 6
<p>The SET and OTL displacements in 29 long-strip differential interferograms and the error analysis of residual tidal displacement produced by bilinear fitting ramp, and the interferogram with largest residual tidal displacement (red box).</p>
Full article ">Figure 7
<p>OTL displacement (<b>a1</b>), SET (<b>a2</b>) and their superposition displacement (<b>a3</b>) based on the differential interferogram acquisition on 6 September 2018 and 12 October 2018; (<b>a4</b>) residual tidal displacement generated by the bilinear ramp fitting method; and (<b>a5</b>) residual tidal displacement from the multi-frame mosaic bilinear ramp fitting method.</p>
Full article ">Figure 8
<p>(<b>a</b><b>1</b>−<b>a3</b>) The long-strip differential interferograms after the correction of (<b>b</b><b>1</b>−<b>b3</b>) atmospheric delay error, (<b>c</b><b>1</b>−<b>c3</b>) SET and OTL, and (<b>d</b><b>1</b>−<b>d3</b>) bilinear ramp elimination. (<b>e</b><b>1</b>−<b>e3</b>) The SET and OTL displacements, whose spatial trend variations are similar to the long-strip differential interferograms after the atmospheric delay is corrected.</p>
Full article ">Figure 9
<p>The comparison of the displacement residuals in the long-strip differential interferograms after the OTL and SET displacements corrected using the proposed tidal method (interferograms-atmospheric-SET-OTL-bilinear fitting ramp) and traditional bilinear ramp fitting method (interferograms-atmospheric-bilinear fitting ramp).</p>
Full article ">Figure 10
<p>The comparison of the displacement residuals in a complex coastline area of the long-strip differential interferograms after the tidal displacements was corrected using the proposed tidal method and FES2014b+osu.usawest model.</p>
Full article ">
20 pages, 5636 KiB  
Article
Classification of Electronic Devices Using a Frequency-Swept Harmonic Radar Approach
by Handan Ilbegi, Halil Ibrahim Turan, Imam Samil Yetik and Harun Taha Hayvaci
Remote Sens. 2022, 14(12), 2953; https://doi.org/10.3390/rs14122953 - 20 Jun 2022
Cited by 5 | Viewed by 2363
Abstract
A new method to classify electronic devices using a Frequency-Swept Harmonic Radar (FSHR) approach is proposed in this paper. The FSHR approach enables us to utilize the frequency diversity of the harmonic responses of the electronic circuits. Unlike previous studies, a frequency-swept signal [...] Read more.
A new method to classify electronic devices using a Frequency-Swept Harmonic Radar (FSHR) approach is proposed in this paper. The FSHR approach enables us to utilize the frequency diversity of the harmonic responses of the electronic circuits. Unlike previous studies, a frequency-swept signal with a constant power is transmitted to Electronic Circuits Under Test (ECUTs). The harmonic response to a frequency-swept transmitted signal is found to be distinguishable for different types of ECUTs. Statistical and Fourier features of the harmonic responses are derived for classification. Later, the harmonic characteristics of the ECUTs are depicted in 3D harmonic and feature spaces for classification. Three-dimensional harmonic and feature spaces are composed of the first three harmonics of the re-radiated signal and the statistical or Fourier features, respectively. We extensively evaluate the performance of our novel method through Monte Carlo simulations in the presence of noise. Full article
(This article belongs to the Special Issue Nonlinear Junction Detection and Harmonic Radar)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Schematics of the ECUTs.</p>
Full article ">Figure 2
<p>Illustration of the Frequency-Swept Harmonic Radar (FSHR) approach.</p>
Full article ">Figure 3
<p>Received powers at first three harmonics of the cascaded amplifier circuits.</p>
Full article ">Figure 4
<p>Received powers at first three harmonics of the common emitter amplifier circuits.</p>
Full article ">Figure 5
<p>Received powers at first three harmonics of the sawtooth oscillator circuits.</p>
Full article ">Figure 6
<p>Scatter plots of variance, skewness, and kurtosis data in 3D harmonic space, SNR = 2 dB.</p>
Full article ">Figure 7
<p>Scatter plots of variance, skewness, and kurtosis data in 3D harmonic space, SNR = 14 dB.</p>
Full article ">Figure 8
<p>Classification performance of the statistical features for varying SNR values.</p>
Full article ">Figure 9
<p>Statistical feature space representation of second and third harmonics data, SNR = 2 dB.</p>
Full article ">Figure 10
<p>Classification performance of the harmonics in statistical feature space for varying SNR values.</p>
Full article ">Figure 11
<p>Classification performance presented in [<a href="#B9-remotesensing-14-02953" class="html-bibr">9</a>]: (<b>a</b>) Classification performance via statistical features [<a href="#B9-remotesensing-14-02953" class="html-bibr">9</a>] and (<b>b</b>) classification performance via each harmonic [<a href="#B9-remotesensing-14-02953" class="html-bibr">9</a>].</p>
Full article ">Figure 12
<p>Scatter plots of Fourier features in 3D harmonic space, SNR = 2 dB.</p>
Full article ">Figure 13
<p>Scatter plots of Fourier features in 3D harmonic space, SNR = 14 dB.</p>
Full article ">Figure 14
<p>Classification performance of the Fourier features in harmonic space for varying SNR values.</p>
Full article ">Figure 15
<p>Fourier feature space representation of second and third harmonics data, SNR = 2 dB.</p>
Full article ">Figure 16
<p>Fourier feature space representation of second and third harmonics data, SNR = 14 dB.</p>
Full article ">Figure 17
<p>Classification performance of the harmonics in Fourier feature space for varying SNR values.</p>
Full article ">Figure 18
<p>Classification performance via Fourier features in [<a href="#B9-remotesensing-14-02953" class="html-bibr">9</a>]: (<b>a</b>) classification performance via Fourier features for each frequency level in harmonic space [<a href="#B9-remotesensing-14-02953" class="html-bibr">9</a>] and (<b>b</b>) classification performance via Fourier features for each harmonic [<a href="#B9-remotesensing-14-02953" class="html-bibr">9</a>].</p>
Full article ">
18 pages, 9908 KiB  
Article
Spectroscopic and Petrographic Investigations of Lunar Mg-Suite Meteorite Northwest Africa 8687
by Lang Qin, Xing Wu, Liying Huang, Yang Liu and Yongliao Zou
Remote Sens. 2022, 14(12), 2952; https://doi.org/10.3390/rs14122952 - 20 Jun 2022
Cited by 2 | Viewed by 3002
Abstract
Magnesian suite (Mg-suite) rocks represent plutonic materials from the lunar crust, and their global distribution can provide critical information for the early magmatic differentiation and crustal asymmetries of the Moon. Visible and near-infrared (VNIR) spectrometers mounted on orbiters and rovers have been proven [...] Read more.
Magnesian suite (Mg-suite) rocks represent plutonic materials from the lunar crust, and their global distribution can provide critical information for the early magmatic differentiation and crustal asymmetries of the Moon. Visible and near-infrared (VNIR) spectrometers mounted on orbiters and rovers have been proven to be powerful approaches for planetary mineral mapping, which are instrumental in diagnosing Mg-suite rocks. However, due to the scarcity and diversity of Mg-suite samples, laboratory measurements with variable proportions of minerals are imperative for spectral characterization. In this study, spectroscopic investigation and petrographic study were conducted on lunar Mg-suite meteorite Northwest Africa 8687. We classify the sample as a pink spinel-bearing anorthositic norite through spectral and petrographic characteristics. The ground-truth information of the Mg-suite rock is provided for future exploration. Meanwhile, the results imply that the VNIR technique has the potential to identify highland rock types by mineral modal abundance and could further be applied in extraterrestrial samples for primary examination due to its advantage of being fast and non-destructive. Full article
(This article belongs to the Special Issue Planetary Geologic Mapping and Remote Sensing)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Photograph of lunar meteorite Northwest Africa 8687 and the location of point spectral measurements (<b>A</b>–<b>E</b>) in this study. The frame shows the rough location of the polished thin section for the petrographic study.</p>
Full article ">Figure 2
<p>(<b>a</b>) Visible and near-infrared (VNIR) reflectance spectra and (<b>b</b>) spectra after continuum removal of endmembers used for spectral unmixing. Spectral ID: clinopyroxene (C1DL84A), orthopyroxene (C1LR209), plagioclase (C1LR223), maskelynite (C1LR222), ilmenite (C1LR222), and olivine (C2PO81).</p>
Full article ">Figure 3
<p>(<b>a</b>) backscattered electron (BSE) image of a polished section of NWA 8687. (<b>b</b>) Distribution of main minerals in NWA 8687 based on X-ray element maps.</p>
Full article ">Figure 4
<p>Backscattered electron (BSE) images of (<b>a</b>) an olivine phenocryst. The frame shows the location of figure (<b>c</b>). (<b>b</b>) A Mg-spinel in the matrix. (<b>c</b>) Fine-grained matrix, and (<b>d</b>) some mineral relicts are visible in the impact melt vein.</p>
Full article ">Figure 5
<p>Pyroxene (<b>a</b>) and olivine (<b>b</b>) composition for NWA 8687, compared with a counterpart from pairs NWA 5744 and NWA 10401 [<a href="#B15-remotesensing-14-02952" class="html-bibr">15</a>,<a href="#B32-remotesensing-14-02952" class="html-bibr">32</a>]. Di—diopside, Hd—hedenbergite, En—enstatite, Fs—ferrosilite, and Fo—forsterite.</p>
Full article ">Figure 6
<p>Spinel composition in NWA 8687 relative to compositional range from M<sup>3</sup> (after Ref. [<a href="#B21-remotesensing-14-02952" class="html-bibr">21</a>]); also compared with spinel in pair stone NWA 10401 [<a href="#B15-remotesensing-14-02952" class="html-bibr">15</a>], as well as in mare basalts [<a href="#B45-remotesensing-14-02952" class="html-bibr">45</a>] and pink spinel troctolites [<a href="#B46-remotesensing-14-02952" class="html-bibr">46</a>,<a href="#B47-remotesensing-14-02952" class="html-bibr">47</a>,<a href="#B48-remotesensing-14-02952" class="html-bibr">48</a>,<a href="#B49-remotesensing-14-02952" class="html-bibr">49</a>,<a href="#B50-remotesensing-14-02952" class="html-bibr">50</a>,<a href="#B51-remotesensing-14-02952" class="html-bibr">51</a>,<a href="#B52-remotesensing-14-02952" class="html-bibr">52</a>,<a href="#B53-remotesensing-14-02952" class="html-bibr">53</a>,<a href="#B54-remotesensing-14-02952" class="html-bibr">54</a>].</p>
Full article ">Figure 7
<p>(<b>a</b>) Visible and near-infrared (VNIR) reflectance spectra of NWA 8687 chip surface spots marked in <a href="#remotesensing-14-02952-f001" class="html-fig">Figure 1</a>. (<b>b</b>) VNIR spectra after continuum removal.</p>
Full article ">Figure 8
<p>Comparisons of NWA 8687 and the best-matched RELAB reflectance spectra of returned lunar sample 60019, 214, and lunar meteorite ALHA 81005 using SAM algorithm. (<b>a</b>) Reflectance. (<b>b</b>) Spectra normalized at 750 nm, and (<b>c</b>) spectra after continuum removal. Spectrum ID for sample 60019.214 and ALHA 81005 is S3LS07 and CFLM34, respectively.</p>
Full article ">Figure 9
<p>(<b>a</b>) Plot of the Band I center versus the Band II center. The laboratory pyroxene spectra are from Ref. [<a href="#B58-remotesensing-14-02952" class="html-bibr">58</a>]. (<b>b</b>) Plot of the Band Area Ratio versus the Band I center. The olivine–orthopyroxene mixing line is represented by the solid line. Different spectral regions are defined by Refs. [<a href="#B37-remotesensing-14-02952" class="html-bibr">37</a>,<a href="#B59-remotesensing-14-02952" class="html-bibr">59</a>].</p>
Full article ">Figure 10
<p>Comparisons between the measured point spectra (<b>A</b>–<b>E</b>) (in blue) and their modeled spectra (in red) of the NWA 8687 chip.</p>
Full article ">Figure 11
<p>Mn versus Fe<sup>2+</sup> atoms per formula unit (afu) in pyroxene (<b>a</b>) and olivine (<b>b</b>) of NWA 8687, compared with pairs NWA 5744 [<a href="#B32-remotesensing-14-02952" class="html-bibr">32</a>] and NWA 10401 [<a href="#B15-remotesensing-14-02952" class="html-bibr">15</a>]. Apollo data and best-fit trend lines are from Ref. [<a href="#B60-remotesensing-14-02952" class="html-bibr">60</a>] (pyroxene) and Ref. [<a href="#B61-remotesensing-14-02952" class="html-bibr">61</a>] (olivine).</p>
Full article ">Figure 12
<p>(<b>a</b>) Mg# in olivine versus An in plagioclase for NWA 8687 and pairs NWA 5744 and NWA 10401 (after Ref. [<a href="#B13-remotesensing-14-02952" class="html-bibr">13</a>]). (<b>b</b>) Bulk Mg# versus Al<sub>2</sub>O<sub>3</sub> content for NWA 8687 and other lunar samples (after Refs. [<a href="#B15-remotesensing-14-02952" class="html-bibr">15</a>,<a href="#B16-remotesensing-14-02952" class="html-bibr">16</a>,<a href="#B32-remotesensing-14-02952" class="html-bibr">32</a>] and references therein). Lines show mixtures among four endmembers: a hypothetical magnesian peridotite, Mg-suite rocks, ferroan anorthosite, and mare basalts. Data source: NWA 5744 [<a href="#B32-remotesensing-14-02952" class="html-bibr">32</a>], NWA 10401 [<a href="#B15-remotesensing-14-02952" class="html-bibr">15</a>], ALH 81005 [<a href="#B16-remotesensing-14-02952" class="html-bibr">16</a>], and other lunar meteorites [<a href="#B33-remotesensing-14-02952" class="html-bibr">33</a>]; Apollo Mg-suite [<a href="#B68-remotesensing-14-02952" class="html-bibr">68</a>], FAN, and mare basalts [<a href="#B9-remotesensing-14-02952" class="html-bibr">9</a>].</p>
Full article ">Figure 13
<p>Classification of lunar meteorite NWA 8687 and comparison with pairs NWA 5744 and NWA 10401. Modal abundance (in vol.%) obtained by petrography (star) and by spectroscopy (blue circle) is expressed in the plot. This figure is modified from <a href="#remotesensing-14-02952-f002" class="html-fig">Figure 2</a> in Ref. [69].</p>
Full article ">
25 pages, 9355 KiB  
Article
Turbulence Detection in the Atmospheric Boundary Layer Using Coherent Doppler Wind Lidar and Microwave Radiometer
by Pu Jiang, Jinlong Yuan, Kenan Wu, Lu Wang and Haiyun Xia
Remote Sens. 2022, 14(12), 2951; https://doi.org/10.3390/rs14122951 - 20 Jun 2022
Cited by 15 | Viewed by 3718
Abstract
The refractive index structure constant (Cn2) is a key parameter used in describing the influence of turbulence on laser transmissions in the atmosphere. Three different methods for estimating Cn2 were analyzed in detail. A new method that [...] Read more.
The refractive index structure constant (Cn2) is a key parameter used in describing the influence of turbulence on laser transmissions in the atmosphere. Three different methods for estimating Cn2 were analyzed in detail. A new method that uses a combination of these methods for continuous Cn2 profiling with both high temporal and spatial resolution is proposed and demonstrated. Under the assumption of the Kolmogorov “2/3 law”, the Cn2 profile can be calculated by using the wind field and turbulent kinetic energy dissipation rate (TKEDR) measured by coherent Doppler wind lidar (CDWL) and other meteorological parameters derived from a microwave radiometer (MWR). In a horizontal experiment, a comparison between the results from our new method and measurements made by a large aperture scintillometer (LAS) is conducted. The correlation coefficient, mean error, and standard deviation between them in a six-day observation are 0.8073, 8.18 × 10−16 m−2/3 and 1.27 × 10−15 m−2/3, respectively. In the vertical direction, the continuous profiling results of Cn2 and other turbulence parameters with high resolution in the atmospheric boundary layer (ABL) are retrieved. In addition, the limitation and uncertainty of this method under different circumstances were analyzed, which shows that the relative error of Cn2 estimation normally does not exceed 30% under the convective boundary layer (CBL). Full article
(This article belongs to the Special Issue Lidar for Advanced Classification and Retrieval of Aerosols)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Illustration of two experiment sites in the satellite map and the instruments layout in the horizontal experiment (USTC, Hefei City).</p>
Full article ">Figure 2
<p>The horizontal wind speed (<b>a</b>), wind direction (<b>b</b>), vertical wind speed (<b>c</b>), wind shear (<b>d</b>), log<sub>10</sub>(TKEDR) (<b>e</b>), temperature and its gradient (<b>f</b>), variances in vertical wind fluctuation and flux of potential temperature (<b>g</b>), log<sub>10</sub><math display="inline"><semantics> <mrow> <msubsup> <mi>C</mi> <mi>n</mi> <mn>2</mn> </msubsup> </mrow> </semantics></math> estimated by different methods (<b>h</b>), Brunt–Väisälä frequency squared, log<sub>10</sub><math display="inline"><semantics> <mrow> <msubsup> <mi>C</mi> <mi>n</mi> <mn>2</mn> </msubsup> </mrow> </semantics></math> and its shaded area error bar (<b>i</b>) retrieved from CDWL, wind tower, and LAS in the observations from 26 September to 1 October 2020, local time.</p>
Full article ">Figure 3
<p>The integral scale <span class="html-italic">L<sub>v</sub></span> (<b>a</b>), relative error of the estimation of TKEDR and <math display="inline"><semantics> <mrow> <msubsup> <mi>C</mi> <mi>n</mi> <mn>2</mn> </msubsup> </mrow> </semantics></math> (<b>b</b>), and comparative statistical analysis of LAS and CDWL observation results (<b>c</b>) from 26 September to 1 October 2020, local time.</p>
Full article ">Figure 4
<p>The CNR (<b>a</b>), horizontal wind speed (<b>b</b>), wind direction (<b>c</b>), vertical wind speed (<b>e</b>), wind shear (<b>f</b>), and log<sub>10</sub>(TKEDR) (<b>g</b>) derived from CDWL. The temperature (<b>d</b>) and pressure (<b>h</b>) retrieved from MWR and the barometric formula in the observations from 6 September to 7 September 2019, local time.</p>
Full article ">Figure 5
<p>The results of the variances of vertical wind fluctuation (<b>a</b>), instantaneous flux of potential temperature (<b>b</b>), log<sub>10</sub><math display="inline"><semantics> <mrow> <msubsup> <mi>C</mi> <mi>n</mi> <mn>2</mn> </msubsup> </mrow> </semantics></math> estimated by <math display="inline"><semantics> <mrow> <msubsup> <mi>C</mi> <mrow> <mi>n</mi> <mo> </mo> <mi>T</mi> <mi>a</mi> <mi>t</mi> <mi>a</mi> <mi>r</mi> <mi>s</mi> <mi>k</mi> <mi>i</mi> <mi>i</mi> <mo>,</mo> <mn>1961</mn> </mrow> <mn>2</mn> </msubsup> </mrow> </semantics></math> (<b>c</b>), <math display="inline"><semantics> <mrow> <msubsup> <mi>C</mi> <mrow> <mi>n</mi> <mo> </mo> <mi>T</mi> <mi>a</mi> <mi>t</mi> <mi>a</mi> <mi>r</mi> <mi>s</mi> <mi>k</mi> <mi>i</mi> <mi>i</mi> <mo>,</mo> <mn>1971</mn> </mrow> <mn>2</mn> </msubsup> </mrow> </semantics></math> (<b>d</b>), and <math display="inline"><semantics> <mrow> <msubsup> <mi>C</mi> <mrow> <mi>n</mi> <mo> </mo> <mi>L</mi> <mi>u</mi> <mi>c</mi> <mi>e</mi> <mo>,</mo> <mn>2020</mn> </mrow> <mn>2</mn> </msubsup> </mrow> </semantics></math> (<b>e</b>) in the observations from 6 September to 7 September 2019, local time.</p>
Full article ">Figure 6
<p>The potential temperature (<b>a</b>), Brunt–Väisälä frequency squared (<b>b</b>), gradient Richardson number (<b>c</b>), log<sub>10</sub><math display="inline"><semantics> <mrow> <msubsup> <mi>C</mi> <mi>n</mi> <mn>2</mn> </msubsup> </mrow> </semantics></math> (<b>d</b>) profiles, perturbations of temperature (<b>e</b>), horizontal wind speed (<b>f</b>), vertical wind speed (<b>g</b>), and averaged wind speed perturbations within 1–2 km (<b>h</b>) derived from CDWL and MWR in the observations on 6–7 September 2019, local time.</p>
Full article ">Figure 7
<p>The results of potential temperature (<b>a</b>,<b>b</b>), potential temperature gradient (<b>e</b>,<b>f</b>), Brunt–Väisälä frequency squared (<b>c</b>,<b>d</b>), and gradient Richardson number (<b>g</b>,<b>h</b>) profiles derived from CDWL, MWR, and the barometric formula at different times on 6–7 September 2019, local time.</p>
Full article ">Figure 8
<p>The relative error of the estimation of TKEDR (<b>a</b>), combined <math display="inline"><semantics> <mrow> <msubsup> <mi>C</mi> <mi>n</mi> <mn>2</mn> </msubsup> </mrow> </semantics></math> (<b>b</b>), mixing layer height (MLH, (<b>a</b>,<b>b</b>)), and integral scale <span class="html-italic">L<sub>v</sub></span> (<b>c</b>) calculated from CDWL and MWR in the observations from 6 September to 7 September 2019, local time.</p>
Full article ">Figure 9
<p>Profiles of log<sub>10</sub><math display="inline"><semantics> <mrow> <msubsup> <mi>C</mi> <mi>n</mi> <mn>2</mn> </msubsup> </mrow> </semantics></math> in raw data that are distance averaged and show a shaded area error bar, calculated from the HAP model and time averaged integral scale <span class="html-italic">L<sub>v</sub></span> (<b>a7</b>–<b>e7</b>) during the period of 14:59:30–15:11:49 (<b>a1</b>–<b>a6</b>) and 20:59:22–21:11:41 (<b>b1</b>–<b>b6</b>) on 6 September and 02:59:02–03:11:21 (<b>c1</b>–<b>c6</b>), 09:01:17–09:13:36 (<b>d1</b>–<b>d6</b>), and 15:01:19–15:13:40 (<b>e1</b>–<b>e6</b>) on 7 September 2019, local time.</p>
Full article ">
21 pages, 19246 KiB  
Article
Damage Properties of the Block-Stone Embankment in the Qinghai–Tibet Highway Using Ground-Penetrating Radar Imagery
by Shunshun Qi, Guoyu Li, Dun Chen, Mingtang Chai, Yu Zhou, Qingsong Du, Yapeng Cao, Liyun Tang and Hailiang Jia
Remote Sens. 2022, 14(12), 2950; https://doi.org/10.3390/rs14122950 - 20 Jun 2022
Cited by 9 | Viewed by 4913
Abstract
The block-stone embankment is a special type of embankment widely used to protect the stability of the underlying warm and ice-rich permafrost. Under the influence of multiple factors, certain damages will still occur in the block-stone embankment after a period of operation, which [...] Read more.
The block-stone embankment is a special type of embankment widely used to protect the stability of the underlying warm and ice-rich permafrost. Under the influence of multiple factors, certain damages will still occur in the block-stone embankment after a period of operation, which may weaken or destroy its cooling function, introducing more serious damages to the Qinghai–Tibet Highway (QTH). Ground-penetrating radar (GPR), a nondestructive testing technique, was adopted to investigate the damage properties of the damaged block-stone embankment. GPR imagery, together with the other data and methods (structural characteristics, field survey data, GPR parameters, etc.), indicated four categories of damage: (i) loosening of the upper sand-gravel layer; (ii) loosening of the block-stone layer; (iii) settlement of the block-stone layer; and (iv) dense filling of the block-stones layer. The first two conditions were widely distributed, whereas the settlement and dense filling of the block-stone layer were less so, and the other combined damages also occurred frequently. The close correlation between the different damages indicated a causal relationship. A preliminary discussion of these observations about the influences on the formation of the damage of the block-stone embankment is included. The findings provide some points of reference for the future construction and maintenance of block-stone embankments in permafrost regions. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>(<b>a</b>) Proactive temperature-controlling measures and the surrounding environment along the QTH; (<b>b</b>) field photographs, from left to right: thermosyphon embankment; block-stone embankment, and duct-ventilated embankment.</p>
Full article ">Figure 2
<p>Study area: (<b>a</b>) spatial distribution of permafrost region in the QTP (the permafrost map is from [<a href="#B61-remotesensing-14-02950" class="html-bibr">61</a>,<a href="#B62-remotesensing-14-02950" class="html-bibr">62</a>]); (<b>b</b>) the location of k3024 to K3025 in the QTH; (<b>c</b>) field photographs of k3024 to k3025, from top to bottom: “wavy” pavement; local settlement of embankment and the block-stone embankment.</p>
Full article ">Figure 3
<p>Schematic diagram of the block-stone embankment: (<b>a</b>) block-stone embankment structure; (<b>b</b>,<b>c</b>) field photographs of block-stone layer.</p>
Full article ">Figure 4
<p>Schematic diagram of the working principle of a block-stone embankment: (<b>a</b>) cold season; (<b>b</b>) warm season.</p>
Full article ">Figure 5
<p>Schematic diagram of the GPR survey on the highway: (<b>a</b>) highway; (<b>b</b>) the electromagnetic waves emitted by the transmitting antenna are reflected by the different layers and then return to the receiving antenna; (<b>c</b>) strong electromagnetic wave reflection; (<b>d</b>) weak electromagnetic wave reflection.</p>
Full article ">Figure 6
<p>Diagram of the field GPR survey on QTH: (<b>a</b>) GPR survey system; (<b>b</b>) site survey.</p>
Full article ">Figure 7
<p>GPR data processing flow: (<b>a</b>) input data; (<b>b</b>) remove direct waves; (<b>c</b>) exponential gain; (<b>d</b>) average path extraction; (<b>e</b>) 1-D filtering; (<b>f</b>) 2-D filtering.</p>
Full article ">Figure 8
<p>GPR images of the loosening of the upper sand-gravel layers: (<b>a</b>) K3024 <math display="inline"> <semantics> <mo>+</mo> </semantics> </math> 935 to K3024 <math display="inline"> <semantics> <mo>+</mo> </semantics> </math> 956; (<b>b</b>) K3024 <math display="inline"> <semantics> <mo>+</mo> </semantics> </math> 544 to K3024 <math display="inline"> <semantics> <mo>+</mo> </semantics> </math> 565; and (<b>c</b>) K3024 <math display="inline"> <semantics> <mo>+</mo> </semantics> </math> 498 to K3024 <math display="inline"> <semantics> <mo>+</mo> </semantics> </math> 548. The yellow dashed box shows the abnormal area of the sand-gravel layer, and the yellow arrow is used to facilitate reading.</p>
Full article ">Figure 9
<p>GPR images indicating looseness of the block-stone layers: (<b>a</b>) K3024 <math display="inline"> <semantics> <mo>+</mo> </semantics> </math> 676 to K3024 <math display="inline"> <semantics> <mo>+</mo> </semantics> </math> 697; (<b>b</b>) K3024 <math display="inline"> <semantics> <mo>+</mo> </semantics> </math> 488 to K3024 <math display="inline"> <semantics> <mo>+</mo> </semantics> </math> 509. The yellow dashed box shows the abnormal area of the block-stone layer, the red dashed box shows the abnormal area of the ground, and the yellow and red arrow is used to facilitate reading.</p>
Full article ">Figure 10
<p>GPR images indicating the settlement of the block-stone layer: (<b>a</b>) K3024 <math display="inline"> <semantics> <mo>+</mo> </semantics> </math> 842 to K3024 <math display="inline"> <semantics> <mo>+</mo> </semantics> </math> 863; (<b>b</b>) K3024 <math display="inline"> <semantics> <mo>+</mo> </semantics> </math> 887 to K3024 <math display="inline"> <semantics> <mo>+</mo> </semantics> </math> 909. The yellow dashed box shows the abnormal area of the sand-gravel layer, the red dotted line indicates the settlement trend of the block-stone layer, and the yellow and red arrow is used to facilitate reading.</p>
Full article ">Figure 11
<p>GPR images of an abnormally dense block-stone layer: (<b>a</b>) K 3024 <math display="inline"> <semantics> <mo>+</mo> </semantics> </math> 235 to K3024 <math display="inline"> <semantics> <mo>+</mo> </semantics> </math> 256; (<b>b</b>) K 3024 <math display="inline"> <semantics> <mo>+</mo> </semantics> </math> 452 to K3024 <math display="inline"> <semantics> <mo>+</mo> </semantics> </math> 474. The yellow dashed box shows the abnormal area of the block-stone layer, the red dotted line indicates the boundary between the sand-gravel layer and the block-stone layer, and the yellow arrow is used to facilitate reading.</p>
Full article ">Figure 12
<p>Comparison of the results of different processing methods for abnormally dense areas: (<b>a</b>) processed according to the normal flow; (<b>b</b>) after direct wave removal only gain processing was performed; (<b>c</b>) only horizontal signals were removed based on the previous step; (<b>d</b>) only moving-average processing was performed based on the previous step. The yellow arrow is used to increase visibility.</p>
Full article ">Figure 13
<p>For the K3024 to K3025 section: (<b>a</b>) distribution of different types of damage in the block-stone embankment. Note: the unit (five meters as a unit) was considered damaged if the unit was damaged more than half, otherwise, it was a normal unit. (<b>b</b>) overall damage situation of the embankment. Note: The degree of embankment damage was based on the number of damaged species on each unit.</p>
Full article ">Figure 14
<p>Relationship of different types of damage in the block-stone embankment: (<b>a</b>) shows the damaged length and rate of the different types of damage; (<b>b</b>) shows the length and proportion of highway sections with different degrees of damage; (<b>c</b>) shows the relationship between the different types of damage and the abnormality of the lower part of the embankment; (<b>d</b>) shows the relationship between the different types of damage and the loosening of the upper sand-gravel layer.</p>
Full article ">Figure 15
<p>Progressive formation of block-stone embankment damage: (<b>a</b>) shows that the action of various influencing factors causes the thawing of the permafrost beneath the embankment; (<b>b</b>) shows the active layer changes that affected the embankment; (<b>c</b>) shows the worsening damage to the embankment; (<b>d</b>) shows that as the damage progressed, the embankment structure was greatly damaged. The red arrows indicate the influence of the external environment on the embankment (such as sun exposure and vehicle rolling), and the other dotted boxes and arrows are used to facilitate reading.</p>
Full article ">
20 pages, 7878 KiB  
Article
On-Orbit Calibration for Spaceborne Line Array Camera and LiDAR
by Xiangpeng Xu, Sheng Zhuge, Banglei Guan, Bin Lin, Shuwei Gan, Xia Yang and Xiaohu Zhang
Remote Sens. 2022, 14(12), 2949; https://doi.org/10.3390/rs14122949 - 20 Jun 2022
Cited by 5 | Viewed by 2348
Abstract
For a multi-mode Earth observation satellite carrying a line array camera and a multi-beam line array LiDAR, the relative installation attitude of the two sensors is of great significance. In this paper, we propose an on-orbit calibration method for the relative installation attitude [...] Read more.
For a multi-mode Earth observation satellite carrying a line array camera and a multi-beam line array LiDAR, the relative installation attitude of the two sensors is of great significance. In this paper, we propose an on-orbit calibration method for the relative installation attitude of the camera and the LiDAR with no need for the calibration field and additional satellite attitude maneuvers. Firstly, the on-orbit joint calibration model of the relative installation attitude of the two sensors is established. However, there may exist a multi-solution problem in the solving of the above model constrained by non-ground control points. Thus, an alternate iterative method by solving the pseudo-absolute attitude matrix of each sensor in turn is proposed. The numerical validation and simulation experiments results show that the relative positioning error of the line array camera and the LiDAR in the horizontal direction of the ground can be limited to 0.8 m after correction by the method in this paper. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Line array camera coordinate system.</p>
Full article ">Figure 2
<p>The geometric meaning of the matrices of the imaging model of the line array camera.</p>
Full article ">Figure 3
<p>LiDAR coordinate system.</p>
Full article ">Figure 4
<p>Earth observation by LiDAR.</p>
Full article ">Figure 5
<p>The geometric meaning of the matrices of the observation model of the LiDAR.</p>
Full article ">Figure 6
<p>Earth observation by line array camera and LiDAR.</p>
Full article ">Figure 7
<p>Flowchart of numerical validation.</p>
Full article ">Figure 8
<p>Flowchart of simulation experiment.</p>
Full article ">Figure 9
<p>Simulated imaging area.</p>
Full article ">Figure 10
<p>Appliances for the hardware-in-loop experiment.</p>
Full article ">Figure 11
<p>Simulated images of the (<b>a</b>) line array and (<b>b</b>) LiDAR.</p>
Full article ">Figure 12
<p>Distribution of the point pairs for calibration in the camera image and LiDAR image.</p>
Full article ">Figure 13
<p>Distribution of the point pairs for the accuracy verification in the camera image and LiDAR image.</p>
Full article ">
18 pages, 7402 KiB  
Article
Spatial and Temporal Variability of Key Bio-Temperature Indicators and Their Effects on Vegetation Dynamics in the Great Lakes Region of Central Asia
by Xuan Gao and Dongsheng Zhao
Remote Sens. 2022, 14(12), 2948; https://doi.org/10.3390/rs14122948 - 20 Jun 2022
Viewed by 2314
Abstract
Dryland ecosystems are fragile to climate change due to harsh environmental conditions. Climate change affects vegetation growth primarily by altering some key bio-temperature thresholds. Key bio-temperatures are closely related to vegetation growth, and slight changes could produce substantial effects on ecosystem structure and [...] Read more.
Dryland ecosystems are fragile to climate change due to harsh environmental conditions. Climate change affects vegetation growth primarily by altering some key bio-temperature thresholds. Key bio-temperatures are closely related to vegetation growth, and slight changes could produce substantial effects on ecosystem structure and function. Therefore, this study selected the number of days with daily mean temperature above 0 °C (DT0), 5 °C (DT5), 10 °C (DT10), 20 °C (DT20), the start of growing season (SGS), the end of growing season (EGS), and the length of growing season (LGS) as bio-temperature indicators to analyze the response of vegetation dynamics to climate change in the Great Lakes Region of Central Asia (GLRCA) for the period 1982–2014. On the regional scale, DT0, DT5, DT10, and DT20 exhibited an overall increasing trend. Spatially, most of the study area showed that the negative correlation between DT0, DT5, DT10, DT20 with the annual Normalized Difference Vegetation Index (NDVI) increased with increasing bio-temperature thresholds. In particular, more than 88.3% of the study area showed a negative correlation between annual NDVI and DT20, as increased DT20 exacerbated ecosystem drought. Moreover, SGS exhibited a significantly advanced trend at a rate of −0.261 days/year for the regional scale, while EGS experienced a significantly delayed trend at a rate of 0.164 days/year. Because of changes in SGS and EGS, LGS across the GLRCA was extended at a rate of 0.425 days/year, which was mainly attributed to advanced SGS. In addition, our study revealed that about 53.6% of the study area showed a negative correlation between annual NDVI and LGS, especially in the north, indicating a negative effect of climate warming on vegetation growth in the drylands. Overall, the results of this study will help predict the response of vegetation to future climate change in the GLRCA, and support decision-making for implementing effective ecosystem management in arid and semi-arid regions. Full article
(This article belongs to the Topic Climate Change and Environmental Sustainability)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>(<b>a</b>) The land cover types (Source: MCD12Q1, <a href="https://modis.gsfc.nasa.gov/data/" target="_blank">https://modis.gsfc.nasa.gov/data/</a>, accessed on 20 September 2021), and (<b>b</b>) the climate classification in the study area.</p>
Full article ">Figure 2
<p>Spatial distribution of the annual means for (<b>a</b>) GLDAS temperature, (<b>b</b>) CRU temperature over the GLRCA from 1982–2014.</p>
Full article ">Figure 3
<p>Spatial distribution of trends in (<b>a</b>) annual mean temperature, (<b>b</b>) annual NDVI, and (<b>c</b>) correlation between annual mean temperature and annual NDVI over the GLRCA from 1982–2014.</p>
Full article ">Figure 4
<p>Spatial distribution of (<b>a</b>) the annual means for maximum temperature, (<b>b</b>) the trend for maximum temperature, and (<b>c</b>) the uncertainty in the trend of maximum temperature in the study area from 1982–2014.</p>
Full article ">Figure 5
<p>Interannual variations in (<b>a</b>) DT<sub>0</sub>, (<b>b</b>) DT<sub>5</sub>, (<b>c</b>) DT<sub>10</sub>, and (<b>d</b>) DT<sub>20</sub> at a regional scale in the GLRCA during the period 1982–2014.</p>
Full article ">Figure 6
<p>Spatial distribution of trends in (<b>a</b>) DT<sub>0</sub>, (<b>b</b>) DT<sub>5</sub>, (<b>c</b>) DT<sub>10</sub>, and (<b>d</b>) DT<sub>20</sub> during the period 1982–2014.</p>
Full article ">Figure 7
<p>Spatial distribution of uncertainty in the trends of (<b>a</b>) DT<sub>0</sub>, (<b>b</b>) DT<sub>5</sub>, (<b>c</b>) DT<sub>10</sub>, and (<b>d</b>) DT<sub>20</sub> during the period 1982–2014.</p>
Full article ">Figure 8
<p>Spatial distribution of correlations between annual NDVI and (<b>a</b>) DT<sub>0</sub>, (<b>b</b>) DT<sub>5</sub>, (<b>c</b>) DT<sub>10</sub>, and (<b>d</b>) DT<sub>20</sub> for the period 1982–2014.</p>
Full article ">Figure 9
<p>Interannual variations in (<b>a</b>) SGS, (<b>b</b>) EGS, and (<b>c</b>) LGS at a regional scale in the GLRCA during the period 1982–2014.</p>
Full article ">Figure 10
<p>Spatial distribution of trends in (<b>a</b>) SGS, (<b>b</b>) EGS, and (<b>c</b>) LGS for the period 1982–2014.</p>
Full article ">Figure 11
<p>Spatial distribution of uncertainty in the trends of (<b>a</b>) SGS, (<b>b</b>) EGS, and (<b>c</b>) LGS during the period 1982–2014.</p>
Full article ">Figure 12
<p>Spatial distribution of correlations between annual NDVI and (<b>a</b>) SGS, (<b>b</b>) EGS, and (<b>c</b>) LGS from 1982–2014.</p>
Full article ">
15 pages, 15759 KiB  
Article
Modeling Potential Impacts on Regional Climate Due to Land Surface Changes across Mongolia Plateau
by Guangshuai Li, Lingxue Yu, Tingxiang Liu, Yue Jiao and Jiaxin Yu
Remote Sens. 2022, 14(12), 2947; https://doi.org/10.3390/rs14122947 - 20 Jun 2022
Cited by 9 | Viewed by 2328
Abstract
Although desertification has greatly increased across the Mongolian Plateau during the last decades of the 20th century, recent satellite records documented increasing vegetation growth since the 21st century in some areas of the Mongolian Plateau. Compared to the study of desertification, the opposite [...] Read more.
Although desertification has greatly increased across the Mongolian Plateau during the last decades of the 20th century, recent satellite records documented increasing vegetation growth since the 21st century in some areas of the Mongolian Plateau. Compared to the study of desertification, the opposite characteristics of land use and vegetation cover changes and their different effects on regional land–atmosphere interaction factors still lack enough attention across this vulnerable region. Using long-term time-series multi-source satellite records and regional climate model, this study investigated the climate feedback to the observed land surface changes from the 1990s to the 2010s in the Mongolia Plateau. Model simulation suggests that vegetation greening induced a local cooling effect, while the warming effect is mainly located in the vegetation degradation area. For the typical vegetation greening area in the southeast of Inner Mongolia, latent heat flux increased over 2 W/m2 along with the decrease of sensible heat flux over 2 W/m2, resulting in a total evapotranspiration increase by 0.1~0.2 mm/d and soil moisture decreased by 0.01~0.03 mm/d. For the typical vegetation degradation area in the east of Mongolia and mid-east of Inner Mongolia, the latent heat flux decreased over 2 W/m2 along with the increase of sensible heat flux over 2 W/m2 obviously, while changes in moisture cycling were spatially more associated with variations of precipitation. It means that precipitation still plays an important role in soil moisture for most areas, and some areas would be at potential risk of drought with the asynchronous increase of evapotranspiration and precipitation. Full article
(This article belongs to the Special Issue Remote Sensing for Advancing Nature-Based Climate Solutions)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Location of study area, simulation domain (D01), and the land use pattern.</p>
Full article ">Figure 2
<p>Changes in growing season LAI from the 1990s to the 2010s.</p>
Full article ">Figure 3
<p>Relationship between simulated T2 and observed T2 (<b>a</b>), simulated precipitation and observed precipitation (<b>b</b>), from 2011 to 2018.</p>
Full article ">Figure 4
<p>Model-simulated differences in T2 (<b>a</b>) and TSK (<b>b</b>) due to vegetation changes in the growing season between scenarios of SL1990s and SL2010s.</p>
Full article ">Figure 5
<p>Model-simulated differences in the growing season net radiation (<b>a</b>), latent heat flux (<b>b</b>), and sensible heat flux (<b>c</b>) due to vegetation changes between scenarios of SL1990s and SL2010s.</p>
Full article ">Figure 6
<p>Model-simulated differences in ECAN (<b>a</b>), ESOIL (<b>b</b>), and ETRAN (<b>c</b>) due to vegetation changes in the growing season between scenarios of SL1990s and SL2010s.</p>
Full article ">Figure 7
<p>Model-simulated differences in precipitation (<b>a</b>), and soil moisture (<b>b</b>) due to vegetation changes in the growing season between scenarios of SL1990s and SL2010s.</p>
Full article ">
22 pages, 7365 KiB  
Article
A Dual-Generator Translation Network Fusing Texture and Structure Features for SAR and Optical Image Matching
by Han Nie, Zhitao Fu, Bo-Hui Tang, Ziqian Li, Sijing Chen and Leiguang Wang
Remote Sens. 2022, 14(12), 2946; https://doi.org/10.3390/rs14122946 - 20 Jun 2022
Cited by 8 | Viewed by 3549
Abstract
The matching problem for heterologous remote sensing images can be simplified to the matching problem for pseudo homologous remote sensing images via image translation to improve the matching performance. Among such applications, the translation of synthetic aperture radar (SAR) and optical images is [...] Read more.
The matching problem for heterologous remote sensing images can be simplified to the matching problem for pseudo homologous remote sensing images via image translation to improve the matching performance. Among such applications, the translation of synthetic aperture radar (SAR) and optical images is the current focus of research. However, the existing methods for SAR-to-optical translation have two main drawbacks. First, single generators usually sacrifice either structure or texture features to balance the model performance and complexity, which often results in textural or structural distortion; second, due to large nonlinear radiation distortions (NRDs) in SAR images, there are still visual differences between the pseudo-optical images generated by current generative adversarial networks (GANs) and real optical images. Therefore, we propose a dual-generator translation network for fusing structure and texture features. On the one hand, the proposed network has dual generators, a texture generator, and a structure generator, with good cross-coupling to obtain high-accuracy structure and texture features; on the other hand, frequency-domain and spatial-domain loss functions are introduced to reduce the differences between pseudo-optical images and real optical images. Extensive quantitative and qualitative experiments show that our method achieves state-of-the-art performance on publicly available optical and SAR datasets. Our method improves the peak signal-to-noise ratio (PSNR) by 21.0%, the chromatic feature similarity (FSIMc) by 6.9%, and the structural similarity (SSIM) by 161.7% in terms of the average metric values on all test images compared with the next best results. In addition, we present a before-and-after translation comparison experiment to show that our method improves the average keypoint repeatability by approximately 111.7% and the matching accuracy by approximately 5.25%. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>High-quality image translation results obtained with our method. The pseudo-optical image is the image generated from the SAR image through our method.</p>
Full article ">Figure 1 Cont.
<p>High-quality image translation results obtained with our method. The pseudo-optical image is the image generated from the SAR image through our method.</p>
Full article ">Figure 2
<p>The generators and discriminator of our network. <b>Generators</b>: The SAR-to-optical translation process is divided between two generators, i.e., a structure generator and a texture generator, which borrow each other’s depth features, and the Bi-GFF and CFA modules are used to refine and fuse the features from these structure and texture reconstruction branches to form the final pseudo-optical image. <b>Discriminator</b>: The texture branch guides texture generation, and the structure branch guides structure generation.</p>
Full article ">Figure 3
<p>Structural diagram of the Bi-GFF module, in which deep fusion of texture and structure features is performed.</p>
Full article ">Figure 4
<p>Structural diagram of the CFA module, which effectively models long-term spatial dependence through multiscale information.</p>
Full article ">Figure 5
<p>Some examples of training data. First row: SAR image blocks. Second row: optical image blocks.</p>
Full article ">Figure 6
<p>Some of the image blocks were selected as test samples. First row: SAR image blocks. Second row: optical image blocks.</p>
Full article ">Figure 7
<p>A visual comparison of different SAR-to-optical methods. The size of all images is 256 × 256.</p>
Full article ">Figure 7 Cont.
<p>A visual comparison of different SAR-to-optical methods. The size of all images is 256 × 256.</p>
Full article ">Figure 8
<p>A visual comparison of the results from the ablation experiment. (1) SEN-1 SAR image. (2) Guo [<a href="#B38-remotesensing-14-02946" class="html-bibr">38</a>]. Adapted with permission from Ref. [<a href="#B38-remotesensing-14-02946" class="html-bibr">38</a>]. 2021, Xiefan Guo. (3) Ours (+<span class="html-italic">MSE</span> Loss). (4) Ours (+<span class="html-italic">MSE</span> loss +<span class="html-italic">FFL</span> Loss). (5) SEN-2 optical image.</p>
Full article ">Figure 8 Cont.
<p>A visual comparison of the results from the ablation experiment. (1) SEN-1 SAR image. (2) Guo [<a href="#B38-remotesensing-14-02946" class="html-bibr">38</a>]. Adapted with permission from Ref. [<a href="#B38-remotesensing-14-02946" class="html-bibr">38</a>]. 2021, Xiefan Guo. (3) Ours (+<span class="html-italic">MSE</span> Loss). (4) Ours (+<span class="html-italic">MSE</span> loss +<span class="html-italic">FFL</span> Loss). (5) SEN-2 optical image.</p>
Full article ">Figure 9
<p>A visual comparison of the results of the matching experiment using RIFT and LoFTR.</p>
Full article ">Figure 10
<p>Test 1: qualitative similarity comparison of the structural features of the real optical image and the pseudo-optical images obtained using different SAR-to-optical image translation methods.</p>
Full article ">
21 pages, 11037 KiB  
Article
Positioning of Quadruped Robot Based on Tightly Coupled LiDAR Vision Inertial Odometer
by Fangzheng Gao, Wenjun Tang, Jiacai Huang and Haiyang Chen
Remote Sens. 2022, 14(12), 2945; https://doi.org/10.3390/rs14122945 - 20 Jun 2022
Cited by 8 | Viewed by 3238
Abstract
Quadruped robots, an important class of unmanned aerial vehicles, have broad potential for applications in education, service, industry, military, and other fields. Their independent positioning plays a key role for completing assigned tasks in a complex environment. However, positioning based on global navigation [...] Read more.
Quadruped robots, an important class of unmanned aerial vehicles, have broad potential for applications in education, service, industry, military, and other fields. Their independent positioning plays a key role for completing assigned tasks in a complex environment. However, positioning based on global navigation satellite systems (GNSS) may result in GNSS jamming and quadruped robots not operating properly in environments sheltered by buildings. In this paper, a tightly coupled LiDAR vision inertial odometer (LVIO) is proposed to address the positioning inaccuracy of quadruped robots, which have poor mileage information obtained though legs and feet structures only. With this optimization method, the point cloud data obtained by 3D LiDAR, the image feature information obtained by binocular vision, and the IMU inertial data are combined to improve the precise indoor and outdoor positioning of a quadruped robot. This method reduces the errors caused by the uniform motion model in laser odometer as well as the image blur caused by rapid movements of the robot, which can lead to error-matching in a dynamic scene; at the same time, it alleviates the impact of drift on inertial measurements. Finally, the quadruped robot in the laboratory is used to build a physical platform for verification. The experimental results show that the designed LVIO effectively realizes the positioning of four groups of robots with high precision and strong robustness, both indoors or outdoors, which verify the feasibility and effectiveness of the proposed method. Full article
(This article belongs to the Special Issue UAV Positioning: From Ground to Sky)
Show Figures

Figure 1

Figure 1
<p>The moving track of quadruped robot.</p>
Full article ">Figure 2
<p>Feature tracking by optical flow method for binocular camera (where the green points represent the characteristic points of tracking, the red points represent the low tracking times, and the blue points represent the high tracking times).</p>
Full article ">Figure 3
<p>Camera-IMU model.</p>
Full article ">Figure 4
<p>Motion distortion compensation and timestamp synchronization of laser frame point cloud.</p>
Full article ">Figure 5
<p>Point-to-plane ICP algorithm.</p>
Full article ">Figure 6
<p>Point-to-plane ICP algorithm.</p>
Full article ">Figure 7
<p>Schematic diagram of sliding window marginalization. <span class="html-fig-inline" id="remotesensing-14-02945-i001"> <img alt="Remotesensing 14 02945 i001" src="/remotesensing/remotesensing-14-02945/article_deploy/html/images/remotesensing-14-02945-i001.png"/></span> represents the old key frame, <span class="html-fig-inline" id="remotesensing-14-02945-i002"> <img alt="Remotesensing 14 02945 i002" src="/remotesensing/remotesensing-14-02945/article_deploy/html/images/remotesensing-14-02945-i002.png"/></span> represents the latest frame, <span class="html-fig-inline" id="remotesensing-14-02945-i003"> <img alt="Remotesensing 14 02945 i003" src="/remotesensing/remotesensing-14-02945/article_deploy/html/images/remotesensing-14-02945-i003.png"/></span> represents IMU constraint, <span class="html-fig-inline" id="remotesensing-14-02945-i004"> <img alt="Remotesensing 14 02945 i004" src="/remotesensing/remotesensing-14-02945/article_deploy/html/images/remotesensing-14-02945-i004.png"/></span> represents visual feature, <span class="html-fig-inline" id="remotesensing-14-02945-i005"> <img alt="Remotesensing 14 02945 i005" src="/remotesensing/remotesensing-14-02945/article_deploy/html/images/remotesensing-14-02945-i005.png"/></span> represents fixed state, and <span class="html-fig-inline" id="remotesensing-14-02945-i006"> <img alt="Remotesensing 14 02945 i006" src="/remotesensing/remotesensing-14-02945/article_deploy/html/images/remotesensing-14-02945-i006.png"/></span> represents estimated state.</p>
Full article ">Figure 8
<p>Loopback detection and feature matching between loopback candidate frames and tightly coupled relocation.</p>
Full article ">Figure 9
<p>Feature matching between loopback candidate frames.</p>
Full article ">Figure 10
<p>LiDAR vision inertial odometer algorithm flow chart.</p>
Full article ">Figure 11
<p>Experimental platform, (<b>a</b>) VLP-16 LiDAR, (<b>b</b>) quadruped robot, and (<b>c</b>) RealSense D455 camera.</p>
Full article ">Figure 12
<p>The indoor and outdoor experimental environment: (<b>a</b>,<b>b</b>) the indoor experimental environment and (<b>c</b>,<b>d</b>) the outdoor experimental environment.</p>
Full article ">Figure 13
<p>Joint calibration process of binocular camera and IMU.</p>
Full article ">Figure 14
<p>Calibration reprojection error of binocular camera. (<b>a</b>) Reprojection error of the left-eye camera and (<b>b</b>) reprojection error of the right-eye camera.</p>
Full article ">Figure 15
<p>The LiDAR–camera joint calibration process: (<b>a</b>) the real environment of the calibration and (<b>b</b>) the LiDAR point cloud of the calibration.</p>
Full article ">Figure 16
<p>Odometer operating trajectory in indoor environment.</p>
Full article ">Figure 17
<p>The track error visualizations of the LVIO algorithm and LVI-SAM algorithm: (<b>a</b>) the trajectory error of the LVIO algorithm and (<b>b</b>) the trajectory error of LVI-SAM algorithm.</p>
Full article ">Figure 18
<p>Specific error curve: (<b>a</b>) mean error, median error, RMSE error, and STD error of the proposed method; (<b>b</b>) mean error, median error, RMSE error, and STD error of LVI-SAM algorithm; (<b>c</b>) the error trajectory in the x, y, and z directions; and (<b>d</b>) the error trajectory in pitch, roll, and yaw direction angle.</p>
Full article ">Figure 19
<p>The lack of texture in an interior white wall.</p>
Full article ">Figure 20
<p>LVIO loop experiment running track in indoor environment: (<b>a</b>) before adding loopback detection and (<b>b</b>) after adding loopback detection.</p>
Full article ">Figure 21
<p>Outdoor LVIO operation trajectory (the red line).</p>
Full article ">
14 pages, 859 KiB  
Technical Note
Ant Colony Pheromone Mechanism-Based Passive Localization Using UAV Swarm
by Yongkun Zhou, Dan Song, Bowen Ding, Bin Rao, Man Su and Wei Wang
Remote Sens. 2022, 14(12), 2944; https://doi.org/10.3390/rs14122944 - 20 Jun 2022
Cited by 6 | Viewed by 2144
Abstract
The problem of passive localization using an unmanned aerial vehicle (UAV) swarm is studied. For multi-UAV localization systems with limited communication and observation range, the challenge is how to obtain accurate target state consistency estimates through local UAV communication. In this paper, an [...] Read more.
The problem of passive localization using an unmanned aerial vehicle (UAV) swarm is studied. For multi-UAV localization systems with limited communication and observation range, the challenge is how to obtain accurate target state consistency estimates through local UAV communication. In this paper, an ant colony pheromone mechanism-based passive localization method using a UAV swarm is proposed. Different from traditional distributed fusion localization algorithms, the proposed method makes use of local interactions among individuals to process the observation data with UAVs, which greatly reduces the cost of the system. First, the UAVs that have detected the radiation source target estimate the rough target position based on the pseudo-linear estimation (PLE). Then, the ant colony pheromone mechanism is introduced to further improve localization accuracy. The ant colony pheromone mechanism consists of two stages: pheromone injection and pheromone transmission. In the pheromone injection mechanism, each UAV uses the maximum likelihood (ML) algorithm with the current observed target bearing information to correct the initial target position estimate. Then, the UAV swarm weights and fuses the target position information between individuals based on the pheromone transmission mechanism. Numerical results demonstrate that the accuracy of the proposed method is better than that of traditional localization algorithms and close to the Cramer–Rao lower bound (CRLB) for small measurement noise. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>UAV swarm passive localization geometry.</p>
Full article ">Figure 2
<p>Initial distribution of UAV swarm passive localization.</p>
Full article ">Figure 3
<p>Distribution of UAV swarm passive localization after pheromone transmission.</p>
Full article ">Figure 4
<p>The flowchart of the proposed algorithm.</p>
Full article ">Figure 5
<p>The simulation scenario of a single-fixed target.</p>
Full article ">Figure 6
<p>RMSE of position estimation versus bearing standard deviation.</p>
Full article ">Figure 7
<p>Bias of position estimation versus bearing standard deviation.</p>
Full article ">Figure 8
<p>RMSE of position estimation versus bearing standard deviation.</p>
Full article ">Figure 9
<p>Bias of position estimation versus bearing standard deviation.</p>
Full article ">Figure 10
<p>RMSE of position estimation versus communication radius.</p>
Full article ">Figure 11
<p>Bias of position estimation versus communication radius.</p>
Full article ">Figure 12
<p>The simulation scenario of a single-unfixed target.</p>
Full article ">Figure 13
<p>RMSE of position estimation versus target distance.</p>
Full article ">Figure 14
<p>Bias of position estimation versus target distance.</p>
Full article ">
16 pages, 8460 KiB  
Article
Tight Integration of GNSS and Static Level for High Accuracy Dilapidated House Deformation Monitoring
by Jian Yang, Weiming Tang, Wei Xuan and Ruijie Xi
Remote Sens. 2022, 14(12), 2943; https://doi.org/10.3390/rs14122943 - 20 Jun 2022
Cited by 1 | Viewed by 1913
Abstract
Global Navigation Satellite System (GNSS) can provide high-precision three-dimensional real-time or quasi-real-time changes of monitoring points automatically in house monitoring applications. However, due to the signal sheltering problem, large observation noise and multipath effects in urban observing environment with dense buildings, ambiguity resolution [...] Read more.
Global Navigation Satellite System (GNSS) can provide high-precision three-dimensional real-time or quasi-real-time changes of monitoring points automatically in house monitoring applications. However, due to the signal sheltering problem, large observation noise and multipath effects in urban observing environment with dense buildings, ambiguity resolution would be hard, and GNSS accuracy cannot always achieve millimeter level to satisfy the requirement of house monitoring. Static level is a precision instrument for measuring elevation difference and its variations, with a precision up to sub-millimeter level. It could be integrated with GNSS to improve the positioning accuracy in height direction. However, the existing integration of GNSS and static level is mostly on a respective results level. In this study, we proposed a method of integrating GNSS and static level observations tightly to enhance the GNSS positioning performance. The hardware design and integration mathematic model in data processing were introduced, and a group of experiments were carried out to verify the performance in positioning with and without the static level observation constraints. It found that the vertical monitoring measurement results of static level can achieve less than 1 mm. The GNSS ambiguity resolution performance can be improved by incorporating the measurement of static level into GNSS positioning equation as external constraints, and the precision of GNSS float solutions was significantly improved. Finally, the static level constraint can further improve the accuracy of the fixed solution from about 2 cm to better than 2 mm in vertical direction, which is even better than the accuracy in horizontal directions with about 3–6 mm with the static level constraint. The tight combination data processing algorithm can significantly improve the working efficiency, accuracy, and reliability of the application of dangerous house monitoring. Full article
(This article belongs to the Section Earth Observation Data)
Show Figures

Figure 1

Figure 1
<p>The principle of the static level. The original stage (<b>a</b>) and the new statement when the vertical displacement occurs (<b>b</b>).</p>
Full article ">Figure 2
<p>The mechanical design drawings of the GNSS antenna and static level integration equipment.</p>
Full article ">Figure 3
<p>The appearance diagram of the GNSS antenna and static level integration equipment.</p>
Full article ">Figure 4
<p>System structure drawing.</p>
Full article ">Figure 5
<p>The left-hand coordinate system.</p>
Full article ">Figure 6
<p>(<b>a</b>) Observing environment; (<b>b</b>) site settings; (<b>c</b>) height adjustment device.</p>
Full article ">Figure 7
<p>Height distance variations of B (<b>a</b>) and C (<b>b</b>) toward A provided by the static level.</p>
Full article ">Figure 8
<p>Interactive interface of the visual quality inspection software.</p>
Full article ">Figure 9
<p>SNR time series of a data case.</p>
Full article ">Figure 10
<p>Skyplot of site A.</p>
Full article ">Figure 11
<p>Flow chart of data processing.</p>
Full article ">Figure 12
<p>Ratio time series of baseline A–B (<b>a</b>) and A–C (<b>b</b>) with and without the static level constraints.</p>
Full article ">Figure 13
<p>Widelane ambiguity float solutions of baseline A–B without (<b>a</b>) and with (<b>b</b>) static levelling constraint.</p>
Full article ">Figure 14
<p>L1/B1 ambiguity float solutions of baseline A–B without (<b>a</b>) and with (<b>b</b>) static levelling constraint.</p>
Full article ">Figure 15
<p>Widelane ambiguity float solutions of baseline A–B without (<b>a</b>) and with (<b>b</b>) static levelling constraint.</p>
Full article ">Figure 16
<p>L1/B1 ambiguity float solutions of baseline A–C without (<b>a</b>) and with (<b>b</b>) static levelling constraint.</p>
Full article ">
29 pages, 31768 KiB  
Article
Risk Assessment of Debris Flow in a Mountain-Basin Area, Western China
by Yanyan Zhou, Dongxia Yue, Geng Liang, Shuangying Li, Yan Zhao, Zengzu Chao and Xingmin Meng
Remote Sens. 2022, 14(12), 2942; https://doi.org/10.3390/rs14122942 - 20 Jun 2022
Cited by 17 | Viewed by 3550
Abstract
Debris flow risk comprehensively reflects the natural and social properties of debris flow disasters and is composed of the risk of the disaster-causing body and the vulnerability of the carrier. The Bailong River Basin (BRB) is a typical mountainous environment where regional debris [...] Read more.
Debris flow risk comprehensively reflects the natural and social properties of debris flow disasters and is composed of the risk of the disaster-causing body and the vulnerability of the carrier. The Bailong River Basin (BRB) is a typical mountainous environment where regional debris flow disasters occur frequently, seriously threatening the lives of residents, infrastructure, and regional ecological security. However, there are few studies on the risk assessment of mountainous debris flow disasters in the BRB. By considering a complete catchment, based on remote sensing and GIS methods, we selected 17 influencing factors, such as area, average slope, lithology, NPP, average annual precipitation, landslide density, river density, fault density, etc. and applied a machine learning algorithm to establish a hazard assessment model. The analysis shows that the Extra Trees model is the most effective for debris flow hazard assessments, with an accuracy rate of 88%. Based on socio-economic data and debris flow disaster survey data, we established a vulnerability assessment model by applying the Contributing Weight Superposition method. We used the product of debris flow hazard and vulnerability to construct a debris flow risk assessment model. The catchments at a very high-risk were distributed mainly in the urban area of Wudu District and the northern part of Tanchang County, that is, areas with relatively dense economic activities and a high disaster frequency. These findings indicate that the assessment results provide scientific support for planning measures to prevent or reduce debris flow hazards. The proposed assessment methods can also be used to provide relevant guidance for a regional risk assessment of debris flows in the BRB and other regions. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>The location of Bailong River Basin and boundary extraction results of 1986 catchments.</p>
Full article ">Figure 2
<p>Technical flow chart of this study.</p>
Full article ">Figure 3
<p>Distribution of factors related to geomorphological conditions.</p>
Full article ">Figure 4
<p>Distribution of factors related to geological structure.</p>
Full article ">Figure 4 Cont.
<p>Distribution of factors related to geological structure.</p>
Full article ">Figure 5
<p>Distribution of factors related to vegetation and soil.</p>
Full article ">Figure 5 Cont.
<p>Distribution of factors related to vegetation and soil.</p>
Full article ">Figure 6
<p>Distribution of factors related to rain conditions.</p>
Full article ">Figure 7
<p>Spatial distribution of factors characterising exposure.</p>
Full article ">Figure 8
<p>Spatial distribution of factors characterising coping capacity.</p>
Full article ">Figure 9
<p>Spatial distribution of factors characterising resilience.</p>
Full article ">Figure 10
<p>Ranking of Model Accuracy Scores.</p>
Full article ">Figure 11
<p>Ten-fold cross-validation ROC curve for the Extra Trees classifier model.</p>
Full article ">Figure 12
<p>Spatial distribution of debris flow hazards in the BRB.</p>
Full article ">Figure 13
<p>RDA ordination diagram of risk and environmental variables in the BRB.</p>
Full article ">Figure 14
<p>Spatial distribution of debris flow vulnerability in the BRB.</p>
Full article ">Figure 15
<p>Spatial distribution of debris flow hazard and landscape photos in the BRB.</p>
Full article ">Figure 16
<p>Ranking of the importance of influencing factors.</p>
Full article ">
14 pages, 5736 KiB  
Technical Note
Improving WRF-Fire Wildfire Simulation Accuracy Using SAR and Time Series of Satellite-Based Vegetation Indices
by Yaron Michael, Gilad Kozokaro, Steve Brenner and Itamar M. Lensky
Remote Sens. 2022, 14(12), 2941; https://doi.org/10.3390/rs14122941 - 20 Jun 2022
Cited by 4 | Viewed by 3426
Abstract
Wildfire simulations depend on fuel representation. Present fuel models are mainly based on the density and properties of different vegetation types. This study aims to improve the accuracy of WRF-Fire wildfire simulations, by using synthetic-aperture radar (SAR) data to estimate the fuel load [...] Read more.
Wildfire simulations depend on fuel representation. Present fuel models are mainly based on the density and properties of different vegetation types. This study aims to improve the accuracy of WRF-Fire wildfire simulations, by using synthetic-aperture radar (SAR) data to estimate the fuel load and the trend of vegetation index to estimate the dryness of woody vegetation. We updated the chaparral and timber standard woody fuel classes in the WRF-Fire fuel settings. We used the ESA global above-ground biomass (AGB) based on SAR data to estimate the fuel load, and the Landsat normalized difference vegetation index (NDVI) trends of woody vegetation to estimate the fuel moisture content. These fuel sub-parameters represent the dynamic changes and spatial variability of woody fuel. We simulated two wildfires in Israel while using three different fuel models: the original 13 Anderson Fire Behavior fuel model, and two modified fuel models introducing AGB alone, and AGB and dryness. The updated fuel model (the basic fuel model plus the AGB and dryness) improved the simulation results significantly, i.e., the Jaccard similarity coefficient increased by 283% on average. Our results demonstrate the potential of combining satellite SAR data and Landsat NDVI trends to improve WRF-Fire wildfire simulations. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Flowchart summarizing the steps of this study.</p>
Full article ">Figure 2
<p>WRF-Fire model workflow: inputs, two-way interactions, and outputs.</p>
Full article ">Figure 3
<p>(<b>a</b>) The multi-scale WRF setup in this study (red boxes) D01 and D02 with horizontal grid spacings of 6000 m and 2000 m, respectively. The sub-grid resolution of the fire sub-model is of 200 m. (<b>b</b>) Locations of the two wildfires in central Israel (case A and B) are marked by the blue stars, and the locations of the meteorological stations by the green circles.</p>
Full article ">Figure 4
<p>The two case studies: (<b>a</b>–<b>d</b>) Bet Shemesh (case A in <a href="#remotesensing-14-02941-f002" class="html-fig">Figure 2</a>), and (<b>e</b>–<b>h</b>) Modiin (case B in <a href="#remotesensing-14-02941-f002" class="html-fig">Figure 2</a>). (<b>a</b>,<b>e</b>) Maps of the actual (black line) and predicted wildfire perimeters (light blue, dark blue, and yellow lines, for fuel maps I, II, and III in <a href="#remotesensing-14-02941-t002" class="html-table">Table 2</a>); (<b>b</b>,<b>f</b>) topography at the two case studies; (<b>c</b>,<b>g</b>) land-cover fuel maps: the dots in the pixels represent the negative trend of woody vegetation, the letters represent low (L), medium (M), and high (H) levels of AGB for chaparral (L = 120, m = 154, H = 180 t ha<sup>−1</sup>) and timber (L = 110, M = 120, G = 130 t ha<sup>−1</sup>); (<b>d</b>,<b>h</b>) 2 m air temperature and wind vectors in the zone of the wildfire at 12:00 using the WRF forecasting model. The star represents the location of ignition. (Additional details of the WRF-Fire results from the three fuel map simulations for three observed times are shown in <a href="#remotesensing-14-02941-f0A1" class="html-fig">Figure A1</a> in <a href="#app1-remotesensing-14-02941" class="html-app">Appendix A</a>).</p>
Full article ">Figure A1
<p>The WRF-Fire simulation results for the three fuel maps are shown for 9:00 (<b>a</b>,<b>d</b>), 10:00 (<b>b</b>,<b>e</b>), and 11:00 AM (<b>c</b>,<b>f</b>) local time.</p>
Full article ">
21 pages, 1819 KiB  
Review
Remote Sensing of Surface and Subsurface Soil Organic Carbon in Tidal Wetlands: A Review and Ideas for Future Research
by Rajneesh Sharma, Deepak R. Mishra, Matthew R. Levi and Lori A. Sutter
Remote Sens. 2022, 14(12), 2940; https://doi.org/10.3390/rs14122940 - 20 Jun 2022
Cited by 13 | Viewed by 5156
Abstract
Tidal wetlands, widely considered the most extensive reservoir of soil organic carbon (SOC), can benefit from remote sensing studies enabling spatiotemporal estimation and mapping of SOC stock. We found that a majority of the remote-sensing-based SOC mapping efforts have been focused on upland [...] Read more.
Tidal wetlands, widely considered the most extensive reservoir of soil organic carbon (SOC), can benefit from remote sensing studies enabling spatiotemporal estimation and mapping of SOC stock. We found that a majority of the remote-sensing-based SOC mapping efforts have been focused on upland ecosystems, not on tidal wetlands. We present a comprehensive review detailing the types of remote sensing models and methods used, standard input variables, results, and limitations for the handful of studies on tidal wetland SOC. Based on that synthesis, we pose several unexplored research questions and methods that are critical for moving tidal wetland SOC science forward. Among these, the applicability of machine learning and deep learning models for predicting surface SOC and the modeling requirements for SOC in subsurface soils (soils without a remote sensing signal, i.e., a soil depth greater than 5 cm) are the most important. We did not find any remote sensing study aimed at modeling subsurface SOC in tidal wetlands. Since tidal wetlands store a significant amount of SOC at greater depths, we hypothesized that surface SOC could be an important covariable along with other biophysical and climate variables for predicting subsurface SOC. Preliminary results using field data from tidal wetlands in the southeastern United States and machine learning model output from mangrove ecosystems in India revealed a strong nonlinear but significant relationship (r2 = 0.68 and 0.20, respectively, p < 2.2 × 10−16 for both) between surface and subsurface SOC at different depths. We investigated the applicability of the Soil Survey Geographic Database (SSURGO) for tidal wetlands by comparing the data with SOC data from the Smithsonian’s Coastal Blue Carbon Network collected during the same decade and found that the SSURGO data consistently over-reported SOC stock in tidal wetlands. We concluded that a novel machine learning framework that utilizes remote sensing data and derived products, the standard covariables reported in the limited literature, and more importantly, other new and potentially informative covariables specific to tidal wetlands such as tidal inundation frequency and height, vegetation species, and soil algal biomass could improve remote-sensing-based tidal wetland SOC studies. Full article
(This article belongs to the Special Issue Remote Sensing for Soil Organic Carbon Mapping and Monitoring)
Show Figures

Figure 1

Figure 1
<p>Soil reflectance spectrum using a hyperspectral sensor in tidal wetlands at Sapelo Island, GA. The sensor was a handheld sensor, and the reading was taken with the proximity (approx. 10–15 cm from the soil surface) of soil (Bohicket soil series) under field conditions.</p>
Full article ">Figure 2
<p>Study locations starting top right and moving clockwise: (<b>a</b>) the state of Odisha, India; (<b>b</b>) Bhitarkanika National Park in Odisha, India; (<b>c</b>) distribution of SOC in upper 5 cm of the soil surface in mangrove wetlands of Bhitarkanika National Park, Odisha, India.</p>
Full article ">Figure 3
<p>Distribution of 261 point locations of SOC measurements in the southeastern coastal United States. These points, collected between 2007 and 2018, were retrieved from the Coastal Blue Carbon Network (CBCN) database.</p>
Full article ">Figure 4
<p>The Spearman’s coefficients between surface SOC percent (0–5 cm) and SOC percent in various depths of soil for the tidal wetlands in the southeastern coastal United States (shown in blue) and mangrove wetlands in the Odisha state of India (shown in orange) [<a href="#B69-remotesensing-14-02940" class="html-bibr">69</a>,<a href="#B113-remotesensing-14-02940" class="html-bibr">113</a>]. Data for Bhitarkanika National Park mangrove site were used with permission from Chakraborty et al. [<a href="#B69-remotesensing-14-02940" class="html-bibr">69</a>].</p>
Full article ">Figure 5
<p>Model diagram for the suggested subsurface SOC modeling technique in wetland ecosystems (ML: machine learning; DL: deep learning; MLR: multiple linear regression; GAM: generalized additive models).</p>
Full article ">
22 pages, 2760 KiB  
Article
A Mask-Guided Transformer Network with Topic Token for Remote Sensing Image Captioning
by Zihao Ren, Shuiping Gou, Zhang Guo, Shasha Mao and Ruimin Li
Remote Sens. 2022, 14(12), 2939; https://doi.org/10.3390/rs14122939 - 20 Jun 2022
Cited by 22 | Viewed by 6228
Abstract
Remote sensing image captioning aims to describe the content of images using natural language. In contrast with natural images, the scale, distribution, and number of objects generally vary in remote sensing images, making it hard to capture global semantic information and the relationships [...] Read more.
Remote sensing image captioning aims to describe the content of images using natural language. In contrast with natural images, the scale, distribution, and number of objects generally vary in remote sensing images, making it hard to capture global semantic information and the relationships between objects at different scales. In this paper, in order to improve the accuracy and diversity of captioning, a mask-guided Transformer network with a topic token is proposed. Multi-head attention is introduced to extract features and capture the relationships between objects. On this basis, a topic token is added into the encoder, which represents the scene topic and serves as a prior in the decoder to help us focus better on global semantic information. Moreover, a new Mask-Cross-Entropy strategy is designed in order to improve the diversity of the generated captions, which randomly replaces some input words with a special word (named [Mask]) in the training stage, with the aim of enhancing the model’s learning ability and forcing exploration of uncommon word relations. Experiments on three data sets show that the proposed method can generate captions with high accuracy and diversity, and the experimental results illustrate that the proposed method can outperform state-of-the-art models. Furthermore, the CIDEr score on the RSICD data set increased from 275.49 to 298.39. Full article
(This article belongs to the Topic Big Data and Artificial Intelligence)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>An example in the RSICD data set. An image is shown on the left and its corresponding descriptions are on the right.</p>
Full article ">Figure 2
<p>The structure of proposed mask-guided Transformer with topic token. For the sake of clarity, the AddNorm and embedding operations are not shown.</p>
Full article ">Figure 3
<p>Some remote sensing images and corresponding topics.</p>
Full article ">Figure 4
<p>The structure of topic encoder.</p>
Full article ">Figure 5
<p>Two examples of extracting topics from data sets. The total topic names are also shown.</p>
Full article ">Figure 6
<p>Some images and two image-captioning pairs in Sydney data set.</p>
Full article ">Figure 7
<p>Some images and two image-captioning pairs in RSICD data sets. It can be observed that the large scene has various objects and the captions are abundant.</p>
Full article ">Figure 8
<p>Sentence comparison under different training strategies. The GT is one of the ground-truth sentences, CE with SC and Mask-CE are two different sentences generated in two corresponding training stages. Some common words are colored in red, and the diverse words are colored in blue.</p>
Full article ">Figure 9
<p>Visualization of cross-modal multi-head attention in decoder.</p>
Full article ">
18 pages, 3758 KiB  
Article
Atmospheric Effects and Precursors of Rainfall over the Swiss Plateau
by Wenyue Wang and Klemens Hocke
Remote Sens. 2022, 14(12), 2938; https://doi.org/10.3390/rs14122938 - 20 Jun 2022
Cited by 5 | Viewed by 2299
Abstract
In this study, we investigate the characteristics of atmospheric parameters before, during, and after rain events in Bern, Switzerland. Ground-based microwave radiometer data of the TROpospheric WAter RAdiometer (TROWARA) with a time resolution of 7 s, observations of a weather station, and the [...] Read more.
In this study, we investigate the characteristics of atmospheric parameters before, during, and after rain events in Bern, Switzerland. Ground-based microwave radiometer data of the TROpospheric WAter RAdiometer (TROWARA) with a time resolution of 7 s, observations of a weather station, and the composite analysis method are used to derive the temporal evolution of rain events and to identify possible rainfall precursors during a 10-year period (1199 available rain events). A rainfall climatology is developed using parameters integrated water vapor (IWV), integrated liquid water (ILW), rain rate, infrared brightness temperature (TIR), temperature, pressure, relative humidity, wind speed, and air density. It was found that the IWV is reduced by about 2.2 mm at the end of rain compared to the beginning. IWV and TIR rapidly increase to a peak at the onset of the rainfall. Precursors of rainfall are that the temperature reaches its maximum around 30 to 60 min before rain, while the pressure and relative humidity are minimal. IWV fluctuates the most before rain (obtained with a 10 min bandpass). In 60% of rain events, the air density decreases 2 to 6 h before the onset of rain. The seasonality and the duration of rain events as well as the diurnal cycle of atmospheric parameters are also considered. Thus, a prediction of rainfall is possible with a true detection rate of 60% by using the air density as a precursor. Further improvements in the nowcasting of rainfall are possible by using a combination of various atmospheric parameters which are monitored by a weather station and a ground-based microwave radiometer. Full article
Show Figures

Figure 1

Figure 1
<p>Composites of (<b>a</b>) integrated water vapor (IWV, green solid line) and (<b>b</b>) integrated liquid water (ILW, green dashed line) 8 h before and 16 h after rain, (<b>c</b>) the IWV composites in summer (<span class="html-italic">N</span> = 345) and winter (<span class="html-italic">N</span> = 273), as well as (<b>d</b>) the IWV composites for long (<span class="html-italic">N</span> = 428)- and short (<span class="html-italic">N</span> = 497)-duration rain events. The subplot is from 5 min before rain to 5 min after rain. IWV and ILW acquired by the TROpospheric WAter RAdiometer (TROWARA) in Bern. The shaded area shows the standard deviation of the mean (error of the mean <math display="inline"><semantics> <mrow> <mi>σ</mi> <mo>/</mo> <msqrt> <mi>n</mi> </msqrt> </mrow> </semantics></math>). Time <span class="html-italic">t</span> is the duration of rainfall. Short horizontal lines mark the beginning and end of rainfall.</p>
Full article ">Figure 2
<p>Composites of (<b>a</b>) IWV (green solid line), T<math display="inline"><semantics> <msub> <mrow/> <mi>IR</mi> </msub> </semantics></math> (red solid line), rain rate (blue area) 8 h before and 16 h during rain, considering the rain events (<b>b</b>) in summer and winter, as well as (<b>c</b>) with long and short duration.</p>
Full article ">Figure 3
<p>Composites of (<b>a</b>) temperature (red solid line), pressure (green solid line), relative humidity (blue solid line), and wind speed (black solid line) 8 h before and 16 h during rain, considering the rain events (<b>b</b>) in summer and winter, as well as (<b>c</b>) with long and short duration. Subplot (<b>a</b>) is from 60 min before rain to the onset of rain.</p>
Full article ">Figure 4
<p>(<b>a</b>) IWV fluctuation series 16 h before and after a rain event between 02:11 UT and 07:05 UT on 8 January 2011 (red). The blue vertical line represent the occurrence/duration time of rain. (<b>b</b>) Composite of IWV fluctuation series of 759 rain events 16 h before and during rain (green).</p>
Full article ">Figure 5
<p>Composite of air density 16 h before and during rain. Red dots represent the 6 h, 2 h before rain, and the onset of the rain. The subplot is from 60 min before rain to the onset of rain. A total of 784 rain events is shown.</p>
Full article ">Figure 6
<p>Time series of TROWARA and meteorological data from 1 January at 12:00 UT to 2 January at 24:00 UT, 2016. (<b>a</b>) Time series of IWV (green solid line) and ILW (red area) provided by TROWARA, and the precipitation (black area) observed by the rain gauge at the ExWi weather station. The two purple vertical lines represent the onset of light rain and moderate rain, respectively. The purple horizontal line is ILW = 0.4 mm. (<b>b</b>) Time series of surface pressure (green solid line) and wind speed (black solid line) observed by the ExWi weather station. (<b>c</b>) Time series of surface temperature (red solid line) and relative humidity (blue dashed line) observed by the ExWi weather station, as well as air density (black solid line).</p>
Full article ">Figure 7
<p>Time series of TROWARA and meteorological data from 1–31 January 2016. (<b>a</b>) Time series of IWV (green solid line), ILW (red area), and the precipitation (black area). The purple horizontal line is ILW = 0.4 mm. The subplot is from 21:00 UT on 2 January to 6:00 UT on 3 January. (<b>b</b>) Time series of surface pressure (green solid line) and wind speed (black solid line). (<b>c</b>) Time series of surface temperature (red solid line), relative humidity (blue dashed line), and air density (black solid line).</p>
Full article ">Figure 8
<p>Time series of TROWARA and meteorological data from 24–29 January 2016. Between the two purple vertical lines is a long clear sky period that lasts from 25–28 January. (<b>a</b>) Time series of IWV (green line) and ILW (red area).The purple horizontal line is ILW = 0.4 mm. (<b>b</b>) Time series of surface pressure (green line) and wind speed (black line). (<b>c</b>) Time series of surface temperature (red line), relative humidity (blue dotted line), and air density (black line).</p>
Full article ">
11 pages, 4116 KiB  
Article
A Martian Analogues Library (MAL) Applicable for Tianwen-1 MarSCoDe-LIBS Data Interpretation
by Changqing Liu, Zhongchen Wu, Xiaohui Fu, Ping Liu, Yanqing Xin, Ayang Xiao, Hongchun Bai, Shangke Tian, Sheng Wan, Yiheng Liu, Enming Ju, Guobin Jin, Xuejin Lu, Xiaobin Qi and Zongcheng Ling
Remote Sens. 2022, 14(12), 2937; https://doi.org/10.3390/rs14122937 - 20 Jun 2022
Cited by 4 | Viewed by 2700
Abstract
China’s first Mars exploration mission, named Tianwen-1, landed on Mars on 15 May 2021. The Mars Surface Composition Detector (MarSCoDe) payload onboard the Zhurong rover applied the laser-induced breakdown spectroscopy (LIBS) technique to acquire chemical compositions of Martian rocks and soils. The quantitative [...] Read more.
China’s first Mars exploration mission, named Tianwen-1, landed on Mars on 15 May 2021. The Mars Surface Composition Detector (MarSCoDe) payload onboard the Zhurong rover applied the laser-induced breakdown spectroscopy (LIBS) technique to acquire chemical compositions of Martian rocks and soils. The quantitative interpretation of MarSCoDe-LIBS spectra needs to establish a LIBS spectral database that requires plenty of terrestrial geological standards. In this work, we selected 316 terrestrial standards including igneous rocks, sedimentary rocks, metamorphic rocks, and ores, whose chemical compositions, rock types, and chemical weathering characteristics were comparable to those of Martian materials from previous orbital and in situ detections. These rocks were crushed, ground, and sieved into powders less than <38 μm and pressed into pellets to minimize heterogeneity at the scale of laser spot. The chemical compositions of these standards were independently measured by X-ray fluorescence (XRF). Subsequently, the LIBS spectra of MAL standards were acquired using an established LIBS system at Shandong University (SDU-LIBS). In order to evaluate the performance of these standards in LIBS spectral interpretation, we established multivariate models using partial least squares (PLS) and least absolute shrinkage and selection (LASSO) algorithms to predict the abundance of major elements based on SDU-LIBS spectra. The root mean squared error (RMSE) values of these models are comparable to those of the published models for MarSCoDe, ChemCam, and SuperCam, suggesting these PLS and LASSO models work well. From our research, we can conclude that these 316 MAL targets are good candidates to acquire geochemistry information based on the LIBS technique. These targets could be regarded as geological standards to build a LIBS database using a prototype of MarSCoDe in the near future, which is critical to obtain accurate chemical compositions of Martian rocks and soils based on MarSCoDe-LIBS spectral data. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Scheme for MAL standards preparation and images of feldspathic quartz sandstone (SDU-001-SEC): (<b>a</b>) Original rock; (<b>b</b>) fragments; (<b>c</b>) powders &lt; 38 μm; (<b>d</b>) pellet.</p>
Full article ">Figure 2
<p>Pie graph of number and mass (kg) of MAL standards.</p>
Full article ">Figure 3
<p>Example of a LIBS spectrum from feldspathic quartz sandstone (SDU-001-SEC). The emission lines of eleven elements (Si, Fe, Mg, K, Na, Ti, Ca, Al, O, H, and C) were identified based on the NIST database [<a href="#B22-remotesensing-14-02937" class="html-bibr">22</a>]. The black, red, and blue segments are corresponding to the three spectral channels.</p>
Full article ">Figure 4
<p>The distribution of oxides of major elements in 408 standards for ChemCam (gray) and 316 standards for MarSCoDe.</p>
Full article ">Figure 5
<p>Geochemical characteristics of MAL standards: (<b>a</b>) Total alkali-silica plot of all 316 MAL standards, with the TES-derived data and GRS results, as well as compositions of rocks and soils at the Pathfinder landing site, Gusev crater, and Meridiani Planum from McSween et al. (2009) [<a href="#B29-remotesensing-14-02937" class="html-bibr">29</a>]. The polyline in the region of typical igneous rocks (cyan) is the alkaline–subalkaline boundary from Irvine et al. (1971) and Sautter et al. (2015) [<a href="#B27-remotesensing-14-02937" class="html-bibr">27</a>,<a href="#B28-remotesensing-14-02937" class="html-bibr">28</a>]; (<b>b</b>) ternary diagram of Al<sub>2</sub>O<sub>3</sub>-(CaO + Na<sub>2</sub>O)-K<sub>2</sub>O with CIA values on the vertical axis; (<b>c</b>) ternary diagram of Al<sub>2</sub>O<sub>3</sub>-(CaO + Na<sub>2</sub>O + K<sub>2</sub>O)-(FeO<sub>T</sub> + MgO). ChemCam results from Yellowknife Bay in the Gale crater (green polygon) are from McLennan et al. (2014) [<a href="#B30-remotesensing-14-02937" class="html-bibr">30</a>].</p>
Full article ">Figure 6
<p>RMSEs of test sets in 50 PLS and 50 LASSO models. Dash lines indicate the average RMSE of 50 RLS models (black line) and 50 LASSO models (red line).</p>
Full article ">
30 pages, 5790 KiB  
Article
Evaluation and Hydrological Application of Four Gridded Precipitation Datasets over a Large Southeastern Tibetan Plateau Basin
by Yueguan Zhang, Qin Ju, Leilei Zhang, Chong-Yu Xu and Xide Lai
Remote Sens. 2022, 14(12), 2936; https://doi.org/10.3390/rs14122936 - 19 Jun 2022
Cited by 13 | Viewed by 2510
Abstract
Reliable precipitation is crucial for hydrological studies over Tibetan Plateau (TP) basins with sparsely distributed rainfall gauges. In this study, four widely used precipitation products, including the Asian Precipitation Highly Resolved Observational Data Integration Towards Evaluation of the water resources (APHRODITE), the High [...] Read more.
Reliable precipitation is crucial for hydrological studies over Tibetan Plateau (TP) basins with sparsely distributed rainfall gauges. In this study, four widely used precipitation products, including the Asian Precipitation Highly Resolved Observational Data Integration Towards Evaluation of the water resources (APHRODITE), the High Asia Reanalysis (HAR), and the satellite-based precipitation estimates from Global Precipitation Measurement (GPM) and Tropical Rainfall Measurement Mission (TRMM), were comprehensively evaluated by combining statistical analysis and hydrological simulation over the Upper Brahmaputra (UB) River Basin of TP during 2001–2013. In respect to the statistical assessment, the overall performances of GPM and HAR are comparable to each other, and both are superior to the other two datasets. For hydrological assessment, both daily and monthly GPM-based streamflow simulations perform the best not only at the UB outlet with very good results, but they also illustrate satisfactory results at Yangcun and Lhasa hydrological stations within the UB. Runoff simulation using HAR only performs well at the UB outlet, whereas it shows poor results at both Yangcun and Lhasa stations. The simulated results based on APHRODITE and TRMM show poor performances at UB. Generally, the GPM shows an encouraging potential for hydro-meteorological investigation over UB, although with some bias in flood simulation. Full article
(This article belongs to the Section Remote Sensing in Geology, Geomorphology and Hydrology)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Location and topography of the Upper Brahmaputra River Basin. The black solid dots denote the meteorological stations.</p>
Full article ">Figure 2
<p>Spatial distributions of annual mean precipitation for reference precipitation (PCP_Sun) and four gridded precipitation datasets among 2001–2013 over the UB (mm/year).</p>
Full article ">Figure 3
<p>Mean monthly precipitation from four products (<b>a</b>) HAR, (<b>b</b>) APHRODITE, (<b>c</b>) TRMM and (<b>d</b>) GPM versus reference precipitation during 2001–2013 over the UB.</p>
Full article ">Figure 4
<p>Box plots of grid-based statistical metrics at daily scale for years 2001–2013: (<b>a</b>) RB, (<b>b</b>) RMSE, (<b>c</b>) CC, (<b>d</b>) POD, (<b>e</b>) FAR, and (<b>f</b>) CSI. The number at the base of the box plot is the median value of each statistical metric. The short dash line in box plots denotes the mean value of statistical metric. The green horizontal line indicates the optimum value. The symbol “+” in box plots represents the outlier of statistical metric value.</p>
Full article ">Figure 5
<p>Scatterplots of daily basin-wide precipitation between gridded precipitation and reference precipitation during rainy season ((<b>a</b>) HAR, (<b>c</b>) APHRODITE, (<b>e</b>) TRMM, (<b>g</b>) GPM), and non-rainy season ((<b>b</b>) HAR, (<b>d</b>) APHRODITE, (<b>f</b>) TRMM, (<b>h</b>) GPM). The solid line shows the linear least-squares result, and the diagonal dash line shows the 1:1 line.</p>
Full article ">Figure 6
<p>Daily flow simulation at Nuxia (<b>a</b>) and Lhasa hydrological stations (<b>b</b>) during the calibration period (1990–2000).</p>
Full article ">Figure 7
<p>Daily flow simulation driven by (<b>a</b>) reference precipitationand four precipitation datasets: (<b>b</b>) HAR, (<b>c</b>) APHRODITE, (<b>d</b>) TRMM and (<b>e</b>) GPM, at Nuxia hydrological station during years 2001–2013. Qobs represents the observed streamflow. Qcal(PCP_Sun), Qcal(HAR), Qcal(APHRODITE), Qcal(TRMM) and Qcal(GPM) represent the calculated streamflow using the corresponding individual precipitation dataset.</p>
Full article ">Figure 8
<p>Daily flow simulation driven by (<b>a</b>) reference precipitation and four precipitation datasets: (<b>b</b>) HAR, (<b>c</b>) APHRODITE, (<b>d</b>) TRMM and (<b>e</b>) GPM, at Lhasa hydrological station during years 2001–2013.</p>
Full article ">Figure 9
<p>Monthly flow simulation at (<b>a</b>) Nuxia, (<b>b</b>) Yangcun, and (<b>c</b>) Lhasa hydrological stations during calibration period (1990–2000).</p>
Full article ">Figure 10
<p>Monthly flow simulation driven by (<b>a</b>) reference precipitation and four precipitation datasets: (<b>b</b>) HAR, (<b>c</b>) APHRODITE, (<b>d</b>) TRMM and (<b>e</b>) GPM, at Nuxia hydrological station during years 2001–2013.</p>
Full article ">Figure 11
<p>Monthly flow simulation driven by (<b>a</b>) reference precipitation and four precipitation datasets: (<b>b</b>) HAR, (<b>c</b>) APHRODITE, (<b>d</b>) TRMM and (<b>e</b>) GPM, at Yangcun hydrological station during years 2001–2013.</p>
Full article ">Figure 12
<p>Monthly flow simulation driven by (<b>a</b>) reference precipitation and four precipitation datasets: (<b>b</b>) HAR, (<b>c</b>) APHRODITE, (<b>d</b>) TRMM and (<b>e</b>) GPM, at Lhasa hydrological station during years 2001–2013.</p>
Full article ">Figure 13
<p>Flow duration curves of the observed and simulated daily discharge at (<b>a</b>) Nuxia and (<b>b</b>) Lhasa hydrological stations.</p>
Full article ">Figure 14
<p>Relationship between the statistical metrics and the altitude for GPM: (<b>a</b>) POD, (<b>b</b>) FAR, and (<b>c</b>) CSI. Black lines show the linear trends estimated with the least-squares-error method. CC, correlation coefficient.</p>
Full article ">Figure 15
<p>Box plots of 15 CMA stations-based statistical metrics over the UB at the grid scale: (<b>a</b>) RB, (<b>b</b>) RMSE, (<b>c</b>) CC, (<b>d</b>) POD, (<b>e</b>) FAR, and (<b>f</b>) CSI at daily scale for years 2001–2013. The number at the base of the box plot is the median value of each statistical metric. The short dash line in box plots denotes the mean value of statistical metric. The green horizontal line indicates the optimum value.</p>
Full article ">
16 pages, 3684 KiB  
Article
Space-Time Adaptive Processing Clutter-Suppression Algorithm Based on Beam Reshaping for High-Frequency Surface Wave Radar
by Jiaming Li, Qiang Yang, Xin Zhang, Xiaowei Ji and Dezhu Xiao
Remote Sens. 2022, 14(12), 2935; https://doi.org/10.3390/rs14122935 - 19 Jun 2022
Cited by 3 | Viewed by 2695
Abstract
In high-frequency surface wave radar (HFSWR) systems, clutter is a common phenomenon that causes objects to be submerged. Space-time adaptive processing (STAP), which uses two-dimensional data to increase the degrees of freedom, has recently become a crucial tool for clutter suppression in advanced [...] Read more.
In high-frequency surface wave radar (HFSWR) systems, clutter is a common phenomenon that causes objects to be submerged. Space-time adaptive processing (STAP), which uses two-dimensional data to increase the degrees of freedom, has recently become a crucial tool for clutter suppression in advanced HFSWR systems. However, in STAP, the pattern is distorted if a clutter component is contained in the main lobe, which leads to errors in estimating the target angle and Doppler frequency. To solve the main-lobe distortion problem, this study developed a clutter-suppression method based on beam reshaping (BR). In this method, clutter components were estimated and maximally suppressed in the side lobe while ensuring that the main lobe remained intact. The results of the proposed algorithm were evaluated by comparison with those of standard STAP and sparse-representation STAP (SR-STAP). Among the tested algorithms, the proposed BR algorithm had the best suppression performance and the most accurate main-lobe peak response, thereby preserving the target angle and Doppler frequency information. The BR algorithm can assist with target detection and tracking despite a background with ionospheric clutter. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>(<b>a</b>) Correlation between target and clutter components in the angle domain. (<b>b</b>) Correlation between target and clutter components in the Doppler frequency domain. (<b>c</b>) Spatial gain in the target direction via the correlation between the clutter component and target for different CNRs. (<b>d</b>) Partial enlargement of (<b>c</b>).</p>
Full article ">Figure 2
<p>(<b>a</b>) Angle–Doppler map of simulation data; (<b>b</b>) angle domain curve at the target Doppler frequency; and (<b>c</b>) Doppler frequency curve in the target angle.</p>
Full article ">Figure 2 Cont.
<p>(<b>a</b>) Angle–Doppler map of simulation data; (<b>b</b>) angle domain curve at the target Doppler frequency; and (<b>c</b>) Doppler frequency curve in the target angle.</p>
Full article ">Figure 3
<p>Comparison of response patterns for the weight vector in (<b>a</b>) the angle domain and (<b>b</b>) the Doppler frequency domain for different algorithms.</p>
Full article ">Figure 4
<p>Result of the angle offset varying with (<b>a</b>) the energy of the main-lobe clutter and (<b>b</b>) the angle of the main-lobe clutter.</p>
Full article ">Figure 5
<p>(<b>a</b>) Range–Doppler map of the measured data; (<b>b</b>) angle–Doppler map of the measured data; (<b>c</b>) angle curve at the Doppler frequency of the target; and (<b>d</b>) the Doppler frequency curve at the angle of the target.</p>
Full article ">Figure 6
<p>(<b>a</b>,<b>b</b>) Angle–Doppler maps of the No. 94 and No. 95 range bins after sparse representation. Comparison of response patterns for the weight vector in (<b>c</b>) the angle domain and (<b>d</b>) Doppler frequency domain for different algorithms for the measured data.</p>
Full article ">
13 pages, 6727 KiB  
Technical Note
Eye-Safe Aerosol and Cloud Lidar Based on Free-Space Intracavity Upconversion Detection
by Wenjie Yue, Tao Chen, Wei Kong, Xin Chen, Genghua Huang and Rong Shu
Remote Sens. 2022, 14(12), 2934; https://doi.org/10.3390/rs14122934 - 19 Jun 2022
Cited by 6 | Viewed by 2635
Abstract
We report an eye-safe aerosol and cloud lidar with an Erbium-doped fiber laser (EDFL) and a free-space intracavity upconversion detector as the transmitter and receiver, respectively. The EDFL was home-made, which could produce linearly-polarized pulses at a repetition rate of 15 kHz with [...] Read more.
We report an eye-safe aerosol and cloud lidar with an Erbium-doped fiber laser (EDFL) and a free-space intracavity upconversion detector as the transmitter and receiver, respectively. The EDFL was home-made, which could produce linearly-polarized pulses at a repetition rate of 15 kHz with pulse energies of ~70 μJ and pulse durations of ~7 ns centered at 1550 nm. The echo photons were upconverted to ~631 nm via the sum frequency generation process in a bow-tie cavity, where a Nd:YVO4 and a PPLN crystal served as the pump and nonlinear frequency conversion devices, respectively. The upconverted visible photons were recorded by a photomultiplier tube and their timestamps were registered by a customized time-to-digital converter for distance-resolved measurement. Reflected signals peaked at ~6.8 km from a hard target were measured with a distance resolution of 0.6 m for an integral duration of 10 s. Atmospheric backscattered signals, with a range of ~6 km, were also detectable for longer integral durations. The evolution of aerosols and clouds were recorded by this lidar in a preliminary experiment with a continuous measuring time of over 18 h. Clear boundary and fine structures of clouds were identified with a spatial resolution of 9.6 m during the measurement, showing its great potential for practical aerosol and cloud monitoring. Full article
(This article belongs to the Special Issue Lidar for Advanced Classification and Retrieval of Aerosols)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>(<b>a</b>) Schematic diagram of noncollinear phase-matching in PPLN. (<b>b</b>) Comparison of normalized efficiencies with respect to incident angles in three different cases.</p>
Full article ">Figure 2
<p>(<b>a</b>) The layout of the designed resonant cavity of the UPD, (<b>b</b>) the beam size evolution along the cavity.</p>
Full article ">Figure 3
<p>Experimental setup of the eye-safe aerosol lidar based on free-space intracavity upconversion.</p>
Full article ">Figure 4
<p>Image of the back-illuminated field stop of the receiver (yellow circle) and the laser transmitter (pink circle) captured from the infrared camera at the focal plane of the Newtonian telescope with focal length of 1.2 m.</p>
Full article ">Figure 5
<p>(<b>a</b>) The bird view of the lidar transmit–receive path in hard target ranging test, the pink line and circles denoted the path and the detected buildings, respectively. (<b>b</b>) The view along the lidar pointing direction from the laboratory, the red arrow indicated the farthest target of the ranging test.</p>
Full article ">Figure 6
<p>The photon histogram acquired in the hard target ranging test with and without the external pump at 808 nm.</p>
Full article ">Figure 7
<p>(<b>a</b>) The original signal and its fitting results for overlap factor measurement. Red dashed curve is the fitted signal evolution from the measured far-field signal (<b>b</b>) The overlap factor retrieved by comparing the measured and fitted signal.</p>
Full article ">Figure 7 Cont.
<p>(<b>a</b>) The original signal and its fitting results for overlap factor measurement. Red dashed curve is the fitted signal evolution from the measured far-field signal (<b>b</b>) The overlap factor retrieved by comparing the measured and fitted signal.</p>
Full article ">Figure 8
<p>The acquired aerosol and cloud signal with spatial resolution of 9.8 m in different integral durations of 5 min and 60 min, respectively. (<b>a</b>) Signals at midnight; (<b>b</b>) Signals at noon.</p>
Full article ">Figure 8 Cont.
<p>The acquired aerosol and cloud signal with spatial resolution of 9.8 m in different integral durations of 5 min and 60 min, respectively. (<b>a</b>) Signals at midnight; (<b>b</b>) Signals at noon.</p>
Full article ">Figure 9
<p>Temporal and spatial evolution of the range squared signal from 25–26 December 2021 with a temporal resolution of 1 min and a spatial resolution of 9.8 m.</p>
Full article ">Figure 10
<p>Temporal and spatial evolution of the range squared signal on 25 December 2021 with a temporal resolution of 5 s and a spatial resolution of 9.8 m.</p>
Full article ">
18 pages, 5600 KiB  
Article
Influence of Spatial Resolution on Satellite-Based PM2.5 Estimation: Implications for Health Assessment
by Heming Bai, Yuli Shi, Myeongsu Seong, Wenkang Gao and Yuanhui Li
Remote Sens. 2022, 14(12), 2933; https://doi.org/10.3390/rs14122933 - 19 Jun 2022
Cited by 8 | Viewed by 2733
Abstract
Satellite-based PM2.5 estimation has been widely used to assess health impact associated with PM2.5 exposure and might be affected by spatial resolutions of satellite input data, e.g., aerosol optical depth (AOD). Here, based on Multi-Angle Implementation of Atmospheric Correction (MAIAC) AOD [...] Read more.
Satellite-based PM2.5 estimation has been widely used to assess health impact associated with PM2.5 exposure and might be affected by spatial resolutions of satellite input data, e.g., aerosol optical depth (AOD). Here, based on Multi-Angle Implementation of Atmospheric Correction (MAIAC) AOD in 2020 over the Yangtze River Delta (YRD) and three PM2.5 retrieval models, i.e., the mixed effects model (ME), the land-use regression model (LUR) and the Random Forest model (RF), we compare these model performances at different spatial resolutions (1, 3, 5 and 10 km). The PM2.5 estimations are further used to investigate the impact of spatial resolution on health assessment. Our cross-validated results show that the model performance is not sensitive to spatial resolution change for the ME and LUR models. By contrast, the RF model can create a more accurate PM2.5 prediction with a finer AOD spatial resolution. Additionally, we find that annual population-weighted mean (PWM) PM2.5 concentration and attributable mortality strongly depend on spatial resolution, with larger values estimated from coarser resolution. Specifically, compared to PWM PM2.5 at 1 km resolution, the estimation at 10 km resolution increases by 7.8%, 22.9%, and 9.7% for ME, LUR, and RF models, respectively. The corresponding increases in mortality are 7.3%, 18.3%, and 8.4%. Our results also show that PWM PM2.5 at 10 km resolution from the three models fails to meet the national air quality standard, whereas the estimations at 1, 3 and 5 km resolutions generally meet the standard. These findings suggest that satellite-based health assessment should consider the spatial resolution effect. Full article
Show Figures

Figure 1

Figure 1
<p>Panel (<b>a</b>) presents the locations of PM<sub>2.5</sub> monitoring stations (magenta dots) and DEM at 1 km spatial resolution in the Yangtze River Delta. Panel (<b>b</b>) shows the spatial distribution of the natural logarithm of population count in 2020.</p>
Full article ">Figure 2
<p>(<b>a</b>) The correlation coefficients between PM<sub>2.5</sub> concentration and AOD at the resolutions of 1, 3, 5 and 10 km over the Yangtze River Delta for different seasons. Based on the spatial AOD–PM<sub>2.5</sub> correlation for each day in summer, the box-and-whisker in panel (<b>b</b>) shows 10th, 25th, 50th, 75th and 90th percentile values of the correlation for different resolutions. Note that the spatial AOD–PM<sub>2.5</sub> correlation for a given day is excluded in panel (<b>b</b>) if the number of AOD–PM<sub>2.5</sub> pairs is less than 10.</p>
Full article ">Figure 3
<p>The spatial distributions of AOD and PM<sub>2.5</sub> concentration on 15 August 2020 at the resolutions of (<b>a</b>) 1, (<b>b</b>) 3, (<b>c</b>) 5 and (<b>d</b>) 10 km. The colored dots and background color represent PM<sub>2.5</sub> concentration and AOD, respectively. The inset in each panel is the scatterplot between PM<sub>2.5</sub> and AOD.</p>
Full article ">Figure 4
<p>Scatter plots of 10-fold cross validation (CV) results for the mixed effects (ME) model at the resolutions of (<b>a</b>) 1, (<b>b</b>) 3, (<b>c</b>) 5 and (<b>d</b>) 10 km. The color bar represents the counts of samples. R<sup>2</sup>, RMSE and N are the coefficient of determination, root-mean-square error and number of samples, respectively. The dashed line stands for the 1:1 line.</p>
Full article ">Figure 5
<p>Scatter plots of 10-fold cross validation (CV) results for the land-use regression (LUR) model at the resolutions of (<b>a</b>) 1, (<b>b</b>) 3, (<b>c</b>) 5 and (<b>d</b>) 10 km. The color bar represents the counts of samples. R<sup>2</sup>, RMSE and N are the coefficient of determination, root-mean-square error and number of samples, respectively. The dashed line stands for the 1:1 line.</p>
Full article ">Figure 6
<p>Scatter plots of 10-fold cross validation (CV) results for the Random Forest (RF) model at the resolutions of (<b>a</b>) 1, (<b>b</b>) 3, (<b>c</b>) 5 and (<b>d</b>) 10 km. The color bar represents the counts of samples. R<sup>2</sup>, RMSE and N are the coefficient of determination, root-mean-square error and number of samples, respectively. The dashed line stands for the 1:1 line.</p>
Full article ">Figure 7
<p>Spatial distribution of annual PM<sub>2.5</sub> in 2020 over the YRD region from our (<b>a</b>) 1 km RF model and (<b>b</b>) CHAP dataset (<a href="https://weijing-rs.github.io/product.html" target="_blank">https://weijing-rs.github.io/product.html</a> (accessed on 19 June 2021)).</p>
Full article ">Figure 8
<p>(<b>a</b>) Annual population-weighted mean (PWM) PM<sub>2.5</sub> and (<b>b</b>) difference in attributable mortality as a function of spatial resolution. Red, green and blue lines in each panel represent results estimated by the mixed effects (ME), the land-use regression (LUR) and the Random Forest (RF) models, respectively.</p>
Full article ">Figure 9
<p>(<b>a</b>) Spatial distribution of annual PM<sub>2.5</sub> at 1 km and (<b>b</b>–<b>d</b>) its difference from coarser resolutions. PM<sub>2.5</sub> retrievals in this figure are based on the Random Forest (RF) model. Annual PM<sub>2.5</sub> estimations at 1 km are first aggregated to arithmetic mean values at a coarse resolution. Then, these aggregated values are compared to that from the PM<sub>2.5</sub> retrieval model at the coarse resolution.</p>
Full article ">Figure 10
<p>Annual population-weighted mean (PWM) PM<sub>2.5</sub> as a function of spatial resolution. Panels (<b>a</b>,<b>b</b>) are based on PM<sub>2.5</sub> retrievals from the mixed effects (ME) and Random Forest (RF) models. Red, green and blue lines in each panel stand for sub-region 1 (northeast YRD), sub-region 2 (northwest YRD) and sub-region 3 (southern YRD), respectively. See <a href="#remotesensing-14-02933-f001" class="html-fig">Figure 1</a>a for different sub-regions.</p>
Full article ">Figure 11
<p>Spatial distribution of (top panels) annual PM<sub>2.5</sub> estimated by the Random Forest (RF) model and (bottom panels) population in northeast YRD (sub-region 1 in <a href="#remotesensing-14-02933-f001" class="html-fig">Figure 1</a>a) for the resolutions of (<b>a</b>,<b>e</b>) 1, (<b>b</b>,<b>f</b>) 3, (<b>c</b>,<b>g</b>) 5 and (<b>d</b>,<b>h</b>) 10 km. Population data in bottom panels are normalized by using standardized z-scores.</p>
Full article ">
20 pages, 5789 KiB  
Article
Variability in the Spatiotemporal Distribution Patterns of Greater Amberjack in Response to Environmental Factors in the Taiwan Strait Using Remote Sensing Data
by Mubarak Mammel, Muhamad Naimullah, Ali Haghi Vayghan, Jhen Hsu, Ming-An Lee, Jun-Hong Wu, Yi-Chen Wang and Kuo-Wei Lan
Remote Sens. 2022, 14(12), 2932; https://doi.org/10.3390/rs14122932 - 19 Jun 2022
Cited by 11 | Viewed by 2872
Abstract
The environmental characteristics of the Taiwan Strait (TS) have been linked to variations in the abundance and distribution of greater amberjack (Seriola dumerili) populations. Greater amberjack is a commercially and ecologically valuable species in ecosystems, and its spatial distribution patterns are [...] Read more.
The environmental characteristics of the Taiwan Strait (TS) have been linked to variations in the abundance and distribution of greater amberjack (Seriola dumerili) populations. Greater amberjack is a commercially and ecologically valuable species in ecosystems, and its spatial distribution patterns are pivotal to fisheries management and conservation. However, the relationship between the catch rates of S. dumerili and the environmental changes and their impact on fish communities remains undetermined in the TS. The goal of this study was to determine the spatiotemporal distribution pattern of S. dumerili with environmental characteristics in the TS from south to north (20°N–29°N and 115°E–127°E), applying generalized additive models (GAMs) and spatiotemporal fisheries data from logbooks and voyage data recorders from Taiwanese fishing vessels (2014–2017) as well as satellite-derived remote sensing environmental data. We used the generalized linear model (GLM) and GAM to analyze the effect of environmental factors and catch rates. The predictive performance of the two statistical models was quantitatively assessed by using the root mean square difference. Results reveal that the GAM outperforms the GLM model in terms of the functional relationship of the GAM for generating a reliable predictive tool. The model selection process was based on the significance of model terms, increase in deviance explained, decrease in residual factor, and reduction in Akaike’s information criterion. We then developed a species distribution model based on the best GAMs. The deviance explained indicated that sea surface temperature, linked to high catch rates, was the key factor influencing S. dumerili distributions, whereas mixed layer depth was the least relevant factor. The model predicted a relatively high S. dumerili catch rate in the northwestern region of the TS in summer, with the area extending to the East China Sea. The target species is strongly influenced by biophysical environmental conditions, and potential fishing areas are located throughout the waters of the TS. The findings of this study showed how S. dumerili populations respond to environmental variables and predict species distributions. Data on the habitat preferences and distribution patterns of S. dumerili are essential for understanding the environmental conditions of the TS, which can inform future priorities for conservation planning and management. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Spatial distribution pattern of <span class="html-italic">S. dumerili</span> caught using angling gear from 2014 to 2017 in the Taiwan Strait (TS). On the right side is a highlighted map of topographic features of bathymetry representing the isobaths, as shown by the orange line.</p>
Full article ">Figure 2
<p>Seasonal variations and distribution patterns of the fishing locations (in 0.1 spatial grids) of greater amberjack from 2014 to 2017 in the TS.</p>
Full article ">Figure 3
<p>Monthly average trends in the time series of (<b>a</b>) latitudinal (Lat.) and (<b>b</b>) longitudinal (Lon.) gravitational centers of observed catch rates (<span class="html-italic">G</span>) compared with the monthly average values of (<b>c</b>) sea surface temperature (SST), (<b>d</b>) sea surface salinity (SSS), (<b>e</b>) sea surface height (SSH), (<b>f</b>) chlorophyll-<span class="html-italic">a</span> (Chl-<span class="html-italic">a</span>), (<b>g</b>) mixed layer depth (MLD), and (<b>h</b>) eddy kinetic energy (EKE) in the TS. On the right side is the correlation between latitude and longitude in relation to all the environmental factors.</p>
Full article ">Figure 4
<p>Seasonal average observed catch rates of <span class="html-italic">S. dumerili</span> overlaid with environmental factors: SST (<b>a</b>–<b>d</b>), SSS (<b>e</b>–<b>h</b>), SSH (<b>i</b>–<b>l</b>), Chl-<span class="html-italic">a</span> (<b>m</b>–<b>p</b>), MLD (<b>q</b>–<b>t</b>), and EKE (<b>u</b>–<b>x</b>).</p>
Full article ">Figure 5
<p>The boxplot graph depicts the GAM and GLM prediction performance based on the root mean square difference (RMSD) value; the model improves if the RMSD value is close to zero.</p>
Full article ">Figure 6
<p>Estimated environmental variable effects derived from the optimal generalized additive model (GAM) analysis of the catch rates of <span class="html-italic">S. dumerili</span> in the TS: (<b>a</b>) normal quantile–quantile plots, (<b>b</b>) latitude, (<b>c</b>) longitude, (<b>d</b>) SST, (<b>e</b>) SSS, (<b>f</b>) SSH, (<b>g</b>) Chl-<span class="html-italic">a</span>, (<b>h</b>) MLD, and (<b>i</b>) EKE. The solid and black dotted lines indicate the fitted GAM function and 95 percent confidence intervals in (<b>a</b>–<b>i</b>). On the <span class="html-italic">x</span>-axis, the rug plot depicts the relative density of data points, while the <span class="html-italic">y</span>-axis depicts the results of smoothing the fitted values. Moreover, <span class="html-italic">s</span>(<span class="html-italic">x<sub>n</sub></span>) denotes each model covariate’s spline smoothing function <span class="html-italic">x<sub>n</sub></span>.</p>
Full article ">Figure 7
<p>Three-dimensional partial dependence plots representing the environmental variables with Log (Catch Rates) in the interaction between the (<b>a</b>) SST and SSS; (<b>b</b>) SST and SSH; (<b>c</b>) SST and Chl-<span class="html-italic">a</span>; (<b>d</b>) SST and MLD; (<b>e</b>) SST and EKE in the model.</p>
Full article ">Figure 8
<p>Seasonal average spatiotemporal distribution pattern of observed catch rates for <span class="html-italic">S. dumerili</span> overlaid with predicted catch rates from selected GAMs in the 2014–2017 period in the TS: (<b>a</b>) spring, (<b>b</b>) summer, (<b>c</b>) autumn, and (<b>d</b>) winter.</p>
Full article ">
17 pages, 28369 KiB  
Article
Wide and Deep Fourier Neural Network for Hyperspectral Remote Sensing Image Classification
by Jiangbo Xi, Okan K. Ersoy, Ming Cong, Chaoying Zhao, Wei Qu and Tianjun Wu
Remote Sens. 2022, 14(12), 2931; https://doi.org/10.3390/rs14122931 - 19 Jun 2022
Cited by 18 | Viewed by 3182
Abstract
Hyperspectral remote sensing image (HSI) classification is very useful in different applications, and recently, deep learning has been applied for HSI classification successfully. However, the number of training samples is usually limited, causing difficulty in use of very deep learning models. We propose [...] Read more.
Hyperspectral remote sensing image (HSI) classification is very useful in different applications, and recently, deep learning has been applied for HSI classification successfully. However, the number of training samples is usually limited, causing difficulty in use of very deep learning models. We propose a wide and deep Fourier network to learn features efficiently by using pruned features extracted in the frequency domain. It is composed of multiple wide Fourier layers to extract hierarchical features layer-by-layer efficiently. Each wide Fourier layer includes a large number of Fourier transforms to extract features in the frequency domain from a local spatial area using sliding windows with given strides.These extracted features are pruned to retain important features and reduce computations. The weights in the final fully connected layers are computed using least squares. The transform amplitudes are used for nonlinear processing with pruned features. The proposed method was evaluated with HSI datasets including Pavia University, KSC, and Salinas datasets. The overall accuracies (OAs) of the proposed method can reach 99.77%, 99.97%, and 99.95%, respectively. The average accuracies (AAs) can achieve 99.55%, 99.95%, and 99.95%, respectively. The Kappa coefficients are as high as 99.69%, 99.96%, and 99.94%, respectively. The experimental results show that the proposed method achieved excellent performance among other compared methods. The proposed method can be used for applications including classification, and image segmentation tasks, and has the ability to be implemented with lightweight embedded computing platforms. The future work is to improve the method to make it available for use in applications including object detection, time serial data prediction, and fast implementation. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Architecture of the wide and deep Fourier neural network for hyperspcetral image classification.</p>
Full article ">Figure 2
<p>Architecture of the wide Fourier neural layer.</p>
Full article ">Figure 3
<p>Classification results of Pavia University data (the unit for both horizontal and vertical axes is: pixel).</p>
Full article ">Figure 4
<p>Classification results of KSC data (the unit for both horizontal and vertical axes is: pixel).</p>
Full article ">Figure 5
<p>Classification results of Salinas data (the unit for both horizontal and vertical axes is: pixel).</p>
Full article ">Figure 6
<p>Visualization of the DFT layers of the WD-FNet. (<b>a</b>–<b>c</b>) Input training vectors of class 1 from Pavia University, KSC, and Salinas datasets, respectively. (<b>d</b>–<b>o</b>) are transformed outputs from the last sliding windows in layer 1, 2, 3, and 4 for Pavia University, KSC, and Salinas datasets, respectively.</p>
Full article ">
19 pages, 21421 KiB  
Article
Land Subsidence Monitoring and Dynamic Prediction of Reclaimed Islands with Multi-Temporal InSAR Techniques in Xiamen and Zhangzhou Cities, China
by Guangrong Li, Chaoying Zhao, Baohang Wang, Xiaojie Liu and Hengyi Chen
Remote Sens. 2022, 14(12), 2930; https://doi.org/10.3390/rs14122930 - 19 Jun 2022
Cited by 14 | Viewed by 2783
Abstract
Artificial islands and land reclamation are one of the most important ways to expand urban space in coastal cities. Long-term consolidation of reclaimed material and compaction of marine sediments can cause ground subsidence, which may threaten the buildings and infrastructure on the reclaimed [...] Read more.
Artificial islands and land reclamation are one of the most important ways to expand urban space in coastal cities. Long-term consolidation of reclaimed material and compaction of marine sediments can cause ground subsidence, which may threaten the buildings and infrastructure on the reclaimed lands. Therefore, it is crucial to monitor the land subsidence and predict the future deformation trend to mitigate the damage and take measures for the land reclamation and any infrastructure. In this paper, a total of 125 SAR images acquired by the C-band Sentinel-1A satellite between June 2017 and September 2021 are collected. The small baseline subsets (SBAS) SAR interferometry (InSAR) method is first conducted to detect the land deformation in Xiamen and Zhangzhou cities of Fujian Province, China, and the distributed scatterers (DS)-InSAR method is used to recover the complete deformation history of some typical areas including Xiamen Airport in Dadeng Island and Shuangyu Island. Then, the sequential estimation and the geotechnical model are jointly applied to demonstrate the current and future evolution of land subsidence of the constructed roads on Shuangyu Island. The results show that the maximum cumulative deformation reaches 425 mm of Xiamen Xiang’an Airport and 626 mm of Shuangyu Island, and the maximum deformation is predicted to be as large as 1.1 m by 2026 of Shuangyu Island. This research will provide important guidelines for the design and construction of Xiamen Xiang’an Airport and Shuangyu Island to prevent and control land subsidence. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Study area location and the coverage of synthetic aperture radar (SAR) image. The background is the Landsat-8 image acquired on 30 January 2021.</p>
Full article ">Figure 2
<p>Optical images of the study area. (<b>a</b>,<b>c</b>) are images of Dadeng Island acquired on 27 July 2014 and 11 February 2017, respectively. (<b>b</b>,<b>d</b>) are images of Shuangyu Island acquired on 20 December 2010 and 26 January 2017, respectively. Two regions A and B in (<b>c</b>) will be enlarged and discussed in the discussion section.</p>
Full article ">Figure 3
<p>Temporal and perpendicular baseline combination of SAR interferograms. The blue lines indicate interferograms generated between archived SAR images in the first group by the SBAS-InSAR method, and the red lines show the new interferograms generated between the newly received SAR images and older archived SAR images for the sequential estimation.</p>
Full article ">Figure 4
<p>Annual vertical deformation rate map over the whole coastal areas of Xiamen and Zhangzhou cities, China.</p>
Full article ">Figure 5
<p>(<b>a</b>) Annual vertical deformation rate map of Dadeng Island. Cross-sections of the annual vertical deformation rate of Shuangyu Island along four profiles, whose positions are marked in <a href="#remotesensing-14-02930-f005" class="html-fig">Figure 5</a>a, (<b>b</b>) Profile A–A’; (<b>c</b>) profile B–B’; (<b>d</b>) profile C–C’; (<b>e</b>) profile D–D’.</p>
Full article ">Figure 6
<p>Cumulative vertical deformation time series maps of Xiamen Xiang’an Airport.</p>
Full article ">Figure 7
<p>(<b>a</b>) Annual vertical deformation rate map of Shuangyu Island. Cross-sections of the annual vertical deformation rate of Shuangyu Island along four profiles, whose positions are marked in <a href="#remotesensing-14-02930-f007" class="html-fig">Figure 7</a>a; (<b>b</b>) profile E–E’; (<b>c</b>) profile F–F’; (<b>d</b>) profile G–G’; (<b>e</b>) profile H–H’.</p>
Full article ">Figure 8
<p>(<b>a</b>) Cumulative vertical deformation map of Shuangyu Island, and (<b>b</b>) of the roads in Shuangyu Island, where S1–S4 are the locations for deformation time series analysis in <a href="#remotesensing-14-02930-f009" class="html-fig">Figure 9</a>.</p>
Full article ">Figure 9
<p>The geotechnical model predicted vertical deformation time series at four points S1 to S4 (<b>a</b>–<b>d</b>), whose positions are marked in <a href="#remotesensing-14-02930-f008" class="html-fig">Figure 8</a>b.</p>
Full article ">Figure 10
<p>The difference between the prediction by archived deformation time sereis and the estimation by the sequential method at four points S1 to S4 (<b>a</b>–<b>d</b>).</p>
Full article ">Figure 11
<p>Cumulative vertical deformation time series maps. Cumulative deformation time series in black frame (estimated with respect to the first SAR time acquisition on 13 June 2017) from 2021 to 2026, which is obtained by the prediction model in this study.</p>
Full article ">Figure 12
<p>Four remote sensing images of Shuangyu Island were acquired on (<b>a</b>) 20 December 2010, (<b>b</b>) 1 January 2012, (<b>c</b>) 12 January 2014, and (<b>d</b>) 27 July 2021. Four remote sensing images of Dadeng Island were acquired on (<b>e</b>) 27 July 2014, (<b>f</b>) 16 June 2015, (<b>g</b>) 24 June 2016, and (<b>h</b>) 14 January 2021. The white lines indicate the first stage of land reclamation, and the yellow lines indicate the second stage.</p>
Full article ">Figure 13
<p>The vertical deformation rate maps and optical remote sensing images of regions A and B in <a href="#remotesensing-14-02930-f002" class="html-fig">Figure 2</a>c. Remote sensing image acquired on (<b>a</b>) 18 June 2019; (<b>b</b>) 26 July 2019; (<b>e</b>) 30 April 2020; and (<b>f</b>) 26 October 2020. (<b>c</b>,<b>d</b>) are annual vertical deformation rate maps.</p>
Full article ">Figure 14
<p>Deformation time series in the vertical direction at four points P1 to P4 (<b>a</b>–<b>d</b>) of regions A and B.</p>
Full article ">Figure 15
<p>Deformation time series in the vertical direction at two points P5 (<b>a</b>) and P6 (<b>b</b>) of Xiamen Xiang’an Airport, whose positions are marked in <a href="#remotesensing-14-02930-f006" class="html-fig">Figure 6</a>.</p>
Full article ">Figure 16
<p>Reclamation elevation of Shuangyu Island [<a href="#B34-remotesensing-14-02930" class="html-bibr">34</a>]. Background is the cumulative subsidence of Shuangyu Island.</p>
Full article ">
23 pages, 16367 KiB  
Article
The Different Impacts of Climate Variability and Human Activities on NPP in the Guangdong–Hong Kong–Macao Greater Bay Area
by Yanyan Wu, Zhaohui Luo and Zhifeng Wu
Remote Sens. 2022, 14(12), 2929; https://doi.org/10.3390/rs14122929 - 19 Jun 2022
Cited by 11 | Viewed by 2902
Abstract
As two main drivers of vegetation dynamics, climate variability and human activities greatly influence net primary productivity (NPP) variability by altering the hydrothermal conditions and biogeochemical cycles. Therefore, studying NPP variability and its drivers is crucial to understanding the patterns and mechanisms that [...] Read more.
As two main drivers of vegetation dynamics, climate variability and human activities greatly influence net primary productivity (NPP) variability by altering the hydrothermal conditions and biogeochemical cycles. Therefore, studying NPP variability and its drivers is crucial to understanding the patterns and mechanisms that sustain regional ecosystem structures and functions under ongoing climate variability and human activities. In this study, three indexes, namely the potential NPP (NPPp), actual NPP (NPPa), and human-induced NPP (NPPh), and their variability from 2000 to 2020 in the Guangdong–Hong Kong–Macao Greater Bay Area (GBA) were estimated and analyzed. Six main scenarios were generated based on change trends in the three indexes over the past 21 years, and the different relative impacts of climate variability and human activities on NPPa variability were quantitatively analyzed and identified. The results showed that the NPPp, NPPa, and NPPh had heterogeneous spatial distributions, and the average NPPp and NPPa values over the whole study area increased at rates of 3.63 and 6.94 gC·m2·yr−1 from 2000 to 2020, respectively, while the NPPh decreased at a rate of −4.43 gC·m2·yr−1. Climate variability and the combined effects of climate variability and human activities were the major driving factors of the NPPa increases, accounting for more than 72% of the total pixels, while the combined effects of the two factors caused the NPPa values to increase by 32–54% of the area in all cities expect Macao and across all vegetation ecosystems. Human activities often led to decreases in NPPa over more than 16% of the total pixels, and were mainly concentrated in the central cities of the GBA. The results can provide a reference for understanding NPP changes and can offer a theoretical basis for implementing ecosystem restoration, ecological construction, and conservation practices in the GBA. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>The location of the GBA, with a DEM and land covers in 2020 for the study area.</p>
Full article ">Figure 2
<p>Identification of the relative impacts of climate variability and human activities on the NPP<sub>a</sub> according to the slopes of the NPP<sub>a</sub>, NPP<sub>p</sub>, and NPP<sub>h</sub>. Six scenarios are represented in (<b>a</b>) slope NPP<sub>a</sub> &gt; 0 and (<b>b</b>) slope NPP<sub>a</sub> &lt; 0.</p>
Full article ">Figure 3
<p>Validation of the CASA-modeled NPP results for the GBA. Note that the 95% confidence intervals of slopes for the CASA model range from 0.79 to 1.33.</p>
Full article ">Figure 4
<p>Spatial patterns of mean annual (<b>a</b>) NPP<sub>p</sub>, (<b>b</b>) NPP<sub>a</sub>, and (<b>c</b>) NPP<sub>h</sub> from 2000 to 2020.</p>
Full article ">Figure 5
<p>Trends in annual (<b>a</b>) NPP<sub>p</sub>, (<b>b</b>) NPP<sub>a</sub>, and (<b>c</b>) NPP<sub>h</sub> from 2000 to 2020, and slope significance tests for (<b>d</b>) NPP<sub>p</sub>, (<b>e</b>) NPP<sub>a</sub>, and (<b>f</b>) NPP<sub>h</sub> in the study area. DS, decreased significantly (0.01&lt; <span class="html-italic">p</span> ≤ 0.05); IS, increased significantly (0.01 &lt; <span class="html-italic">p</span> ≤ 0.05); DVS, decreased very significantly (<span class="html-italic">p</span> ≤ 0.01); IVS, increased very significantly (<span class="html-italic">p</span> ≤ 0.01); D_NS, decreased or non-significantly (<span class="html-italic">p</span> &gt; 0.05); I_NS, increased non-significantly (<span class="html-italic">p</span> &gt; 0.05).</p>
Full article ">Figure 6
<p>Partial correlations between NPP<sub>a</sub> and (<b>a</b>,<b>d</b>) temperature, (<b>b</b>,<b>e</b>) precipitation, and (<b>c</b>,<b>f</b>) solar radiation in the GBA from 2000 to 2020, with corresponding frequency distributions.</p>
Full article ">Figure 7
<p>Sankey diagram for conversion between different land cover types in the GBA from 2000 to 2020. Note that different color lines represent the direction of land type transfer, and the width of the line represents the area transferred.</p>
Full article ">Figure 8
<p>Spatial patterns of (<b>a</b>) six driving forces of NPP<sub>a</sub> changes from 2000 to 2020 and (<b>b</b>) area proportion. Definitions of the driving forces are given in <a href="#remotesensing-14-02929-f002" class="html-fig">Figure 2</a>.</p>
Full article ">Figure 9
<p>The area of NPP<sub>a</sub> changes due to different driving forces in the 11 cities (<b>a</b>) and land cover types (<b>b</b>) from 2000 to 2020. Definitions of the driving forces are given in <a href="#remotesensing-14-02929-f002" class="html-fig">Figure 2</a>.</p>
Full article ">Figure 10
<p>Trends in annual (<b>a</b>) NPP<sub>a</sub>, (<b>b</b>) NPP<sub>h</sub>, and (<b>c</b>) EVI from 2000 to 2020 in the areas where NPP<sub>a</sub> was negatively influenced by human activities (ADH), and corresponding frequencies of (<b>d</b>) slope NPP<sub>a</sub>, (<b>e</b>) slope NPP<sub>h</sub>, and (<b>f</b>) slope EVI values.</p>
Full article ">Figure A1
<p>Temporal variability of (<b>a</b>) temperature, (<b>b</b>) precipitation, and (<b>c</b>) total solar radiation data in the GBA from 2000 to 2020.</p>
Full article ">
20 pages, 24276 KiB  
Article
Measuring Vegetation Heights and Their Seasonal Changes in the Western Namibian Savanna Using Spaceborne Lidars
by Farid Atmani, Bodo Bookhagen and Taylor Smith
Remote Sens. 2022, 14(12), 2928; https://doi.org/10.3390/rs14122928 - 19 Jun 2022
Cited by 5 | Viewed by 3105
Abstract
The Ice, Cloud, and Land Elevation Satellite-2 (ICESat-2) with its land and vegetation height data product (ATL08), and Global Ecosystem Dynamics Investigation (GEDI) with its terrain elevation and height metrics data product (GEDI Level 2A) missions have great potential to globally map ground [...] Read more.
The Ice, Cloud, and Land Elevation Satellite-2 (ICESat-2) with its land and vegetation height data product (ATL08), and Global Ecosystem Dynamics Investigation (GEDI) with its terrain elevation and height metrics data product (GEDI Level 2A) missions have great potential to globally map ground and canopy heights. Canopy height is a key factor in estimating above-ground biomass and its seasonal changes; these satellite missions can also improve estimated above-ground carbon stocks. This study presents a novel Sparse Vegetation Detection Algorithm (SVDA) which uses ICESat-2 (ATL03, geolocated photons) data to map tree and vegetation heights in a sparsely vegetated savanna ecosystem. The SVDA consists of three main steps: First, noise photons are filtered using the signal confidence flag from ATL03 data and local point statistics. Second, we classify ground photons based on photon height percentiles. Third, tree and grass photons are classified based on the number of neighbors. We validated tree heights with field measurements (n = 55), finding a root-mean-square error (RMSE) of 1.82 m using SVDA, GEDI Level 2A (Geolocated Elevation and Height Metrics product): 1.33 m, and ATL08: 5.59 m. Our results indicate that the SVDA is effective in identifying canopy photons in savanna ecosystems, where ATL08 performs poorly. We further identify seasonal vegetation height changes with an emphasis on vegetation below 3 m; widespread height changes in this class from two wet-dry cycles show maximum seasonal changes of 1 m, possibly related to seasonal grass-height differences. Our study shows the difficulties of vegetation measurements in savanna ecosystems but provides the first estimates of seasonal biomass changes. Full article
Show Figures

Figure 1

Figure 1
<p>Geographic overview of the study area (<b>A</b>). International borders are indicated by black lines, lakes and wetlands are shown in light blue colors, and major rivers are in blue color. (<b>B</b>): Mean NDVI derived from Landsat 8 (2014–2021) for the study area including savanna and mountainous terrain (red polygon) and an area confined to the savanna ecosystem (blue polygon) [<a href="#B31-remotesensing-14-02928" class="html-bibr">31</a>]. Major roads are shown by black lines. (<b>C</b>): Range of the NDVI values between the 10th and 90th percentile showing the variability of vegetation cover. Higher values indicate higher seasonal vegetation changes.</p>
Full article ">Figure 2
<p>ICESat-2 and GEDI ground tracks for the study region (cf. <a href="#remotesensing-14-02928-f001" class="html-fig">Figure 1</a>). (<b>A</b>) Canopy height derived from the SVDA processed ATL03 data described in this study. Red crosses indicate the 55 field-based tree height measurements of different species for validation. A detailed map of the tree-height field measurements is shown in <a href="#app1-remotesensing-14-02928" class="html-app">Supplementary Figure S2</a>. Red polygon outlines study area including more densely vegetated mountainous terrain and blue polygon is limited to the low-slope, sparsely vegetated savanna ecosystem. (<b>B</b>) Canopy height measurements taken from GEDI product L2A (version 2).</p>
Full article ">Figure 3
<p>(<b>A</b>) Sentinel 1 median polarization ratio (VV/VH) averaged from October 2018 to April 2021 extracted for ATL03 and GEDI locations (only GEDI locations are shown here). Polarization ratios indicate the amount of depolarization usually associated with scattering on vegetation, especially trees. (<b>B</b>) Sentinel 2 NDVI differences between the rainy and dry season of 2020. The full S1 polarization ratio (VV/VH) coverage is shown in <a href="#app1-remotesensing-14-02928" class="html-app">Figure S1</a>.</p>
Full article ">Figure 4
<p>Flowchart of the Sparse Vegetation Detection Algorithm (SVDA) using ICESat-2 ATL03 data.</p>
Full article ">Figure 5
<p>Characteristic example of ground-photon classification with gray dots showing all extracted signal photons from the ATL03 product (strong beam, gt1l, daytime acquisition with ID: ATL03_20190116100224_02890214_004_01). (<b>A</b>) Preliminary ground photons in 30 m steps are in red after extracting photons within the 25–75th height percentiles. This corresponds to step 1 and 2 as described in <a href="#sec3dot6dot1-remotesensing-14-02928" class="html-sec">Section 3.6.1</a>. (<b>B</b>) Final ground photons are selected based on additional filtering steps and detrending of height (steps 3 and 4 in <a href="#sec3dot6dot1-remotesensing-14-02928" class="html-sec">Section 3.6.1</a>).</p>
Full article ">Figure 6
<p>Characteristic filtering steps of a canopy classification exemplified on the strong beam, gt1l, daytime acquisition with ID: ATL03_20190116100224_02890214_004_01: (<b>A</b>) Ground photons in red derived from steps described in <a href="#sec3dot6dot1-remotesensing-14-02928" class="html-sec">Section 3.6.1</a> and shown in <a href="#remotesensing-14-02928-f005" class="html-fig">Figure 5</a>. (<b>B</b>) Preliminary canopy photons (green) as described in the first step in <a href="#sec3dot6dot2-remotesensing-14-02928" class="html-sec">Section 3.6.2</a>. (<b>C</b>) The second step described in <a href="#sec3dot6dot2-remotesensing-14-02928" class="html-sec">Section 3.6.2</a> using height and neighbor filtering distinguishes between top of canopy (green), canopy (blue) photons, and the remaining noise (gray) photons. (<b>D</b>) Final photons classification with a cubic spline interpolation of the ground photons (black line).</p>
Full article ">Figure 7
<p>Signal photons extraction using the signal confidence flag for (<b>A</b>) strong beam during day-time acquisition, (<b>B</b>) strong beam during night-time acquisition, (<b>C</b>) weak beam during day-time acquisition, and (<b>D</b>) weak beam during night-time acquisition.</p>
Full article ">Figure 8
<p>Copernicus DEM elevation and ATL03 SVDA based on ATL03 ground height measurements difference as described in <a href="#sec3dot6dot1-remotesensing-14-02928" class="html-sec">Section 3.6.1</a> (<b>A</b>). (<b>B</b>): Copernicus DEM elevation and ATL08 ground height difference. We note that two ATL08 ground height points resulted in a difference more than 40 m compared to the ground elevation from Copernicus. (<b>C</b>): Copernicus DEM elevation and GEDI difference. rRMSE shows the relative RMSE (RMSE/mean of elevation difference). All elevation points are shown in this figure; see <a href="#app1-remotesensing-14-02928" class="html-app">Supplementary Figure S8</a> for a comparison only of DEM slopes less than 5 degrees.</p>
Full article ">Figure 9
<p>Field measurements (<span class="html-italic">n</span> = 55) intersection with ATL03 SVDA based on ATL03 canopy height measurements as described in <a href="#sec3dot6dot1-remotesensing-14-02928" class="html-sec">Section 3.6.1</a> and <a href="#sec3dot6dot2-remotesensing-14-02928" class="html-sec">Section 3.6.2</a> within a buffer of 10 m (<b>A</b>). (<b>B</b>): Field measurements intersection with ATL08 within a buffer of 30 m. We note that one canopy height field measurements of 5.4 m resulted in an ATL08 canopy height of 31.24 m and we do not show this outlier to keep axes constant between the plots. (<b>C</b>): Field measurements intersection with GEDI within a buffer of 30 m. Blue line shows the weighted least squares regression, where the inverse of canopy height difference between GEDI/ICESat-2 (ATL03 SVDA and ATL08) and field measurements were used as weights. rRMSE shows the relative RMSE (RMSE/mean of the height differences between ATL and field measurements).</p>
Full article ">Figure 10
<p>ATL03 SVDA and ATL08 canopy height relationship (<b>A</b>,<b>B</b>) for the study areas shown red and blue polygons in <a href="#remotesensing-14-02928-f001" class="html-fig">Figure 1</a>B,C and <a href="#remotesensing-14-02928-f002" class="html-fig">Figure 2</a>. (<b>C</b>) The full distribution of canopy height differences, showing height differences of several meters in places and generally higher values of ATL08 estimates. rRMSE shows the relative RMSE (RMSE/mean of canopy height difference).</p>
Full article ">Figure 11
<p>GEDI L2A version 2 and ATL03 SVDA canopy height relationship (<b>A</b>) and the distribution of ATL03 SVDA minus GEDI canopy heights (<b>B</b>). Comparison is based on 600 overlapping measurements within a buffer of 5 m and was carried out on the entire study area shown by the red polygon in <a href="#remotesensing-14-02928-f001" class="html-fig">Figure 1</a>B,C and <a href="#remotesensing-14-02928-f002" class="html-fig">Figure 2</a>. GEDI canopy heights are generally lower than ATL03 SVDA canopy heights. Red line shows the weighted least squares regression, where the inverse of canopy height difference between ATL03 SVDA and GEDI measurements were used as weights.</p>
Full article ">Figure 12
<p>GEDI and ATL08 canopy height relationship (<b>A</b>) and the distribution of ATL08 minus GEDI canopy heights (<b>B</b>). Comparison is based on 91 overlapping measurements within a buffer of 5 m and was carried out on the entire study area shown by the red polygon in <a href="#remotesensing-14-02928-f001" class="html-fig">Figure 1</a>B,C and <a href="#remotesensing-14-02928-f002" class="html-fig">Figure 2</a>. Red line shows the weighted least squares-regression, where the inverse of canopy height difference between ATL08 and GEDI measurements were used as weights.</p>
Full article ">Figure 13
<p>Top of canopy seasonal changes from ATL03 SVDA (<b>A</b>–<b>C</b>), ATL08 (<b>D</b>–<b>F</b>), and GEDI (<b>G</b>–<b>I</b>) data between the dry and the rainy seasons of 2019 and 2020. First and second column show the kernel density estimation (KDE) plots of the canopy height and NDVI difference distribution using a gaussian kernel with 100 equal steps for the x (NDVI difference from –0.05 to 0.35) and y axes (canopy height from 3 to 6.2 m for ATL03 SVDA and GEDI and canopy height from 3 to 8.5 m for ATL08) data. Left column: dry seasons, center column: rainy seasons. Third column shows the distributions indicating the small (but statistically significant) differences.</p>
Full article ">Figure 14
<p>Seasonal vegetation height changes from ATL03 SVDA data between the dry and the rainy seasons of 2019 and 2020 for 231,534 (dry season) and 233,947 (wet season) points for the entire study area shown in the red polygon in <a href="#remotesensing-14-02928-f001" class="html-fig">Figure 1</a>B,C and <a href="#remotesensing-14-02928-f002" class="html-fig">Figure 2</a>. (<b>A</b>) is the kernel density estimation of the joint distribution between NDVI difference and ATL03 SVDA dry season vegetation heights with a gaussian smoothing kernel and a step size of 100 in the x and y directions (with similar parameters in (<b>B</b>)). Left column: dry seasons, center column: rainy seasons. The differences are visible at all percentiles, but most strongly for the 30% of the data at higher vegetation heights. The distribution of the dry and rainy season vegetation heights (<b>C</b>) show the largest differences for higher vegetation height, which is expected in the seasonal savanna ecosystem. A non-parametric KS tests indicate that the seasonal vegetation heights are drawn from different distributions (<span class="html-italic">p</span> = 0.0).</p>
Full article ">Figure 15
<p>Canopy height from spaceborne lasers and Sentinel-1 VV/VH median relationships. (<b>A</b>) ATL03 SVDA canopy height and VV/VH median relationship, (<b>B</b>) ATL08 canopy height and VV/VH median relationship, and (<b>C</b>) GEDI canopy height and VV/VH median relationship. We note the weak, but statistically significant linear relation between VV/VH and canopy height measurements.</p>
Full article ">
24 pages, 4691 KiB  
Article
UAV Remote Sensing for High-Throughput Phenotyping and for Yield Prediction of Miscanthus by Machine Learning Techniques
by Giorgio Impollonia, Michele Croci, Andrea Ferrarini, Jason Brook, Enrico Martani, Henri Blandinières, Andrea Marcone, Danny Awty-Carroll, Chris Ashman, Jason Kam, Andreas Kiesel, Luisa M. Trindade, Mirco Boschetti, John Clifton-Brown and Stefano Amaducci
Remote Sens. 2022, 14(12), 2927; https://doi.org/10.3390/rs14122927 - 19 Jun 2022
Cited by 18 | Viewed by 5220
Abstract
Miscanthus holds a great potential in the frame of the bioeconomy, and yield prediction can help improve Miscanthus’ logistic supply chain. Breeding programs in several countries are attempting to produce high-yielding Miscanthus hybrids better adapted to different climates and end-uses. Multispectral images acquired [...] Read more.
Miscanthus holds a great potential in the frame of the bioeconomy, and yield prediction can help improve Miscanthus’ logistic supply chain. Breeding programs in several countries are attempting to produce high-yielding Miscanthus hybrids better adapted to different climates and end-uses. Multispectral images acquired from unmanned aerial vehicles (UAVs) in Italy and in the UK in 2021 and 2022 were used to investigate the feasibility of high-throughput phenotyping (HTP) of novel Miscanthus hybrids for yield prediction and crop traits estimation. An intercalibration procedure was performed using simulated data from the PROSAIL model to link vegetation indices (VIs) derived from two different multispectral sensors. The random forest algorithm estimated with good accuracy yield traits (light interception, plant height, green leaf biomass, and standing biomass) using 15 VIs time series, and predicted yield using peak descriptors derived from these VIs time series with root mean square error of 2.3 Mg DM ha−1. The study demonstrates the potential of UAVs’ multispectral images in HTP applications and in yield prediction, providing important information needed to increase sustainable biomass production. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Field experiment locations: PAC 1 is situated in Piacenza (North-West Italy) and TWS in Aberystwyth (Mid-West Wales).</p>
Full article ">Figure 2
<p>Linear regressions of each VI simulated by the PROSAIL model, between the two sensors (MicaSense and SlantRange). The color code represents the point count distribution scaled to maximum of 1 in each hexagon. The blue dashed line represents the 1:1 relationship and the black line represents the linear regression.</p>
Full article ">Figure 3
<p>Importance of the RF models variables for crop trait estimations and for yield prediction, expressed as the drop-out loss of model performance (RMSE) for each variable related to the drop-out loss of the full model (dotted line). * The RMSE values are in (%), (cm), (Mg DM ha<sup>−1</sup>), (Mg DM ha<sup>−1</sup>), and (Mg DM ha<sup>−1</sup>), respectively, for the light interception, plant height, green leaf biomass, standing biomass, and yield.</p>
Full article ">Figure 4
<p>Frequency distribution of <span class="html-italic">Miscanthus</span> traits at the two locations in PAC 1 and TWS 1: (<b>a</b>) light interception (%), (<b>b</b>) plant height (cm), (<b>c</b>) green leaf biomass (Mg DM ha<sup>−1</sup>), and (<b>d</b>) standing biomass (Mg DM ha<sup>−1</sup>).</p>
Full article ">Figure 5
<p>Estimated vs. measured crop traits on the ground of four <span class="html-italic">Miscanthus</span> hybrids’ growth at PAC 1 and TWS 1: (<b>a</b>) light interception (%), (<b>b</b>) plant height (cm), (<b>c</b>) green leaf biomass (Mg DM ha<sup>−1</sup>), and (<b>d</b>) standing biomass (Mg DM ha<sup>−1</sup>). The blue dashed line represents the 1:1 relationship and the black line represents the linear regression.</p>
Full article ">Figure 6
<p>NRMSE values of the RF models obtained for each crop trait assessed in four hybrids grown at the two locations in PAC 1 and TWS 1. Note: lower values indicate higher estimation accuracies.</p>
Full article ">Figure 7
<p>Frequency distribution of yield (Mg DM ha<sup>−1</sup>) for the different <span class="html-italic">Miscanthus</span> hybrids at the two locations, PAC 1 and TWS 1.</p>
Full article ">Figure 8
<p>The (<b>a</b>) predicted vs. measured yield (Mg DM ha<sup>−1</sup>), (<b>b</b>) NRMSE of the RF model for yield prediction, and (<b>c</b>) boxplot of the modified days of the year (DOY<sub>M</sub>) of the peak derived from the complete time series of the 5 VIs (greenWDRVI, GNDVI, MTVI2, NDVI, and WDRVI,) of different <span class="html-italic">Miscanthus</span> hybrids, at two locations, PAC 1 and TWS 1.</p>
Full article ">Figure 9
<p>(<b>a</b>) Time series of the five VIs (GNDVI, greenWDRVI, MTVI2, NDVI, and WDRVI) fitted via the generalized additive model (GAM) throughout the growing season of <span class="html-italic">Miscanthus</span> in PAC 1. Modified days of the year (DOY<sub>M</sub>) were calculated by adding 365 to the DOY of the corresponding year from January on. (<b>b</b>) Variation of the peak of the VIs derived from the complete time series of the VIs as compared to the peak of the VIs derived from the partial time series of the VIs. The DOY<sub>M</sub> of the UAV flights performed during the season in PAC 1 are reported in the x-axis. In the y-axis, the peak differences between the peak derived to the end of the season in PAC 1 (397 DOY<sub>M</sub>) and the peak derived from partial time series fitted until the DOY<sub>M</sub> of the UAV flight, are reported.</p>
Full article ">Figure 10
<p>(<b>a</b>) Time series of the five VIs (GNDVI, greenWDRVI, MTVI2, NDVI, and WDRVI) fitted via the generalized additive model (GAM) throughout the growing season of <span class="html-italic">Miscanthus</span> in TWS 1. Modified days of the year (DOY<sub>M</sub>) were calculated by adding 365 to the DOY of the corresponding year from January on. (<b>b</b>) Variation of the peak of the VIs derived from the complete time series of the VIs as compared to the peak of the VIs derived from the partial time series of the VIs. In the x-axis, the DOY<sub>M</sub> of the UAV flights performed during the season in TWS 1 are reported. In the y-axis, the peak differences between the peak derived to the end of the season in TWS 1 (421 DOY<sub>M</sub>) and the peak derived from partial time series fitted until the DOY<sub>M</sub> of the UAV flight are reported.</p>
Full article ">Figure 11
<p>NRMSE values of the RF model for yield prediction trained with the peak of the VIs, derived from the complete time series and tested with the peak of the VIs derived from the partial time series fitted until the modified days of the year (DOY<sub>M</sub>) of the UAV flight, at the two locations PAC 1 and TWS 1.</p>
Full article ">
17 pages, 8705 KiB  
Article
Smartphone-Based Unconstrained Step Detection Fusing a Variable Sliding Window and an Adaptive Threshold
by Ying Xu, Guofeng Li, Zeyu Li, Hao Yu, Jianhui Cui, Jin Wang and Yu Chen
Remote Sens. 2022, 14(12), 2926; https://doi.org/10.3390/rs14122926 - 19 Jun 2022
Cited by 4 | Viewed by 2970
Abstract
Step detection for smartphones plays an important role in the pedestrian dead reckoning (PDR) for indoor positioning. Aiming at the problem of low step detection accuracy of smartphones in complex unconstrained states in PDR, smartphone-based unconstrained step detection method fusing a variable sliding [...] Read more.
Step detection for smartphones plays an important role in the pedestrian dead reckoning (PDR) for indoor positioning. Aiming at the problem of low step detection accuracy of smartphones in complex unconstrained states in PDR, smartphone-based unconstrained step detection method fusing a variable sliding window and an adaptive threshold is proposed. In this method, the dynamic updating algorithm of a peak threshold is developed, and the minimum peak value filtered after a sliding window filter is used as the adaptive peak threshold, which solves the problem that the peak threshold of different motion states is difficult to update adaptively. Then, a variable sliding window collaborative time threshold method is proposed, which solves the problem that the adjacent windows cannot be contacted, and the initial peak and the end peak are difficult to accurately identify. To evaluate the performance of the proposed unconstrained step detection algorithm, 50 experiments in constrained and unconstrained states are conducted by 25 volunteers holding 21 different types of smartphones. Experimental results show: The average step counting accuracy of the proposed unconstrained step detection algorithm is over 98%. Compared with the open source program Stepcount, the average step counting accuracy of the proposed algorithm is improved by 10.0%. The smartphone-based unconstrained step detection fusing a variable sliding window and an adaptive threshold has a strong ability to adapt to complex unconstrained states, and the average step counting accuracy rate is only 0.6% lower than that of constrained states. This algorithm has a wide audience and is friendly for different genders and smartphones with different prices. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Triaxial acceleration and overall acceleration when walking at will.</p>
Full article ">Figure 2
<p>Raw overall acceleration signal.</p>
Full article ">Figure 3
<p>FIR low-pass filter and sliding window filter.</p>
Full article ">Figure 4
<p>Single step action decomposition.</p>
Full article ">Figure 5
<p>Spectrum of pedestrian stroll walking, normal walking and running: (<b>a</b>) Stroll walking (<b>b</b>) Normal walking (<b>c</b>) Running.</p>
Full article ">Figure 6
<p>Flow chart of the unconstrained step detection algorithm.</p>
Full article ">Figure 7
<p>Pseudo-peak elimination by adaptive threshold in running state.</p>
Full article ">Figure 8
<p>Pseudo-peaks in front of the window and at the end of the window.</p>
Full article ">Figure 9
<p>Flow chart of pseudo-peak elimination method with variable sliding window cooperative time threshold.</p>
Full article ">Figure 10
<p>Schematic diagram of the variable sliding window: (<b>a</b>) The <span class="html-italic">i</span> (<span class="html-italic">i</span> ≥ 1) window contains a peak (<b>b</b>) The <span class="html-italic">i</span> window does not contain a true peak. Note: the green arrows indicate the starting point of new windows.</p>
Full article ">Figure 11
<p>Brands and models of 25 smartphones.</p>
Full article ">Figure 12
<p>Constrained experiment process.</p>
Full article ">Figure 13
<p>Unconstrained experiment process: (<b>a</b>) Schematic diagram of the second group of experimental process (<b>b</b>) Schematic diagram of the fourth group of experimental process.</p>
Full article ">Figure 14
<p>Step counting accuracy in constrained state.</p>
Full article ">Figure 15
<p>Step counting accuracy in unconstrained state.</p>
Full article ">Figure 16
<p>Step counting accuracy rate of smartphones with different prices in unconstrained state.</p>
Full article ">
19 pages, 4850 KiB  
Article
Attention-Unet-Based Near-Real-Time Precipitation Estimation from Fengyun-4A Satellite Imageries
by Yanbo Gao, Jiping Guan, Fuhan Zhang, Xiaodong Wang and Zhiyong Long
Remote Sens. 2022, 14(12), 2925; https://doi.org/10.3390/rs14122925 - 18 Jun 2022
Cited by 26 | Viewed by 3341
Abstract
Reliable near-real-time precipitation estimation is crucial for scientific research and resistance to natural disasters such as floods. Compared with ground-based precipitation measurements, satellite-based precipitation measurements have great advantages, but precipitation estimation based on satellite is still a challenging issue. In this paper, we [...] Read more.
Reliable near-real-time precipitation estimation is crucial for scientific research and resistance to natural disasters such as floods. Compared with ground-based precipitation measurements, satellite-based precipitation measurements have great advantages, but precipitation estimation based on satellite is still a challenging issue. In this paper, we propose a deep learning model named Attention-Unet for precipitation estimation. The model utilizes the high temporal, spatial and spectral resolution data of the FY4A satellite to improve the accuracy of precipitation estimation. To evaluate the effectiveness of the proposed model, we compare it with operational near-real-time satellite-based precipitation products and deep learning models which proved to be effective in precipitation estimation. We use classification metrics such as Probability of detection (POD), False Alarm Ratio (FAR), Critical success index (CSI), and regression metrics including Root Mean Square Error (RMSE) and Pearson correlation coefficient (CC) to evaluate the performance of precipitation identification and precipitation amounts estimation, respectively. Furthermore, we select an extreme precipitation event to validate the generalization ability of our proposed model. Statistics and visualizations of the experimental results show the proposed model has better performance than operational precipitation products and baseline deep learning models in both precipitation identification and precipitation amounts estimation. Therefore, the proposed model has the potential to serve as a more accurate and reliable satellite-based precipitation estimation product. This study suggests that applying an appropriate deep learning algorithm may provide an opportunity to improve the quality of satellite-based precipitation products. Full article
(This article belongs to the Topic Advanced Research in Precipitation Measurements)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Visualized structure of Attention-Unet network.</p>
Full article ">Figure 2
<p>Visualized structure of Unet network.</p>
Full article ">Figure 3
<p>POD, CSI, and FAR of CMORPH, FY4A-QPE and Attention-Unet model over the southeastern coast of China during the validation period (June–August 2020). (<b>a</b>–<b>c</b>) POD; (<b>d</b>–<b>f</b>) CSI; and (<b>g</b>–<b>i</b>) FAR. (<b>a</b>) POD: FY4A-QPE; (<b>b</b>) POD: CMORPH; (<b>c</b>) POD: Attention-Unet; (<b>d</b>) CSI: FY4A-QPE; (<b>e</b>) CSI: CMORPH; (<b>f</b>) CSI: Attention-Unet; (<b>g</b>) FAR: FY4A-QPE; (<b>h</b>) FAR: CMORPH; (<b>i</b>) FAR: Attention-Unet.</p>
Full article ">Figure 4
<p>The Pearson correlation coefficient (CC) and Root Mean Square Error (RMSE) values for the CMORPH, FY4A-QPE and Attention-Unet model over the southeastern coast of China during the validation period (June–August 2020). (<b>a</b>) RMSE: FY4A-QPE; (<b>b</b>) RMSE: CMORPH; (<b>c</b>) RMSE: Attention-Unet; (<b>d</b>) CC: FY4A-QPE; (<b>e</b>) CC: CMORPH; (<b>f</b>) CC: Attention-Unet.</p>
Full article ">Figure 5
<p>Visualization of precipitation identification performance of CMORPH, FY4A-QPE and Attention-Unet model over the Henan province at 0430 UTC 19 July 2021. (<b>a</b>) FY4A-QPE; (<b>b</b>) CMORPH; (<b>c</b>) Attention-Unet.</p>
Full article ">Figure 6
<p>POD, CSI, and FAR of Unet, PERSIANN-CNN and Attention-Unet model over the southeastern coast of China during the validation period (June–August 2020). (<b>a</b>–<b>c</b>) POD; (<b>d</b>–<b>f</b>) FAR; and (<b>g</b>–<b>i</b>) CSI. (<b>a</b>) POD: Unet; (<b>b</b>) POD: PERSIANN-CNN; (<b>c</b>) POD: Attention-Unet; (<b>d</b>) FAR: Unet; (<b>e</b>) FAR: PERSIANN-CNN; (<b>f</b>) FAR: Attention-Unet; (<b>g</b>) CSI: Unet; (<b>h</b>) CSI: PERSIANN-CNN; (<b>i</b>) CSI: Attention-Unet.</p>
Full article ">Figure 7
<p>The Pearson correlation coefficient (CC) and Root Mean Square Error (RMSE) values of Unet, PERSIANN-CNN and Attention-Unet model over the southeastern coast of China during the validation period (June–August 2020). (<b>a</b>–<b>c</b>) RMSE and (<b>d</b>–<b>f</b>) CC. (<b>a</b>) RMSE: Unet; (<b>b</b>) RMSE: PERSIANN-CNN; (<b>c</b>) RMSE: Attention-Unet; (<b>d</b>) CC: Unet; (<b>e</b>) CC: PERSIANN-CNN; (<b>f</b>) CC: Attention-Unet.</p>
Full article ">Figure 8
<p>Visualization of precipitation identification performance of Unet, PERSIANN-CNN and Attention-Unet model over the Henan province at 0430 UTC 19 July 2021. (<b>a</b>) Unet; (<b>b</b>) PERSIANN-CNN; (<b>c</b>) Attention-Unet.</p>
Full article ">
Previous Issue
Next Issue
Back to TopTop