Nothing Special   »   [go: up one dir, main page]

Next Issue
Volume 15, July-1
Previous Issue
Volume 15, 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 15, Issue 12 (June-2 2023) – 256 articles

Cover Story (view full-size image): This review article, with its principal basis in French research, describes data, tools and methods that use remote sensing (RS) to support the spatial predictions of soil properties, and discusses their pros and cons. The review demonstrates that RS data are frequently used in soil mapping, (i) by considering them as a substitute for analytical measurements (left part of the graph), or (ii) by considering them as covariates related to the controlling factors of soil formation and evolution used in digital soil mapping (DSM) approaches (right part of the graph). It further highlights the great potential of RS imagery to improve DSM, providing an overview of the primary challenges and prospects related to DSM and future RS sensors. The discussion opens up broad prospects for the use of RS for DSM and natural resource monitoring. 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:
36 pages, 9387 KiB  
Article
Solid Angle Geometry-Based Modeling of Volume Scattering with Application in the Adaptive Decomposition of GF-3 Data of Sea Ice in Antarctica
by Dong Li, He Lu and Yunhua Zhang
Remote Sens. 2023, 15(12), 3208; https://doi.org/10.3390/rs15123208 - 20 Jun 2023
Viewed by 2659
Abstract
Over the last two decades, spaceborne polarimetric synthetic aperture radar (PolSAR) has been widely used to penetrate sea ice surfaces to achieve fully polarimetric high-resolution imaging at all times of day and in a range of weather conditions. Model-based polarimetric decomposition is a [...] Read more.
Over the last two decades, spaceborne polarimetric synthetic aperture radar (PolSAR) has been widely used to penetrate sea ice surfaces to achieve fully polarimetric high-resolution imaging at all times of day and in a range of weather conditions. Model-based polarimetric decomposition is a powerful tool used to extract useful physical and geometric information about sea ice from the matrix datasets acquired by PolSAR. The volume scattering of sea ice is usually modeled as the incoherent average of scatterings of a large volume of oriented ellipsoid particles that are uniformly distributed in 3D space. This uniform spatial distribution is often approximated as a uniform orientation distribution (UOD), i.e., the particles are uniformly oriented in all directions. This is achieved in the existing literature by ensuring the canting angle φ and tilt angle τ of particles uniformly distributed in their respective ranges and introducing a factor cosτ in the ensemble average. However, we find this implementation of UOD is not always effective, while a real UOD can be realized by distributing the solid angles of particles uniformly in 3D space. By deriving the total solid angle of the canting-tilt cell spanned by particles and combining the differential relationship between solid angle and Euler angles φ and τ, a complete expression of the joint probability density function pφ,τ that can always ensure the uniform orientation of particles of sea ice is realized. By ensemble integrating the coherency matrix of φ,τ-oriented particle with pφ,τ, a generalized modeling of the volume coherency matrix of 3D uniformly oriented spheroid particles is obtained, which covers factors such as radar observation geometry, particle shape, canting geometry, tilt geometry and transmission effect in a multiplicative way. The existing volume scattering models of sea ice constitute special cases. The performance of the model in the characterization of the volume behaviors was investigated via simulations on a volume of oblate and prolate particles with the differential reflectivity ZDR, polarimetric entropy H and scattering α angle as descriptors. Based on the model, several interesting orientation geometries were also studied, including the aligned orientation, complement tilt geometry and reflection symmetry, among which the complement tilt geometry is specifically highlighted. It involves three volume models that correspond to the horizontal tilt, vertical tilt and random tilt of particles within sea ice, respectively. To match the models to PolSAR data for adaptive decomposition, two selection strategies are provided. One is based on ZDR, and the other is based on the maximum power fitting. The scattering power that reduces the rank of coherency matrix by exactly one without violating the physical realizability condition is obtained to make full use of the polarimetric scattering information. Both the models and decomposition were finally validated on the Gaofen-3 PolSAR data of a young ice area in Prydz Bay, Antarctica. The adaptive decomposition result demonstrates not only the dominant vertical tilt preference of brine inclusions within sea ice, but also the subordinate random tilt preference and non-negligible horizontal tilt preference, which are consistent with the geometric selection mechanism that the c-axes of polycrystallines within sea ice would gradually align with depth. The experiment also indicates that, compared to the strategy based on ZDR, the maximum power fitting is preferable because it is entirely driven by the model and data and is independent of any empirical thresholds. Such soft thresholding enables this strategy to adaptively estimate the negative ZDR offset introduced by the transmission effect, which provides a novel inversion of the refractive index of sea ice based on polarimetric model-based decomposition. Full article
Show Figures

Figure 1

Figure 1
<p>Volume scattering mechanism induced by the 3D-oriented brine inclusions within sea ice, and the Cartesian, polarization, ellipsoidal and solid angle geometry in modeling of the polarimetric backscattering and angular distribution of 3D-oriented particles.</p>
Full article ">Figure 2
<p>The complementary relationship among the geometries of horizontal tilt, vertical tilt and random tilt when (<bold>a</bold>) <inline-formula><mml:math id="mm954"><mml:semantics><mml:mrow><mml:msub><mml:mrow><mml:mi>τ</mml:mi></mml:mrow><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:msub><mml:mo>&gt;</mml:mo><mml:mn>0</mml:mn></mml:mrow></mml:semantics></mml:math></inline-formula>, (<bold>b</bold>) <inline-formula><mml:math id="mm427"><mml:semantics><mml:mrow><mml:msub><mml:mrow><mml:mi>τ</mml:mi></mml:mrow><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:msub><mml:mo>&lt;</mml:mo><mml:mn>0</mml:mn></mml:mrow></mml:semantics></mml:math></inline-formula> and (<bold>c</bold>) <inline-formula><mml:math id="mm428"><mml:semantics><mml:mrow><mml:msub><mml:mrow><mml:mi>τ</mml:mi></mml:mrow><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mn>0</mml:mn></mml:mrow></mml:semantics></mml:math></inline-formula>.</p>
Full article ">Figure 3
<p>Polarimetric descriptors, including (<bold>a</bold>,<bold>d</bold>,<bold>g</bold>,<bold>j</bold>,<bold>m</bold>) the differential reflectivity <inline-formula><mml:math id="mm955"><mml:semantics><mml:mrow><mml:msub><mml:mrow><mml:mi>Z</mml:mi></mml:mrow><mml:mrow><mml:mi>D</mml:mi><mml:mi>R</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:semantics></mml:math></inline-formula>, (<bold>b</bold>,<bold>e</bold>,<bold>h</bold>,<bold>k</bold>,<bold>n</bold>) the polarimetric entropy <inline-formula><mml:math id="mm531"><mml:semantics><mml:mrow><mml:mi>H</mml:mi></mml:mrow></mml:semantics></mml:math></inline-formula> and (<bold>c</bold>,<bold>f</bold>,<bold>i</bold>,<bold>l</bold>,<bold>o</bold>) the mean alpha angle α simulated on a volume of reflection-symmetric (i.e., the canting width <inline-formula><mml:math id="mm532"><mml:semantics><mml:mrow><mml:mo>∆</mml:mo><mml:mi>φ</mml:mi><mml:mo>=</mml:mo><mml:mi mathvariant="sans-serif">π</mml:mi></mml:mrow></mml:semantics></mml:math></inline-formula>) oblate particles (<inline-formula><mml:math id="mm533"><mml:semantics><mml:mrow><mml:msub><mml:mrow><mml:mi>ρ</mml:mi></mml:mrow><mml:mrow><mml:mi mathvariant="sans-serif">Σ</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mn>11</mml:mn><mml:mo>,</mml:mo><mml:msub><mml:mrow><mml:mi>ρ</mml:mi></mml:mrow><mml:mrow><mml:mi mathvariant="sans-serif">Δ</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mo>−</mml:mo><mml:mn>9</mml:mn></mml:mrow></mml:semantics></mml:math></inline-formula>) with incidence <inline-formula><mml:math id="mm534"><mml:semantics><mml:mrow><mml:mi>θ</mml:mi></mml:mrow></mml:semantics></mml:math></inline-formula> and tilt width <inline-formula><mml:math id="mm535"><mml:semantics><mml:mrow><mml:mo>∆</mml:mo><mml:mi>τ</mml:mi></mml:mrow></mml:semantics></mml:math></inline-formula> varying from <inline-formula><mml:math id="mm536"><mml:semantics><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:semantics></mml:math></inline-formula> to <inline-formula><mml:math id="mm537"><mml:semantics><mml:mrow><mml:mfrac><mml:mrow><mml:mi mathvariant="sans-serif">π</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:mfrac></mml:mrow></mml:semantics></mml:math></inline-formula>, and mean tilt <inline-formula><mml:math id="mm538"><mml:semantics><mml:mrow><mml:msub><mml:mrow><mml:mi>τ</mml:mi></mml:mrow><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:semantics></mml:math></inline-formula> taking the values of (<bold>a</bold>–<bold>c</bold>) <inline-formula><mml:math id="mm539"><mml:semantics><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:semantics></mml:math></inline-formula>, (<bold>d</bold>–<bold>f</bold>) <inline-formula><mml:math id="mm540"><mml:semantics><mml:mrow><mml:mfrac><mml:mrow><mml:mi mathvariant="sans-serif">π</mml:mi></mml:mrow><mml:mrow><mml:mn>6</mml:mn></mml:mrow></mml:mfrac></mml:mrow></mml:semantics></mml:math></inline-formula>, (<bold>g</bold>–<bold>i</bold>) <inline-formula><mml:math id="mm541"><mml:semantics><mml:mrow><mml:mfrac><mml:mrow><mml:mi mathvariant="sans-serif">π</mml:mi></mml:mrow><mml:mrow><mml:mn>4</mml:mn></mml:mrow></mml:mfrac></mml:mrow></mml:semantics></mml:math></inline-formula>, (<bold>j</bold>–<bold>l</bold>) <inline-formula><mml:math id="mm542"><mml:semantics><mml:mrow><mml:mfrac><mml:mrow><mml:mi mathvariant="sans-serif">π</mml:mi></mml:mrow><mml:mrow><mml:mn>3</mml:mn></mml:mrow></mml:mfrac></mml:mrow></mml:semantics></mml:math></inline-formula> and (<bold>m</bold>–<bold>o</bold>) <inline-formula><mml:math id="mm543"><mml:semantics><mml:mrow><mml:mfrac><mml:mrow><mml:mi mathvariant="sans-serif">π</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:mfrac></mml:mrow></mml:semantics></mml:math></inline-formula>, respectively.</p>
Full article ">Figure 4
<p>Polarimetric descriptors, including (<bold>a</bold>,<bold>d</bold>,<bold>g</bold>,<bold>j</bold>,<bold>m</bold>) the differential reflectivity <inline-formula><mml:math id="mm956"><mml:semantics><mml:mrow><mml:msub><mml:mrow><mml:mi>Z</mml:mi></mml:mrow><mml:mrow><mml:mi>D</mml:mi><mml:mi>R</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:semantics></mml:math></inline-formula>, (<bold>b</bold>,<bold>e</bold>,<bold>h</bold>,<bold>k</bold>,<bold>n</bold>) the polarimetric entropy <inline-formula><mml:math id="mm544"><mml:semantics><mml:mrow><mml:mi>H</mml:mi></mml:mrow></mml:semantics></mml:math></inline-formula> and (<bold>c</bold>,<bold>f</bold>,<bold>i</bold>,<bold>l</bold>,<bold>o</bold>) the mean alpha angle α simulated on a volume of reflection-symmetric (i.e., the canting width <inline-formula><mml:math id="mm545"><mml:semantics><mml:mrow><mml:mo>∆</mml:mo><mml:mi>φ</mml:mi><mml:mo>=</mml:mo><mml:mi mathvariant="sans-serif">π</mml:mi></mml:mrow></mml:semantics></mml:math></inline-formula>) prolate particles (<inline-formula><mml:math id="mm546"><mml:semantics><mml:mrow><mml:msub><mml:mrow><mml:mi>ρ</mml:mi></mml:mrow><mml:mrow><mml:mi mathvariant="sans-serif">Σ</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mn>11</mml:mn><mml:mo>,</mml:mo><mml:msub><mml:mrow><mml:mi>ρ</mml:mi></mml:mrow><mml:mrow><mml:mi mathvariant="sans-serif">Δ</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mn>9</mml:mn></mml:mrow></mml:semantics></mml:math></inline-formula>) with incidence <inline-formula><mml:math id="mm547"><mml:semantics><mml:mrow><mml:mi>θ</mml:mi></mml:mrow></mml:semantics></mml:math></inline-formula> and tilt width <inline-formula><mml:math id="mm548"><mml:semantics><mml:mrow><mml:mo>∆</mml:mo><mml:mi>τ</mml:mi></mml:mrow></mml:semantics></mml:math></inline-formula> ranging from <inline-formula><mml:math id="mm549"><mml:semantics><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:semantics></mml:math></inline-formula> to <inline-formula><mml:math id="mm550"><mml:semantics><mml:mrow><mml:mfrac><mml:mrow><mml:mi mathvariant="sans-serif">π</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:mfrac></mml:mrow></mml:semantics></mml:math></inline-formula>, and mean tilt <inline-formula><mml:math id="mm551"><mml:semantics><mml:mrow><mml:msub><mml:mrow><mml:mi>τ</mml:mi></mml:mrow><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:semantics></mml:math></inline-formula> taking the values of (<bold>a</bold>–<bold>c</bold>) <inline-formula><mml:math id="mm552"><mml:semantics><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:semantics></mml:math></inline-formula>, (<bold>d</bold>–<bold>f</bold>) <inline-formula><mml:math id="mm553"><mml:semantics><mml:mrow><mml:mfrac><mml:mrow><mml:mi mathvariant="sans-serif">π</mml:mi></mml:mrow><mml:mrow><mml:mn>6</mml:mn></mml:mrow></mml:mfrac></mml:mrow></mml:semantics></mml:math></inline-formula>, (<bold>g</bold>–<bold>i</bold>) <inline-formula><mml:math id="mm554"><mml:semantics><mml:mrow><mml:mfrac><mml:mrow><mml:mi mathvariant="sans-serif">π</mml:mi></mml:mrow><mml:mrow><mml:mn>4</mml:mn></mml:mrow></mml:mfrac></mml:mrow></mml:semantics></mml:math></inline-formula>, (<bold>j</bold>–<bold>l</bold>) <inline-formula><mml:math id="mm555"><mml:semantics><mml:mrow><mml:mfrac><mml:mrow><mml:mi mathvariant="sans-serif">π</mml:mi></mml:mrow><mml:mrow><mml:mn>3</mml:mn></mml:mrow></mml:mfrac></mml:mrow></mml:semantics></mml:math></inline-formula> and (<bold>m</bold>–<bold>o</bold>) <inline-formula><mml:math id="mm556"><mml:semantics><mml:mrow><mml:mfrac><mml:mrow><mml:mi mathvariant="sans-serif">π</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:mfrac></mml:mrow></mml:semantics></mml:math></inline-formula>, respectively.</p>
Full article ">Figure 5
<p>Tilt geometries of particles within the sea ice. (<bold>a</bold>) Particles prefer horizontal tilt; (<bold>b</bold>) particles prefer vertical tilt; (<bold>c</bold>) particles show no preference but random tilt.</p>
Full article ">Figure 6
<p>Experimental sea ice dataset and information. (<bold>a</bold>) Pauli RGB composite image; (<bold>b</bold>) ice atlas chart of Antarctica produced by the National Ice Center on 2 May 2019. The red star indicates the location of the study area.</p>
Full article ">Figure 7
<p>Volume powers extracted by (<bold>a</bold>) the maximum power fitting <inline-formula><mml:math id="mm957"><mml:semantics><mml:mrow><mml:mfenced open="[" close="]" separators="|"><mml:mrow><mml:msubsup><mml:mrow><mml:mi>P</mml:mi></mml:mrow><mml:mrow><mml:mi>V</mml:mi></mml:mrow><mml:mrow><mml:mi mathvariant="normal">m</mml:mi><mml:mi mathvariant="normal">a</mml:mi><mml:mi mathvariant="normal">x</mml:mi></mml:mrow></mml:msubsup></mml:mrow></mml:mfenced></mml:mrow></mml:semantics></mml:math></inline-formula> and (<bold>b</bold>) <inline-formula><mml:math id="mm886"><mml:semantics><mml:mrow><mml:msub><mml:mrow><mml:mi>Z</mml:mi></mml:mrow><mml:mrow><mml:mi>D</mml:mi><mml:mi>R</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:semantics></mml:math></inline-formula>-based strategy <inline-formula><mml:math id="mm887"><mml:semantics><mml:mrow><mml:mfenced open="[" close="]" separators="|"><mml:mrow><mml:msubsup><mml:mrow><mml:mi>P</mml:mi></mml:mrow><mml:mrow><mml:mi>V</mml:mi></mml:mrow><mml:mrow><mml:mi mathvariant="normal">z</mml:mi><mml:mi mathvariant="normal">d</mml:mi><mml:mi mathvariant="normal">r</mml:mi></mml:mrow></mml:msubsup></mml:mrow></mml:mfenced></mml:mrow></mml:semantics></mml:math></inline-formula>, as well as (<bold>c</bold>) their difference <inline-formula><mml:math id="mm888"><mml:semantics><mml:mrow><mml:mfenced open="[" close="]" separators="|"><mml:mrow><mml:msubsup><mml:mrow><mml:mi>P</mml:mi></mml:mrow><mml:mrow><mml:mi>V</mml:mi></mml:mrow><mml:mrow><mml:mo>∆</mml:mo></mml:mrow></mml:msubsup></mml:mrow></mml:mfenced></mml:mrow></mml:semantics></mml:math></inline-formula> and the (<bold>a</bold>,<bold>d</bold>,<bold>g</bold>) horizontal component, (<bold>b</bold>,<bold>e</bold>,<bold>h</bold>) vertical component and (<bold>c</bold>,<bold>f</bold>,<bold>i</bold>) random component of (<bold>d</bold>–<bold>f</bold>) <inline-formula><mml:math id="mm889"><mml:semantics><mml:mrow><mml:mfenced open="[" close="]" separators="|"><mml:mrow><mml:msubsup><mml:mrow><mml:mi>P</mml:mi></mml:mrow><mml:mrow><mml:mi>V</mml:mi></mml:mrow><mml:mrow><mml:mi mathvariant="normal">m</mml:mi><mml:mi mathvariant="normal">a</mml:mi><mml:mi mathvariant="normal">x</mml:mi></mml:mrow></mml:msubsup></mml:mrow></mml:mfenced></mml:mrow></mml:semantics></mml:math></inline-formula> and (<bold>g</bold>–<bold>i</bold>) <inline-formula><mml:math id="mm890"><mml:semantics><mml:mrow><mml:mfenced open="[" close="]" separators="|"><mml:mrow><mml:msubsup><mml:mrow><mml:mi>P</mml:mi></mml:mrow><mml:mrow><mml:mi>V</mml:mi></mml:mrow><mml:mrow><mml:mi mathvariant="normal">z</mml:mi><mml:mi mathvariant="normal">d</mml:mi><mml:mi mathvariant="normal">r</mml:mi></mml:mrow></mml:msubsup></mml:mrow></mml:mfenced></mml:mrow></mml:semantics></mml:math></inline-formula>.</p>
Full article ">Figure 8
<p>Quantitative statistics of the selection of the three tilt scattering models in (<bold>a</bold>) the maximum power fitting and (<bold>b</bold>) <inline-formula><mml:math id="mm958"><mml:semantics><mml:mrow><mml:msub><mml:mrow><mml:mi>Z</mml:mi></mml:mrow><mml:mrow><mml:mi>D</mml:mi><mml:mi>R</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:semantics></mml:math></inline-formula>-based strategy.</p>
Full article ">Figure 9
<p><inline-formula><mml:math id="mm959"><mml:semantics><mml:mrow><mml:msub><mml:mrow><mml:mi>Z</mml:mi></mml:mrow><mml:mrow><mml:mi>D</mml:mi><mml:mi>R</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:semantics></mml:math></inline-formula> histogram of pixels attributed to the three tilt geometries in (<bold>a</bold>) the maximum power fitting and (<bold>b</bold>) <inline-formula><mml:math id="mm925"><mml:semantics><mml:mrow><mml:msub><mml:mrow><mml:mi>Z</mml:mi></mml:mrow><mml:mrow><mml:mi>D</mml:mi><mml:mi>R</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:semantics></mml:math></inline-formula>-based strategy.</p>
Full article ">
25 pages, 12883 KiB  
Article
Remote Sensing Monitoring and Analysis of Spatiotemporal Changes in China’s Anthropogenic Carbon Emissions Based on XCO2 Data
by Yanjun Wang, Mengjie Wang, Fei Teng and Yiye Ji
Remote Sens. 2023, 15(12), 3207; https://doi.org/10.3390/rs15123207 - 20 Jun 2023
Cited by 9 | Viewed by 2248
Abstract
The monitoring and analysis of the spatiotemporal distribution of anthropogenic carbon emissions is an important part of realizing China’s regional “dual carbon” goals; that is, the aim is for carbon emissions to peak in 2030 an to achieve carbon neutrality by 2060, as [...] Read more.
The monitoring and analysis of the spatiotemporal distribution of anthropogenic carbon emissions is an important part of realizing China’s regional “dual carbon” goals; that is, the aim is for carbon emissions to peak in 2030 an to achieve carbon neutrality by 2060, as well as achieving sustainable development of the ecological environment. The column-averaged CO2 dry air mole fraction (XCO2) of greenhouse gas remote sensing satellites has been widely used to monitor anthropogenic carbon emissions. However, selecting a reasonable background region to eliminate the influence of uncertainty factors is still an important challenge to monitor anthropogenic carbon emissions by using XCO2. Aiming at the problems of the imprecise selection of background regions, this study proposes to enhance the anthropogenic carbon emission signal in the XCO2 by using the regional comparison method based on the idea of zoning. First, this study determines the background region based on the Open-Data Inventory for Anthropogenic Carbon dioxide (ODIAC) dataset and potential temperature data. Second, the average value of the XCO2 in the background area was extracted and taken as the XCO2 background. On this basis, the XCO2 anomaly (XCO2ano) was obtained by regional comparison method. Finally, the spatiotemporal variation characteristics and trends of XCO2ano were analyzed, and the correlations between the number of residential areas and fossil fuel emissions were calculated. The results of the satellite observation data experiments over China from 2010 to 2020 show that the XCO2ano and anthropogenic carbon emissions have similar spatial distribution patterns. The XCO2ano in China changed significantly and was in a positive growth trend as a whole. The XCO2ano values have a certain positive correlation with the number of residential areas and observations of fossil fuel emissions. The purpose of this research is to enhance the anthropogenic carbon emission signals in satellite observation XCO2 data by combining ODIAC data and potential temperature data, achieve the remote sensing monitoring and analysis of spatiotemporal changes in anthropogenic carbon emissions over China, and provide technical support for the policies and paths of regional carbon emission reductions and ecological environmental protection. Full article
Show Figures

Figure 1

Figure 1
<p>Partitioning based on potential temperature data. Beijing–Tianjin–Hebei (BTH); Pearl River Delta (PRD); Yangtze River Delta (YRD).</p>
Full article ">Figure 2
<p>Research flow chart.</p>
Full article ">Figure 3
<p>Time series of the XCO<sub>2</sub>. (<b>a</b>) is the time series of XCO<sub>2</sub> in the background area and emission area of China, and (<b>b</b>) is the time series of XCO<sub>2</sub> in the background area and emission area of China excluding the Tibetan Plateau with an altitude higher than 3000 m. X-axis: January 2010 represents 1, increasing by 1 for every month of growth, and finally December 2020 represents 132.</p>
Full article ">Figure 4
<p>XCO<sub>2ano</sub> and ODIAC data fossil fuel emissions spatial distribution and land use data. (<b>a</b>) is the average value of the XCO<sub>2ano</sub> from 2010 to 2020, (<b>b</b>) is the average value of ODIAC fossil fuel emissions from 2010 to 2019, and (<b>c</b>) is the 300 m spatial resolution land use data.</p>
Full article ">Figure 5
<p>Spatial distribution of average XCO<sub>2ano</sub>. (<b>a</b>–<b>c</b>) are the spatial distribution of average XCO<sub>2ano</sub> in China in 2010, 2015 and 2020, respectively.</p>
Full article ">Figure 6
<p>Spatial distribution of the average XCO<sub>2ano</sub>. (<b>a</b>–<b>d</b>) are the average values of XCO<sub>2ano</sub> in spring, summer, autumn, and winter from 2010 to 2020, respectively.</p>
Full article ">Figure 7
<p>Seasonal variation in the mean value of the XCO<sub>2ano</sub> in five regions and China.</p>
Full article ">Figure 8
<p>Spatial distribution of the CV of the XCO<sub>2ano</sub> in China from 2010 to 2020.</p>
Full article ">Figure 9
<p>Spatial distribution of the SKEW of XCO<sub>2ano</sub> in China from 2010 to 2020.</p>
Full article ">Figure 10
<p>The correlation analysis between the average XCO<sub>2ano</sub> and the number of residential areas, in which (<b>a</b>) is the correlation analysis result of all study areas and (<b>b</b>) is the correlation analysis result of Area III.</p>
Full article ">Figure 11
<p>Correlation analysis between the average XCO<sub>2ano</sub> and fossil fuel emissions, in which (<b>a</b>) is the correlation analysis result of all study areas and (<b>b</b>) is the correlation analysis result of Area V.</p>
Full article ">Figure 12
<p>XCO<sub>2ano</sub> distribution and average fossil fuel carbon emissions for different background region selection methods. (<b>a</b>) is the XCO<sub>2ano</sub> obtained with the Chinese median XCO<sub>2</sub> is taken as the background value, (<b>b</b>) is the XCO<sub>2ano</sub> obtained with the regional median XCO<sub>2</sub> as the background value.</p>
Full article ">Figure 13
<p>Correlation analysis of XCO<sub>2ano</sub> and fossil fuel emissions. (<b>a</b>,<b>b</b>) are the correlations between the XCO<sub>2ano</sub> and fossil fuel emissions using the Chinese median XCO<sub>2</sub> and the regional median XCO<sub>2ano</sub> as the background value, respectively.</p>
Full article ">Figure A1
<p>Spatial distribution of monthly average XCO<sub>2ano</sub> values from 2010 to 2020 with the background without the high altitudes.</p>
Full article ">Figure A2
<p>Spatial distribution of XCO<sub>2ano</sub> in China from 2011 to 2019 (<b>a</b>–<b>h</b>). The results showed that the spatial distribution pattern of XCO<sub>2ano</sub> in China in 2011 and 2012 was similar to that in 2010, and the high-value area was mainly distributed in southern and northwestern China. Since 2013, the high-value area of China XCO<sub>2ano</sub> is mainly located in eastern China. At the same time, there are inter-annual differences. The high value of XCO<sub>2ano</sub> in 2012 and 2014 is significantly less than that in other years.</p>
Full article ">Figure A2 Cont.
<p>Spatial distribution of XCO<sub>2ano</sub> in China from 2011 to 2019 (<b>a</b>–<b>h</b>). The results showed that the spatial distribution pattern of XCO<sub>2ano</sub> in China in 2011 and 2012 was similar to that in 2010, and the high-value area was mainly distributed in southern and northwestern China. Since 2013, the high-value area of China XCO<sub>2ano</sub> is mainly located in eastern China. At the same time, there are inter-annual differences. The high value of XCO<sub>2ano</sub> in 2012 and 2014 is significantly less than that in other years.</p>
Full article ">Figure A3
<p>Spatial distribution of XCO<sub>2ano</sub> in China from 2010 to 2020 (<b>a</b>–<b>k</b>) with the background without the high altitudes. The spatial distribution pattern of XCO<sub>2ano</sub> in China from 2010 to 2020 is generally similar, showing a distribution pattern of high in the east and low in the west.</p>
Full article ">Figure A3 Cont.
<p>Spatial distribution of XCO<sub>2ano</sub> in China from 2010 to 2020 (<b>a</b>–<b>k</b>) with the background without the high altitudes. The spatial distribution pattern of XCO<sub>2ano</sub> in China from 2010 to 2020 is generally similar, showing a distribution pattern of high in the east and low in the west.</p>
Full article ">Figure A3 Cont.
<p>Spatial distribution of XCO<sub>2ano</sub> in China from 2010 to 2020 (<b>a</b>–<b>k</b>) with the background without the high altitudes. The spatial distribution pattern of XCO<sub>2ano</sub> in China from 2010 to 2020 is generally similar, showing a distribution pattern of high in the east and low in the west.</p>
Full article ">Figure A4
<p>Spatial distribution of the average XCO<sub>2ano</sub> in four seasons with the background without the high altitudes. The XCO<sub>2ano</sub> distribution patterns were similar in spring (<b>a</b>) and autumn (<b>c</b>), with XCO<sub>2ano</sub> significantly lower in summer (<b>b</b>) than in the remaining three seasons, and higher in winter (<b>d</b>).</p>
Full article ">Figure A5
<p>Spatial distribution of the CV of the XCO<sub>2ano</sub> in China from 2010 to 2020 with the background without the high altitudes.</p>
Full article ">Figure A6
<p>Spatial distribution of the SKEW of XCO<sub>2ano</sub> in China from 2010 to 2020 with the background without the high altitudes.</p>
Full article ">Figure A7
<p>The correlation analysis between average XCO<sub>2ano</sub> and the number of residential areas in different areas. There is a high correlation between the number of residential area and the average XCO<sub>2ano</sub> in areas I (<b>a</b>), II (<b>b</b>), and IV (<b>c</b>), while the lowest correlation between the two is found in area V (<b>d</b>).</p>
Full article ">Figure A8
<p>The correlation analysis between XCO<sub>2ano</sub> and fossil fuel emissions in different area. There is a high correlation between the FFCO<sub>2</sub> emissions and the average XCO<sub>2ano</sub> in areas II (<b>b</b>), III (<b>c</b>), and IV (<b>d</b>), while the lowest correlation between the two is found in area I (<b>a</b>).</p>
Full article ">Figure A8 Cont.
<p>The correlation analysis between XCO<sub>2ano</sub> and fossil fuel emissions in different area. There is a high correlation between the FFCO<sub>2</sub> emissions and the average XCO<sub>2ano</sub> in areas II (<b>b</b>), III (<b>c</b>), and IV (<b>d</b>), while the lowest correlation between the two is found in area I (<b>a</b>).</p>
Full article ">
20 pages, 7456 KiB  
Article
Unsupervised Machine Learning for GNSS Reflectometry Inland Water Body Detection
by Stylianos Kossieris, Milad Asgarimehr and Jens Wickert
Remote Sens. 2023, 15(12), 3206; https://doi.org/10.3390/rs15123206 - 20 Jun 2023
Cited by 5 | Viewed by 2118
Abstract
Inland water bodies, wetlands and their dynamics have a key role in a variety of scientific, economic, and social applications. They are significant in identifying climate change, water resource management, agricultural productivity, and the modeling of land–atmosphere exchange. Changes in the extent and [...] Read more.
Inland water bodies, wetlands and their dynamics have a key role in a variety of scientific, economic, and social applications. They are significant in identifying climate change, water resource management, agricultural productivity, and the modeling of land–atmosphere exchange. Changes in the extent and position of water bodies are crucial to the ecosystems. Mapping water bodies at a global scale is a challenging task due to the global variety of terrains and water surface. However, the sensitivity of spaceborne Global Navigation Satellite System Reflectometry (GNSS-R) to different land surface properties offers the potential to detect and monitor inland water bodies. The extensive dataset available in the Cyclone Global Navigation Satellite System (CYGNSS), launched in December 2016, is used in our investigation. Although the main mission of CYGNSS was to measure the ocean wind speed in hurricanes and tropical cyclones, we show its capability of detecting and mapping inland water bodies. Both bistatic radar cross section (BRCS) and signal-to-noise ratio (SNR) can be used to detect, identify, and map the changes in the position and extent of inland waterbodies. We exploit the potential of unsupervised machine learning algorithms, more specifically the clustering methods, K-Means, Agglomerative, and Density-based Spatial Clustering of Applications with Noise (DBSCAN), for the detection of inland waterbodies. The results are evaluated based on the Copernicus land cover classification gridded maps, at 300 m spatial resolution. The outcomes demonstrate that CYGNSS data can identify and monitor inland waterbodies and their tributaries at high temporal resolution. K-Means has the highest Accuracy (93.5%) compared to the DBSCAN (90.3%) and Agglomerative (91.6%) methods. However, the DBSCAN method has the highest Recall (83.1%) as compared to Agglomerative (82.7%) and K-Means (79.2%). The current study offers valuable insights and analysis for further investigations in the field of GNSS-R and machine learning. Full article
(This article belongs to the Special Issue Applications of GNSS Reflectometry for Earth Observation III)
Show Figures

Figure 1

Figure 1
<p>GNSS Reflectometry using receivers on satellites (image credit: GGOS, <a href="https://ggos.org/item/gnss-reflectometry/" target="_blank">https://ggos.org/item/gnss-reflectometry/</a>, accessed on 21 May 2023).</p>
Full article ">Figure 2
<p>(<b>a</b>) Track over river Rio Negro and (<b>b</b>) the maximum values of DDMs bin BRCS of the track.</p>
Full article ">Figure 3
<p>(<b>a</b>) Track over river Rio Uapes and (<b>b</b>) the maximum values of DDMs bin BRCS of the track.</p>
Full article ">Figure 4
<p>DDM SNR measurements over river Rio Negro track.</p>
Full article ">Figure 5
<p>DDM SNR measurements over river Rio Uapes track.</p>
Full article ">Figure 6
<p>(<b>a</b>) Track over the main part of the Amazon river and (<b>b</b>) the DDM SNR measurements of the track.</p>
Full article ">Figure 7
<p>(<b>a</b>) Track over the Congo river and (<b>b</b>) the DDM SNR measurements of the track.</p>
Full article ">Figure 8
<p>Evaluation results using outcomes of K-Means clustering as predicted values.</p>
Full article ">Figure 9
<p>Confusion matrix using K-Means and false positive classes.</p>
Full article ">Figure 10
<p>Small tributaries detected by CYGNSS observations but not by Copernicus land cover mask.</p>
Full article ">Figure 11
<p>Confusion matrix using Agglomerative and false positive classes.</p>
Full article ">Figure 12
<p>Confusion matrix using DBSCAN and false positive classes.</p>
Full article ">Figure 13
<p>Map over the Congo basin using 1-month observations (1–31 August 2018).</p>
Full article ">Figure 14
<p>False clustering because of track-based SNR measurements.</p>
Full article ">Figure 15
<p>SNR values not related to surface properties (Track 1056, on 3 August 2018).</p>
Full article ">Figure 16
<p>False clustering over large inland water bodies, such as Tanganyika.</p>
Full article ">
21 pages, 10793 KiB  
Article
Augmented GBM Nonlinear Model to Address Spectral Variability for Hyperspectral Unmixing
by Linghong Meng, Danfeng Liu, Liguo Wang, Jón Atli Benediktsson, Xiaohan Yue and Yuetao Pan
Remote Sens. 2023, 15(12), 3205; https://doi.org/10.3390/rs15123205 - 20 Jun 2023
Cited by 2 | Viewed by 1993
Abstract
Spectral unmixing (SU) is a significant preprocessing task for handling hyperspectral images (HSI), but its process is affected by nonlinearity and spectral variability (SV). Currently, SV is considered within the framework of linear mixing models (LMM), which ignores the nonlinear effects in the [...] Read more.
Spectral unmixing (SU) is a significant preprocessing task for handling hyperspectral images (HSI), but its process is affected by nonlinearity and spectral variability (SV). Currently, SV is considered within the framework of linear mixing models (LMM), which ignores the nonlinear effects in the scene. To address that issue, we consider the effects of SV on SU while investigating the nonlinear effects of hyperspectral images. Furthermore, an augmented generalized bilinear model is proposed to address spectral variability (abbreviated AGBM-SV). First, AGBM-SV adopts a generalized bilinear model (GBM) as the basic framework to address the nonlinear effects caused by second-order scattering. Secondly, scaling factors and spectral variability dictionaries are introduced to model the variability issues caused by the illumination conditions, material intrinsic variability, and other environmental factors. Then, a data-driven learning strategy is employed to set sparse and orthogonal bases for the abundance and spectral variability dictionaries according to the distribution characteristics of real materials. Finally, the alternating direction method of multipliers (ADMM) optimization method is used to split and solve the objective function, enabling the AGBM-SV algorithm to estimate the abundance and learn the spectral variability dictionary more effectively. The experimental results demonstrate the comparative superiority of the AGBM-SV method in both qualitative and quantitative perspectives, which can effectively solve the problem of spectral variability in nonlinear mixing scenes and to improve unmixing accuracy. Full article
(This article belongs to the Special Issue Self-Supervised Learning in Remote Sensing)
Show Figures

Figure 1

Figure 1
<p>The holistic diagram of the proposed AGBM-SV.</p>
Full article ">Figure 2
<p>Convergence analyses of AGBM-SV were experimentally performed on the Synthetic dataset, the Urban dataset and the Cuprite dataset.</p>
Full article ">Figure 3
<p>(<b>a</b>) A false color representation of the synthetic dataset. (<b>b</b>) The five endmembers used for data simulation.</p>
Full article ">Figure 4
<p>Sensitivity analysis of 5 regularization parameters using the proposed AGBM-SV algorithm on synthetic dataset (<math display="inline"><semantics> <mrow> <mi>α</mi> <mo>,</mo> <mi>β</mi> <mo>,</mo> <mi>γ</mi> <mo>,</mo> <mi>η</mi> </mrow> </semantics></math> and the number of basis vectors (D) of <math display="inline"><semantics> <mstyle mathvariant="bold" mathsize="normal"> <mi>W</mi> </mstyle> </semantics></math>).</p>
Full article ">Figure 5
<p>Visualization of the aRMSE evaluation metric between different algorithms for the Synthetic dataset and Urban dataset, and the OA evaluation metric for the Cuprite dataset. (<b>a</b>) Synthetic; (<b>b</b>) Urban; (<b>c</b>) Cuprite.</p>
Full article ">Figure 6
<p>Reference abundance maps and the estimated abundance maps obtained by using all the tested methods for the synthetic dataset. The first row represents the classification results obtained by different methods based on estimated abundance values.</p>
Full article ">Figure 7
<p>(<b>a</b>) A false color representation of the urban dataset. (<b>b</b>) Four endmembers extracted by VCA in spectral unmixing.</p>
Full article ">Figure 8
<p>Reference abundance maps and the estimated abundance maps obtained by using all of the tested methods for the Urban dataset.</p>
Full article ">Figure 9
<p>(<b>a</b>) A false color image display of the cuprite dataset. (<b>b</b>) The fourteen endmembers extracted by VCA in spectral unmixing.</p>
Full article ">Figure 10
<p>The abundance maps comparison between the proposed method and the state-of-the-art methods. The first row represents the classification results obtained by different methods based on estimated abundance values.</p>
Full article ">
11 pages, 2620 KiB  
Communication
Enhanced Radar Detection in the Presence of Specular Reflection Using a Single Transmitting Antenna and Three Receiving Antennas
by Yong Yang and Xue-Song Wang
Remote Sens. 2023, 15(12), 3204; https://doi.org/10.3390/rs15123204 - 20 Jun 2023
Cited by 2 | Viewed by 1379
Abstract
Radar target echoes undergo fading in the presence of specular reflection, which is adverse to radar detection. To address this problem, this paper proposes a radar detection method that uses a single transmitting antenna and three receiving antennas. The proposed method uses the [...] Read more.
Radar target echoes undergo fading in the presence of specular reflection, which is adverse to radar detection. To address this problem, this paper proposes a radar detection method that uses a single transmitting antenna and three receiving antennas. The proposed method uses the maximum absolute value of the difference in the radar received signal power among the three receiving antennas as the test statistic. First, the target echo in the presence of specular reflection is analyzed. Then, selection of the required number of radar antennas and the heights at which they must be situated are discussed. Subsequently, analytical expressions of the radar detection probability and the false alarm probability are derived. Finally, simulation results are presented, which show that the proposed method improves radar detection performance in the presence of specular reflection. Full article
(This article belongs to the Special Issue Theory and Applications of MIMO Radar)
Show Figures

Figure 1

Figure 1
<p>Schematic diagram of radar specular reflection.</p>
Full article ">Figure 2
<p><math display="inline"><semantics> <mrow> <msub> <mi>ϕ</mi> <mi>l</mi> </msub> </mrow> </semantics></math> versus radar antenna height and target location.</p>
Full article ">Figure 3
<p>Received target echo powers received by two antennas with different heights.</p>
Full article ">Figure 4
<p>Target echo powers received by three antennas with different heights.</p>
Full article ">Figure 5
<p>Schematic diagram of the configuration of the three antennas.</p>
Full article ">Figure 6
<p>False alarm probability versus detection threshold.</p>
Full article ">Figure 7
<p>Theoretical and simulated detection probabilities for radar using a single transmitting antenna and three receiving antennas.</p>
Full article ">Figure 8
<p>Radar detection probabilities under various <math display="inline"><semantics> <mrow> <msub> <mi>ϕ</mi> <mi>l</mi> </msub> </mrow> </semantics></math>.</p>
Full article ">
21 pages, 11097 KiB  
Article
Implementing Cloud Computing for the Digital Mapping of Agricultural Soil Properties from High Resolution UAV Multispectral Imagery
by Samuel Pizarro, Narcisa G. Pricope, Deyanira Figueroa, Carlos Carbajal, Miriam Quispe, Jesús Vera, Lidiana Alejandro, Lino Achallma, Izamar Gonzalez, Wilian Salazar, Hildo Loayza, Juancarlos Cruz and Carlos I. Arbizu
Remote Sens. 2023, 15(12), 3203; https://doi.org/10.3390/rs15123203 - 20 Jun 2023
Cited by 5 | Viewed by 4669
Abstract
The spatial heterogeneity of soil properties has a significant impact on crop growth, making it difficult to adopt site-specific crop management practices. Traditional laboratory-based analyses are costly, and data extrapolation for mapping soil properties using high-resolution imagery becomes a computationally expensive procedure, taking [...] Read more.
The spatial heterogeneity of soil properties has a significant impact on crop growth, making it difficult to adopt site-specific crop management practices. Traditional laboratory-based analyses are costly, and data extrapolation for mapping soil properties using high-resolution imagery becomes a computationally expensive procedure, taking days or weeks to obtain accurate results using a desktop workstation. To overcome these challenges, cloud-based solutions such as Google Earth Engine (GEE) have been used to analyze complex data with machine learning algorithms. In this study, we explored the feasibility of designing and implementing a digital soil mapping approach in the GEE platform using high-resolution reflectance imagery derived from a thermal infrared and multispectral camera Altum (MicaSense, Seattle, WA, USA). We compared a suite of multispectral-derived soil and vegetation indices with in situ measurements of physical-chemical soil properties in agricultural lands in the Peruvian Mantaro Valley. The prediction ability of several machine learning algorithms (CART, XGBoost, and Random Forest) was evaluated using R2, to select the best predicted maps (R2 > 0.80), for ten soil properties, including Lime, Clay, Sand, N, P, K, OM, Al, EC, and pH, using multispectral imagery and derived products such as spectral indices and a digital surface model (DSM). Our results indicate that the predictions based on spectral indices, most notably, SRI, GNDWI, NDWI, and ExG, in combination with CART and RF algorithms are superior to those based on individual spectral bands. Additionally, the DSM improves the model prediction accuracy, especially for K and Al. We demonstrate that high-resolution multispectral imagery processed in the GEE platform has the potential to develop soil properties prediction models essential in establishing adaptive soil monitoring programs for agricultural regions. Full article
Show Figures

Figure 1

Figure 1
<p>Location of the study area, Santa Ana, Junin (Peru).</p>
Full article ">Figure 2
<p>Representation of the methodological framework used in this study.</p>
Full article ">Figure 3
<p>(<b>a</b>) Matrice 300 UAV integrated with the Altum sensor serving as the imaging platform used in this study, (<b>b</b>) Altum camera, (<b>c</b>) flight plan for the study image, and (<b>d</b>) GCP.</p>
Full article ">Figure 4
<p>Bagging/random forest example (<b>left</b>), boosting/XG-boost example (<b>right</b>).</p>
Full article ">Figure 5
<p>Correlation coefficients between measured soil physical chemical properties and predictors. r—Pearson’s correlation coefficient; Significant at 5% probability; X—non-significant.</p>
Full article ">Figure 6
<p>Prediction maps for Lime and Clay (<b>left</b>) and relative importance of predictors (<b>right</b>).</p>
Full article ">Figure 7
<p>Prediction maps for sand, N, P and K (<b>left</b>) and relative importance of predictors (<b>right</b>).</p>
Full article ">Figure 8
<p>Prediction maps for OM, Al, EC and pH (<b>left</b>) and relative importance of predictors (<b>right</b>).</p>
Full article ">Figure 8 Cont.
<p>Prediction maps for OM, Al, EC and pH (<b>left</b>) and relative importance of predictors (<b>right</b>).</p>
Full article ">
20 pages, 7275 KiB  
Article
Suppression of Mainlobe Jammers with Quadratic Element Pulse Coding in MIMO Radar
by Yiqun Zhang, Guisheng Liao, Lan Lan, Jingwei Xu and Xuepan Zhang
Remote Sens. 2023, 15(12), 3202; https://doi.org/10.3390/rs15123202 - 20 Jun 2023
Cited by 4 | Viewed by 1261
Abstract
The problem of suppressing mainlobe deceptive jammers, which spoof radar systems by generating multiple false targets, has attracted widespread attention. To tackle this problem, in this paper, the multiple-input multiple-output (MIMO) radar system was utilized by applying a quadratic element phase code (QEPC) [...] Read more.
The problem of suppressing mainlobe deceptive jammers, which spoof radar systems by generating multiple false targets, has attracted widespread attention. To tackle this problem, in this paper, the multiple-input multiple-output (MIMO) radar system was utilized by applying a quadratic element phase code (QEPC) to the transmitted pulses of different elements. In the receiver, by utilizing the spatial frequency and Doppler frequency offset generated after decoding, the jammers were equivalently distributed in the sidelobes of the joint Doppler-transmit-receive domain and were distinguishable from the true target. Then, further spatial frequency compensation and Doppler compensation were performed to align the true target to the zero point in the transmit spatial and Doppler domains. Moreover, by designing appropriate coding coefficients, the jammers were suppressed by data-independent Doppler-transmit-receive three-dimensional beamforming. However, the beamforming performance was sensitive to angular estimation mismatches, resulting in performance degradation of jammer suppression. To this end, a center-boundary null-broadening control (CBNBC) approach was used to broaden the nulls in the equivalent beampattern by generating multiple artificial jammers with preset powers around the nulls. Thus, the false targets (FTs) with deviations were sufficiently suppressed in the broadened notches. Numerical simulations and theoretical analysis demonstrated the performance of the developed jammer suppression method. Full article
Show Figures

Figure 1

Figure 1
<p>Signal model in the QEPC-MIMO.</p>
Full article ">Figure 2
<p>Receive signal processing procedures.</p>
Full article ">Figure 3
<p>Scenario for airborne self-protection jamming.</p>
Full article ">Figure 4
<p>Illumination of true and FTs with the QEPC-MIMO.</p>
Full article ">Figure 5
<p>Distribution of the FTs and the true target in Doppler-transmit-receive space.</p>
Full article ">Figure 6
<p>Processing flow chart.</p>
Full article ">Figure 7
<p>The square null region with artificial jammers.</p>
Full article ">Figure 8
<p>Spectrum Distributions (<b>a</b>) Capon in the spatial domain. (<b>b</b>) 3-D plot in the Transmit–Receive domain. (<b>c</b>) 2-D plot in Range-Doppler domain.</p>
Full article ">Figure 9
<p>Beampatterns. (<b>a</b>) Doppler-transmit-receive 3D beampattern in the QEPC-MIMO radar. (<b>b</b>) Transmit-receive slice of the 3D beampattern in the QEPC-MIMO radar.</p>
Full article ">Figure 10
<p>Comparison of output results. (<b>a</b>) Comparison in one true target (<b>b</b>) Comparison in multiple true targets (<b>c</b>) Output SJNR performance with respect to input SNR.</p>
Full article ">Figure 10 Cont.
<p>Comparison of output results. (<b>a</b>) Comparison in one true target (<b>b</b>) Comparison in multiple true targets (<b>c</b>) Output SJNR performance with respect to input SNR.</p>
Full article ">Figure 11
<p>Suppression nulls of FTs. (<b>a</b>) FDA-MIMO. (<b>b</b>) EPC-MIMO. (<b>c</b>) <math display="inline"><semantics> <mrow> <msub> <mi mathvariant="italic">λ</mi> <mn>1</mn> </msub> </mrow> </semantics></math> versus nulls number in the QEPC-MIMO radar. (<b>d</b>) <math display="inline"><semantics> <mrow> <msub> <mi mathvariant="italic">λ</mi> <mn>2</mn> </msub> </mrow> </semantics></math> versus nulls number in the QEPC-MIMO radar.</p>
Full article ">Figure 12
<p>Capon in Spectrum Distributions based on the CBNBC-QEPC-MIMO radar.</p>
Full article ">Figure 13
<p>Beampatterns with CBNBC method. (<b>a</b>) Transmit-receive beampattern. (<b>b</b>) Equivalent transmit beampattern.</p>
Full article ">Figure 14
<p>Output results of CBNBC.</p>
Full article ">
28 pages, 25351 KiB  
Article
Temporal and Spatial Variations of Potential and Actual Evapotranspiration and the Driving Mechanism over Equatorial Africa Using Satellite and Reanalysis-Based Observation
by Isaac Kwesi Nooni, Faustin Katchele Ogou, Jiao Lu, Francis Mawuli Nakoty, Abdoul Aziz Saidou Chaibou, Birhanu Asmerom Habtemicheal, Linda Sarpong and Zhongfang Jin
Remote Sens. 2023, 15(12), 3201; https://doi.org/10.3390/rs15123201 - 20 Jun 2023
Cited by 4 | Viewed by 2025
Abstract
This study investigated the actual evapotranspiration (AET) and potential evapotranspiration (PET) seasonality, trends, abrupt changes, and driving mechanisms with global sea surface temperature (SST) and atmospheric circulation patterns over Equatorial Africa (EQA) during 1980–2020. The spatiotemporal characteristics of mean ET were computed based [...] Read more.
This study investigated the actual evapotranspiration (AET) and potential evapotranspiration (PET) seasonality, trends, abrupt changes, and driving mechanisms with global sea surface temperature (SST) and atmospheric circulation patterns over Equatorial Africa (EQA) during 1980–2020. The spatiotemporal characteristics of mean ET were computed based on a 40-year average at annual and seasonal scales. The Mann-Kendall statistical test, the Sen slope test, and the Bayesian test were used to analyze trends and detect abrupt changes. The results showed that the mean annual PET (AET) for 1980–2020 was 110 (70) mm. Seasonal mean PET (AET) values were 112 (72) in summer, 110 (85) in autumn, 109 (84) in winter, and 110 (58) in spring. The MK test showed an increasing (decreasing) rate, and the Sen slope identified upward (downward) at a rate of 0.35 (0.05) mm yr−10. The PET and AET abrupt change points were observed to happen in 1995 and 2000. Both dry and wet regions showed observed weak (strong) correlation coefficient values of 0.3 (0.8) between PET/AET and climate factors, but significant spatiotemporal differences existed. Generally, air temperature, soil moisture, and relative humidity best explain ET dynamics rather than precipitation and wind speed. The regional atmospheric circulation patterns are directly linked to ET but vary significantly in space and time. From a policy perspective, these findings may have implications for future water resource management. Full article
Show Figures

Figure 1

Figure 1
<p>The study area map of Equatorial Africa. (<b>a</b>) Elevation (m); (<b>b</b>) annual precipitation (mm) and (<b>c</b>) temperature (°C).</p>
Full article ">Figure 2
<p>Spatial pattern means PET (<b>a1</b>–<b>a4</b>) and AET (<b>b1</b>–<b>b4</b>) for annual (<b>a1</b>,<b>b1</b>), summer (<b>a2</b>,<b>b2</b>), autumn (<b>a3</b>,<b>b3</b>), winter (<b>a4</b>,<b>b4</b>), and spring (<b>a5</b>,<b>b5</b>) from 1980 to 2020. The unit is mm.</p>
Full article ">Figure 3
<p>The annual cycles of PET (red color) and AET (blue color) averaged from 1980 to 2020.</p>
Full article ">Figure 4
<p>Linear trends in annual PET (<b>a1</b>) and AET (<b>b1</b>) and seasonal PET (<b>a2</b>–<b>a5</b>) and AET (<b>b2</b>–<b>b5</b>) during 1980–2020 across EQA. A positive sign denotes an upward trend, and a negative value indicates a significant downward trend at <span class="html-italic">p</span> &lt; 0.05.</p>
Full article ">Figure 5
<p>Spatial patterns of annual and seasonal linear trends in PET and AET during 1980–2020 across EQA. Annual trends in PET (<b>a1</b>) and annual AET (<b>b1</b>). Seasonal PET (<b>a2</b>–<b>a5</b>) and AET (<b>b2</b>–<b>b5</b>). The hatched lines denote a linear trend significant at <span class="html-italic">p</span> &lt; 0.05.</p>
Full article ">Figure 6
<p>Abrupt climate changes in annual mean PET (top panel) and AET (bottom panel) during 1980–2020 across EQ. Africa using the Bayesian test. The blue dashed line indicates the year the change occurred with a probability of less than 0.05.</p>
Full article ">Figure 7
<p>The correlation coefficient between PET, AET, and the corresponding precipitation at the annual scale (<b>a1</b>,<b>b1</b>), summer (<b>a2</b>,<b>b2</b>), autumn (<b>a3</b>,<b>b3</b>), winter (<b>a4</b>,<b>b4</b>), and spring (<b>a5</b>,<b>b5</b>) during 1908 to 2020. The hatched area indicates a trend passing the 0.05 significance test.</p>
Full article ">Figure 8
<p>The correlation coefficient between PET, AET and the corresponding mean air temperature at the annual scale (<b>a1</b>,<b>b1</b>), summer (<b>a2</b>,<b>b2</b>), autumn (<b>a3</b>,<b>b3</b>), winter (<b>a4</b>,<b>b4</b>), and spring (<b>a5</b>,<b>b5</b>) during 1980 to 2020. The hatched area indicates a trend passing the 0.05 significance test.</p>
Full article ">Figure 9
<p>The correlation coefficient between PET, AET, and the corresponding wind speed (10 ms<sup>−1</sup>) at the annual scale (<b>a1</b>,<b>b1</b>), summer (<b>a2</b>,<b>b2</b>), autumn (<b>a3</b>,<b>b3</b>), winter (<b>a4</b>,<b>b4</b>), and spring (<b>a5</b>,<b>b5</b>) during 1980 to 2020. The hatched area indicates a trend passing the 0.05 significance test.</p>
Full article ">Figure 10
<p>Heatmap plots comparing temporal correlation coefficients (r) of mean evaporation (actual and potential) during 1980–2020. The correlation coefficient is significant at the 5% level.</p>
Full article ">Figure 11
<p>Spatial distributions of correlation coefficients of PET and AET and SST at the annual scale (<b>a1</b>,<b>b1</b>), summer (<b>a2</b>,<b>b2</b>), autumn (<b>a3</b>,<b>b3</b>), winter (<b>a4</b>,<b>b4</b>), and spring (<b>a5</b>,<b>b5</b>) during 1908 to 2020. The hatched area indicates a trend passing the 0.05 significance test.</p>
Full article ">Figure 12
<p>Climatology of wind circulation at 850 pha over equatorial Africa from 1980 to 2020. annual scale (<b>a</b>), summer (<b>b</b>), autumn (<b>c</b>), winter (<b>d</b>), and spring (<b>e</b>).</p>
Full article ">Figure A1
<p>Correlation coefficient between PET, AET, and the corresponding soil moisture (root zone) at the annual scale (<b>a1</b>,<b>b1</b>), summer (<b>a2</b>,<b>b2</b>), autumn (<b>a3</b>,<b>b3</b>), winter (<b>a4</b>,<b>b4</b>), and spring (<b>a5</b>,<b>b5</b>) during 1908 to 2020. The hatched area indicates a trend that passes the 0.05 significance test.</p>
Full article ">Figure A2
<p>Correlation coefficient between PET, AET, and the corresponding relative humidity at annual scale (<b>a1</b>,<b>b1</b>), Summer (<b>a2</b>,<b>b2</b>), Autumn (<b>a3</b>,<b>b3</b>), Winter (<b>a4</b>,<b>b4</b>) and Spring (<b>a5</b>,<b>b5</b>) during 1908 to 2020. The hatched area indicates a trend that passes the 0.05 significance test.</p>
Full article ">Figure A3
<p>Correlation coefficient between PET, AET, and the corresponding global SST at the annual scale (<b>a1</b>,<b>b1</b>), summer (<b>a2</b>,<b>b2</b>), autumn (<b>a3</b>,<b>b3</b>), winter (<b>a4</b>,<b>b4</b>), and spring (<b>a5</b>,<b>b5</b>) during 1908 to 2020. The hatched area indicates a trend that passes the 0.05 significance test.</p>
Full article ">
27 pages, 8444 KiB  
Article
An Assessment of Negative Samples and Model Structures in Landslide Susceptibility Characterization Based on Bayesian Network Models
by Sahand Khabiri, Matthew M. Crawford, Hudson J. Koch, William C. Haneberg and Yichuan Zhu
Remote Sens. 2023, 15(12), 3200; https://doi.org/10.3390/rs15123200 - 20 Jun 2023
Cited by 12 | Viewed by 2121
Abstract
Landslide susceptibility mapping (LSM) characterizes landslide potential, which is essential for assessing landslide risk and developing mitigation strategies. Despite the significant progress in LSM research over the past two decades, several long-standing issues, such as uncertainties related to training samples and model selection, [...] Read more.
Landslide susceptibility mapping (LSM) characterizes landslide potential, which is essential for assessing landslide risk and developing mitigation strategies. Despite the significant progress in LSM research over the past two decades, several long-standing issues, such as uncertainties related to training samples and model selection, remain inadequately addressed in the literature. In this study, we employed a physically based susceptibility model, PISA-m, to generate four different non-landslide data scenarios and combine them with mapped landslides from Magoffin County, Kentucky, for model training. We utilized two Bayesian network model structures, Naïve Bayes (NB) and Tree-Augmented Naïve Bayes (TAN), to produce LSMs based on regional geomorphic conditions. After internal validation, we evaluated the robustness and reliability of the models using an independent landslide inventory from Owsley County, Kentucky. The results revealed considerable differences between the most effective model in internal validation (AUC = 0.969), which used non-landslide samples extracted exclusively from low susceptibility areas predicted by PISA-m, and the models’ unsatisfactory performance in external validation, as manifested by the identification of only 79.1% of landslide initiation points as high susceptibility areas. The obtained results from both internal and external validation highlighted the potential overfitting problem, which has largely been overlooked by previous studies. Additionally, our findings also indicate that TAN models consistently outperformed NB models when training datasets were the same due to the ability to account for variables’ dependencies by the former. Full article
(This article belongs to the Special Issue Machine Learning and Remote Sensing for Geohazards)
Show Figures

Figure 1

Figure 1
<p>Index maps showing the location of Kentucky and, within Kentucky, the locations of Magoffin and Owsley Counties.</p>
Full article ">Figure 2
<p>Methodological flow chart of the study.</p>
Full article ">Figure 3
<p>Six mapped landslides and selected initiation points. Contours are PISA-m classified areas of landslide susceptibility, overlaying the hillshade of Magoffin County, Kentucky.</p>
Full article ">Figure 4
<p>(<b>a</b>) Magoffin PISA-m susceptibility map; (<b>b</b>) probability density function (pdf) of probability of sliding.</p>
Full article ">Figure 5
<p>Landslide inventory dataset combined with different negative sample scenarios for Magoffin County: (<b>a</b>) Data-1; (<b>b</b>) Data-2; (<b>c</b>) Data-3; (<b>d</b>) Data-4.</p>
Full article ">Figure 6
<p>The geomorphic variables used in model training: (<b>a</b>) elevation, (<b>b</b>) slope, (<b>c</b>) relief, (<b>d</b>) profile curvature, (<b>e</b>) plan curvature, (<b>f</b>) distance to river, (<b>g</b>) distance to road, (<b>h</b>) annual rainfall, (<b>i</b>) geology, and (<b>j</b>) percentage of clay in unit soil profile.</p>
Full article ">Figure 6 Cont.
<p>The geomorphic variables used in model training: (<b>a</b>) elevation, (<b>b</b>) slope, (<b>c</b>) relief, (<b>d</b>) profile curvature, (<b>e</b>) plan curvature, (<b>f</b>) distance to river, (<b>g</b>) distance to road, (<b>h</b>) annual rainfall, (<b>i</b>) geology, and (<b>j</b>) percentage of clay in unit soil profile.</p>
Full article ">Figure 7
<p>Correlation matrices for multicollinearity analysis of datasets: (<b>a</b>) Data-1; (<b>b</b>) Data-2; (<b>c</b>) Data-3; (<b>d</b>) Data-4.</p>
Full article ">Figure 8
<p>Results of validation metrics for eight model–dataset combinations: (<b>a</b>) mean squared error (MSE); (<b>b</b>) area under ROC curve (AUC); (<b>c</b>) accuracy score; (<b>d</b>) precision; (<b>e</b>) recall.</p>
Full article ">Figure 9
<p>Mean of ROC curves for eight model–dataset combinations.</p>
Full article ">Figure 10
<p>Magoffin landslide susceptibility mapping (LSM) based on predictions from four trained models. (<b>c</b>) shows the LSM based on trained model NB-Data-1, with (<b>a</b>,<b>b</b>,<b>d</b>) as zoomed-in views of local areas. (<b>g</b>) shows the LSM based on trained model NB-Data-3, with (<b>e</b>,<b>f</b>,<b>h</b>) as zoomed-in views of local areas. (<b>k</b>) shows the LSM based on trained model TAN-Data-1, with (<b>i</b>,<b>j</b>,<b>l</b>) as zoomed-in views of local areas. (<b>o</b>) shows the LSM based on trained model TAN-Data-3, with (<b>m</b>,<b>n</b>,<b>p</b>) as zoomed-in views of local areas.</p>
Full article ">Figure 11
<p>Probability density function of Magoffin LSM calculated from (<b>a</b>) NB-Data-1; (<b>b</b>) TAN-Data-1; (<b>c</b>) NB-Data-2; (<b>d</b>) TAN-Data-2; (<b>e</b>) NB-Data-3; (<b>f</b>) TAN-Data-3; (<b>g</b>) NB-Data-4; (<b>h</b>) TAN-Data-4.</p>
Full article ">Figure 12
<p>Summarized probability of each susceptibility class from eight trained models.</p>
Full article ">Figure 13
<p>Landslide inventory dataset of Owsley County, KY, USA.</p>
Full article ">Figure 14
<p>The geomorphic variables used in model training for Owsley County, KY, USA: (<b>a</b>) elevation, (<b>b</b>) slope, (<b>c</b>) relief, (<b>d</b>) profile curvature, (<b>e</b>) plan curvature, (<b>f</b>) distance to river, (<b>g</b>) distance to road, (<b>h</b>) annual rainfall, (<b>i</b>) geology, and (<b>j</b>) percentage of clay in unit soil profile.</p>
Full article ">Figure 14 Cont.
<p>The geomorphic variables used in model training for Owsley County, KY, USA: (<b>a</b>) elevation, (<b>b</b>) slope, (<b>c</b>) relief, (<b>d</b>) profile curvature, (<b>e</b>) plan curvature, (<b>f</b>) distance to river, (<b>g</b>) distance to road, (<b>h</b>) annual rainfall, (<b>i</b>) geology, and (<b>j</b>) percentage of clay in unit soil profile.</p>
Full article ">Figure 15
<p>Susceptibility classifications for Owsley County computed using eight trained models: (<b>a</b>) percentages of mapped landslide initiation points and total county area classified as low susceptibility; (<b>b</b>) percentages of mapped landslide initiation points and total county area classified as low-moderate susceptibility; (<b>c</b>) percentages of mapped landslide initiation points and total county area classified as moderate-high susceptibility; (<b>d</b>) percentages of mapped landslide initiation points and total county area classified as high susceptibility.</p>
Full article ">Figure 16
<p>Local susceptibility classifications for Owsley County computed using TAN-Data1 and TAN-Data2 models. Subplots (<b>a</b>,<b>b</b>) represent the same area with different susceptibility characterizations based on predictions from the TAN-Data-1 and TAN-Data-2 models, respectively. Similarly, subplots (<b>c</b>,<b>d</b>) display another comparative analysis of the same area, showing susceptibility predictions from the TAN-Data-1 and TAN-Data-2 models.</p>
Full article ">
15 pages, 6296 KiB  
Technical Note
Region-Based Sea Ice Mapping Using Compact Polarimetric Synthetic Aperture Radar Imagery with Learned Features and Contextual Information
by Saeid Taleghanidoozdoozan, Linlin Xu and David A. Clausi
Remote Sens. 2023, 15(12), 3199; https://doi.org/10.3390/rs15123199 - 20 Jun 2023
Viewed by 1469
Abstract
Operational sea ice maps are usually generated manually using dual-polarization (DP) synthetic aperture radar (SAR) satellite imagery, but there is strong interest in automating this process. Recently launched satellites offer compact polarimetry (CP) imagery that provides more comprehensive polarimetric information compared to DP, [...] Read more.
Operational sea ice maps are usually generated manually using dual-polarization (DP) synthetic aperture radar (SAR) satellite imagery, but there is strong interest in automating this process. Recently launched satellites offer compact polarimetry (CP) imagery that provides more comprehensive polarimetric information compared to DP, which compels the use of CP for automated classification of SAR sea ice imagery. Existing sea ice scene classification algorithms using CP imagery rely on handcrafted features, while neural networks offer the potential of features that are more discriminating. We have developed a new and effective sea ice classification algorithm that leverages the nature of CP data. First, a residual-based convolutional neural network (ResCNN) is implemented to classify each pixel. In parallel, an unsupervised segmentation is performed to generate regions based on CP statistical properties. Regions are assigned a single class label by majority voting using the ResCNN output. For testing, quad-polarimetric (QP) SAR sea ice scenes from the RADARSAT Constellation Mission (RCM) are used, and QP, DP, CP, and reconstructed QP modes are compared for classification accuracy, while also comparing them to other classification approaches. Using CP achieves an overall accuracy of 96.86%, which is comparable to QP (97.16%), and higher than reconstructed QP and DP data by about 2% and 10%, respectively. The implemented algorithm using CP imagery provides an improved option for automated sea ice mapping. Full article
(This article belongs to the Section Ocean Remote Sensing)
Show Figures

Figure 1

Figure 1
<p>Flowchart of the main steps of the proposed classification method.</p>
Full article ">Figure 2
<p>(<b>a</b>) The locations of the two RADARSAT-2 fine QP scenes in the Barrow Strait. The first element of the covariance matrix of Scenes (<b>b</b>) 56 and (<b>c</b>) 58, along with the overlaid labeled pixels of open water/new ice class (blue), young ice (violet), first-year ice (yellow), and multi-year ice (red). (<b>d</b>–<b>h</b>) are the segmentation images obtained using IRGS-based methods.</p>
Full article ">Figure 3
<p>(<b>a</b>–<b>d</b>) Pixel-level sea ice maps generated by the ResCNN feature learning classifier using DP (DP–ResCNN), derived CP (CP–ResCNN), RQP (RQP–ResCNN), and QP SAR data (QP–ResCNN). (<b>e</b>–<b>h</b>) segmentation combined with ResCNN classification results along with their overall accuracy(OA). It should be noted that the QP results have been rescaled to match CP data size for presentation purposes.</p>
Full article ">Figure 4
<p>Sea ice maps indicating OA for baseline approaches. (<b>a</b>) is pixel-based and (<b>b</b>) is region-based, using method by Ghanbari et al. [<a href="#B2-remotesensing-15-03199" class="html-bibr">2</a>], while (<b>c</b>) is pixel-based and (<b>d</b>) is region-based, using method by Leigh et al. [<a href="#B28-remotesensing-15-03199" class="html-bibr">28</a>].</p>
Full article ">
20 pages, 9264 KiB  
Article
Shoreline Analysis and Extraction Tool (SAET): A New Tool for the Automatic Extraction of Satellite-Derived Shorelines with Subpixel Accuracy
by Jesús Palomar-Vázquez, Josep E. Pardo-Pascual, Jaime Almonacid-Caballer and Carlos Cabezas-Rabadán
Remote Sens. 2023, 15(12), 3198; https://doi.org/10.3390/rs15123198 - 20 Jun 2023
Cited by 9 | Viewed by 4391
Abstract
SAET (Shoreline Analysis and Extraction Tool) is a novel open-source tool to enable the completely automatic detection of shoreline position changes using the optical imagery acquired by the Sentinel-2 and Landsat 8 and 9 satellites. SAET has been developed within the ECFAS (European [...] Read more.
SAET (Shoreline Analysis and Extraction Tool) is a novel open-source tool to enable the completely automatic detection of shoreline position changes using the optical imagery acquired by the Sentinel-2 and Landsat 8 and 9 satellites. SAET has been developed within the ECFAS (European Coastal Flood Awareness System) project, which is intended to be the first European service for coastal flood forecasting, management, and recovery analysis. The tool is developed to characterise the shoreline response associated with punctual events such as coastal storms as well as any other phenomenon. For a given beach segment, SAET facilitates the selection of the satellite images closest in time to the analysed events that offer an adequate cloud coverage level for analysing the shoreline change. Subsequently, the tool automatically downloads the images from their official repositories, pre-processes them and extracts the shoreline position with sub-pixel accuracy. In order to do so, an initial approximate definition of the shoreline is carried out at the pixel level using a water index thresholding, followed by an accurate extraction operating on the shortwave infrared bands to produce a sub-pixel line in vector formats (points and lines). The tool offers different settings to be adapted to the different coastal environments and beach typologies. Its main advantages refer to its autonomy, its efficiency in extracting complete satellite scenes, its flexibility in adapting to different environments and conditions, and its high subpixel accuracy. This work presents an accuracy assessment on a long Mediterranean sandy beach of SDSs extracted from L8 and S2 imagery against coincident alongshore reference lines, showing an accuracy of about 3 m RMSE. At the same time, the work shows an example of the usage of SAET for characterising the response to Storm Gloria (January 2020) on the beaches of Valencia (E Spain). SAET provides an efficient and completely automatic workflow that leads to accurate SDSs while only relying on publicly available information. The tool appears to be a useful extraction tool for beach monitoring, both for public administrations and individual users. Full article
Show Figures

Figure 1

Figure 1
<p>SAET workflow.</p>
Full article ">Figure 2
<p>Example of the results obtained in the searching mode (S2 images). On the left, the first list shows the found images and their availability, while the second list is an ordered-by-date list including the central date and showing only the images that are available (online). On the right, the HTML file displays an overview of the found images. The yellow box highlights the information for the closest image to the analysed date. The parameters to obtain this result are the following (see <a href="#remotesensing-15-03198-t001" class="html-table">Table 1</a>): --rm=os --fp=NONE --sd=20220801 --cd=20220815 --ed=20220830 --mc=10 --lp=NONE --ll=NONE --sp=S2MSI1C --sl=30SYJ.</p>
Full article ">Figure 3
<p>For the coast of Cornwall (U.K.), the four first images show (<b>a</b>) the RGB Sentinel-2, the indices (<b>b</b>) AWEInsh, (<b>c</b>) AWEIsh, and (<b>d</b>) MNDWI using the 0 threshold for the date 26 February 2022. The next four images show (<b>e</b>) the Sentinel-2 and the indices in the same order (<b>f</b>–<b>h</b>) for the date 8 April 2017.</p>
Full article ">Figure 4
<p>S2 (<b>a</b>) and L9 (<b>b</b>) images of St. Ives Bay (Cornwall, U.K.) were acquired the same day with only 10 min of difference. The classification, mainly focused on water and clouds (high confidence) classes, shows several problems of misclassification.</p>
Full article ">Figure 5
<p>Location map of the study area between Valencia Port and the Cape of Cullera (Spain). The sections considered in the shoreline assessment appeared in black (7 May 2018), green (14 September 2015), yellow (10 July 2018) and blue (26 January 2020). The SIMAR point was employed for defining the oceanographic data at the moment of the assessment. Source: EO Browser, <a href="https://apps.sentinel-hub.com/eo-browser/" target="_blank">https://apps.sentinel-hub.com/eo-browser/</a> (accessed on 10 April 2023), Sinergise Ltd. (Ljubljana, Slovenia).</p>
Full article ">Figure 6
<p>Errors of the SDSs along the study site. Positive values represent a seaward bias, while negative values indicate a landward bias.</p>
Full article ">Figure 7
<p>Significant wave height (m) and availability of satellite imagery during the occurrence of Storm Gloria (January 2020). Data from the SIMAR station 2,081,111 (0.25°O, 39.25°N) by the Spanish Ports Authority (<a href="https://www.puertos.es/es-es/oceanografia/Paginas/portus.aspx" target="_blank">https://www.puertos.es/es-es/oceanografia/Paginas/portus.aspx</a>, accessed on 10 April 2023). The pre-storm and post-storm images employed for the shoreline extraction are shown in green and red colour respectively. Likewise, coinciding with the post-storm image, the very high-resolution image used for the accuracy assessment also appears.</p>
Full article ">Figure 8
<p>The study area (in red) is covered by both a Sentinel-2 tile (30SYJ, in blue) and a Landsat scene (198,033, in purple). The choice of one or the other product will depend on which images are closest to the peak of the storm.</p>
Full article ">Figure 9
<p>Ordered-by-date list of the available satellite images. This list offers some information about each image, such as the identifier, the cloud coverage percentage, and the difference in days with respect to the peak of the storm. In this example, Sentinel-2 numbers 0 (26 January) and 1 (16 January) are selected to run SAET in the ‘downloading and processing’ mode.</p>
Full article ">Figure 10
<p>Detail of the extracted shorelines associated with Storm Gloria. The pre-storm (in green) and post-storm (in red) SDSs are displayed together with the reference line from the VHR image (in yellow). Coordinate system: WGS84/UTM 30N. Background image: Pléiades-1A © CNES (2020), distributed by Airbus DS, provided under COPERNICUS by the European Union and ESA, all rights reserved.</p>
Full article ">Figure 11
<p>Shoreline change measured with SAET as response to Storm Gloria along 28 km of the Valencian coast. Coordinate system: WGS84/UTM 30N. Background image: Pléiades-1A © CNES (2020), distributed by Airbus DS, provided under COPERNICUS by the European Union and ESA, all rights reserved.</p>
Full article ">Figure 12
<p>Tidal range along the coast of Europe. Data obtained from the Global Tidal Range dataset computed by ESRI (<a href="https://www.arcgis.com/home/item.html?id=d5354dea41b14f0689860bf4b2cf5e8a" target="_blank">https://www.arcgis.com/home/item.html?id=d5354dea41b14f0689860bf4b2cf5e8a</a>, accessed on 2 June 2023) using the FES2014 tide model produced by Noveltis, Legos and CLS and distributed by Aviso+, with support from Cnes (<a href="https://www.aviso.altimetry.fr/" target="_blank">https://www.aviso.altimetry.fr/</a>, accessed on 2 June 2023).</p>
Full article ">
20 pages, 11704 KiB  
Article
Winter Wheat Drought Risk Assessment by Coupling Improved Moisture-Sensitive Crop Model and Gridded Vulnerability Curve
by Haibo Yang, Zenglan Li, Qingying Du and Zheng Duan
Remote Sens. 2023, 15(12), 3197; https://doi.org/10.3390/rs15123197 - 20 Jun 2023
Viewed by 1474
Abstract
The crop drought risk assessment is an important basis for mitigating the effects of drought on crops. The study of drought using crop growth models is an integral part of agricultural drought risk research. The current Decision Support System for Agrotechnology Transfer (DSSAT) [...] Read more.
The crop drought risk assessment is an important basis for mitigating the effects of drought on crops. The study of drought using crop growth models is an integral part of agricultural drought risk research. The current Decision Support System for Agrotechnology Transfer (DSSAT) model is not sufficiently sensitive to moisture parameters when performing simulations, and most studies that conduct different scenario simulations to assess crop drought vulnerability are based on the site-scale. In this paper, we improved the moisture sensitivity of the Crop Environment Resource Synthesis System (CERES)-Wheat to improve the simulation accuracy of winter wheat under water stress, and then we assessed the drought intensity in the Beijing-Tianjin-Hebei region and constructed a gridded vulnerability curve. The grid vulnerability curves (1 km × 1 km) were quantitatively characterized using key points, and the drought risk distribution and zoning of winter wheat were evaluated under different return periods. The results show that the stress mechanism of coupled water and photosynthetic behavior improved the CERES-Wheat model. The accuracy of the modified model improved in terms of the above-ground biomass and yield compared with that before the modification, with increases of 20.39% and 11.45% in accuracy, respectively. The drought hazard intensity index of winter wheat in the study area from 1970 to 2019 exhibited a trend of high in the southwest and low in the southeast. The range of the multi-year average drought hazard intensity across the region was 0.29–0.61. There were some differences in the shape and characteristic covariates of the drought vulnerability curves among the different sub-zones. In terms of the cumulative loss rates, almost the entire region had a cumulative drought loss rate of 49.00–54.00%. Overall, the drought risk index decreased from west to east and from north to south under different return periods. This quantitative evaluation of the drought hazard intensity index provides a reference for agricultural drought risk evaluation. Full article
(This article belongs to the Special Issue Crop Quantitative Monitoring with Remote Sensing)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>The study area and distribution of winter wheat.</p>
Full article ">Figure 2
<p>Meteorological geographical data for the wheat growing periods from 1970 to 2019 partitioning in the Beijing-Tianjin-Hebei region: (<b>a</b>) average annual precipitation, (<b>b</b>) average annual sunshine hours, (<b>c</b>) average annual maximum temperature, (<b>d</b>) average annual minimum temperature and (<b>e</b>) DEM.</p>
Full article ">Figure 3
<p>Simulated and measured values of above-ground biomass throughout the growing period before and after model modification.</p>
Full article ">Figure 4
<p>(<b>a</b>) Original model and (<b>b</b>) modified model simulated yield values vs. measured values.</p>
Full article ">Figure 5
<p>Comparison of simulated and measured values for the (<b>a</b>) anthesis and (<b>b</b>) maturity stages after localization of genetic parameters.</p>
Full article ">Figure 6
<p>The spatial distribution of the DHI for winter wheat in a time series: (<b>a</b>) 1970, (<b>b</b>) 1975, (<b>c</b>) 1980, (<b>d</b>) 1985, (<b>e</b>) 1990, (<b>f</b>) 1995, (<b>g</b>) 2000, (<b>h</b>) 2005, (<b>i</b>) 2010, (<b>j</b>) 2015 and (<b>k</b>) 2018.</p>
Full article ">Figure 7
<p>Winter wheat zoning in the Beijing-Tianjin-Hebei region.</p>
Full article ">Figure 8
<p>Drought vulnerability curves for winter wheat in different sub-zones ((<b>a</b>–<b>l</b>) represent Nos. 1–12).</p>
Full article ">Figure 9
<p>Characteristic drought vulnerability curve parameters for winter wheat in different sub-zones ((<b>a</b>–<b>l</b>) represent Nos.1–12).</p>
Full article ">Figure 10
<p>Distribution of yield loss rates for fixed risk levels of winter wheat drought in the Beijing-Tianjin-Hebei region for return periods of (<b>a</b>) 2, (<b>b</b>) 5, (<b>c</b>) 10, (<b>d</b>) 20, (<b>e</b>) 25 and (<b>f</b>) 50 years.</p>
Full article ">Figure 11
<p>Distribution of winter wheat drought risk zones in the Beijing-Tianjin-Hebei region for return periods of (<b>a</b>) 2, (<b>b</b>) 5, (<b>c</b>) 10, (<b>d</b>) 20, (<b>e</b>) 25 and (<b>f</b>) 50 years.</p>
Full article ">Figure 12
<p>Global sensitivity indices of the (<b>a</b>) maximum leaf area index, (<b>b</b>) above-ground biomass and (<b>c</b>) yield to 10 soil parameters simulated using the modified model and the original model for the 2015 wheat growing season, where 1 represents the improved model and 2 represents the original model.</p>
Full article ">
19 pages, 7596 KiB  
Article
Assessing the Magnitude of the Amazonian Forest Blowdowns and Post-Disturbance Recovery Using Landsat-8 and Time Series of PlanetScope Satellite Constellation Data
by Dazhou Ping, Ricardo Dalagnol, Lênio Soares Galvão, Bruce Nelson, Fabien Wagner, David M. Schultz and Polyanna da C. Bispo
Remote Sens. 2023, 15(12), 3196; https://doi.org/10.3390/rs15123196 - 20 Jun 2023
Cited by 4 | Viewed by 7703
Abstract
Blowdown events are a major natural disturbance in the central Amazon Forest, but their impact and subsequent vegetation recovery have been poorly understood. This study aimed to track post-disturbance regeneration after blowdown events in the Amazon Forest. We analyzed 45 blowdown sites identified [...] Read more.
Blowdown events are a major natural disturbance in the central Amazon Forest, but their impact and subsequent vegetation recovery have been poorly understood. This study aimed to track post-disturbance regeneration after blowdown events in the Amazon Forest. We analyzed 45 blowdown sites identified after September 2020 at Amazonas, Mato Grosso, and Colombia jurisdictions using Landsat-8 and PlanetScope NICFI satellite imagery. Non-photosynthetic vegetation (NPV), green vegetation (GV), and shade fractions were calculated for each image and sensor using spectral mixture analysis in Google Earth Engine. The results showed that PlanetScope NICFI data provided more regular and higher-spatial-resolution observations of blowdown areas than Landsat-8, allowing for more accurate characterization of post-disturbance vegetation recovery. Specifically, NICFI data indicated that just four months after the blowdown event, nearly half of ΔNPV, which represents the difference between the NPV after blowdown and the NPV before blowdown, had disappeared. ΔNPV and GV values recovered to pre-blowdown levels after approximately 15 months of regeneration. Our findings highlight that the precise timing of blowdown detection has huge implications on quantification of the magnitude of damage. Landsat data may miss important changes in signal due to the difficulty of obtaining regular monthly observations. These findings provide valuable insights into vegetation recovery dynamics following blowdown events. Full article
(This article belongs to the Special Issue Remote Sensing of the Amazon Region)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Study sites with related blowdown events found in the Amazon Forest. The background image corresponds to the areas covered by rainforests and other land covers in the Amazon Forest, according to the Tropical Moist Forests product [<a href="#B42-remotesensing-15-03196" class="html-bibr">42</a>].</p>
Full article ">Figure 2
<p>Spectral reflectance of Landsat-8 OLI and PlanetScope NICFI endmembers.</p>
Full article ">Figure 3
<p>Examples of true-color composites of blowdown-disturbed areas, as indicated by red lines. Images (<b>A1</b>,<b>B1</b>,<b>C1</b>) refer to PlanetScope NICFI, whereas images (<b>A2</b>,<b>B2</b>,<b>C2</b>) correspond to OLI Landsat-8.</p>
Full article ">Figure 4
<p>Examples of PlanetScope NICFI basemap (<b>A</b>), blowdown mapping considering a ΔNPV threshold greater than 0.3 for PlanetScope NICFI (<b>B</b>), and Landsat-8 (<b>C</b>). The longitude and latitude of this blowdown are 63.75°W, 2.31°S; image data: October 2020.</p>
Full article ">Figure 5
<p>Boxplots of (<b>A</b>) ΔNPV<sub>first</sub> (the difference between the average NPV of the first post-blowdown image acquisition and of the pre-blowdown images) and (<b>B</b>) max ΔNPV (the max ΔNPV value among all months post-blowdown). The notch in the boxplots represents the 95% confidence interval around the median.</p>
Full article ">Figure 6
<p>Changes in ΔNPV over time for all 45 blowdowns based on (<b>A</b>) PlanetScope NICFI and (<b>B</b>) Landsat-8.</p>
Full article ">Figure 7
<p>Changes in ΔNPV over time based on PlanetScope NICFI. (Blue line: mean ΔNPV values; grey curve: each blowdown event; green shades: 95% CI.)</p>
Full article ">Figure 8
<p>Changes in GV ratio (GV/pre-blowdown mean GV) over time based on PlanetScope NICFI. (Blue line: mean values; gray curve: each blowdown event; green shades: 95% CI.) There are negative GV values because the SMA did not apply non-negativity constraints.</p>
Full article ">Figure 9
<p>Time series at monthly scale of NPV fractions from PlanetScope NICFI. (The longitude and latitude of this blowdown are 72.42°W, 1.22°S; image of February 2021 is absent due to clouds).</p>
Full article ">
18 pages, 4427 KiB  
Article
Characteristics of Regional GPS Crustal Deformation before the 2021 Yunnan Yangbi Ms 6.4 Earthquake and Its Implications for Determining Potential Areas of Future Strong Earthquakes
by Chenglong Dai, Weijun Gan, Zhangjun Li, Shiming Liang, Genru Xiao, Keliang Zhang and Ling Zhang
Remote Sens. 2023, 15(12), 3195; https://doi.org/10.3390/rs15123195 - 20 Jun 2023
Viewed by 1528
Abstract
The 2021 Yangbi Ms 6.4 earthquake in Yunnan, China, occurred in an area where the Global Positioning System (GPS) geodetic observations are particularly intensive. Based on a detailed retrospective analysis of the GPS observations of about 133 stations distributed in the proximately 400 [...] Read more.
The 2021 Yangbi Ms 6.4 earthquake in Yunnan, China, occurred in an area where the Global Positioning System (GPS) geodetic observations are particularly intensive. Based on a detailed retrospective analysis of the GPS observations of about 133 stations distributed in the proximately 400 km × 400 km region that contains the area affected by the earthquake., we obtain a high-resolution GPS velocity field and strain rate field and then derive the present-day slip rates of major faults in the region with the commonly used half-space elastic dislocation model and constraints from the GPS velocity field. Furthermore, by calculating the seismic moment accumulation and release and deficit rates in the main fault segments and combining with the distribution characteristics of small earthquakes, we evaluate the regional seismic risk. The results show that (1) there was a localized prominent strain accumulation rate around the seismogenic area of the impending Yangbi Ms 6.4 earthquake, although this was not the only area with a prominent strain rate in the whole region. (2) The seismogenic area of the earthquake was just located where the strain direction was deflected, which, together with the localized outstanding maximum shear strain and dilatation rates, provides us with important hints to determine the potential areas of future strong earthquakes. (3) Of all the seismogenic fault segments with relatively high potentials, judged using the elapsed time of historical earthquakes and effective strain accumulation rate, the middle section of the Weixi–Qiaohou fault has a higher earthquake risk than the southern section, the Midu–Binchuan section of the Chenghai fault has a higher risk than the Yongsheng section and the Jianchuan section of the Jianchuan–Qiaohou–Lijiang–Xiaojinhe fault has a higher risk than the Lijiang section. Full article
(This article belongs to the Special Issue Monitoring Subtle Ground Deformation of Geohazards from Space)
Show Figures

Figure 1

Figure 1
<p>Schematic diagram of the distribution of active faults and moderately strong earthquakes in northwest Yunnan, where the Yangbi Ms 6.4 earthquake of 2021 was located. The green star denotes epicenter and the dense blue circles denote foreshocks and aftershocks before 24 May (earthquake catalog is from Yang et al. [<a href="#B8-remotesensing-15-03195" class="html-bibr">8</a>]). The colored dots are historical earthquakes with M ≥ 5 from 1900 to the present, where the radius of the dots distinguishes the magnitude and the color distinguishes the approximate age. Seismic data were obtained from Xu et al. [<a href="#B9-remotesensing-15-03195" class="html-bibr">9</a>] and the China Earthquake Network Center (<a href="http://www.cenc.ac.cn/" target="_blank">http://www.cenc.ac.cn/</a> (accessed on 1 January 2021)). NTHF: Nantinghe fault; WTADF: Wanding–Anding fault; NJF: Nujiang fault; LCJF: Lancangjiang fault; WX-QHF: Weixi–Qiaohou fault; LJXJHF: Lijiang–Xiaojinhe fault; NH-CHF: Nanhua–Chuxiong fault; RRF: Red River fault; HQ-ERF: Heqing–Eryuan fault; CHF: Chenghai fault; NJ-WSF: Nanjian–Weishan fault.</p>
Full article ">Figure 2
<p>GPS crustal movement horizontal velocity field under the Eurasian reference frame. The blue arrows are the velocity field of the ‘China Crustal Movement Observation Network’ (Phase I and II) calculated by Wang and Shen [<a href="#B6-remotesensing-15-03195" class="html-bibr">6</a>]. The red arrows are the velocity field of the self-built observation station, and the error ellipses represent 70% confidence.</p>
Full article ">Figure 3
<p>Smoothing distance of strain rate field.</p>
Full article ">Figure 4
<p>Strain rate field and rotation rate field. (<b>a</b>) The maximum shear strain rate (background color), principal strain rate; (<b>b</b>) dilatational strain rate (background color, the negative indicates compression and the positive denotes extension), principal strain rate; (<b>c</b>) rotation rate (the negative indicates clockwise and the positive indicates anti-clockwise).</p>
Full article ">Figure 5
<p>Three-dimensional fault model. F1 Nantinghe fault; F2 Wanding fault; F3 Nujiang fault; F4 Lancangjiang fault; F5 Lanping–Yunlong fault F6 Yongping fault; F7 Weixi–Qiaohou fault; F8 south section of Jinshajiang fault (Jidala section); F9 Deqin–Zhongdian–Daju fault; F10 Lijiang–Xiaojinhe fault; F11 Yuanmou fault; F12 Chuxiong–Nanhua fault; F13 Honghe fault; F14 Chenghai fault; F15 Heqing–Eryuan fault; F16 Jianchuan–Qiaohou fault; F17 Nanjian–Weishan fault.</p>
Full article ">Figure 6
<p>(<b>a</b>) Residual value of model fitting; (<b>b</b>) modeled strike—slip rates. Red lines denote sinistral movements and blue lines represent dextral movements. (<b>c</b>) Modeled dip—slip rates. Red lines indicate tension rate and blue lines indicate extrusion rate. The faults abbreviations are shown in <a href="#remotesensing-15-03195-f001" class="html-fig">Figure 1</a>.</p>
Full article ">Figure 7
<p>The seismic moment deficit of the main faults, where the color of the fault segment indicates the seismic moment deficit, the thickness indicates the sum of the total sliding rates of the fault and the number in parentheses indicates the magnitude of the equivalent moment of the earthquake if the earthquake released all moments.</p>
Full article ">Figure 8
<p>Seismic risk distribution map. The red area satisfies three characteristics and the blue area satisfies two characteristics. The small red dot indicates small earthquakes from April 2018 to August 2020.</p>
Full article ">Figure 9
<p>Block division, fitting residuals and distribution of fault locking.</p>
Full article ">
17 pages, 2862 KiB  
Article
Satellite Sensed Data-Dose Response Functions: A Totally New Approach for Estimating Materials’ Deterioration from Space
by Georgios Kouremadas, John Christodoulakis, Costas Varotsos and Yong Xue
Remote Sens. 2023, 15(12), 3194; https://doi.org/10.3390/rs15123194 - 20 Jun 2023
Viewed by 2470
Abstract
When construction materials are exposed to the atmospheric environment, they are subject to deterioration, which varies according to the time period of exposure and the location. A tool named Dose–Response Functions (DRFs) has been developed to estimate this deterioration. DRFs use specific air [...] Read more.
When construction materials are exposed to the atmospheric environment, they are subject to deterioration, which varies according to the time period of exposure and the location. A tool named Dose–Response Functions (DRFs) has been developed to estimate this deterioration. DRFs use specific air pollutants and climatic parameters as input data. Existing DRFs in the literature use only ground-based measurements as input data. This fact constitutes a limitation for the application of this tool because it is too expensive to establish and maintain such a large network of ground-based stations for pollution monitoring. In this study, we present the development of new DRFs using only satellite data as an input named Satellite Sensed Data Dose-Response Functions (SSD-DRFs). Due to the global coverage provided by satellites, this new tool for monitoring the corrosion/soiling of materials overcomes the previous limitation because it can be applied to any area of interest. To develop SSD-DRFs, we used measurements from MODIS (Moderate Resolution Imaging Spectroradiometer) and AIRS (Atmospheric Infrared Sounder) on board Aqua and OMI (Ozone Monitoring Instrument) on Aura. According to the obtained results, SSD-DRFs were developed for the case of carbon steel, zinc, limestone and modern glass materials. SSD-DRFs are shown to produce more reliable corrosion/soiling estimates than “traditional” DRFs using ground-based data. Furthermore, research into the development of the SSD-DRFs revealed that the different corrosion mechanisms taking place on the surface of a material do not act additively with each other but rather synergistically. Full article
Show Figures

Figure 1

Figure 1
<p>Thin film of water formation on the material’s surface.</p>
Full article ">Figure 2
<p>The experimentally observed mass loss of carbon steel (red bars), the estimated mass loss of carbon steel using SSD-DRF (blue bars) and the estimated mass loss of carbon steel using G-DRF (yellow bars) for case (<b>a</b>) Athens, (<b>b</b>) Bottrop, (<b>c</b>) Kopisty, (<b>d</b>) Madrid, (<b>e</b>) Prague and (<b>f</b>) Toledo after five different one-year exposure periods (2005–2006, 2008–2009, 2011–2012, 2014–2015, 2017–2018). The missing yellow bars in the case of Athens are due to the lack of the necessary ground-based data necessary for the application of G-DRF.</p>
Full article ">Figure 3
<p>As in <a href="#remotesensing-15-03194-f002" class="html-fig">Figure 2</a>, but for the case of zinc. (<b>a</b>) Athens, (<b>b</b>) Bottrop, (<b>c</b>) Kopisty, (<b>d</b>) Madrid, (<b>e</b>) Prague and (<b>f</b>) Toledo.</p>
Full article ">Figure 4
<p>As in <a href="#remotesensing-15-03194-f002" class="html-fig">Figure 2</a>, but for the case of Limestone. (<b>a</b>) Athens, (<b>b</b>) Bottrop, (<b>c</b>) Kopisty, (<b>d</b>) Madrid, (<b>e</b>) Prague and (<b>f</b>) Toledo.</p>
Full article ">Figure 5
<p>As in <a href="#remotesensing-15-03194-f002" class="html-fig">Figure 2</a>, but for the case of modern glass haze. (<b>a</b>) Athens, (<b>b</b>) Bottrop, (<b>c</b>) Kopisty, (<b>d</b>) Madrid, (<b>e</b>) Prague and (<b>f</b>) Toledo. The gaps at Athens and Bottrop at the exposure periods of 2005–2006 and 2008–2009, respectively, are due to a lack of experimental data.</p>
Full article ">Figure 6
<p>The relative differences between corrosion/soiling estimates calculated using SSD-DRFs and G-DRFs and experimentally obtained data for the case of (<b>a</b>) carbon steel mass loss, (<b>b</b>) zinc mass loss, (<b>c</b>) limestone recession and (<b>d</b>) modern glass haze.</p>
Full article ">
21 pages, 3709 KiB  
Article
Exploring Deep Learning Models on GPR Data: A Comparative Study of AlexNet and VGG on a Dataset from Archaeological Sites
by Merope Manataki, Nikos Papadopoulos, Nikolaos Schetakis and Alessio Di Iorio
Remote Sens. 2023, 15(12), 3193; https://doi.org/10.3390/rs15123193 - 20 Jun 2023
Cited by 3 | Viewed by 2607
Abstract
This comparative study evaluates the performance of three popular deep learning architectures, AlexNet, VGG-16, and VGG-19, on a custom-made dataset of GPR C-scans collected from several archaeological sites. The introduced dataset has 15,000 training images and 3750 test images assigned to three classes: [...] Read more.
This comparative study evaluates the performance of three popular deep learning architectures, AlexNet, VGG-16, and VGG-19, on a custom-made dataset of GPR C-scans collected from several archaeological sites. The introduced dataset has 15,000 training images and 3750 test images assigned to three classes: Anomaly, Noise, and Structure. The aim is to assess the performance of the selected architectures applied to the custom dataset and examine the potential gains of using deeper and more complex architectures. Further, this study aims to improve the training dataset using augmentation techniques. For the comparisons, learning curves, confusion matrices, precision, recall, and f1-score metrics are employed. The Grad-CAM technique is also used to gain insights into the models’ learning. The results suggest that using more convolutional layers improves overall performance. Further, augmentation techniques can also be used to increase the dataset volume without causing overfitting. In more detail, the best-obtained model was trained using VGG-19 architecture and the modified dataset, where the training samples were raised to 60,000 images through augmentation techniques. This model reached a classification accuracy of 94.12% on an evaluation set with 170 unseen data. Full article
(This article belongs to the Special Issue Application of Remote Sensing in Cultural Heritage Research II)
Show Figures

Figure 1

Figure 1
<p>Random training set samples for the classes Anomaly, Noise, and Structure. On the top row are 100 samples from dataset-1 per class, while on the bottom row are 100 samples from dataset-2, where the volume is increased using data augmentation techniques to produce synthetic data.</p>
Full article ">Figure 2
<p>Schematic description of AlexNet, VGG-16, and VGG-19 architectures as they were implemented for the purposes of the study. The filter size, filter number, stride, and pad size are given on each convolutional layer. Similarly, the kernel size, stride, and pad size are given on the pooling layers. A pad of 1 indicates no changes in dimension, while a pad of 2 indicates dimensionality reduction. On the bottom right of each layer, the output dimensions are noted. Last, the feature extraction and classification stages are pointed out in each architecture.</p>
Full article ">Figure 3
<p>The resulting loss and accuracy learning curves for all 9 models grouped by the architecture. Blue color represents the results obtained from models AlexNet-1, VGG16-1, and VGG19-1 that were trained with dataset-1 (~15,000 training samples) and orange represents the results obtained from models AlexNet-2, VGG16-2, and VGG19-2 trained with dataset-2 (~60,000 training samples), while green color represents the results obtained from AlexNet-3, VGG-16-3, and VGG19-3 trained with dataset-1 using image augmentation techniques to replace training samples. The dashed line indicates the accuracy and loss calculated on the training set, while the solid line is the validation accuracy and validation loss calculated on the test set for each dataset.</p>
Full article ">Figure 4
<p>The confusion matrices were calculated on the evaluation set for all 9 models. On top of each matrix, the total accuracy score is presented.</p>
Full article ">Figure 5
<p>Selective Grad-CAM results of the evaluation set for correct predictions under the Anomaly class (<b>a</b>), the Noise class (<b>b</b>), and the Structure class (<b>c</b>). For all three cases, the first column is the input to trained models, and the following columns are the generated heatmaps for each model overlaid on the input image. Warm colors indicate the most important regions for each model’s correct prediction, while cooler colors suggest little to no contribution to the classification prediction.</p>
Full article ">Figure 6
<p>Compilation of Grad-CAM results for the samples where models made the most mistakes. On the top row are the inputs, where each class is mentioned at the bottom. The generated heatmaps for all 9 models are overlaid on each input. Warm colors indicate the highest impact, while cooler colors indicate little to no impact on the classification prediction. The classification prediction of the highest two percentages is presented at the bottom of each sample, with the first one being the classification outcome. Red indicates the wrong class, while blue indicates the correct class.</p>
Full article ">
17 pages, 5528 KiB  
Article
Applying Reconstructed Daily Water Storage and Modified Wetness Index to Flood Monitoring: A Case Study in the Yangtze River Basin
by Cuiyu Xiao, Yulong Zhong, Yunlong Wu, Hongbing Bai, Wanqiu Li, Dingcheng Wu, Changqing Wang and Baoming Tian
Remote Sens. 2023, 15(12), 3192; https://doi.org/10.3390/rs15123192 - 20 Jun 2023
Cited by 3 | Viewed by 2481
Abstract
The terrestrial water storage anomaly (TWSA) observed by the Gravity Recovery and Climate Experiment (GRACE) satellite and its successor GRACE Follow-On (GRACE-FO) provides a new means for monitoring floods. However, due to the coarse temporal resolution of GRACE/GRACE-FO, the understanding of flood occurrence [...] Read more.
The terrestrial water storage anomaly (TWSA) observed by the Gravity Recovery and Climate Experiment (GRACE) satellite and its successor GRACE Follow-On (GRACE-FO) provides a new means for monitoring floods. However, due to the coarse temporal resolution of GRACE/GRACE-FO, the understanding of flood occurrence mechanisms and the monitoring of short-term floods are limited. This study utilizes a statistical model to reconstruct daily TWS by combining monthly GRACE observations with daily temperature and precipitation data. The reconstructed daily TWSA is utilized to monitor the catastrophic flood event that occurred in the middle and lower reaches of the Yangtze River basin in 2020. Furthermore, the study compares the reconstructed daily TWSA with the vertical displacements of eight Global Navigation Satellite System (GNSS) stations at grid scale. A modified wetness index (MWI) and a normalized daily flood potential index (NDFPI) are introduced and compared with in situ daily streamflow to assess their potential for flood monitoring and early warning. The results show that terrestrial water storage (TWS) in the study area increases from early June, reaching a peak on 19 July, and then receding till September. The reconstructed TWSA better captures the changes in water storage on a daily scale compared to monthly GRACE data. The MWI and NDFPI based on the reconstructed daily TWSA both exceed the 90th percentile 7 days earlier than the in situ streamflow, demonstrating their potential for daily flood monitoring. Collectively, these findings suggest that the reconstructed TWSA can serve as an effective tool for flood monitoring and early warning. Full article
(This article belongs to the Special Issue GRACE for Earth System Mass Change: Monitoring and Measurement)
Show Figures

Figure 1

Figure 1
<p>Map of the study region. The red line represents the study area boundary defined in this study. The red triangle represents the Datong hydrometric station, and the yellow pentagrams indicate GNSS stations. The sky-blue and black lines represent the Three Gorges Reservoir region and the YRB, respectively. The left part of the study area is the sub-study area.</p>
Full article ">Figure 2
<p>Time series of reconstructed monthly/daily TWSA, detrended monthly TWSA from CSRM and detrended daily TWSA from ITSG−Grace2018 solutions in the study area from 2003 to 2020.</p>
Full article ">Figure 3
<p>Comparison of reconstructed daily TWSA, soil moisture, and the vertical displacement of eight GNSS stations during 2020. It should be noted that the axis of vertical displacements from GNSS are inverted for comparison.</p>
Full article ">Figure 4
<p>(<b>a</b>) Time series of reconstructed daily TWSA in the sub−study area and daily streamflow from Datong station from 2003 to 2020; (<b>b</b>) scatter graph. The X-axis represents streamflow, and the Y-axis represents TWSA.</p>
Full article ">Figure 5
<p>Spatial distribution of the reconstructed daily TWSA during the flooding period from 1 June to 3 September 2020.</p>
Full article ">Figure 6
<p>Spatial distribution of TWSA derived from CSRM from May to August 2020.</p>
Full article ">Figure 7
<p>Time series of reconstructed daily TWSA, soil moisture, daily precipitation, and accumulated precipitation of the study area during 1 June to 31 August 2020.</p>
Full article ">Figure 8
<p>Comparison of MWI, NDFPI, and WI (NDFPI is moved up by 0.1 for a better view) of the sub-study area and daily streamflow at Datong station in 2020. The dashed line in the figure represents the thresholds of the MWI, NDFPI, WI, and daily streamflow for the 90th percentile floods from 2003 to 2020.</p>
Full article ">Figure 9
<p>Comparison of NDFPI, WI, and MWI (NDFPI moves up by 0.1 for a better view) of the sub-study area and daily streamflow at Datong station in (<b>a</b>) 2016 for the 90th percentile floods, (<b>b</b>) 2020 for the 95th percentile floods, and (<b>c</b>) 2016 for the 95th percentile floods.</p>
Full article ">Figure 10
<p>Time series of deseasonalized and detrended reconstructed monthly/daily TWSA, monthly TWSA from CSRM, and daily TWSA from ITSG−Grace2018 solutions in the study area from 2003 to 2020.</p>
Full article ">
17 pages, 10813 KiB  
Article
A Framework for Retrieving Soil Organic Matter by Coupling Multi-Temporal Remote Sensing Images and Variable Selection in the Sanjiang Plain, China
by Haiyi Ma, Changkun Wang, Jie Liu, Xinyi Wang, Fangfang Zhang, Ziran Yuan, Chengshuo Yao and Xianzhang Pan
Remote Sens. 2023, 15(12), 3191; https://doi.org/10.3390/rs15123191 - 20 Jun 2023
Cited by 5 | Viewed by 1805
Abstract
Soil organic matter (SOM) is an important soil property for agricultural production. Rising grain demand has increased the intensity of cultivated land development in the Sanjiang Plain of China, and there is a strong demand for SOM monitoring in this region. Therefore, Baoqing [...] Read more.
Soil organic matter (SOM) is an important soil property for agricultural production. Rising grain demand has increased the intensity of cultivated land development in the Sanjiang Plain of China, and there is a strong demand for SOM monitoring in this region. Therefore, Baoqing County of the Sanjiang Plain, an important grain production area, was considered the study area. In the study, we proposed a framework for high-accuracy SOM retrieval by coupling multi-temporal remote sensing (RS) images and variable selection algorithms. A total of 73 surface soil samples (0–20 cm) were collected in 2010, and Landsat 5 images acquired during the bare soil period (April, May, and June) were selected from 2008 to 2011. Three variable selection algorithms, namely, Genetic Algorithm, Random Frog and Competitive Adaptive Reweighted Sampling (CARS), were combined with Partial Least Squares Regression (PLSR) to build SOM retrieval models on the spectral bands and indices of the images. The results using a single-date image showed that the combination of variable selection algorithms and PLSR outperformed using PLSR alone, and CARS showed the best performance (R2 = 0.34, RMSE = 15.66 g/kg) among all the algorithms. Therefore, only CARS was applied to SOM retrieval in the different year interval groups. To investigate the effect of the image acquisition time, all images were divided into various year interval groups, and the resulting images were then stacked. The results using multi-temporal images showed that the SOM retrieval accuracy improved as the year interval lengthened. The optimal result (R2 = 0.59, RMSE = 11.81 g/kg) was obtained from the 2008–2011 group, wherein the difference indices derived from the images of 2009, 2010, and 2011 dominated the selected spectral variables. Moreover, the spatial prediction of SOM based on the optimal model was consistent with the distribution of SOM. Our study suggested that the proposed framework that couples stacked multi-temporal RS images with variable selection algorithms has potential for SOM retrieval. Full article
(This article belongs to the Special Issue Remote Sensing for Soil Mapping and Monitoring)
Show Figures

Figure 1

Figure 1
<p>Overview of the study area and distribution of soil samples.</p>
Full article ">Figure 2
<p>Retrieval results of SOM for each single-date image. None: taking all variables as input variables without variable selection; GA: Genetic Algorithm; RF: Random Frog; CARS: Competitive Adaptive Reweighted Sampling.</p>
Full article ">Figure 3
<p>Scatter plots between measured vs. predicted SOM by the optimal model based on the multi-temporal images. Note: the dotted line is trend line and the solid line is 1:1 reference line.</p>
Full article ">Figure 4
<p>The variables selected by CARS on multi-temporal images in the 2008–2011 group.</p>
Full article ">Figure 5
<p>The SOM distribution in cultivated soils predicted by the model using the images in the 2008–2011 group.</p>
Full article ">
17 pages, 17478 KiB  
Article
Estimating the Evolution of a Post-Little Ice Age Deglaciated Alpine Valley through the DEM of Difference (DoD)
by Roberto Sergio Azzoni, Manuela Pelfini and Andrea Zerboni
Remote Sens. 2023, 15(12), 3190; https://doi.org/10.3390/rs15123190 - 19 Jun 2023
Cited by 2 | Viewed by 1624
Abstract
Since the end of the Little Ice Age (LIA, ~1830), the accelerated glaciers’ shrinkage along mid-latitude high mountain areas promoted a quick readjustment of geomorphological processes with the onset of the paraglacial dynamic, making proglacial areas among the most sensitive Earth’s landscapes to [...] Read more.
Since the end of the Little Ice Age (LIA, ~1830), the accelerated glaciers’ shrinkage along mid-latitude high mountain areas promoted a quick readjustment of geomorphological processes with the onset of the paraglacial dynamic, making proglacial areas among the most sensitive Earth’s landscapes to ongoing climate change. A potentially useful remote-sensing method for investigating such dynamic areas is the DEM (Digital Elevation Model) of Difference (DoD) technique, which quantifies volumetric changes in a territory between successive topographic surveys. After a detailed geomorphological analysis and comparison with historical maps of the Martello Valley (central Italian Alps), we applied the DoD for reconstructing post-LIA deglaciation dynamics and reported on the surface effects of freshly-onset paraglacial processes. The head of the valley is still glacierized, with three main ice bodies resulting from the huge reduction of a single glacier present at the apogee of the LIA. Aftermath: the glaciers lose 60% of their initial surface area, largely modifying local landforms and expanding the surface of the proglacial areas. The DoD analysis of the 2006–2015 timeframe (based on registered DEM derived from LiDAR—Light Detection and Ranging—data) highlights deep surface elevation changes ranging from +38 ± 4.01 m along the foot of rock walls, where gravitative processes increased their intensity, to −47 ± 4.01 m where the melting of buried ice caused collapses of the proglacial surface. This approach permits estimating the volume of sediments mobilized and reworked by paraglacial processes. Here, in less than 10 years, −23,675 ± 1165 m3 of sediment were removed along the proglacial area and transported down valley, highlighting the dynamicity of proglacial areas. Full article
(This article belongs to the Topic Cryosphere: Changes, Impacts and Adaptation)
Show Figures

Figure 1

Figure 1
<p>Location map of the Martello Valley in the eastern-central Italian Alps and its glacierized basin highlighted in red. The 2016 glacier limits of the whole Ortles-Cevedale area are represented in light blue (data from Paul et al. [<a href="#B34-remotesensing-15-03190" class="html-bibr">34</a>]).</p>
Full article ">Figure 2
<p>Evolution of the glacierized area (light blue) and proglacial area (yellow) at the head of Martello Valley. Glacier limits reconstructed after manual digitalization of available topographic maps and orthophotos. Toponyms in the last box refer to the New Italian Glacier Inventory [<a href="#B26-remotesensing-15-03190" class="html-bibr">26</a>].</p>
Full article ">Figure 3
<p>Simplified geomorphological sketch of glacial landforms and deposits of the Upper Martello Valley; see the text for discussion on landform distribution. The positions of some examples of landforms presented in <a href="#remotesensing-15-03190-f004" class="html-fig">Figure 4</a> are reported in black boxes.</p>
Full article ">Figure 4
<p>Some examples of landforms observed along the upper Martello Valley: (<b>A</b>) lateral moraine of the Forcola/Furkele Glacier that leans against the slope; (<b>B</b>) lateral moraine of the Cevedale/Zufall Glacier; (<b>C</b>) fluted moraine in the proglacial plain of the Forcola/Furkele; (<b>D</b>) sheepback rocks at the valley floor; (<b>E</b>) anastomosed pattern along the Plima stream; (<b>F</b>) braided pattern along the Plima stream.</p>
Full article ">Figure 5
<p>Glacier elevation changes in the timeframe 2005–2016.</p>
Full article ">Figure 6
<p>The 2016–2005 DEM of the difference of the proglacial area of Martello Valley with the main moraine ridges highlighted. The elevation profiles of four transects are reported, with their localization indicated on the map. The y-axis range of AA′ and BB′ is different from the CC′ and DD′ elevation profiles for better readability of the data.</p>
Full article ">Figure 7
<p>Detail of the 2016–2005 DEM of the difference in the proglacial area of the Cevedale/Zufall Glacier, reporting the main hydrographic features. The positions of some examples of processes presented in <a href="#remotesensing-15-03190-f008" class="html-fig">Figure 8</a> and <a href="#remotesensing-15-03190-f009" class="html-fig">Figure 9</a> are reported in black boxes.</p>
Full article ">Figure 8
<p>Main geomorphic processes that occur in the proglacial plain of Zufall/Cevedale Glacier, promoting wide topographic variation: (<b>A</b>) moraine degradation, (<b>B</b>) slope processes (the line indicates the moraine edge), and (<b>C</b>) fluvial reworking of glacial deposits.</p>
Full article ">Figure 9
<p>(<b>A</b>–<b>C</b>) Formation and evolution of a proglacial lake in the deglaciated area of Cevedale/Zufall Glacier in the 2000–2020 timeframe from aerial orthophotos (the dashed blue line indicates the 2020 lake limit). (<b>D</b>,<b>E</b>) field images of the proglacial lakes enclosed by an Eighties moraine ridge marked with the purple line.</p>
Full article ">
19 pages, 19495 KiB  
Article
A Self-Supervised Learning Approach for Extracting China Physical Urban Boundaries Based on Multi-Source Data
by Yuan Tao, Wanzeng Liu, Jun Chen, Jingxiang Gao, Ran Li, Jiaxin Ren and Xiuli Zhu
Remote Sens. 2023, 15(12), 3189; https://doi.org/10.3390/rs15123189 - 19 Jun 2023
Cited by 5 | Viewed by 1930
Abstract
Physical urban boundaries (PUBs) are basic geographic information data for defining the spatial extent of urban landscapes with non-agricultural land and non-agricultural economic activities. Accurately mapping PUBs provides a spatiotemporal database for urban dynamic monitoring, territorial spatial planning, and ecological environment protection. However, [...] Read more.
Physical urban boundaries (PUBs) are basic geographic information data for defining the spatial extent of urban landscapes with non-agricultural land and non-agricultural economic activities. Accurately mapping PUBs provides a spatiotemporal database for urban dynamic monitoring, territorial spatial planning, and ecological environment protection. However, traditional extraction methods often have problems, such as subjective parameter settings and inconsistent cartographic scales, making it difficult to identify PUBs objectively and accurately. To address these problems, we proposed a self-supervised learning approach for PUB extraction. First, we used nighttime light and OpenStreetMap road data to map the initial urban boundary for data preparation. Then, we designed a pretext task of self-supervised learning based on an unsupervised mutation detection algorithm to automatically mine supervised information in unlabeled data, which can avoid subjective human interference. Finally, a downstream task was designed as a supervised learning task in Google Earth Engine to classify urban and non-urban areas using impervious surface density and nighttime light data, which can solve the scale inconsistency problem. Based on the proposed method, we produced a 30 m resolution China PUB dataset containing six years (i.e., 1995, 2000, 2005, 2010, 2015, and 2020). Our PUBs show good agreement with existing products and accurately describe the spatial extent of urban areas, effectively distinguishing urban and non-urban areas. Moreover, we found that the gap between the national per capita GDP and the urban per capita GDP is gradually decreasing, but regional coordinated development and intensive development still need to be strengthened. Full article
Show Figures

Figure 1

Figure 1
<p>Flowchart of producing CPUB (NTL means the nighttime light data; OSM means OpenStreetMap; ISA means the impervious surface).</p>
Full article ">Figure 2
<p>Schematic diagram explaining the scale inconsistency problem caused by directly connecting mutation points for different strategies (black lines: sampling lines; blue lines: mapped boundaries; red dots: detected mutation points). (<b>a</b>) Connection of mutation points obtained by radial lines; (<b>b</b>) connection of mutation points obtained by a fan-shaped area-constrained method; (<b>c</b>) connection of mutation points obtained by a grid.</p>
Full article ">Figure 3
<p>The spatial extent of CPUB urban areas in the Beijing–Tianjin region and Hefei, Anhui Province, in 2020 (Subfigures (<b>a</b>–<b>f</b>) show the local details of the urban-scale boundaries in the Beijing–Tianjin region. Subfigures (<b>g</b>–<b>l</b>) show the local details of the building-scale boundaries in Hefei. Basemap: High-resolution Google Earth images).</p>
Full article ">Figure 4
<p>PUBs at multiple periods and the initial urban boundaries of Chengdu (Sichuan Province), Hefei (Anhui Province), Wuhan (Hubei Province), and Xi’an (Shaanxi Province) (Red region means the urban extent. Basemap: Landsat; the natural color band combination strategy. Initial UB means the initial urban boundary).</p>
Full article ">Figure 5
<p>Comparison of urban boundaries derived from CPUB and the GUB in Hefei and Xi’an (Subfigures (<b>a</b>–<b>f</b>) show the local details of Hefei. Subfigures (<b>g</b>–<b>l</b>) show the local details of Xi’an. Basemap: Google Earth images).</p>
Full article ">Figure 6
<p>The size rank relationship of Chinese cities (an area greater than 2 km<sup>2</sup>) from 1995 to 2020.</p>
Full article ">Figure 7
<p>Zipf’s distribution of Chinese cities from 1995 to 2020.</p>
Full article ">Figure 8
<p>The urban areas of China in the longitude and latitude directions from 1995 to 2020 (interval: 10 km; basemap: OpenStreetMap).</p>
Full article ">Figure 9
<p>Impervious surface areas and urban areas and their ratios in different regions from 1995 to 2020.</p>
Full article ">Figure 10
<p>Scatter distribution of urban area and GDP at the province level and its fitting line.</p>
Full article ">Figure 11
<p>Distribution of urban intensification in each province of China (blue line: national per capita GDP; red line: province-level average per urban capita GDP).</p>
Full article ">
20 pages, 10226 KiB  
Article
A Novel Optimization Strategy of Sidelobe Suppression for Pulse Compression Weather Radar
by Jiaqi Hu, Xichao Dong, Weiming Tian, Cheng Hu, Kai Feng and Jun Lu
Remote Sens. 2023, 15(12), 3188; https://doi.org/10.3390/rs15123188 - 19 Jun 2023
Cited by 6 | Viewed by 2051
Abstract
The solid-state transmitters are widely adopted for weather radars, where pulse compression is operated to provide the required sensitivity and range resolution. Therefore, effective sidelobe suppression strategies must be employed, especially for weather observation. Currently, many methods can suppress the sidelobe to a [...] Read more.
The solid-state transmitters are widely adopted for weather radars, where pulse compression is operated to provide the required sensitivity and range resolution. Therefore, effective sidelobe suppression strategies must be employed, especially for weather observation. Currently, many methods can suppress the sidelobe to a very low level in the case of point targets or uniformly distributed targets. However, in strong convection weather process, the weather echo amplitude lies in a wide dynamic range and the main lobe of weak target is prone to being contaminated by the sidelobe of strong target, causing the degradation of weather fundamental data estimation, even generating artifacts and affecting the quantitative precipitation evaluation. In this paper, we propose a novel strategy which is the further processing of a general pulse compression radar to mitigate the effects of sidelobes. The proposed method is called the predominant component extraction (PCE), in which the re-weighting processing is operated after pulse compression, and then the echo of each bin is optimized and its energy will approach the real targets in each bin. It can improve the estimation of weak signals or even eliminate the artifact at the edge of the scene. Numerical simulation experiments and real-data verifications are implemented to validate the feasibility and superiority. It is noted that the proposed method has no requirement on the transmitted waveform and can be realized only by adding a step after pulse compression in the actual system. Full article
(This article belongs to the Special Issue Processing and Application of Weather Radar Data)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>The schematic diagram of PC-WTF.</p>
Full article ">Figure 2
<p>The schematic diagram of received echoes.</p>
Full article ">Figure 3
<p>The block diagram of the PCE algorithm.</p>
Full article ">Figure 4
<p>The pulse compression results of point target.</p>
Full article ">Figure 5
<p>The input Z with different gradients. (<b>a</b>) Gradient of Z is 15 dBZ/km; (<b>b</b>) gradient of Z is 25 dBZ/km; (<b>c</b>) gradient of Z is 40 dBZ/km.</p>
Full article ">Figure 6
<p>The reflectivity calculated from unoptimized echoes. (<b>a</b>) Gradient of Z is 15 dBZ/km; (<b>b</b>) gradient of Z is 25 dBZ/km; (<b>c</b>) gradient of Z is 40 dBZ/km.</p>
Full article ">Figure 7
<p>The mean velocity calculated from unoptimized echoes. (<b>a</b>) Gradient of Z is 15 dBZ/km; (<b>b</b>) gradient of Z is 25 dBZ/km; (<b>c</b>) gradient of Z is 40 dBZ/km.</p>
Full article ">Figure 8
<p>The spectral width calculated from unoptimized echoes. (<b>a</b>) Gradient of Z is 15 dBZ/km; (<b>b</b>) gradient of Z is 25 dBZ/km; (<b>c</b>) gradient of Z is 40 dBZ/km.</p>
Full article ">Figure 9
<p>The comparison results of reflectivity. (<b>a</b>) Gradient of Z is 15 dBZ/km; (<b>b</b>) gradient of Z is 25 dBZ/km; (<b>c</b>) gradient of Z is 40 dBZ/km.</p>
Full article ">Figure 10
<p>The comparison results of mean velocity. (<b>a</b>) Gradient of Z is 15 dBZ/km; (<b>b</b>) gradient of Z is 25 dBZ/km; (<b>c</b>) gradient of Z is 40 dBZ/km.</p>
Full article ">Figure 11
<p>The comparison results of spectral width. (<b>a</b>) Gradient of Z is 15 dBZ/km; (<b>b</b>) gradient of Z is 25 dBZ/km; (<b>c</b>) gradient of Z is 40 dBZ/km.</p>
Full article ">Figure 12
<p>The RMSE of ZVW changing with the gradient of reflectivity. (<b>a</b>) Reflectivity; (<b>b</b>) mean velocity; (<b>c</b>) spectral width.</p>
Full article ">Figure 13
<p>The AR of ZVW changing with the gradient of reflectivity. (<b>a</b>) Reflectivity; (<b>b</b>) mean velocity; (<b>c</b>) spectral width.</p>
Full article ">Figure 14
<p>The SNR loss changing with reflectivity gradient.</p>
Full article ">Figure 15
<p>The ZVW of the supercell. ((<b>a</b>) Reflectivity; (<b>b</b>) velocity; (<b>c</b>) spectral width).</p>
Full article ">Figure 16
<p>The results after optimization by PCE. ((<b>a</b>) Reflectivity; (<b>b</b>) velocity; (<b>c</b>) spectral width).</p>
Full article ">Figure 17
<p>Local reflectivity before and after optimization. ((<b>a</b>) Not optimized; (<b>b</b>) optimized).</p>
Full article ">Figure 18
<p>Selected region for QPE.</p>
Full article ">Figure 19
<p>The comparison results of precipitations. (<b>a</b>) Position A; (<b>b</b>) position B.</p>
Full article ">Figure 20
<p>The estimation errors of precipitations. (<b>a</b>) Position A; (<b>b</b>) position B.</p>
Full article ">
14 pages, 2514 KiB  
Communication
Sensitivity of Grassland Coverage to Climate across Environmental Gradients on the Qinghai-Tibet Plateau
by Rihan Wu, Guozheng Hu, Hasbagan Ganjurjav and Qingzhu Gao
Remote Sens. 2023, 15(12), 3187; https://doi.org/10.3390/rs15123187 - 19 Jun 2023
Cited by 3 | Viewed by 1433
Abstract
Grassland cover is strongly influenced by climate change. The response of grassland cover to climate change becomes complex with background climate. There have been some advances in research on the sensitivity of grassland vegetation to climate change around the world, but the differences [...] Read more.
Grassland cover is strongly influenced by climate change. The response of grassland cover to climate change becomes complex with background climate. There have been some advances in research on the sensitivity of grassland vegetation to climate change around the world, but the differences in climate sensitivity among grassland types are still unclear in alpine grassland. Therefore, we applied MODIS NDVI data and trend analysis methods to quantify the spatial and temporal variation of grassland vegetation cover on the Qinghai-Tibet Plateau. Then, we used multiple regression models to analyze the sensitivity of fractional vegetation cover (FVC) to climatic factors (Temperature, Precipitation, Solar radiation, Palmer drought severity index) and summarized the potential mechanisms of vegetation sensitivity to different climatic gradients. Our results showed (1) a significant increasing trend in alpine desert FVC from 2000–2018 (1.12 × 10−3/a, R2 = 0.56, p < 0.001) but no significant trend in other grassland types. (2) FVC sensitivity to climatic factors varied among grassland types, especially in the alpine desert, which had over 60% of the area with positive sensitivity to temperature, precipitation and PDSI. (3) The sensitivity of grassland FVC to heat factors decreases with rising ambient temperature while the sensitivity to moisture increases. Similarly, the sensitivity to moisture decreases while the sensitivity to thermal factors increases along the moisture gradient. Furthermore, the results suggest that future climate warming will promote grassland in cold and wet areas of the Qinghai-Tibet Plateau and may suppress vegetation in warmer areas. In contrast, the response of the alpine desert to future climate is more stable. Studying the impact of climate variation at a regional scale could enhance the adaptability of vegetation in future global climates. Full article
Show Figures

Figure 1

Figure 1
<p>Qinghai-Tibet Plateau Grassland Classification.</p>
Full article ">Figure 2
<p>Spatial distribution of (<b>a</b>) mean FVC and (<b>b</b>) FVC trends in grasslands during 2000–2018.</p>
Full article ">Figure 3
<p>Interannual FVC variation in different grassland types during 2000–2018.TP: total grassland of the QTP; AM, alpine meadow; AS, alpine steppe; AD, alpine desert.</p>
Full article ">Figure 4
<p>Proportions of areas displaying differing degrees of sensitivity for grassland FVC as related to mean annual temperature (<b>a</b>,<b>b</b>), mean annual Srad (<b>c</b>,<b>d</b>), annual precipitation (<b>e</b>,<b>f</b>) and mean annual PDSI (<b>g</b>,<b>h</b>). Insets in the upper right corners of the graphs represent the spatial distribution of significant sensitivity as follows; red, positive and blue, negative. TP: total grasslands of the QTP; AM, alpine meadow; AS, alpine steppe; AD, alpine desert; SN, significant negative; NSN, not significant negative; NSP, not significant positive and SP, significant positive.</p>
Full article ">Figure 5
<p>Distribution of sensitivity of grasslands to temperature (<b>a</b>–<b>d</b>), solar radiation (<b>e</b>–<b>h</b>), precipitation (<b>i</b>–<b>l</b>) and drought level (<b>m</b>–<b>p</b>) with gradients of temperature, Srad, precipitation and PDSI. (NS, not significant, TP: total grassland of the QTP (in green); AM, alpine meadow (in red); AS, alpine steppe (in blue); AD, alpine desert (in yellow)).</p>
Full article ">
31 pages, 2896 KiB  
Article
Accuracy Assessment of Eleven Medium Resolution Global and Regional Land Cover Land Use Products: A Case Study over the Conterminous United States
by Zhixin Wang and Giorgos Mountrakis
Remote Sens. 2023, 15(12), 3186; https://doi.org/10.3390/rs15123186 - 19 Jun 2023
Cited by 11 | Viewed by 1883
Abstract
Land cover land use (LCLU) products provide essential information for numerous environmental and human studies. Here, we assess the accuracy of eleven global and regional products over the conterminous U.S. using 25,000 high-confidence randomly distributed samples. Results show that in general, the National [...] Read more.
Land cover land use (LCLU) products provide essential information for numerous environmental and human studies. Here, we assess the accuracy of eleven global and regional products over the conterminous U.S. using 25,000 high-confidence randomly distributed samples. Results show that in general, the National Land Cover Database (NLCD) and the Land Change Monitoring, Assessment and Projection (LCMAP) outperform other multi-class products, both in terms of higher individual class accuracy and with accuracy variability across classes. More specifically, F1 accuracy comparisons between the best performing USGS and non-USGS products indicate: (i) similar performance for the water class, (ii) USGS product outperformance in the developed (+1.3%), grass/shrub (+3.2%) and tree cover (+4.2%) classes, and (iii) non-USGS product (WorldCover) gains in the cropland (+5.1%) class. The NLCD and LCMAP also outperformed specialized single-class products, such as the Hansen Global Forest Change, the Cropland Data Layer and the Global Artificial Impervious Areas, while offering comparable results to the Global Surface Water Dynamics product. Spatial visualizations also allowed accuracy comparisons across different geographic areas. In general, the NLCD and LCMAP have disagreements mainly in the middle and southeastern part of conterminous U.S. while Esri, WorldCover and Dynamic World have most errors in the western U.S. Comparisons were also undertaken on a subset of the reference data, called spatial edge samples, that identifies samples surrounded by neighboring samples of different class labels, thus excluding easy-to-classify homogenous areas. There, the WorldCover product offers higher accuracies for the highly dynamic grass/shrub (+4.4%) and cropland (+8.1%) classes when compared to the NLCD and LCMAP products. An important conclusion while looking at these challenging samples is that except for the tree class (78%), the best performing products per class range in accuracy between 55% and 70%, which suggests that there is substantial room for improvement. Full article
(This article belongs to the Section Earth Observation Data)
Show Figures

Figure 1

Figure 1
<p>Workflow for statistical and spatial accuracy assessments. Six 30 m products are GlobeLand30, FROM-GLC, GSWD, GAIA, HGFC and CDL.</p>
Full article ">Figure 2
<p>Five climatic zones for the conterminous United States.</p>
Full article ">Figure 3
<p>F1 score of HGFC with thresholds for year 2000 (<b>left</b>) and 2012 (<b>right</b>).</p>
Full article ">Figure 4
<p>F1 score of GSWD with different thresholds for years from 2000 to 2015.</p>
Full article ">Figure 5
<p>Spatial agreement percentage maps of LCMAP, NLCD, ESRI, World Cover and Dynamic World for recent years using developed, cropland, grass/shrub, tree cover and water classes.</p>
Full article ">Figure 6
<p>The average annual number of non-fill non-cloudy observations at each CONUS ARD 30 m pixel location over the 36-year study period (1982–2017) for Landsat 7 ETM+. Adapted from <a href="#remotesensing-15-03186-f006" class="html-fig">Figure 6</a> of [<a href="#B78-remotesensing-15-03186" class="html-bibr">78</a>].</p>
Full article ">
26 pages, 11502 KiB  
Article
A Global Structure and Adaptive Weight Aware ICP Algorithm for Image Registration
by Lin Cao, Shengbin Zhuang, Shu Tian, Zongmin Zhao, Chong Fu, Yanan Guo and Dongfeng Wang
Remote Sens. 2023, 15(12), 3185; https://doi.org/10.3390/rs15123185 - 19 Jun 2023
Cited by 4 | Viewed by 2588
Abstract
As an important technology in 3D vision, point-cloud registration has broad development prospects in the fields of space-based remote sensing, photogrammetry, robotics, and so on. Of the available algorithms, the Iterative Closest Point (ICP) algorithm has been used as the classic algorithm for [...] Read more.
As an important technology in 3D vision, point-cloud registration has broad development prospects in the fields of space-based remote sensing, photogrammetry, robotics, and so on. Of the available algorithms, the Iterative Closest Point (ICP) algorithm has been used as the classic algorithm for solving point cloud registration. However, with the point cloud data being under the influence of noise, outliers, overlapping values, and other issues, the performance of the ICP algorithm will be affected to varying degrees. This paper proposes a global structure and adaptive weight aware ICP algorithm (GSAW-ICP) for image registration. Specifically, we first proposed a global structure mathematical model based on the reconstruction of local surfaces using both the rotation of normal vectors and the change in curvature, so as to better describe the deformation of the object. The model was optimized for the convergence strategy, so that it had a wider convergence domain and a better convergence effect than either of the original point-to-point or point-to-point constrained models. Secondly, for outliers and overlapping values, the GSAW-ICP algorithm was able to assign appropriate weights, so as to optimize both the noise and outlier interference of the overall system. Our proposed algorithm was extensively tested on noisy, anomalous, and real datasets, and the proposed method was proven to have a better performance than other state-of-the-art algorithms. Full article
Show Figures

Figure 1

Figure 1
<p>GSAW-ICP algorithm flow chart.</p>
Full article ">Figure 2
<p>Ideal case of point cloud alignment.</p>
Full article ">Figure 3
<p>Find the target point graphically.</p>
Full article ">Figure 4
<p>Iterative alignment schematic.</p>
Full article ">Figure 5
<p>Schematic diagram of uniform sampling and normal vector sampling.</p>
Full article ">Figure 6
<p>Schematic diagram of initial value iteration.</p>
Full article ">Figure 7
<p>Schematic diagram of normal vector projection.</p>
Full article ">Figure 8
<p>Schematic diagram of exception handling.</p>
Full article ">Figure 9
<p>Schematic diagram of large differences in feature constraints.</p>
Full article ">Figure 10
<p>The relationship between the root mean squared error of the registration effect and the signal-to-noise ratio. From (<b>a</b>–<b>d</b>), the number of iterations increases by 10, while (<b>a</b>) has 50 iterations.</p>
Full article ">Figure 11
<p>The relationship between the root mean squared error of the registration effect and the signal-to-noise ratio. From (<b>a</b>–<b>d</b>), the noise interference increases by 5 db in sequence, while (<b>a</b>) has a noise interference of 0 db.</p>
Full article ">Figure 12
<p>Image to be registered.</p>
Full article ">Figure 13
<p>Image registration under different rotation angles.</p>
Full article ">Figure 14
<p>Remote sensing registration effect detection.</p>
Full article ">Figure 15
<p>Point cloud map of the original site 1.</p>
Full article ">Figure 16
<p>Point cloud map of the original site 2.</p>
Full article ">Figure 17
<p>Point cloud map of the original site 3.</p>
Full article ">Figure 18
<p>Raw point cloud 1TruSlicer image.</p>
Full article ">Figure 19
<p>Raw point cloud 2TruSlicer image.</p>
Full article ">Figure 20
<p>Raw point cloud 3TruSlicer image.</p>
Full article ">Figure 21
<p>Original 1.2 site point cloud registration map.</p>
Full article ">Figure 22
<p>Original 1.2 point cloud TruSlicer image registration.</p>
Full article ">Figure 23
<p>Registration of original 2 and 3 site point cloud images.</p>
Full article ">Figure 24
<p>Registration of original 2 and 3 point cloud TruSlicer images.</p>
Full article ">Figure 25
<p>Registration of original 1, 2, 3 site point cloud images.</p>
Full article ">Figure 26
<p>Original 1, 2, 3 point cloud TruSlicer image registration.</p>
Full article ">
27 pages, 13888 KiB  
Article
A New Blind Selection Approach for Lunar Landing Zones Based on Engineering Constraints Using Sliding Window
by Hengxi Liu, Yongzhi Wang, Shibo Wen, Jianzhong Liu, Jiaxiang Wang, Yaqin Cao, Zhiguo Meng and Yuanzhi Zhang
Remote Sens. 2023, 15(12), 3184; https://doi.org/10.3390/rs15123184 - 19 Jun 2023
Cited by 1 | Viewed by 1829
Abstract
Deep space exploration has risen in interest among scientists in recent years, with soft landings being one of the most straightforward ways to acquire knowledge about the Moon. In general, landing mission success depends on the selection of landing zones, and there are [...] Read more.
Deep space exploration has risen in interest among scientists in recent years, with soft landings being one of the most straightforward ways to acquire knowledge about the Moon. In general, landing mission success depends on the selection of landing zones, and there are currently few effective quantitative models that can be used to select suitable landing zones. When automatic landing zones are selected, the grid method used for data partitioning tends to miss potentially suitable landing sites between grids. Therefore, this study proposes a new engineering-constrained approach for landing zone selection using LRO LOLA-based slope data as original data based on the sliding window method, which solves the spatial omission problem of the grid method. Using the threshold ratio, mean, coefficient of variation, Moran’s I, and overall rating, this method quantifies the suitability of each sliding window. The k-means clustering algorithm is adopted to determine the suitability threshold for the overall rating. The results show that 20 of 22 lunar soft landing sites are suitable for landing. Additionally, 43 of 50 landing sites preselected by the experts (suitable landing sites considering a combination of conditions) are suitable for landing, accounting for 90.9% and 86% of the total number, respectively, for a window size of 0.5° × 0.5°. Among them, there are four soft landing sites: Surveyor 3, 6, 7, and Apollo 15, which are not suitable for landing in the evaluation results of the grid method. However, they are suitable for landing in the overall evaluation results of the sliding window method, which significantly reduces the spatial omission problem of the grid method. In addition, four candidate landing regions, including Aristarchus Crater, Marius Hills, Moscoviense Basin, and Orientale Basin, were evaluated for landing suitability using the sliding window method. The suitability of the landing area within the candidate range of small window sizes was 0.90, 0.97, 0.49, and 0.55. This indicates the capacity of the method to analyze an arbitrary range during blind landing zone selection. The results can quantify the slope suitability of the landing zones from an engineering perspective and provide different landing window options. The proposed method for selecting lunar landing zones is clearly superior to the gridding method. It enhances data processing for automatic lunar landing zone selection and progresses the selection process from qualitative to quantitative. Full article
Show Figures

Figure 1

Figure 1
<p>Schematic diagram of slope calculation. A certain pixel “e” is usually combined with eight adjacent pixels to form a slope, a~i stands for the elevation of each pixel.</p>
Full article ">Figure 2
<p>Slope of the whole Moon with a horizontal resolution of 512 pixels per degree (60 m at the equator) and a typical vertical accuracy of 3 to 4 m.</p>
Full article ">Figure 3
<p>The basic theory of sliding window algorithm. ① and ② represent the sliding window cyclic decision condition, ③ represents a single sliding window containing information and processing, ④ demonstrates the overlapping characteristics of sliding windows, ⑤ represents the completion of the sliding cycle of the original data.</p>
Full article ">Figure 4
<p>Spatial units based on the Queen proximity rule. “A” represents the spatial unit being studied, “B” represents the neighboring spatial unit for “A” under the Queen’s proximity rule, which contains the common vertex connection and common adjacent edge connection.</p>
Full article ">Figure 5
<p>The results of the evaluation of the whole Moon slope based on sliding windows combined with multiple indicators are shown in Figure (<b>a</b>–<b>f</b>) with a window size of 1° × 1° and in Figure (<b>e</b>) with a window size of 0.5° × 0.5°. (<b>a</b>) Results for the ratio of the whole Moon slope with a threshold of 20. (<b>b</b>) Results for the ratio of the whole Moon slope with a threshold of 8. (<b>c</b>) Results for the coefficient of variation of the whole Moon slope. (<b>d</b>) Results for the mean of the whole Moon slope. (<b>e</b>) Results for the whole Moon overall rating. (<b>f</b>) Results for the 0.5° × 0.5° whole Moon overall rating. Green snowflakes are expertly preselected landing sites, and red asterisks are available soft landing sites. Considering that the lower mean value reflects the overall flatness of the slope and the lower coefficient of variation reflects the dispersion of slope data, the color bars of results in <a href="#remotesensing-15-03184-f005" class="html-fig">Figure 5</a>c,d were flipped (low value (blue) to high value (yellow) flip to high value (yellow) to low value (yellow) from top to bottom) to ensure consistency in the results and improve the overall presentation.</p>
Full article ">Figure 5 Cont.
<p>The results of the evaluation of the whole Moon slope based on sliding windows combined with multiple indicators are shown in Figure (<b>a</b>–<b>f</b>) with a window size of 1° × 1° and in Figure (<b>e</b>) with a window size of 0.5° × 0.5°. (<b>a</b>) Results for the ratio of the whole Moon slope with a threshold of 20. (<b>b</b>) Results for the ratio of the whole Moon slope with a threshold of 8. (<b>c</b>) Results for the coefficient of variation of the whole Moon slope. (<b>d</b>) Results for the mean of the whole Moon slope. (<b>e</b>) Results for the whole Moon overall rating. (<b>f</b>) Results for the 0.5° × 0.5° whole Moon overall rating. Green snowflakes are expertly preselected landing sites, and red asterisks are available soft landing sites. Considering that the lower mean value reflects the overall flatness of the slope and the lower coefficient of variation reflects the dispersion of slope data, the color bars of results in <a href="#remotesensing-15-03184-f005" class="html-fig">Figure 5</a>c,d were flipped (low value (blue) to high value (yellow) flip to high value (yellow) to low value (yellow) from top to bottom) to ensure consistency in the results and improve the overall presentation.</p>
Full article ">Figure 6
<p>K value clustering deviation map for different K values with different window sizes. The red circles represent the within-cluster sum of squares at the corresponding K values. (<b>a</b>) 1° × 1° window size. (<b>b</b>) 0.5° × 0.5° window size.</p>
Full article ">Figure 7
<p>Comparison of the whole Moon evaluation results of the sliding window and grid methods. (<b>a</b>) The whole Moon evaluation results of the grid algorithm are presented in the form of surface elements. (<b>b</b>) Whole Moon evaluation results of the sliding window algorithm in the form of surface elements. The black box represents the area selected for local comparison. For example, some areas with more landing points were selected to compare the differences between the grid and sliding window methods. The red asterisk represents the soft landing sites.</p>
Full article ">Figure 7 Cont.
<p>Comparison of the whole Moon evaluation results of the sliding window and grid methods. (<b>a</b>) The whole Moon evaluation results of the grid algorithm are presented in the form of surface elements. (<b>b</b>) Whole Moon evaluation results of the sliding window algorithm in the form of surface elements. The black box represents the area selected for local comparison. For example, some areas with more landing points were selected to compare the differences between the grid and sliding window methods. The red asterisk represents the soft landing sites.</p>
Full article ">Figure 8
<p>WAC image maps of five local candidate regions (left column), grid method results (middle column), and sliding window method results (right column). The red asterisks are soft landing sites. (<b>a</b>) Sinus Iridum and Jura Mountains in northwestern Mare Imbrium. (<b>b</b>) Mare Vaporum and vicinity. (<b>c</b>) Mare Crisium and vicinity. (<b>d</b>) Tycho and vicinity. (<b>e</b>) Von Kármán impact crater and vicinity. According to the clustering results, the blue areas ranging from 0.81 to 1 represent suitable landing areas, and the green and red areas represent the quantitative results of less suitable and unsuitable landings, respectively.</p>
Full article ">Figure 8 Cont.
<p>WAC image maps of five local candidate regions (left column), grid method results (middle column), and sliding window method results (right column). The red asterisks are soft landing sites. (<b>a</b>) Sinus Iridum and Jura Mountains in northwestern Mare Imbrium. (<b>b</b>) Mare Vaporum and vicinity. (<b>c</b>) Mare Crisium and vicinity. (<b>d</b>) Tycho and vicinity. (<b>e</b>) Von Kármán impact crater and vicinity. According to the clustering results, the blue areas ranging from 0.81 to 1 represent suitable landing areas, and the green and red areas represent the quantitative results of less suitable and unsuitable landings, respectively.</p>
Full article ">Figure 9
<p>WAC images of the four candidate landing sites and the combined index results at different window sizes: first column (WAC images), second column (evaluation results at 1° × 1° window size), third column (evaluation results at 0.5° × 0.5° window size), and fourth column (evaluation results at 0.03125° × 0.03125° window size). (<b>a</b>) Candidate landing site, Aristarchus Crater. (<b>b</b>) Marius Hills. (<b>c</b>) Moscoviense Basin. (<b>d</b>) Orientale Basin.</p>
Full article ">Figure 9 Cont.
<p>WAC images of the four candidate landing sites and the combined index results at different window sizes: first column (WAC images), second column (evaluation results at 1° × 1° window size), third column (evaluation results at 0.5° × 0.5° window size), and fourth column (evaluation results at 0.03125° × 0.03125° window size). (<b>a</b>) Candidate landing site, Aristarchus Crater. (<b>b</b>) Marius Hills. (<b>c</b>) Moscoviense Basin. (<b>d</b>) Orientale Basin.</p>
Full article ">
17 pages, 5727 KiB  
Article
The Transport Path and Vertical Structure of Dust Storms in East Asia and the Impacts on Cities in Northern China
by Tana Bao, Guilin Xi, Yanling Hao, I-Shin Chang, Jing Wu, Zhichao Xue, Erdemtu Jin, Wenxing Zhang and Yuhai Bao
Remote Sens. 2023, 15(12), 3183; https://doi.org/10.3390/rs15123183 - 19 Jun 2023
Viewed by 1894
Abstract
Dust storm disasters have emerged as a significant environmental challenge in East Asia. However, relying on a single monitoring method to track dust storms presents limitations and can be variable. Therefore, it is necessary to use a combination of ground and remote sensing [...] Read more.
Dust storm disasters have emerged as a significant environmental challenge in East Asia. However, relying on a single monitoring method to track dust storms presents limitations and can be variable. Therefore, it is necessary to use a combination of ground and remote sensing monitoring methods to explore the source and impact range of dust storms in order to fully characterize them. To achieve this, we examined the sources and impact ranges of dust storms in East Asia from 1980 to 2020 using both ground station data and remote sensing data. In addition, we focused on three specific dust storm events in the region. Our results indicate that the central source areas of dust storms are located in southern Mongolia and the Taklamakan Desert in China. Dust storms are mainly transported and spread in the northwestern region, while they are relatively rare in the southeastern region. The HYSPLIT model simulations reveal that the primary source directions of dust storms in East Asia are northwest, west, and north, the region involved includes Kazakhstan, southern Mongolia, and the Taklimakan Desert in China. The vertical structure of the dust storm layer depends on the source of the dust storm and the intensity of the dust storm event. Dust grain stratification probably occurs due to differences in dust storm sources, grain size, and regularity. These findings demonstrate that a combination of ground-based and remote sensing monitoring methods is an effective approach to fully characterize dust storms and can provide more comprehensive information for dust storm studies. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Location of the study area.</p>
Full article ">Figure 2
<p>Spatial and temporal distribution of dust outbreaks.</p>
Full article ">Figure 3
<p>The trends of DO variations in 1980–2020. Note: Green circles represent decreasing trends of DO and red triangles represent increasing trends of DO. Only trends with a statistical significance &gt;95% are presented using colors.</p>
Full article ">Figure 4
<p>Seventy-two hours of backward trajectory (<b>a</b>,<b>c</b>,<b>e</b>,<b>g</b>) and trajectory clustering (<b>b</b>,<b>d</b>,<b>f</b>,<b>h</b>) of dust storms. The light yellow line indicates the airflow backward trajectory of each dust storms, purple, red, green, blue for clustering results.</p>
Full article ">Figure 5
<p>HYSPLIT simulation of dust transport path on May 5, 2015, and dust events station monitored, including flowing dust (FD), blowing dust (BD), and dust storm (DS). The red square is Yumen station.</p>
Full article ">Figure 6
<p>On 5 May 2015, at 20:00 UTC, CALIPSO satellite 532 nm total attenuation backscatter coefficient, depolarization ratio, color ratio vertical profile and running trajectory. (<b>a</b>) Aerosol classification, the red line is satellite orbital paths of the dust storm events; (<b>b</b>) attenuated radar color ratio; (<b>c</b>) polarization ratio and (<b>d</b>) Total attenuated backscatter coefficient in 532 nm channel (Unit: km<sup>−1</sup> ·sr<sup>−1</sup>).</p>
Full article ">Figure 7
<p>HYSPLIT simulation of dust transport path on 12 April 2018, and dust events monitored at the station, including flowing dust (FD), blowing dust (BD), and dust storm (DS). The red square is Erlianhot station.</p>
Full article ">Figure 8
<p>On 12 April 2018, at 5:00 UTC, CALIPSO satellite 532 nm total attenuation backscatter coefficient, depolarization ratio, color ratio vertical profile and running trajectory. (<b>a</b>) Aerosol classification, the red line is satellite orbital paths of the dust storm events; (<b>b</b>) attenuated radar color ratio; (<b>c</b>) polarization ratio and (<b>d</b>) Total attenuated backscatter coefficient in 532 nm channel (Unit: km<sup>−1</sup> … sr<sup>−1</sup>).</p>
Full article ">Figure 9
<p>HYSPLIT simulation of dust transport path on 15 March 2021, and dust events monitored at the station, including flowing dust (FD), blowing dust (BD), and dust storm (DS). The red square is Guaizihu station.</p>
Full article ">Figure 10
<p>On 15 March 2021, at 5:00 UTC, CALIPSO satellite 532 nm total attenuation backscatter coefficient, depolarization ratio, color ratio vertical profile and running trajectory. (<b>a</b>) Aerosol classification, the red line is satellite orbital paths of the dust storm events; (<b>b</b>) attenuated radar color ratio; (<b>c</b>) polarization ratio and (<b>d</b>) Total attenuated backscatter coefficient in 532 nm channel (Unit: km<sup>−1</sup> … sr<sup>−1</sup>).</p>
Full article ">
5 pages, 2752 KiB  
Correction
Correction: Sesnie et al. In-Situ and Remote Sensing Platforms for Mapping Fine-Fuels and Fuel-Types in Sonoran Semi-Desert Grasslands. Remote Sens. 2018, 10, 1358
by Steven E. Sesnie, Holly Eagleston, Lacrecia Johnson and Emily Yurcich
Remote Sens. 2023, 15(12), 3182; https://doi.org/10.3390/rs15123182 - 19 Jun 2023
Viewed by 842
Abstract
Text Correction [...] Full article
(This article belongs to the Section Remote Sensing in Agriculture and Vegetation)
Show Figures

Figure 6

Figure 6
<p>Scatter plot of the square root transformed fine-fuel biomass from plots and ratio of herbaceous to bare ground cover from plot data with a curvilinear fit (solid line) and 95th percentile confidence intervals (dashed lines).</p>
Full article ">Figure 8
<p>Scatter plots comparing (<b>A</b>) predicted and observed plot biomass (kg/ha) from field plots and WV3 (<span class="html-italic">r</span><sup>2</sup> = 0.70, solid line) and OLI (<span class="html-italic">r</span><sup>2</sup> = 0.74, dashed line) models and (<b>B</b>) predicted biomass (fine-fuel) for plot locations by each sensor type and plot (<span class="html-italic">r</span><sup>2</sup> = 0.73, solid line). The dashed line in (<b>B</b>) is forced though the origin.</p>
Full article ">Figure 15
<p>Modeled (<b>A</b>) land cover types from WV3 imagery and (<b>B</b>) corresponding fine-fuel biomass within the study areas in 2015. Grasslands highly invaded by the non-native grass species <span class="html-italic">E. lehmaninana</span> overlapped with areas of high fine-fuel biomass accumulation. Grasslands with a greater number of native grass species were associated with lower fine-fuel accumulations.</p>
Full article ">Figure 16
<p>Fine-fuels and fuel-type data compared within principal fire management units (<span class="html-italic">n</span> = 59) on the Buenos Aires National Wildlife Refuge (BANWR) showed a strong positive relationship between non-native grass cover and average fine-fuel biomass.</p>
Full article ">
27 pages, 23704 KiB  
Article
A Proposal for Automatic Coastline Extraction from Landsat 8 OLI Images Combining Modified Optimum Index Factor (MOIF) and K-Means
by Francesco Giuseppe Figliomeni, Francesca Guastaferro, Claudio Parente and Andrea Vallario
Remote Sens. 2023, 15(12), 3181; https://doi.org/10.3390/rs15123181 - 19 Jun 2023
Cited by 6 | Viewed by 3383
Abstract
The coastal environment is a natural and economic resource of extraordinary value, but it is constantly modifying and susceptible to climate change, human activities and natural hazards. Remote sensing techniques have proved to be excellent for coastal area monitoring, but the main issue [...] Read more.
The coastal environment is a natural and economic resource of extraordinary value, but it is constantly modifying and susceptible to climate change, human activities and natural hazards. Remote sensing techniques have proved to be excellent for coastal area monitoring, but the main issue is to detect the borderline between water bodies (ocean, sea, lake or river) and land. This research aims to define a rapid and accurate methodological approach, based on the k-means algorithm, to classify the remotely sensed images in an unsupervised way to distinguish water body pixels and detect coastline. Landsat 8 Operational Land Imager (OLI) multispectral satellite images were considered. The proposal requires applying the k-means algorithm only to the most appropriate multispectral bands, rather than using the entire dataset. In fact, by using only suitable bands to detect the differences between water and no-water (vegetation and bare soil), more accurate results were obtained. For this scope, a new index based on the optimum index factor (OIF) was applied to identify the three best-performing bands for the purpose. The direct comparison between the automatically extracted coastline and the manually digitized one was used to evaluate the product accuracy. The results were very satisfactory and the combination involving bands B2 (blue), B5 (near infrared), and B6 (short-wave infrared-1) provided the best performance. Full article
(This article belongs to the Special Issue Mapping and Change Analysis Applications with Remote Sensing and GIS)
Show Figures

Figure 1

Figure 1
<p>Study area: on the left, the location of the study area in the Tyrrhenian Sea in equirectangular projection and WGS 84 geographic coordinates (EPSG:4326); on the right, the visualization in RGB true color composition of Landsat 8 OLI images in UTM/WGS 84 plane coordinates expressed in meters (EPSG: 32632).</p>
Full article ">Figure 2
<p>Workflow of the methodological approach adopted in our study.</p>
Full article ">Figure 3
<p>False color visualization (on the <b>left</b>) and result of KM clustering (on the <b>right</b>) applied to bands 2-5-6.</p>
Full article ">Figure 4
<p>False color visualization (on the <b>left</b>) and result of KM clustering (on the <b>right</b>) applied to bands 1-3-6.</p>
Full article ">Figure 5
<p>False color visualization (on the <b>left</b>) and result of KM clustering (on the <b>right</b>) applied to bands 1-2-9.</p>
Full article ">Figure 6
<p>Geolocation of the three examined frames.</p>
Full article ">Figure 7
<p>Comparison between the reference coastline (in black) and the automatically vectorized coastline (in red) resulting from B2, B5, B6 band composition in frame 1.</p>
Full article ">Figure 8
<p>Comparison between the reference coastline (in black) and the automatically vectorized coastline (in red) resulting from B3, B4, B5 band composition in frame 1.</p>
Full article ">Figure 9
<p>Comparison between the reference coastline (in black) and the automatically vectorized coastline (in red) resulting from B3, B4, B6 band composition in frame 1.</p>
Full article ">Figure 10
<p>Comparison between the reference coastline (in black) and the automatically vectorized coastline (in red) resulting from B2, B3, B9 band composition in frame 1.</p>
Full article ">Figure 11
<p>Comparison between the reference coastline (in black) and the automatically vectorized coastline (in red) resulting from B2, B5, B6 band composition in frame 2.</p>
Full article ">Figure 12
<p>Comparison between the reference coastline (in black) and the automatically vectorized coastline (in red) resulting from B3, B4, B5 band composition in frame 2.</p>
Full article ">Figure 13
<p>Comparison between the reference coastline (in black) and the automatically vectorized coastline (in red) resulting from B3, B4, B6 band composition in frame 2.</p>
Full article ">Figure 14
<p>Comparison between the reference coastline (in black) and the automatically vectorized coastline (in red) resulting from B2, B3, B9 band composition in frame 2.</p>
Full article ">Figure 15
<p>Comparison between the reference coastline (in black) and the automatically vectorized coastline (in red) resulting from B2, B5, B6 band composition in frame 3.</p>
Full article ">Figure 16
<p>Comparison between the reference coastline (in black) and the automatically vectorized coastline (in red) resulting from B3, B4, B5 band composition in frame 3.</p>
Full article ">Figure 17
<p>Comparison between the reference coastline (in black) and the automatically vectorized coastline (in red) resulting from B3, B4, B6 band composition in frame 3.</p>
Full article ">Figure 18
<p>Comparison between the reference coastline (in black) and the automatically vectorized coastline (in red) resulting from B2, B3, B9 band composition in frame 3.</p>
Full article ">
21 pages, 10765 KiB  
Article
Morphology Dynamics of Ice Cover in a River Bend Revealed by the UAV-GPR and Sentinel-2
by Chunjiang Li, Zhijun Li, Wenfeng Huang, Baosen Zhang, Yu Deng and Guoyu Li
Remote Sens. 2023, 15(12), 3180; https://doi.org/10.3390/rs15123180 - 19 Jun 2023
Cited by 1 | Viewed by 1272
Abstract
After the formation of the bend ice cover, the ice thickness of the bend is not uniformly distributed, and an open-water area is usually formed downstream of the bend. The spatial and temporal variation of the ice thickness in seven cross sections was [...] Read more.
After the formation of the bend ice cover, the ice thickness of the bend is not uniformly distributed, and an open-water area is usually formed downstream of the bend. The spatial and temporal variation of the ice thickness in seven cross sections was determined via Unmanned Aerial Vehicle Ground Penetrating Radar (UAV-GPR) technology and traditional borehole measurements. The plane morphology change of the open water was observed by Sentinel-2. The results show that the average dielectric permittivity of GPR was 3.231, 3.249, and 3.317 on three surveys (5 January 2022, 16 February 2022, and 25 February 2022) of the Yellow River ice growing period, respectively. The average ice thickness of the three surveys was 0.402 m, 0.509 m, and 0.633 m, respectively. The ice thickness of the concave bank was larger than that of the convex bank. The plane morphology of the open water first shrinks rapidly longitudinally and then shrinks slowly transversely. The vertical boundary of the open water was composed of two arcs, in which the slope of Arc I (close to the water surface) was steeper than that of Arc II, and the hazardous distance of the open-water boundary was 10.3 m. The increased flow mostly affected the slope change of Arc I. Finally, we discuss the variation of hummocky ice and flat ice in GPR images and the physical factors affecting GPR detection accuracy, as well as the ice-thickness variation of concave and convex banks in relation to channel curvature. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>The location of the Shisifenzi bend on the Yellow River, China. Red box is the area where the UAV-GPR measured ice thickness.</p>
Full article ">Figure 2
<p>(<b>a</b>) Three flight paths of the UAV-GPR on 5 January 2022, 16 February 2022, and 25 February 2022) and seven planned cross sections (from CS0 to CS6) on the Shisifenzi Bend (red box in <a href="#remotesensing-15-03180-f002" class="html-fig">Figure 2</a>) of the Yellow River, (<b>b</b>) the UAV-GPR system, (<b>c</b>) the schematic diagram of ice-thickness detection by the UAV-GPR. In (<b>a</b>), points A, B, C, D, E, and F represent the intersection of the flight path with the boundary of open water on 16 February 2022; points G, H, J and K represent the intersection of the flight path with the boundary of open water on 25 February 2022. The white dots in (<b>a</b>) represent the locations of the drilled ice holes. The background in (<b>a</b>) is an orthophoto image taken by DJI Phantom 4 on 27 February 2022.</p>
Full article ">Figure 3
<p>The typical GPR profile on 25 February 2022.</p>
Full article ">Figure 4
<p>Ice thickness measured by ice drill vs that measured by GPR.</p>
Full article ">Figure 5
<p>The spatial distribution of the ice thickness on the (<b>a</b>) 5 January 2022, (<b>b</b>) 16 February 2022, (<b>c</b>) 25 February 2022, and (<b>d</b>) changes in the ice thickness between 16 February 2022 and 25 February 2022. The CS0–CS6 are the observation cross sections laid out. These contour maps were interpolated onto the 10 m-by-10 m grid.</p>
Full article ">Figure 6
<p>Plane morphology changes of the open water of the Shisifenzi bend in the winter of 2021–2022. The red solid boxes in (<b>a</b>,<b>b</b>,<b>h</b>) is the same as the range from (<b>c</b>–<b>g</b>). Both CS7 and CS8 are the cross sections for observing the transverse variation of the open water.</p>
Full article ">Figure 7
<p>(<b>a</b>) The change of relative length, width, and area of the open water over time; (<b>b</b>) air temperature in the winter of 2021–2022.</p>
Full article ">Figure 8
<p>Ice thickness at the boundary of the open water on 16 February 2022. The distance starting point of each image is the boundary of the open water. The (<b>a</b>) A, (<b>b</b>) B, (<b>c</b>) C, (<b>d</b>) D, (<b>e</b>) E, (<b>f</b>) F represent the locations in the <a href="#remotesensing-15-03180-f002" class="html-fig">Figure 2</a>a. The blue dotted line is the demarcation between I and II. The blue solid line is the demarcation between water and revetment.</p>
Full article ">Figure 9
<p>Ice thickness at the boundary of the open water on 25 February 2022. The distance starting point of each image is the boundary of the open water. The (<b>a</b>) G, (<b>b</b>) H, (<b>c</b>) J, (<b>d</b>) K represent the locations in the <a href="#remotesensing-15-03180-f002" class="html-fig">Figure 2</a>a. I and II represent two arcs at the boundary of the ice cover. The blue dotted line is the demarcation between I and II. The blue solid line is the demarcation between water and revetment.</p>
Full article ">Figure 10
<p>The slope of (<b>a</b>) Arc I and (<b>b</b>) Arc II at different locations on 16 February 2022 and on 25 February 2022. (<b>c</b>) the flow changes of Toudaoguai hydrological station.</p>
Full article ">Figure 11
<p>Field and schematic diagram of (<b>a</b>) Hummocky Ice and (<b>b</b>) Flat Ice.</p>
Full article ">Figure 12
<p>The part radar image on the CS3 at (<b>a</b>) 5 January 2022, (<b>b</b>) 16 February 2022 and (<b>c</b>) 25 February 2022. In (<b>a</b>), (<b>b</b>) and (<b>c</b>), the 0 ns in the radar image is based on the UAV flight altitude.</p>
Full article ">Figure 13
<p>The natural ice, crystal type, and air bubble of (<b>a</b>) flat ice (FI) and (<b>b</b>) hummocky ice (HI). The ice samples were cut on 14 February 2022 in the Shisifenzi bend of the Yellow River.</p>
Full article ">Figure 14
<p>(<b>a</b>) the average ice thickness and standard deviation of the concave and convex banks at different surveys and locations; (<b>b</b>) the percentage (<math display="inline"><semantics> <mrow> <mi>M</mi> </mrow> </semantics></math>) of the average ice thickness of the concave bank was greater than that of the convex bank. The negative <math display="inline"><semantics> <mrow> <mi>M</mi> </mrow> </semantics></math> represents that the average ice thickness of the convex bank was larger than that of the concave bank.</p>
Full article ">Figure 15
<p>The Sentinel images in early freeze-up of bend on (<b>a</b>) 24 November 2021, (<b>b</b>) 4 December 2021, (<b>c</b>) 9 December 2021, and (<b>d</b>) 19 December 2021. The CS0–CS6 were the planned observation cross sections. The area enclosed by red dotted line represents the calm flow area in the concave bank. The area enclosed by blue dotted line represents the big accumulation area at the early freeze-up.</p>
Full article ">Figure 16
<p>The discharge of river and water level in Toudaoguai station.</p>
Full article ">
17 pages, 6156 KiB  
Article
First Experience with Zhuhai-1 Hyperspectral Data for Urban Dominant Tree Species Classification in Shenzhen, China
by Haiming Qin, Weimin Wang, Yang Yao, Yuguo Qian, Xiangyun Xiong and Weiqi Zhou
Remote Sens. 2023, 15(12), 3179; https://doi.org/10.3390/rs15123179 - 19 Jun 2023
Cited by 4 | Viewed by 2393
Abstract
An accurate spatial distribution map of the urban dominant tree species is crucial for evaluating the ecosystem service value of urban forests and formulating urban sustainable development strategies. Spaceborne hyperspectral remote sensing has been utilized to distinguish tree species, but these hyperspectral data [...] Read more.
An accurate spatial distribution map of the urban dominant tree species is crucial for evaluating the ecosystem service value of urban forests and formulating urban sustainable development strategies. Spaceborne hyperspectral remote sensing has been utilized to distinguish tree species, but these hyperspectral data have a low spatial resolution (pixel size ≥ 30 m), which limits their ability to differentiate tree species in urban areas characterized by fragmented patches and robust spatial heterogeneity. Zhuhai-1 is a new hyperspectral satellite sensor with a higher spatial resolution of 10 m. This study aimed to evaluate the potential of Zhuhai-1 hyperspectral imagery for classifying the urban dominant tree species. We first extracted 32 reflectance bands and 18 vegetation indices from Zhuhai-1 hyperspectral data. We then used the random forest classifier to differentiate 28 dominant tree species in Shenzhen based on these hyperspectral features. Finally, we analyzed the effects of the classification paradigm, classifier, and species number on the classification accuracy. We found that combining the hyperspectral reflectance bands and vegetation indices could effectively distinguish the 28 dominant tree species in Shenzhen, obtaining an overall accuracy of 76.8%. Sensitivity analysis results indicated that the pixel-based classification paradigm was slightly superior to the object-based paradigm. The random forest classifier proved to be the optimal classifier for distinguishing tree species using Zhuhai-1 hyperspectral imagery. Moreover, reducing the species number could slowly improve the classification accuracy. These findings suggest that Zhuhai-1 hyperspectral data can identify the urban dominant tree species with accuracy and holds potential for application in other cities. Full article
(This article belongs to the Section Urban Remote Sensing)
Show Figures

Figure 1

Figure 1
<p>The study area located in Shenzhen, southern China.</p>
Full article ">Figure 2
<p>Field samples used in this study.</p>
Full article ">Figure 3
<p>Average reflectance according to spectral bands for 28 main tree species.</p>
Full article ">Figure 4
<p>The flowchart for urban dominant tree species classification using Zhuhai-1 hyperspectral data.</p>
Full article ">Figure 5
<p>The variable importance of all hyperspectral features for classification of tree species using the random forest classifier.</p>
Full article ">Figure 6
<p>(<b>a</b>) Pixel-based tree species classification results for the entire study area derived from all hyperspectral features using the RF classifier; (<b>b</b>) detailed pixel-based tree species classification results for the area in the black box in (<b>a</b>); (<b>c</b>) detailed object-based tree species classification results for the area in the black box in (<b>a</b>).</p>
Full article ">Figure 7
<p>The relationship between the species number and classification accuracy.</p>
Full article ">
Previous Issue
Next Issue
Back to TopTop