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

Years

Between: -

Subjects

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Journals

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Article Types

Countries / Regions

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Search Results (2,561)

Search Parameters:
Keywords = remote-sensed reflectance

Order results
Result details
Results per page
Select all
Export citation of selected articles as:
18 pages, 4628 KiB  
Technical Note
Characterizing Changes in Paddy Rice Flooding Time over the Sanjiang Plain Using Moderate-Resolution Imaging Spectroradiometer Time Series
by Xiangyu Ning, Huapeng Li and Ruoqi Liu
Remote Sens. 2024, 16(24), 4683; https://doi.org/10.3390/rs16244683 - 15 Dec 2024
Viewed by 341
Abstract
Rice is a primary food crop, and rice production ensures food security and maintains social stability with great significance. Flooding paddy rice fields as an important step in rice production affects the entire growth process of rice. The selection of flooding time is [...] Read more.
Rice is a primary food crop, and rice production ensures food security and maintains social stability with great significance. Flooding paddy rice fields as an important step in rice production affects the entire growth process of rice. The selection of flooding time is highly correlated with paddy rice yield and water resource utilization. In the background of global warming, early flooding in high-latitude paddy rice planting areas can ensure that rice has sufficient growing time to increase yield. However, overly early flooding may cause waste of water resources due to insufficient heat. Currently, research on flooding timing is relatively lacking, and monitoring of temperature during flooding is particularly deficient. To respond to climate change, it is necessary to explore whether the current flooding schedule meets the actual needs. Based on MODIS surface reflectivity data, we identified the First Flooding Day (FFD) and Peak Flooding Day (PFD) in the Sanjiang Plain. Using MODIS Land Surface Temperature (LST) data and meteorological station-provided air temperature data, we analyzed the corresponding LST and air temperature for PFD from 2008 to 2024. The main conclusions are as follows: (1) both FFD and PFD in the Sanjiang Plain have a trend of advancing year by year, with PFD showing stronger advancement than FFD; (2) the LST and air temperature during flooding in the Sanjiang Plain show a downward trend year by year; and (3) by 2024, the flooding temperature of paddy rice fields in the Sanjiang Plain has generally met the needs for the next step of production. This study first attempts to use high-temporal-resolution remote sensing images to identify the flooding time of paddy fields and achieve timely monitoring of flooding and changes in flooding temperature. Full article
Show Figures

Figure 1

Figure 1
<p>Location of the Sanjiang Plain. The rice data were from You et al. (2021) [<a href="#B33-remotesensing-16-04683" class="html-bibr">33</a>]. DEM data were from SRTM. Water data were from Pekel et al. (2016) [<a href="#B34-remotesensing-16-04683" class="html-bibr">34</a>].</p>
Full article ">Figure 2
<p>Method for composite daily flooding images in potential flooding periods, and the definitions of FFD and PFD. Different pixel colors correspond to different DoY.</p>
Full article ">Figure 3
<p>The spatio-temporal distribution of paddy rice flooding time (FFD and PFD) in Sanjiang from 2008 to 2024. The insert numbers of each subfigure mean the different years from 2008 to 2024.</p>
Full article ">Figure 4
<p>Changes in paddy fields’ FFD and PFD year by year. (<b>a</b>) Broken line of FFD; (<b>b</b>) linear fit of FFD; (<b>c</b>) broken line of PFD; and (<b>d</b>) linear fit of PFD.</p>
Full article ">Figure 5
<p>PFD’s LST changes year by year. (<b>a</b>) PFD’s LST at 13:30. (<b>b</b>) PFD’s LST at 01:30 the next day.</p>
Full article ">Figure 6
<p>PFD’s T<sub>air</sub> changes year by year. (<b>a</b>) PFD’s T<sub>air</sub> at 14:00. (<b>b</b>) PFD’s T<sub>air</sub> at 02:00 the next day.</p>
Full article ">Figure 7
<p>The differences between FFD and PFD year by year. (<b>a</b>) Broken line of the difference. (<b>b</b>) Linear fit of the difference.</p>
Full article ">Figure 8
<p>Average temperature from 1 April to 31 May of the Sanjiang Plain in 2008–2024.</p>
Full article ">
21 pages, 13076 KiB  
Article
A Framework for High-Spatiotemporal-Resolution Soil Moisture Retrieval in China Using Multi-Source Remote Sensing Data
by Zhuangzhuang Feng, Xingming Zheng, Xiaofeng Li, Chunmei Wang, Jinfeng Song, Lei Li, Tianhao Guo and Jia Zheng
Land 2024, 13(12), 2189; https://doi.org/10.3390/land13122189 - 15 Dec 2024
Viewed by 585
Abstract
High-spatiotemporal-resolution and accurate soil moisture (SM) data are crucial for investigating climate, hydrology, and agriculture. Existing SM products do not yet meet the demands for high spatiotemporal resolution. The objective is to develop and evaluate a retrieval framework to derive SM estimates with [...] Read more.
High-spatiotemporal-resolution and accurate soil moisture (SM) data are crucial for investigating climate, hydrology, and agriculture. Existing SM products do not yet meet the demands for high spatiotemporal resolution. The objective is to develop and evaluate a retrieval framework to derive SM estimates with high spatial (100 m) and temporal (<3 days) resolution that can be used on a national scale in China. Therefore, this study integrates multi-source data, including optical remote sensing (RS) data from Sentinel-2 and Landsat-7/8/9, synthetic aperture radar (SAR) data from Sentinel-1, and auxiliary data. Four machine learning and deep learning algorithms are applied, including Random Forest Regression (RFR), Extreme Gradient Boosting (XGBoost), Long Short-Term Memory (LSTM) networks, and Ensemble Learning (EL). The integrated framework (IF) considers three feature scenarios (SC1: optical RS + auxiliary data, SC2: SAR + auxiliary data, SC3: optical RS + SAR + auxiliary data), encompassing a total of 33 features. The results are as follows: (1) The correlation coefficients (r) between auxiliary data (such as sand fraction, r = −0.48; silt fraction, r = 0.47; and evapotranspiration, r = −0.42), SAR features (such as the backscatter coefficients for VV-pol (σvv0), r = 0.47), and optical RS features (such as Shortwave Infrared Band 2 (SWIR2) reflectance data from Sentinel-2 and Landsat-7/8/9, r = −0.39) with observed SM are significant. This indicates that multi-source data can provide complementary information for SM monitoring. (2) Compared to XGBoost and LSTM, RFR and EL demonstrate superior overall performance and are the preferred models for SM prediction. Their R2 for the training and test sets exceed 0.969 and 0.743, respectively, and their ubRMSE are below 0.022 and 0.063 m3/m3, respectively. (3) The SM prediction accuracy is highest for the scenario of optical + SAR + auxiliary data, followed by SAR + auxiliary data, and finally optical + auxiliary data. (4) With an increasing Normalized Difference Vegetation Index (NDVI) and SM values, the trained models exhibit a general decrease in prediction performance and accuracy. (5) In 2021 and 2022, without considering cloud cover, the IF theoretically achieved an SM revisit time of 1–3 days across 95.01% and 96.53% of China’s area, respectively. However, SC1 was able to achieve a revisit time of 1–3 days over 60.73% of China’s area in 2021 and 69.36% in 2022, while the area covered by SC2 and SC3 at this revisit time accounted for less than 1% of China’s total area. This study validates the effectiveness of combining multi-source RS data with auxiliary data in large-scale SM monitoring and provides new methods for improving SM retrieval accuracy and spatiotemporal coverage. Full article
(This article belongs to the Section Land – Observation and Monitoring)
Show Figures

Figure 1

Figure 1
<p>The spatial distribution of the SONTE-China 17 sites within the study area.</p>
Full article ">Figure 2
<p>A framework for estimating SM based on multi-source RS data. *** represents the first priority, ** represents the second priority, and * represents the third priority.</p>
Full article ">Figure 3
<p>The training (<b>top</b>) and test (<b>bottom</b>) results of four models from IF at SONTE-China (17 sites). The red dotted line is the trend line. The gray dotted line represents the error line at 0.06 m<sup>3</sup>/m<sup>3</sup>.</p>
Full article ">Figure 4
<p>The training results of four models at SONTE-China (17 sites). The red dotted line is the trend line. The gray dotted line represents the error line at 0.06 m<sup>3</sup>/m<sup>3</sup>.</p>
Full article ">Figure 5
<p>The test results of four models at SONTE-China (17 sites). The red dotted line is the trend line. The gray dotted line represents the error line at 0.06 m<sup>3</sup>/m<sup>3</sup>.</p>
Full article ">Figure 6
<p>The time series of estimated and observed SM from three scenarios at NQ, JYT, and MQ sites. The blue solid line represents the observed SM at 0–5 cm. The green solid line represents the daily NDVI. The red, green, and purple squares represent the estimated SM for SC1, SC2, and SC3, respectively. The blue bars indicate daily precipitation. The red dashed vertical lines distinguish between the training and test sets.</p>
Full article ">Figure 7
<p>Revisit time between SC1, SC2, SC3, and IF for monitoring SM in China (2021). (<b>a</b>) SC1: Optical RS + auxiliary data only; (<b>b</b>) SC2: SAR + auxiliary data only; (<b>c</b>) SC3: optical RS + SAR + auxiliary data; (<b>d</b>) IF: combined SC3, SC2, and SC1 scenarios.</p>
Full article ">Figure 8
<p>Revisit time between SC1, SC2, SC3, and IF for monitoring SM in China (2022). (<b>a</b>) SC1: Optical RS + auxiliary data only; (<b>b</b>) SC2: SAR + auxiliary data only; (<b>c</b>) SC3: optical RS + SAR + auxiliary data; (<b>d</b>) IF: combined SC3, SC2, and SC1 scenarios.</p>
Full article ">Figure 9
<p>Training (<b>top</b>) and test (<b>bottom</b>) results of three categories using the RFR based on the SC3 dataset at SONTE-China (17 sites). The red dotted line is the trend line. The gray dotted line represents the error line at 0.06 m<sup>3</sup>/m<sup>3</sup>.</p>
Full article ">Figure 10
<p>Performance of different models under various NDVI categories in the training set (<b>left</b>) and test set (<b>right</b>). The colored dot lines represent R<sup>2</sup>, and the bar charts represent ubRMSE.</p>
Full article ">Figure 11
<p>Performance of different models under various SM categories in the training set (<b>left</b>) and test set (<b>right</b>). The bar charts represent ubRMSE, and the red dot line represents the average ubRMSE.</p>
Full article ">Figure 12
<p>Revisit time distribution for multi-source RS monitoring of SM under different scenarios (2021–2022).</p>
Full article ">
19 pages, 10993 KiB  
Article
Observation Angle Effect of Near-Ground Thermal Infrared Remote Sensing on the Temperature Results of Urban Land Surface
by Xu Yuan, Zhi Lv, Kati Laakso, Jialiang Han, Xiao Liu, Qinglin Meng and Sihan Xue
Land 2024, 13(12), 2170; https://doi.org/10.3390/land13122170 - 13 Dec 2024
Viewed by 315
Abstract
During the process of urbanization, a large number of impervious land surfaces are replacing the biologically active surface. Land surface temperature is a key factor reflecting the urban thermal environment and a crucial factor affecting city livability and resident comfort. Therefore, the accurate [...] Read more.
During the process of urbanization, a large number of impervious land surfaces are replacing the biologically active surface. Land surface temperature is a key factor reflecting the urban thermal environment and a crucial factor affecting city livability and resident comfort. Therefore, the accurate measurement of land surface temperature is of great significance. Thermal infrared remote sensing is widely applied to study the urban thermal environment due to its distinctive advantages of high sensitivity, wide coverage, high resolution, and continuous measurement. Low-altitude remote sensing, performed using thermal infrared sensors carried by unmanned aerial vehicles (UAVs), is a common method of land surface observation. However, thermal infrared sensors may experience varying degrees of sway due to wind, affecting the quality of the data. It is still uncertain as to what degree angle changes affect thermal infrared data in urban environments. To investigate this effect, a near-ground remote sensing experiment was conducted to observe three common urban land surfaces, namely, marble tiles, cement tiles and grasses, at observation angles of 15°, 30°, 45°, and 60° using a thermal infrared imager. This is accompanied by synchronous ground temperature measurements conducted by iButton digital thermometers. Our results suggest that the temperature differences between the remote sensing data of the land surface and the corresponding ground truth data increase as a function of the increasing observation angle of the three land surfaces. Furthermore, the differences are minor when the observation angle changes are not more than 15° and the changes are not the same for different land surfaces. Our findings increase the current understanding of the effects of different angles on thermal infrared remote sensing in urban land surface temperature monitoring. Full article
(This article belongs to the Section Land – Observation and Monitoring)
Show Figures

Figure 1

Figure 1
<p>Map of study area locations (Reference system: UTM-WGS 1984, EPSG: 4326): (<b>a</b>) Guangzhou City within Guangdong Province, China; (<b>b</b>) study area in Tianhe district, Guangzhou City; (<b>c</b>) Wushan Campus of South China University of Technology. Area A: plaza in front of the Liwu Science and Technology Building; area B: courtyard of the No. 2 Graduate Student Dormitory Building.</p>
Full article ">Figure 2
<p>Measuring points for ground temperature measurement experiments: (<b>a</b>) cement tiles; (<b>b</b>) marble tiles; (<b>c</b>) grasses. Specific measurement points are marked using red dots.</p>
Full article ">Figure 3
<p>Thermal infrared remote sensing data collection. The observation angle of the infrared thermal imager was adjusted according to the reference lines on the white paper.</p>
Full article ">Figure 4
<p>Visible images of the land surface observation area and the ranges of different thermal infrared observation angles.</p>
Full article ">Figure 5
<p>Thermal infrared images of different urban land surfaces at different observation angles. Area A: (<b>a</b>) Cement tiles at an observation angle of 30°. (<b>b</b>) Cement and marble tiles at an observation angle of 45°. (<b>c</b>) Marble tiles at an observation angle of 60°. Area B: (<b>d</b>) Grasses at an observation angle of 15°. (<b>e</b>) Grasses at an observation angle of 30°.</p>
Full article ">Figure 6
<p>Thermal infrared remote sensing data (<math display="inline"><semantics> <mrow> <msub> <mrow> <mi>T</mi> </mrow> <mrow> <mi>i</mi> <mi>n</mi> <mi>f</mi> </mrow> </msub> </mrow> </semantics></math>) and ground data (<math display="inline"><semantics> <mrow> <msub> <mrow> <mi>T</mi> </mrow> <mrow> <mi>i</mi> <mi>B</mi> <mi>u</mi> </mrow> </msub> <mo>)</mo> </mrow> </semantics></math>, and temperature difference between them. (<b>a</b>) Cement tiles: results at observation angles of 30° and 45°. (<b>b</b>) Marble tiles: results at observation angles of 45° and 60°. (<b>c</b>) Grasses: results at observation angles of 15° and 30°.</p>
Full article ">Figure 6 Cont.
<p>Thermal infrared remote sensing data (<math display="inline"><semantics> <mrow> <msub> <mrow> <mi>T</mi> </mrow> <mrow> <mi>i</mi> <mi>n</mi> <mi>f</mi> </mrow> </msub> </mrow> </semantics></math>) and ground data (<math display="inline"><semantics> <mrow> <msub> <mrow> <mi>T</mi> </mrow> <mrow> <mi>i</mi> <mi>B</mi> <mi>u</mi> </mrow> </msub> <mo>)</mo> </mrow> </semantics></math>, and temperature difference between them. (<b>a</b>) Cement tiles: results at observation angles of 30° and 45°. (<b>b</b>) Marble tiles: results at observation angles of 45° and 60°. (<b>c</b>) Grasses: results at observation angles of 15° and 30°.</p>
Full article ">Figure 7
<p>Box plot of temperature differences between thermal infrared remote sensing data (<math display="inline"><semantics> <mrow> <msub> <mrow> <mi>T</mi> </mrow> <mrow> <mi>i</mi> <mi>n</mi> <mi>f</mi> </mrow> </msub> </mrow> </semantics></math>) and ground data (<math display="inline"><semantics> <mrow> <msub> <mrow> <mi>T</mi> </mrow> <mrow> <mi>i</mi> <mi>B</mi> <mi>u</mi> </mrow> </msub> <mo>)</mo> </mrow> </semantics></math>. (<b>a</b>) Cement tiles: results at observation angles of 30° and 45°. (<b>b</b>) Marble tiles: results at observation angles of 45° and 60°. (<b>c</b>) Grasses: results at observation angles of 15° and 30°.</p>
Full article ">Figure 8
<p>The average temperature value of the thermal infrared remote sensing data (<math display="inline"><semantics> <mrow> <msub> <mrow> <mi>T</mi> </mrow> <mrow> <mi>i</mi> <mi>n</mi> <mi>f</mi> </mrow> </msub> </mrow> </semantics></math>) and the ground data (<math display="inline"><semantics> <mrow> <msub> <mrow> <mi>T</mi> </mrow> <mrow> <mi>i</mi> <mi>B</mi> <mi>u</mi> </mrow> </msub> </mrow> </semantics></math>), shown as a function of the observation angle and material (cement tiles, marble tiles and grasses).</p>
Full article ">Figure 9
<p>Normalized results of the thermal infrared remote sensing data (<math display="inline"><semantics> <mrow> <msub> <mrow> <mi>T</mi> </mrow> <mrow> <mi>i</mi> <mi>n</mi> <mi>f</mi> </mrow> </msub> </mrow> </semantics></math>) and the ground data (<math display="inline"><semantics> <mrow> <msub> <mrow> <mi>T</mi> </mrow> <mrow> <mi>i</mi> <mi>B</mi> <mi>u</mi> </mrow> </msub> </mrow> </semantics></math>), shown as a function of the observation angle and land surface (cement tiles, marble tiles and grasses).</p>
Full article ">Figure 10
<p>The temperature differences between the thermal infrared remote sensing data (<math display="inline"><semantics> <mrow> <msub> <mrow> <mi>T</mi> </mrow> <mrow> <mi>i</mi> <mi>n</mi> <mi>f</mi> </mrow> </msub> </mrow> </semantics></math>) and the ground data (<math display="inline"><semantics> <mrow> <msub> <mrow> <mi>T</mi> </mrow> <mrow> <mi>i</mi> <mi>B</mi> <mi>u</mi> </mrow> </msub> </mrow> </semantics></math>), shown as percentages and as a function of the observation angle and land surface (cement tiles, marble tiles and grasses), and the percentage difference between the results from different angles.</p>
Full article ">
16 pages, 4700 KiB  
Technical Note
Precision and Characteristics of Satellite Spatial Quality Estimators’ Measurement Using an Edge Target Imaged with KOMPSAT-3A
by Donghan Lee, Daesoon Park and Daehoon Yoo
Remote Sens. 2024, 16(24), 4660; https://doi.org/10.3390/rs16244660 - 12 Dec 2024
Viewed by 429
Abstract
After the launch of a high-resolution remote sensing satellite, representative spatial quality estimators (RER, FWHM, MTF50, MTFA) are measured from images taken of ground Edge targets. In this work, the best spatial quality estimator is proposed by quantitatively comparing and analyzing the precision [...] Read more.
After the launch of a high-resolution remote sensing satellite, representative spatial quality estimators (RER, FWHM, MTF50, MTFA) are measured from images taken of ground Edge targets. In this work, the best spatial quality estimator is proposed by quantitatively comparing and analyzing the precision between the Relative Edge Response (RER), the Full Width at Half Maximum (FWHM), the MTF value at the Nyquist frequency (MTF50), and the MTF Area between 0 and the Nyquist frequency (MTFA). While the basic method for the measurement of spatial quality estimators on Edge targets is already well established, this work summarizes and explains the uncertain factors and problems in the measurement procedure that affect the accuracy and precision of spatial quality estimators. It also considers how to improve the precision of spatial quality estimators during the measurement procedure. The contents and results of this work were discussed by various satellite development organizations in the Geo-Spatial Working Group within CEOS WGCV IVOS from 2012 to 2019, and the Edge target Spatial quality Measurement Python code (ESMP) was developed in 2019 to reflect the findings of this workshop. Using 483 Edge targets from worldwide images taken by KOMPSAT-3A, which has been in operation since 2017, the results obtained via ESMP show that the precision levels of RER, FWHM, and MTFA are approximately three to four times higher than that of MTF50 when comparing the Coefficient of Variance (CV) statistics. This is the first statistical comparison of spatial quality estimators using 7 years of ground Edge target imagery from a single satellite of KOMPSAT-3A. Full article
(This article belongs to the Special Issue Remote Sensing Satellites Calibration and Validation)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Plot showing the procedure of using the Edge target to measure the spatial quality (RER, FWHM, MTF50, MTFA) using KOMPSAT-3A Level 0 (raw) across Edge image data imaged at Zuunmod with a roll tilt of 29.31 deg on 21 February 2019.</p>
Full article ">Figure 2
<p>Edge targets in catalog of USGS and CEOS Cal/Val portal site [<a href="#B30-remotesensing-16-04660" class="html-bibr">30</a>,<a href="#B31-remotesensing-16-04660" class="html-bibr">31</a>], imaged by KOMPSAT-3 (GSD, 0.7 m).</p>
Full article ">Figure 3
<p>(<b>a</b>) RER is the slope of Edge at −0.5 to 0.5 pixels on ESF (green line); FWHM is the width (pixels) of LSF at 0.5 on the <span class="html-italic">Y</span>-axis over LSF (blue line); RER and FWHM have different measurement locations and units. (<b>b</b>) Procedure for detection of Edge and calculation of the Edge angle and FitErr: obtain line of Edge (green line) and calculate the Edge angle, with FitErr as the standard deviation of Edge points (red dot) over the line of Edge.</p>
Full article ">Figure 4
<p>First CSAPS fitting (<b>left</b>) and second CSAPS fitting of ESF after removing outlier noise (<b>right</b>). The measurement result corrupted by noise is restored. Baotou Edge target imaged by KOMPSAT-3A (13 January 2019).</p>
Full article ">Figure 5
<p>(<b>left</b>) Julian date vs. RER, FWHM, MTF50, and MTFA. (<b>right</b>) Roll tilt angle vs. RER, FWHM, MTF50, and MTFA (upper right in red: IQR CV value). <span class="html-italic">Y</span>-axis scale was changed to CV value.</p>
Full article ">Figure 6
<p>(<b>left</b>) ΔDN vs. RER, FWHM, MTF50, and MTFA. (<b>right</b>) SNR vs. RER, FWHM, MTF50, and MTFA (upper right in red: IQR CV value). <span class="html-italic">Y</span>-axis scale was changed to CV value.</p>
Full article ">Figure 7
<p>(<b>left</b>) Edge angle vs. RER, FWHM, MTF50, and MTFA, separated by the Edge target. (<b>right</b>) <span class="html-italic">X</span>-axis for Edge targets in the order of Baotou (blue), India (red), Salon (green), and Zuunmod (purple), separated by imaging order vs. RER, FWHM, MTF50, and MTFA (bottom left: IQR average and CV). <span class="html-italic">Y</span>-axis scale was changed to CV value.</p>
Full article ">Figure 8
<p>Histograms used to check z-distributions of RER, FWHM, MTF50, and MTFA (upper right: skew value).</p>
Full article ">Figure 9
<p>Scatter plot of RER vs. FWHM vs. MTF50 vs. MTFA (upper right: (red) Pearson’s values of constraint and IQR, and (blue) CV values of spatial quality estimator). Plot scale was changed to z-distribution.</p>
Full article ">
20 pages, 13662 KiB  
Article
Unmanned Aerial Vehicle (UAV) Hyperspectral Imagery Mining to Identify New Spectral Indices for Predicting the Field-Scale Yield of Spring Maize
by Yue Zhang, Yansong Wang, Hang Hao, Ziqi Li, Yumei Long, Xingyu Zhang and Chenzhen Xia
Sustainability 2024, 16(24), 10916; https://doi.org/10.3390/su162410916 - 12 Dec 2024
Viewed by 593
Abstract
A nondestructive approach for accurate crop yield prediction at the field scale is vital for precision agriculture. Considerable progress has been made in the use of the spectral index (SI) derived from unmanned aerial vehicle (UAV) hyperspectral images to predict crop yields before [...] Read more.
A nondestructive approach for accurate crop yield prediction at the field scale is vital for precision agriculture. Considerable progress has been made in the use of the spectral index (SI) derived from unmanned aerial vehicle (UAV) hyperspectral images to predict crop yields before harvest. However, few studies have explored the most sensitive wavelengths and SIs for crop yield prediction, especially for different nitrogen fertilization levels and soil types. This study aimed to investigate the appropriate wavelengths and their combinations to explore the ability of new SIs derived from UAV hyperspectral images to predict yields during the growing season of spring maize. In this study, the hyperspectral canopy reflectance measurement method, a field-based high-throughput method, was evaluated in three field experiments (Wang-Jia-Qiao (WJQ), San-Ke-Shu (SKS), and Fu-Jia-Jie (FJJ)) since 2009 with different soil types (alluvial soil, black soil, and aeolian sandy soil) and various nitrogen (N) fertilization levels (0, 168, 240, 270, and 312 kg/ha) in Lishu County, Northeast China. The measurements of canopy spectral reflectance and maize yield were conducted at critical growth stages of spring maize, including the jointing, silking, and maturity stages, in 2019 and 2020. The best wavelengths and new SIs, including the difference spectral index, ratio spectral index, and normalized difference spectral index forms, were obtained from the contour maps constructed by the coefficient of determination (R2) from the linear regression models between the yield and all possible SIs screened from the 450 to 950 nm wavelengths. The new SIs and eight selected published SIs were subsequently used to predict maize yield via linear regression models. The results showed that (1) the most sensitive wavelengths were 640–714 nm at WJQ, 450–650 nm and 750–950 nm at SKS, and 450–700 nm and 750–950 nm at FJJ; (2) the new SIs established here were different across the three experimental fields, and their performance in maize yield prediction was generally better than that of the published SIs; and (3) the new SIs presented different responses to various N fertilization levels. This study demonstrates the potential of exploring new spectral characteristics from remote sensing technology for predicting the field-scale crop yield in spring maize cropping systems before harvest. Full article
Show Figures

Figure 1

Figure 1
<p>Location of the study area (<b>a</b>), UAV hyperspectral images (<b>b</b>–<b>d</b>) and the nitrogen application rates (<b>e</b>) of three experimental fields (WJQ, SKS, and FJJ).</p>
Full article ">Figure 2
<p>Mean canopy reflectance spectra curves of spring maize under different N treatments across three growth stages in the three experimental fields. (<b>a</b>): WJQ, (<b>b</b>): SKS, (<b>c</b>): FJJ.</p>
Full article ">Figure 3
<p>Contour maps for the linear model between the difference spectral index (DSI), ratio spectral index (RSI), normalized difference spectral index (NDSI), and maize yield for the WJQ experimental field. (<b>a</b>–<b>c</b>): DSI, RSI, and NDSI forms at the jointing stage; (<b>d</b>–<b>f</b>): DSI, RSI, and NDSI forms at the silking stage; (<b>g</b>–<b>i</b>): DSI, RSI, and NDSI forms at the maturity stage.</p>
Full article ">Figure 4
<p>Contour maps for the linear model between the difference spectral index (DSI), ratio spectral index (RSI), normalized difference spectral index (NDSI), and maize yield for the SKS experimental field. (<b>a</b>–<b>c</b>): DSI, RSI, and NDSI forms at the jointing stage; (<b>d</b>–<b>f</b>): DSI, RSI, and NDSI forms at the silking stage; (<b>g</b>–<b>i</b>): DSI, RSI, and NDSI forms at the maturity stage.</p>
Full article ">Figure 5
<p>Contour maps for the linear model between the difference spectral index (DSI), ratio spectral index (RSI), normalized difference spectral index (NDSI), and maize yield for the FJJ experimental field. (<b>a</b>–<b>c</b>): DSI, RSI, and NDSI forms at the jointing stage; (<b>d</b>–<b>f</b>): DSI, RSI, and NDSI forms at the silking stage; (<b>g</b>–<b>i</b>): DSI, RSI, and NDSI forms at the maturity stage.</p>
Full article ">Figure 6
<p>Scatter plots of the measured yield (kg/ha) versus the yield (kg/ha) predicted by the new SIs: (<b>a</b>) NDSI (690, 710) at WJQ, (<b>b</b>) RSI (906, 546) at SKS, and (<b>c</b>) DSI (698, 922) at FJJ.</p>
Full article ">Figure 7
<p>The response of the maize yield to different N application rates on the three experimental fields.</p>
Full article ">Figure 8
<p>The response of the new SIs to different N treatments on the three experimental fields. (<b>a</b>−<b>c</b>): DSI, RSI, and NDSI forms for WJQ; (<b>d</b>−<b>f</b>): DSI, RSI, and NDSI forms for SKS; and (<b>g</b>−<b>i</b>): DSI, RSI, and NDSI forms for FJJ, respectively.</p>
Full article ">Figure 8 Cont.
<p>The response of the new SIs to different N treatments on the three experimental fields. (<b>a</b>−<b>c</b>): DSI, RSI, and NDSI forms for WJQ; (<b>d</b>−<b>f</b>): DSI, RSI, and NDSI forms for SKS; and (<b>g</b>−<b>i</b>): DSI, RSI, and NDSI forms for FJJ, respectively.</p>
Full article ">Figure 8 Cont.
<p>The response of the new SIs to different N treatments on the three experimental fields. (<b>a</b>−<b>c</b>): DSI, RSI, and NDSI forms for WJQ; (<b>d</b>−<b>f</b>): DSI, RSI, and NDSI forms for SKS; and (<b>g</b>−<b>i</b>): DSI, RSI, and NDSI forms for FJJ, respectively.</p>
Full article ">
23 pages, 71464 KiB  
Article
The Decisive Influence of the Improved Remote Sensing Ecological Index on the Terrestrial Ecosystem in Typical Arid Areas of China
by Long Guo, Chao Xu, Hongqi Wu, Mingjie Shi and Yanmin Fan
Land 2024, 13(12), 2162; https://doi.org/10.3390/land13122162 - 12 Dec 2024
Viewed by 332
Abstract
This study aims to assess the spatiotemporal changes in ecological environment quality (EEQ) in arid regions, using Xinjiang as a case study, from 2000 to 2023, with an improved remote sensing ecological index (IRSEI). Due to the complex ecology of arid [...] Read more.
This study aims to assess the spatiotemporal changes in ecological environment quality (EEQ) in arid regions, using Xinjiang as a case study, from 2000 to 2023, with an improved remote sensing ecological index (IRSEI). Due to the complex ecology of arid regions, the traditional remote sensing ecological index (RSEI) has limitations in capturing ecological dynamics. To address this, we propose an enhanced IRSEI model that replaces normalization with standardization, improving robustness against outliers. Additionally, the kernel normalized difference vegetation index (kNDVI) and normalized difference salinity index (NDSI) are integrated to assess saline areas more effectively. The methodology includes time series analysis, spatial distribution analysis, and statistical evaluations using the difference method, coefficient of variation, and the Hurst index. Results show that the IRSEI more accurately reflects ecological dynamics than the RSEI. Temporal analysis reveals stable overall EEQ, with some areas improving. Spatially, the environment is generally better in the north and in mountainous regions than in the south and plains. Statistical evaluations suggest a positive trend in ecological changes, with improved areas surpassing degraded ones. This study contributes to the monitoring, protection, and management of arid region ecosystems, emphasizing the need for high-resolution data and further analysis. Full article
Show Figures

Figure 1

Figure 1
<p>Geographical location of the study area and land use type distribution map.</p>
Full article ">Figure 2
<p>Technology roadmap for this research. The dashed box part on the left side is the process of constructing IRSEI in the GEE platform, and the solid line part on the right side is the process of performing spatio-temporal change analysis.</p>
Full article ">Figure 3
<p><span class="html-italic">IRSEI</span> and RSEI Pixel changes in the diagonal direction (the two solid black lines L1 and L2).</p>
Full article ">Figure 4
<p>(<b>A</b>–<b>D</b>) Comparative analysis of the <span class="html-italic">IRSEI</span> and the RSEI in the study area in 2020. (<b>A</b>) is a comparison in cropland; (<b>B</b>) is a comparison in built-up land; (<b>C</b>) is a comparison in grassland; (<b>D</b>) is a comparison in unutilized land.</p>
Full article ">Figure 5
<p>2001 and 2020 LC type proportions.</p>
Full article ">Figure 6
<p>Percentage of area in each class of <span class="html-italic">IRSEI</span> in Xinjiang, 2000–2023.</p>
Full article ">Figure 7
<p>Spatial distribution of the <span class="html-italic">IRSEI</span> in Xinjiang for six time periods from 2000 to 2023.</p>
Full article ">Figure 8
<p>Spatial distribution map of Xinjiang’s annual average monthly precipitation and annual average temperature.</p>
Full article ">Figure 9
<p>Evolution of the <span class="html-italic">IRSEI</span> grade in Xinjiang from 2000 to 2023 based on the difference method in six time periods.</p>
Full article ">Figure 10
<p>Spatial distribution of the evolution trend and sustainability of the <span class="html-italic">IRSEI</span> ranking in Xinjiang based on the difference method and hurst index. (<b>A</b>) is the spatial distribution and area share of IRSEI trend changes; (<b>B</b>) is the spatial distribution and area share of IRSEI persistent trends.</p>
Full article ">Figure 11
<p>Stability of spatial distribution of <span class="html-italic">IRSEI</span> in Xinjiang from 2000 to 2023 based on the coefficient of variation.</p>
Full article ">Figure 12
<p>Future trends in EEQ in Xinjiang.</p>
Full article ">Figure 13
<p>Spatial distribution and area proportion of LC in Xinjiang in 2020.</p>
Full article ">
23 pages, 10390 KiB  
Article
The Influence of Spatial Scale Effect on Rock Spectral Reflectance: A Case Study of Huangshan Copper–Nickel Ore District
by Ziwei Wang, Huijie Zhao, Guorui Jia and Feixiang Wang
Remote Sens. 2024, 16(24), 4643; https://doi.org/10.3390/rs16244643 - 11 Dec 2024
Viewed by 308
Abstract
The spectral reflectance measured in situ is often regarded as the “truth”. However, its limited coverage and large spatial heterogeneity often make the ground-based reflectance unable to represent the remote sensing images. Since the spatial scale mismatch between ground-based, airborne, and spaceborne measurements, [...] Read more.
The spectral reflectance measured in situ is often regarded as the “truth”. However, its limited coverage and large spatial heterogeneity often make the ground-based reflectance unable to represent the remote sensing images. Since the spatial scale mismatch between ground-based, airborne, and spaceborne measurements, the applications of geological exploration, metallogenic prognosis and mine monitoring are facing severe challenges. In order to explore the influence of spatial scale effect on rock spectra, spectral reflectance with uncertainty caused by differences in illumination view geometry and spatial heterogeneity is introduced into the Bayesian Maximum Entropy (BME) method. Then, the rock spectra are upscaled from the point-scale to meter-scale and to 10 m-scale, respectively. Finally, the influence of spatial scale effect is evaluated based on the reflectance value, spectral shape, and spectral characteristic parameters. The results indicate that the BME model shows better upscaling accuracy and stability than Ordinary Kriging and Ordinary Least Squares model. The maximum Euclidean Distance of rock spectra caused by spatial resolution change is 6.271, and the Spectral Angle Mapper can reach 0.370. The spectral absorption position, absorption depth, and spectral absorption index are less affected by scale effect. For the area with similar spatial heterogeneity to the Huangshan Copper–Nickel Ore District, when the spatial resolution of the image is greater than 10 m, the rock’s spectrum is less influenced by the change in spatial resolution. Otherwise, the influence of spatial scale effect should be considered in applications. In addition, this work puts forward a set of processes to evaluate the influence of spatial scale effect in the study area and carry out the upscaling. Full article
Show Figures

Figure 1

Figure 1
<p>Overview of the study area and the 40 ground sample points. The bottom image is the Hymap hyperspectral image covering the study area on 22 October 2002.</p>
Full article ">Figure 2
<p>Semivariogram of representative bands of spectral reflectance measured in situ in the study area.</p>
Full article ">Figure 2 Cont.
<p>Semivariogram of representative bands of spectral reflectance measured in situ in the study area.</p>
Full article ">Figure 3
<p>The fluctuation interval in the in situ measured spectral reflectance of rocks when the SZA and SAA is inconsistent with the image data, at some sampling points.</p>
Full article ">Figure 4
<p>The fluctuation range in the in situ measured spectral reflectance of rocks caused by spatial heterogeneity at some sampling points.</p>
Full article ">Figure 4 Cont.
<p>The fluctuation range in the in situ measured spectral reflectance of rocks caused by spatial heterogeneity at some sampling points.</p>
Full article ">Figure 5
<p>Scatter plots and linear fitting results of in situ measured and upscaled spectral reflectance at 10 ground validation points in the absorption bands. “*” for a multiplication sign.</p>
Full article ">Figure 6
<p>Comparison between the in situ measured spectral reflectance and the upscaled spectral reflectance with 3 m spatial resolution at 6 ground verification points.</p>
Full article ">Figure 7
<p>Comparison of the reflectance value of rock spectra with different spatial resolutions at some sampling points.</p>
Full article ">Figure 8
<p>Histogram of maximum ED and maximum SAM distribution caused by spatial resolution change at 40 sampling points. (<b>a</b>) Maximum ED. (<b>b</b>) Maximum SAM.</p>
Full article ">Figure 9
<p>Comparison of the rock’s spectral shape with different spatial resolutions at some sampling points.</p>
Full article ">Figure 10
<p>Trend diagram and STD of rock’s spectral characteristic parameters at 2250 nm with the decrease in spatial resolution.</p>
Full article ">Figure 11
<p>Comparison of the ratios of Hymap data and upscaling results at 20 sampling points.</p>
Full article ">
21 pages, 7882 KiB  
Article
Multi-Scale Gross Ecosystem Product (GEP) Valuation for Wetland Ecosystems: A Case Study of Lishui City
by Zhixin Zhu, Keyue Wu, Shuyue Zhou, Zhe Wang and Weiya Chen
Water 2024, 16(24), 3554; https://doi.org/10.3390/w16243554 - 10 Dec 2024
Viewed by 438
Abstract
Traditional gross ecosystem product (GEP) accounting methods often operate at macro scales, failing to reflect the localized and nuanced values of wetland ecosystems. This study addresses these challenges by introducing a fine-grained classification system based on a localized adaptation of international standards. The [...] Read more.
Traditional gross ecosystem product (GEP) accounting methods often operate at macro scales, failing to reflect the localized and nuanced values of wetland ecosystems. This study addresses these challenges by introducing a fine-grained classification system based on a localized adaptation of international standards. The framework integrates high-precision national land surveys and remote sensing quantitative analysis while incorporating fisheries resource models, climate regulation beneficiary mapping, and visitor interpolation to address data scarcity related to human activities. This approach refines the spatial calculation methods for functional quantity accounting at fine scales. The results demonstrate that the refined classification maintains consistency with traditional methods in total value while adapting to multi-scale accounting, filling gaps at small and medium scales and providing a more accurate representation of localized wetland characteristics. Additionally, the study highlights the dominance of cultural services in GEP, emphasizing the need to balance cultural and regulatory services to ensure fairness in decision-making. Finally, a village-scale decision-support model is proposed, offering actionable guidance for wetland management and sustainable development planning. Full article
(This article belongs to the Special Issue Hydro-Economic Models for Sustainable Water Resources Management)
Show Figures

Figure 1

Figure 1
<p>Layout of study area.</p>
Full article ">Figure 2
<p>Wetland classification mapping process.</p>
Full article ">Figure 3
<p>Spatial calculation methods for functional quantities maps based on interpolation optimization. (<b>a</b>) Fishery suitability mapping and adjustment. (<b>b</b>) Wetland climate regulation beneficiaries analysis and (<b>c</b>) wetland tourism distribution.</p>
Full article ">Figure 4
<p>GEP accounting process: In the ecosystem product amout part, the “·” symbol in the figure indicates data calculations based on coefficients or models specified in the standard, while the “+” symbol denotes the spatial calculations or adjustment coefficients added for optimizing functional quantities in the research design.</p>
Full article ">Figure 5
<p>Lishui wetland classification map.</p>
Full article ">Figure 6
<p>Multi-scale wetland waterbody area statistics.</p>
Full article ">Figure 7
<p>Gross ecosystem product result in multi scale.</p>
Full article ">Figure 8
<p>Wetland area vs. GEP value with Z-score selections.</p>
Full article ">Figure 9
<p>Cluster analysis of regulatory vs. cultural service contributions.</p>
Full article ">Figure 10
<p>Village-level wetland development strategy map.</p>
Full article ">
13 pages, 5892 KiB  
Article
Detecting Sensitive Spectral Bands and Vegetation Indices for Potato Yield Using Handheld Spectroradiometer Data
by Diego Gomez, Pablo Salvador, Juan Fernando Rodrigo and Jorge Gil
Plants 2024, 13(23), 3436; https://doi.org/10.3390/plants13233436 - 7 Dec 2024
Viewed by 587
Abstract
Remote sensing is a valuable tool in precision agriculture due to its spatial and temporal coverage, non-destructive method of data collection, and cost-effectiveness. In this study, we measured the canopy reflectance of potato (Solanum tuberosum L.) crops on a plant-by-plant basis with [...] Read more.
Remote sensing is a valuable tool in precision agriculture due to its spatial and temporal coverage, non-destructive method of data collection, and cost-effectiveness. In this study, we measured the canopy reflectance of potato (Solanum tuberosum L.) crops on a plant-by-plant basis with a handheld spectrometer instrument. Our study pursues two primary objectives: (1) determining the optimal temporal aggregation for measuring canopy signals related to potato yield and (2) identifying the best spectral bands in the 350–2500 nm domain and vegetation indices. The study was conducted over two consecutive years (2020 and 2021) with 60 plants per plot, encompassing six potato varieties and three replicates annually throughout the growth season. Employing correlation analysis and dimensionality reduction, we identified 23 independent features significantly correlated with tuber yield. We used multiple linear regression analysis to model the relationship between the selected features and yield and to compare their influence in the fitted model. We used the Leave-One-Out Cross-Validation (LOOCV) method to assess the validity of the model (RMSE = 702 g and %RMSE = 29.2%). The most significant features included the Gitelson2 and Vogelmann indices. The optimal time period for measurements was determined to be from 56 to 100 days after planting. These findings may contribute to the advancement of precision farming by proposing tailored sensor applications, paving the way for improved agricultural practices and enhanced food security. Full article
(This article belongs to the Special Issue Precision Agriculture in Crop Production)
Show Figures

Figure 1

Figure 1
<p>Annual plot layout for the experimental study, showing three separate plots per year. Each plot consists of six rows, one for each potato variety, with 10 plants per row. This results in a total of 60 plants per plot and 180 plants across all three plots for the year.</p>
Full article ">Figure 2
<p>Location of the fields and plots in the study area. The upper image shows the location of Castilla y Leon (yellow) within Spain, with the province of Segovia (purple). (<b>A</b>) corresponds to the plot location (red) for the year 2020 and (<b>B</b>) for the year 2021. The (<b>C</b>–<b>E</b>) corresponds to photos taken during the sowing, ASD measurement, and harvest, respectively.</p>
Full article ">Figure 3
<p>Summary statistics of tuber weight (kg), total number of tubers (n.tuber), and number of tubers according to their diameter size (&gt;75 mm for large, 40 to 75 mm for medium, and &lt;40 mm for small) for each experimental plot, variety, and given year. (<b>A</b>,<b>C</b>,<b>E</b>) correspond to the plots for the year 2020, and (<b>B</b>,<b>D</b>,<b>F</b>) to those for the year 2021.</p>
Full article ">Figure 4
<p>Comparison of rf results (RMSE) based on bootstrapping, expressed as 95%CI for time aggregation and statistic aggregation.</p>
Full article ">Figure 5
<p>Scatter plot of measured yield and predicted yield for different potato varieties for 2020 and 2021. Points correspond to the held-out test set at each iteration of LOOCV.</p>
Full article ">Figure 6
<p>Standardized residual plot showing characteristics of varieties.</p>
Full article ">Figure 7
<p>Variable importance expressed in standardized units (centered and scaled). Y-axis shows the selected variables after the dimensionality reduction process. X-axis shows, in absolute values, the standardized scores of (red) absolute Pearson correlation coefficients with yield, and (blue) importance scores from MLR model using permutation feature importance and 95% CI error bars.</p>
Full article ">
22 pages, 27970 KiB  
Article
Monthly Prediction of Pine Stress Probability Caused by Pine Shoot Beetle Infestation Using Sentinel-2 Satellite Data
by Wen Jia, Shili Meng, Xianlin Qin, Yong Pang, Honggan Wu, Jia Jin and Yunteng Zhang
Remote Sens. 2024, 16(23), 4590; https://doi.org/10.3390/rs16234590 - 6 Dec 2024
Viewed by 414
Abstract
Due to the significant threat to forest health posed by beetle infestations on pine trees, timely and accurate predictions are crucial for effective forest management. This study developed a pine tree stress probability prediction workflow based on monthly cloud-free Sentinel-2 composite images to [...] Read more.
Due to the significant threat to forest health posed by beetle infestations on pine trees, timely and accurate predictions are crucial for effective forest management. This study developed a pine tree stress probability prediction workflow based on monthly cloud-free Sentinel-2 composite images to address this challenge. First, representative pine tree stress samples were selected by combining long-term forest disturbance data using the Continuous Change Detection and Classification (CCDC) algorithm with high-resolution remote sensing imagery. Monthly cloud-free Sentinel-2 images were then composited using the Multifactor Weighting (MFW) method. Finally, a Random Forest (RF) algorithm was employed to build the pine tree stress probability model and analyze the importance of spectral, topographic, and meteorological features. The model achieved prediction precisions of 0.876, 0.900, and 0.883, and overall accuracies of 89.5%, 91.6%, and 90.2% for January, February, and March 2023, respectively. The results indicate that spectral features, such as band reflectance and vegetation indices, ranked among the top five in importance (i.e., SWIR2, SWIR1, Red band, NDVI, and NBR). They more effectively reflected changes in canopy pigments and leaf moisture content under stress compared with topographic and meteorological features. Additionally, combining long-term stress disturbance data with high-resolution imagery to select training samples improved their spatial and temporal representativeness, enhancing the model’s predictive capability. This approach provides valuable insights for improving forest health monitoring and uncovers opportunities to predict future beetle outbreaks and take preventive measures. Full article
(This article belongs to the Section Forest Remote Sensing)
Show Figures

Figure 1

Figure 1
<p>Location of the study area in Ning’er County, Puer City, Yunnan Province, China, overlaid on a false-color Sentinel-2 image (R, G, B = SWIR1, NIR, Red bands). The yellow dashed line delineates the study area’s boundaries.</p>
Full article ">Figure 2
<p>Field survey of pine stress.</p>
Full article ">Figure 3
<p>Overall technical workflow for predicting monthly pine stress probability.</p>
Full article ">Figure 4
<p>Reference data based on stress disturbance results. (<b>a</b>) The monthly stress disturbance results from 2019 to 2023; (<b>b</b>) An example of reference sample points displayed on the GF-1, GF-2, Sentinel-2, and Landsat-8 imagery; (<b>c</b>) The spatial distribution of non-stress sample points selected through visual interpretation; (<b>d</b>) The spatial distribution of pine stress sample points selected through visual interpretation.</p>
Full article ">Figure 5
<p>Comparison of monthly cloud-free Sentinel-2 composite images and vegetation indices from January to March 2023. The images are displayed as false-color composites (RGB = SWIR1, NIR, Red). A specific site was selected for detailed close-up analysis, showing the imagery, NDVI, and NDWI of the pine stress area affected by beetle infestation.</p>
Full article ">Figure 6
<p>Feature importance ranking for pine stress prediction model.</p>
Full article ">Figure 7
<p>Predicted pine stress probability for January, February, and March 2023 (<b>left</b>) and spatial distribution of areas with probability greater than 80% (<b>right</b>).</p>
Full article ">Figure 8
<p>Site 1: Monthly increase in pine stress level and area from January to March 2023, with Sentinel-2 imagery and stress probability distribution. The dash circles are key focus areas of forest stress.</p>
Full article ">Figure 9
<p>Site 2: Monthly decrease in pine stress level and area from January to March 2023, with Sentinel-2 imagery and stress probability distribution. The dash circles are key focus areas of forest stress.</p>
Full article ">Figure 10
<p>Site 3: Monthly changes (increase and decrease) in pine stress levels and areas from January to March 2023, with Sentinel-2 imagery and stress probability distribution. The dash circles are key focus areas of forest stress.</p>
Full article ">
21 pages, 14898 KiB  
Article
Analysis of Economic Vitality and Development Equilibrium of China’s Three Major Urban Agglomerations Based on Nighttime Light Data
by Saimiao Liu, Wenliang Liu, Yi Zhou, Shixin Wang, Zhenqing Wang, Zhuochen Wang, Yanchao Wang, Xinran Wang, Luoyao Hao and Futao Wang
Remote Sens. 2024, 16(23), 4571; https://doi.org/10.3390/rs16234571 - 6 Dec 2024
Viewed by 421
Abstract
Eliminating poverty, reducing inequality, and achieving balanced development are one of the United Nations Sustainable Development Goals. Objectively and accurately measuring regional economic vitality and development equilibrium is a pressing scientific issue that needs to be addressed in order to achieve common prosperity. [...] Read more.
Eliminating poverty, reducing inequality, and achieving balanced development are one of the United Nations Sustainable Development Goals. Objectively and accurately measuring regional economic vitality and development equilibrium is a pressing scientific issue that needs to be addressed in order to achieve common prosperity. Nighttime light (NTL) remote sensing data have been proven to be a good proxy variable for socio-economic development, and are widely used due to their advantages of convenient access and wide spatial coverage. Based on multi-source data, this study constructs an Economic Development Index (EDI) that comprehensively reflects regional economic vitality from two aspects, economic quality and development potential, combines the Nighttime Light Development Index (NLDI) as the evaluation indicators to measure the economic vitality and development equilibrium, analyzes the economic vitality and development equilibrium of 300 district and county units in China’s three major urban agglomerations from 2000 to 2020 and their temporal and spatial variation characteristics, and discusses the connotation of EDI and its availability. The results show the following: (1) From 2000 to 2020, the average growth rate of EDI in China’s three major urban agglomerations reached 36.32%, while the average decrease rate of NLDI reached 38.75%; both economic vitality and the development equilibrium have been continuously enhanced. Among them, the Yangtze River Delta (YRD) urban agglomeration experienced the fastest economic growth, while the Guangdong–Hong Kong–Macao Greater Bay Area (GBA) exhibited the strongest economic strength. (2) Both economic vitality and the development equilibrium in these three urban agglomerations exhibited distinct spatial agglomeration characteristics, namely center-surrounding distribution, coastal–inland distribution, and radial belt–pole distribution, respectively. (3) Over the past two decades, the economic development of these three urban agglomerations has progressed towards the pattern of regional coordinated development, pole-driven development and urban–rural integrated development. The research results can provide new research perspectives and scientific support for promoting regional balanced development, achieving sustainable development goals, and reducing inequality. Full article
Show Figures

Figure 1

Figure 1
<p>Study areas and their nighttime light remote sensing images in 2020.</p>
Full article ">Figure 2
<p>Schematic diagram for the calculation of the NTL.</p>
Full article ">Figure 3
<p>EDI and NLDI statistics of the three urban agglomerations from 2000 to 2020 ((<b>a</b>,<b>b</b>) UA represents the average value of the three urban agglomerations, (<b>c</b>,<b>d</b>) Distribution of EDI and NLDI Interval Quantity in the three urban agglomerations from 2000 to 2020).</p>
Full article ">Figure 4
<p>Spatial distribution of EDI of three urban agglomerations from 2000 to 2020.</p>
Full article ">Figure 5
<p>Spatial distribution of NLDI of three urban agglomerations from 2000 to 2020.</p>
Full article ">Figure 6
<p>EDI hotspot maps of three major urban agglomerations from 2000 to 2020.</p>
Full article ">Figure 7
<p>NLDI hotspot maps of three major urban agglomerations from 2000 to 2020.</p>
Full article ">Figure 8
<p>Trends of EDI and NLDI of three urban agglomerations from 2000 to 2020.</p>
Full article ">Figure 9
<p>(<b>a</b>) Spatial distribution patterns of NPP-VIIRS-like. (<b>b</b>) Spatial distribution patterns of resampled SDGSAT-1. Comparison of NTL spatial distribution patterns.</p>
Full article ">Figure 10
<p>Comparison between high economic vitality areas and urban built-up areas.</p>
Full article ">
16 pages, 3324 KiB  
Article
Matching Vegetation Indices and Tree Vigor in Pyrenean Silver Fir Stands
by Juan Pablo Crespo-Antia, Antonio Gazol, Manuel Pizarro, Ester González de Andrés, Cristina Valeriano, Álvaro Rubio Cuadrado, Juan Carlos Linares and Jesús Julio Camarero
Remote Sens. 2024, 16(23), 4564; https://doi.org/10.3390/rs16234564 - 5 Dec 2024
Viewed by 442
Abstract
Forest health monitoring is crucial for sustainable management, especially with the challenges posed by climate warming. Remote sensing data provide vegetation indices, such as the Normalized Difference Vegetation Index (NDVI) and Enhanced Vegetation Index (EVI), that are widely used in assessing forest health. [...] Read more.
Forest health monitoring is crucial for sustainable management, especially with the challenges posed by climate warming. Remote sensing data provide vegetation indices, such as the Normalized Difference Vegetation Index (NDVI) and Enhanced Vegetation Index (EVI), that are widely used in assessing forest health. However, studies considering the validation of these data with field assessments of tree vigor are still scarce. To address this issue, we explored the relationships in declining (D) and non-declining (N) silver fir (Abies alba Mill.) stands from the Spanish Pyrenees between changes in canopy (a proxy of vigor), vegetation indices (NDVI, EVI) and climate variables. We compared trends in the NDVI and EVI for the period of 1984–2023 for D and N stands showing high and low crown defoliation levels, respectively. The EVI values allowed for the separation of stands according to their vigor earlier and more clearly than NDVI values, which did not show clear patterns throughout the time series. Significant negative correlations were found between the EVI and stand defoliation (r = −0.57) or mean radial growth (r = 0.81). Late-spring drought reduced the EVI. The EVI series reflected similar spatial patterns in terms of stand defoliation and tree growth, offering complementary information, along with the strengths of remote sensing with respect to its spatial and temporal coverage, for the early detection of forest dieback. This study also contributes to a better understanding of remote sensing indices, which is useful for forest health monitoring. Full article
Show Figures

Figure 1

Figure 1
<p>Spatial distribution of the studied silver fir stands in the Spanish Pyrenees. (<b>a</b>) Map of all study sites, with pink polygons representing declining stands, green polygons indicating non-declining stands, and orange squares showing locations where comparisons of stand defoliation or growth and vegetation indices were conducted. (<b>b</b>) Map of Europe and distribution range of silver fir, highlighting the study region (pink square) in the Pyrenees. (<b>c</b>) Panoramic view of a declining stand in Salvatierra de Escá. The white dashed line in panels (<b>a</b>,<b>b</b>) indicates silver fir distribution obtained from [<a href="#B11-remotesensing-16-04564" class="html-bibr">11</a>]. The color scale in panels (<b>a</b>,<b>b</b>) indicates climatic water deficit (CWD), with blue indicating higher moisture availability and red indicating moisture deficit. The CWD was calculated for the period of 1984–2023.</p>
Full article ">Figure 2
<p>Climatic trends for declining (D, black circles) and non-declining (N, white triangles) silver fir stands. Panel (<b>a</b>) shows trends for vapor pressure deficit (VPD) and panel (<b>b</b>) for climatic water deficit (CWD). Linear regressions are fitted for both D and N stands.</p>
Full article ">Figure 3
<p>Box plots of EVI and NDVI indices showing differences between declining (D) and non-declining (N) silver fir stands from 1984 to 2024 in the Spanish Pyrenees. The panels show standard deviation of NDVI (NDVI SD), 95% decile slope of NDVI (NDVI 95), slope of NDVI (NDVI slope) and Enhanced Vegetation Index (EVI).</p>
Full article ">Figure 4
<p>Long-term trends in EVI for declining (dieback) and non-declining (no dieback) silver fir stands from 1984 to 2024. Error bars represent standard error (SE).</p>
Full article ">Figure 5
<p>Relationships between defoliation, tree growth and EVI in silver fir stands. Panel (<b>a</b>) shows the correlation between the mean percentage of defoliated silver fir stands in 2000 and 2020. Panel (<b>b</b>) shows the correlation between the mean tree-ring width and the mean percentage of defoliated silver fir trees in 2020 for both declining (D, black circles) and non-declining (N, white triangles) stands. Panel (<b>c</b>) shows the correlation between EVI and the mean percentage of defoliated silver fir trees in 2020. The black lines represent linear regressions. The thick gray line in panel (<b>a</b>) represents the 1:1 relationship for reference. Linear regressions are fitted with corresponding equations, showing correlation and <span class="html-italic">p</span> values in each panel.</p>
Full article ">Figure 6
<p>Coefficient estimates for the relationships between average climatic variables (CWD, VPD, Tmin, Tmax, soil moisture, precipitation) and EVI for silver fir stands. The black dots represent the coefficient estimates and the bars represent the standard error associated with the coefficient. The red vertical lines indicate the null value (0) for reference.</p>
Full article ">Figure 7
<p>Monthly correlations calculated between the 1-month SPEI drought index or climatic variables (VPD, CWD, Tmax, Tmin, precipitation, soil moisture) and EVI for declining (D) and non-declining (N) silver fir stands. Correlations were obtained from the previous October to the current (growth year) August. Significant correlations are marked with an asterisk (*) for <span class="html-italic">p</span> &lt; 0.05.</p>
Full article ">
18 pages, 15356 KiB  
Article
Implications of Pulse Frequency in Terrestrial Laser Scanning on Forest Point Cloud Quality and Individual Tree Structural Metrics
by Tom E. Verhelst, Kim Calders, Andrew Burt, Miro Demol, Barbara D’hont, Joanne Nightingale, Louise Terryn and Hans Verbeeck
Remote Sens. 2024, 16(23), 4560; https://doi.org/10.3390/rs16234560 - 5 Dec 2024
Viewed by 582
Abstract
Terrestrial laser scanning (TLS) provides highly detailed 3D information of forest environments but is limited to small spatial scales, as data collection is time consuming compared to other remote sensing techniques. Furthermore, TLS data collection is heavily dependent on wind conditions, as the [...] Read more.
Terrestrial laser scanning (TLS) provides highly detailed 3D information of forest environments but is limited to small spatial scales, as data collection is time consuming compared to other remote sensing techniques. Furthermore, TLS data collection is heavily dependent on wind conditions, as the movement of trees negatively impacts the acquired data. Hardware advancements resulting in faster data acquisition times have the potential to be valuable in upscaling efforts but might impact overall data quality. In this study, we investigated the impact of the pulse repetition rate (PRR), or pulse frequency, which is the number of laser pulses emitted per second by the scanner. Increasing the PRR reduces the scan time required for a single scan but decreases the power (amplitude) of the emitted laser pulses commensurately. This trade-off could potentially impact the quality of the acquired data. We used a RIEGL VZ400i laser scanner to test the impact of different PRR settings on the point cloud quality and derived tree structural metrics from individual tree point clouds (diameter, tree height, crown projected area) as well as quantitative structure models (total branch length, tree volume). We investigated this impact across five field plots of different forest complexity and canopy density for three different PRR settings (300, 600 and 1200 kHz). The scan time for a single scan was 180, 90 and 45 s for 300, 600 and 1200 kHz, respectively. Differences among the raw acquired scans from different PRR replicates were largely removed by several necessary data processing steps, notably the removal of uncertain points with a low reflectance attribute. We found strong agreement between the individual tree structural metrics derived from each of the PRR replicates, independent of the forest complexity. This was the case for both point cloud-based metrics and those derived from quantitative structural models (QSMs). The results demonstrate that the PRR in high-end TLS instruments can be increased for data collection with negligible impact on a selection of derived structural metrics that are commonly used in the context of aboveground biomass estimation. Full article
Show Figures

Figure 1

Figure 1
<p>Overview of the different terrestrial laser scanning (TLS) sites, using pictures captured from the scanner. Left to right: BEoff, BEon, AUS, GAB, UK.</p>
Full article ">Figure 2
<p>Flowchart of the processing workflow that was followed to process the data from every PRR replicate (300, 600 and 1200 kHz) of each forest plot (BEoff, BEon, AUS, GAB, UK). Along the processing chain, the outputs that were compared among the different PRR replicates are indicated (green circles). On the level of individual scans, the differences in four point attribute distributions among the PRR replicates were investigated: reflectance, deviation, range and return number. The differences in the reflectance and deviation attributes were assessed on the raw collected scans, whereas the range and return number distributions were assessed for the data after a noise-filtering step (based on reflectance and deviation). After this, the individual scans were co-registered and individual tree point clouds were segmented from the resulting plot point clouds. The segmented point clouds were then used to investigate differences among the PRR replicates on the level of individual trees, looking at both the point distribution and sparsity along the tree height as well as the derived tree structural metrics. The extracted tree structural metrics were diameter at breast height (DBH), tree height (H), crown projected area (CPA), branch length (BL) and total tree volume (V).</p>
Full article ">Figure 3
<p>Overview of the reflectance and deviation distributions in the raw data collected from the same single upright scan position with different pulse repetition rate (PRR) settings (top to bottom: BEoff, BEon, AUS, GAB, UK). The distributions show the number of points (millions) across the reflectance and distribution ranges. The grey bands indicate the areas from the distributions that were filtered out in the processing steps to remove high uncertainty points.</p>
Full article ">Figure 4
<p>Overview of the differences in the amount of remaining points in a single scan among the different pulse repetition rate replicates (PRR; 300, 600 and 1200 kHz) after the attribute filter step (top to bottom: BEoff, BEon, AUS, GAB, UK). The left side of the figure shows the amount of points in function of the range from the scanner, whereas the right side shows the number of points in each return number class (e.g., return number class 2 indicates the second hit from a laser beam). Both y-axes show the number of points in millions. The bar chart on the range plots indicates the total number of points in every PRR replicate. The return number graphs show the cumulative distribution of the points in function of the return number of the points. Thus, this shows the contribution of each consequent return number class on the total amount of points. This distribution is shown for both the filtered and the filtered + downsampled point clouds. At return number 1, the amount of single returns (no further hits) is indicated separately through the dotted lines.</p>
Full article ">Figure 5
<p>Detailed overview of the final tree point cloud for the tallest tree in both the Beon and GAB plots. These are the tree point clouds after filtering, segmentation and 1 cm downsampling. The tree point clouds are plotted in the middle panel for visual reference. The left panel shows the distribution of the amount of points in the tree point cloud (dashed lines) as well as the mean nearest neighbour (four closest neighbours) distance (solid lines) in function of the height. The right panel shows a zoomed-in subset of a 5 m × 5 m × 5 m voxel from the tree canopies in each of the PRR replicates, to visually assess the occlusion impact of the PRR. The points in the zoomed-in subsets are coloured using the x-values (axis perpendicular to the plane), to better display depth. Overviews of the tallest trees from the other plots can be found in <a href="#app1-remotesensing-16-04560" class="html-app">Appendix A</a> <a href="#remotesensing-16-04560-f0A2" class="html-fig">Figure A2</a>.</p>
Full article ">Figure 6
<p>Overview of the individual tree metrics for all plots; values extracted from the 300 kHz replicate are plotted as a reference versus the values from the 1200 kHz PRR replicate. The overview of the 300 kHz replicate metrics versus the 600 kHz replicate values can be found in <a href="#app1-remotesensing-16-04560" class="html-app">Appendix A</a> <a href="#remotesensing-16-04560-f0A3" class="html-fig">Figure A3</a>. The plots are ordered in rows per structural metric (diameter at breast height DBH [cm], tree height H [m], crown projected area CPA [m<sup>2</sup>]), branch length BL [km] and total tree volume V [m<sup>3</sup>], respectively) and in columns per plot (BEoff, BEon, AUS, GAB, UK, respectively). The orange diagonal 1:1 line is plotted on every subplot. The concordance correlation coefficient (CCC) and its 95% confidence interval CI is given on every subplot, to demonstrate the agreement between the metric values from the 300 and 1200 kHz replicates.</p>
Full article ">Figure A1
<p>Overview of the reflectance and deviation distributions in the raw data of a single tilt scan with different PRR settings. These are the tilt scans from the same scan locations as the upright scans from <a href="#remotesensing-16-04560-f001" class="html-fig">Figure 1</a> in <a href="#sec3dot1-remotesensing-16-04560" class="html-sec">Section 3.1</a>. The deviation values are cut off at 20, but the values continue until 500+. The grey bands indicate the areas from the distributions that were filtered out in the processing, to remove high-uncertainty points.</p>
Full article ">Figure A2
<p>Detailed overview of the tallest tree in the BEoff, AUS and UK plot. The tree point clouds are plotted in the middle panel for visual reference. The left panel shows the distribution of the amount of points in the tree point cloud (dashed lines) as well as the mean nearest neighbour (four closest neighbours) distance (solid lines) in function of the height. The right panel shows a zoomed-in subset of a 5 m × 5 m × 5 m voxel from the tree canopy in each of the PRR replicates, to visually assess the occlusion impact of the PRR. The points in the zoomed-in subsets are coloured using the x-values (axis perpendicular to the plane), to better display depth.</p>
Full article ">Figure A3
<p>Overview of the individual tree metrics for all plots; values extracted from the 300 kHz replicate are plotted as a reference versus the values from the 600 kHz PRR replicate, in the same way as <a href="#remotesensing-16-04560-f005" class="html-fig">Figure 5</a>. The plots are ordered in rows per structural metric (DBH (m), H (m), CPA (m<sup>2</sup>), BL (m) and V (m<sup>3</sup>), respectively) and in columns per plot (BEoff, BEon, AUS, GAB, UK, respectively). The orange diagonal 1:1 line is plotted on every subplot. The concordance correlation coefficient (CCC) and its 95% confidence interval CI is given on every subplot, to demonstrate the agreement between the metric values from the 300 and 600 kHz replicates.</p>
Full article ">
18 pages, 4141 KiB  
Article
The Coupling Coordination Relationship Between Urbanization and the Eco-Environment in Resource-Based Cities, Loess Plateau, China
by Shuaizhi Kang, Xia Jia, Yonghua Zhao, Manya Luo, Huanyuan Wang and Ming Zhao
ISPRS Int. J. Geo-Inf. 2024, 13(12), 437; https://doi.org/10.3390/ijgi13120437 - 4 Dec 2024
Viewed by 505
Abstract
Resource-based cities face numerous sustainability challenges, making the coupled and coordinated relationship between urbanization and the eco-environment critical for sustainable development strategies. The Loess Plateau is an essential energy base and ecologically fragile area in China, holding unique and significant research value. This [...] Read more.
Resource-based cities face numerous sustainability challenges, making the coupled and coordinated relationship between urbanization and the eco-environment critical for sustainable development strategies. The Loess Plateau is an essential energy base and ecologically fragile area in China, holding unique and significant research value. This research employed the Remote Sensing Ecological Index (RSEI) and the Compound Night Light Index (CNLI), based on MODIS and night light data, to investigate the socio-economic development and eco-environmental changes across 25 resource-based cities on the Loess Plateau (LP) in China over the past 20 years. The Coupling Coordination Degree Model (CCDM) and Multi-Scale Geographically Weighted Regression (MGWR) were utilized to assess the relationship between urbanization and ecological factors. The average RSEI values for these cities ranged from 0.4524 to 0.4892 over the 20 years, reflecting an upward trend with a growth rate of 8.13%. Simultaneously, the average CNLI values ranged from 1.5700 to 6.0864, with a change of 4.5164. Over the past two decades, all cities in the study area experienced rapid urbanization and ecological development. The correlation between urbanization and ecological factors strengthened, alongside an increasing spatial heterogeneity. While the coupling coordination relationship in most cities showed improvement, many remained within the low to middle grades. These findings enhance the understanding of the intricate relationships between urbanization and ecology, offering valuable insights for policy-making aimed at creating sustainable and livable resource-based cities. Full article
Show Figures

Figure 1

Figure 1
<p>Development stage and location of resource-based cities in the Loess Plateau area.</p>
Full article ">Figure 2
<p>The RSEI map and area percentage in 2000, 2005, 2010, 2015, and 2020.</p>
Full article ">Figure 3
<p>The RSEI mean value distribution between 2000 and 2020.</p>
Full article ">Figure 4
<p>Resource-based cites in China’s Loess Plateau light images and changes from 2000 to 2020.</p>
Full article ">Figure 5
<p>CNLI of resource-based cites in China’s Loess Plateau and its changes.</p>
Full article ">Figure 6
<p>Spatial distribution of coefficients of multi-scale geographically weighted regression model.</p>
Full article ">Figure 7
<p>The coordination degree between eco-environment and urbanization in China’s Loess Plateau.</p>
Full article ">Figure 8
<p>The temporal changes of the CCD between eco-environment and urbanization in China Loess Plateau about prefecture-level cities.</p>
Full article ">
21 pages, 6350 KiB  
Article
Spatiotemporal Urban Evolution Along the China–Laos Railway in Laos Determined Using Multiple Sources of Remote-Sensed Landscape Indicators and Interpretable Machine Learning
by Dongxue Li, Jin Tang, Qiao Hu, Mingjuan Dong and Soukanh Chithpanya
Land 2024, 13(12), 2094; https://doi.org/10.3390/land13122094 - 4 Dec 2024
Viewed by 563
Abstract
Constructing high-speed railways (HSRs) is critical for developing countries to stimulate economic growth and urbanization. This study focuses on the Lao section of the China–Laos Railway (CLR) and employs explicitly spatial remote sensing images to investigate the urban development surrounding HSR stations. Data-driven [...] Read more.
Constructing high-speed railways (HSRs) is critical for developing countries to stimulate economic growth and urbanization. This study focuses on the Lao section of the China–Laos Railway (CLR) and employs explicitly spatial remote sensing images to investigate the urban development surrounding HSR stations. Data-driven machine learning and causal inference approaches are integrated to quantify the spatial–temporal evolution and discover its driving factors. The results suggest that the CLR has had positive spatial spillover effects on the development of the surrounding urban space. These spillover effects have exhibited a distance attenuation pattern, reflecting obvious development in 2D rather than in 3D urban space. Meanwhile, the distance to stations and adjacent city centers as well as functional urban characteristics, such as land use patterns and industrialization level, have significantly influences the surrounding spatial development. Specifically, in industrial-dominated cities, the surrounding spatial changes have been most significant under the influence of the HSR. Change related to industrial and residential land use has shown significant land expansion patterns and increased utilization efficiency, reflecting that industrialization and urbanization have been the primary drivers of land demand surrounding the HSR. The findings offer valuable insights and references for developing nations to formulate and implement spatial management policies and initiatives related to HSR. Full article
Show Figures

Figure 1

Figure 1
<p>Research workflow.</p>
Full article ">Figure 2
<p>HSR stations’ distribution along the Laos section of the CLR.</p>
Full article ">Figure 3
<p>LEI of density and height around each station at different stages.</p>
Full article ">Figure 4
<p>LEI of NTL around each station at different stages.</p>
Full article ">Figure 5
<p>Visualization of density, height, and NTL changes around the HSR stations over time.</p>
Full article ">Figure 6
<p>Characteristics of density, height, and NTL variation with distance.</p>
Full article ">Figure 7
<p>Parameter tuning of RF model.</p>
Full article ">Figure 8
<p>Relative importance of factors affecting density, height, and NTL value changes.</p>
Full article ">Figure 9
<p>PDP diagrams of factors influencing spatial development. Since the spatial resolution of the remote sensing imagery in this study used was 100 m; the unit of the x-axis in (<b>a</b>,<b>b</b>) is consistent with the physical distance represented by the pixels (100 m). (<b>a</b>,<b>b</b>) explain the effects of continuous variables, (<b>c</b>,<b>d</b>) explain the effects of discrete variables.</p>
Full article ">Figure 9 Cont.
<p>PDP diagrams of factors influencing spatial development. Since the spatial resolution of the remote sensing imagery in this study used was 100 m; the unit of the x-axis in (<b>a</b>,<b>b</b>) is consistent with the physical distance represented by the pixels (100 m). (<b>a</b>,<b>b</b>) explain the effects of continuous variables, (<b>c</b>,<b>d</b>) explain the effects of discrete variables.</p>
Full article ">Figure 10
<p>Threshold effect of spatial density and NTL around Vang Vieng railway station.</p>
Full article ">Figure 11
<p>Variations in land use patterns around stations in different cities.</p>
Full article ">Figure 12
<p>PDP analysis of “Distance to city center” (x<sub>4</sub>) in NTL during the CLR operation phase.</p>
Full article ">Figure 13
<p>Mechanisms of HSR’s impact on surrounding spatial development.</p>
Full article ">
Back to TopTop