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

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

Article Types

Countries / Regions

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Search Results (224)

Search Parameters:
Keywords = 3D GPR

Order results
Result details
Results per page
Select all
Export citation of selected articles as:
21 pages, 18859 KiB  
Article
Polarisation Synthesis Applied to 3D Polarimetric Imaging for Enhanced Buried Object Detection and Identification
by Samuel J. I. Forster, Anthony J. Peyton, Frank J. W. Podd and Nigel Davidson
Remote Sens. 2024, 16(22), 4279; https://doi.org/10.3390/rs16224279 (registering DOI) - 17 Nov 2024
Viewed by 174
Abstract
Detecting sub-surface objects poses significant challenges, partly due to attenuation of the ground medium and cluttered environments. The acquisition polarisation and antenna orientation can also yield significant variation of detection performance. These challenges can be mitigated by developing more versatile systems and algorithms [...] Read more.
Detecting sub-surface objects poses significant challenges, partly due to attenuation of the ground medium and cluttered environments. The acquisition polarisation and antenna orientation can also yield significant variation of detection performance. These challenges can be mitigated by developing more versatile systems and algorithms to enhance detection and identification. In this study, a novel application of a 3D SAR inverse algorithm and polarisation synthesis was applied to ultra-wideband polarimetric data of buried objects. The principle of polarisation synthesis facilitates an adaptable technique which can be used to match the target’s polarisation characteristics, and the application of this revealed hidden structures, enhanced detection, and increased received power when compared to single polarisation results. This study emphasises the significance of polarimetry in ground-penetrating radar (GPR), particularly for target discrimination in high-lift-off applications. The findings offer valuable insights that could drive future research and enhance the performance of these sensing systems. Full article
(This article belongs to the Section Remote Sensing Image Processing)
Show Figures

Figure 1

Figure 1
<p>Imaging geometry in 3D.</p>
Full article ">Figure 2
<p>The polarisation ellipse shown with ellipticity angle <math display="inline"><semantics> <mi>χ</mi> </semantics></math>, orientation angle <math display="inline"><semantics> <mi>φ</mi> </semantics></math>, wave amplitudes <math display="inline"><semantics> <msub> <mi>A</mi> <mi>x</mi> </msub> </semantics></math> and <math display="inline"><semantics> <msub> <mi>A</mi> <mi>y</mi> </msub> </semantics></math>, and magnitude <span class="html-italic">A</span>.</p>
Full article ">Figure 3
<p>Experimental setup showing (<b>a</b>) the positioning system, VNA, and antenna and (<b>b</b>) close-up view of the dual-polarised antenna.</p>
Full article ">Figure 4
<p>The chosen targets in the study showing (<b>a</b>) 5 cm diameter metallic sphere, (<b>b</b>,<b>c</b>) a wire in two orientations, and (<b>d</b>) air-filled cylindrical container.</p>
Full article ">Figure 5
<p>Estimated received power versus range for a 5 cm diameter metallic sphere buried in sand for <math display="inline"><semantics> <mrow> <msub> <mi>σ</mi> <mrow> <mi>c</mi> <mi>o</mi> <mi>n</mi> <mi>d</mi> </mrow> </msub> <mo>=</mo> <msup> <mn>10</mn> <mrow> <mo>−</mo> <mn>3</mn> </mrow> </msup> </mrow> </semantics></math> and <math display="inline"><semantics> <msup> <mn>10</mn> <mrow> <mo>−</mo> <mn>5</mn> </mrow> </msup> </semantics></math> <math display="inline"><semantics> <mrow> <msup> <mi>Sm</mi> <mrow> <mo>−</mo> <mn>1</mn> </mrow> </msup> </mrow> </semantics></math>.</p>
Full article ">Figure 6
<p>SAR images of the metallic sphere with (<b>a</b>) HH; (<b>b</b>) HV; (<b>c</b>) VV polarisations.</p>
Full article ">Figure 7
<p>Comparison of sphere target cross-sections with (HH 1) and without (HH 2) refraction corrections, shown in (<b>a</b>) linear scale and (<b>b</b>) dB scale.</p>
Full article ">Figure 8
<p>SAR images of the straight insulated wire with (<b>a</b>) HH; (<b>b</b>) HV; (<b>c</b>) VV polarisations.</p>
Full article ">Figure 9
<p>SAR images of the curved insulated wire with (<b>a</b>) HH; (<b>b</b>) HV; (<b>c</b>) VV polarisations.</p>
Full article ">Figure 10
<p>SAR images of the air-filled cylinder with (<b>a</b>) HH; (<b>b</b>) HV; (<b>c</b>) VV polarisations.</p>
Full article ">Figure 11
<p>Orientation plot of the metallic sphere.</p>
Full article ">Figure 12
<p>Polarimetric images of the sphere synthesised in linear and circular polarisations. (<b>a</b>) Linear horizontal; (<b>b</b>) linear vertical; (<b>c</b>) linear +45°; (<b>d</b>) linear <math display="inline"><semantics> <mrow> <msup> <mrow> <mo>−</mo> <mn>45</mn> </mrow> <mo>°</mo> </msup> </mrow> </semantics></math>; (<b>e</b>) RHC; (<b>f</b>) LHC.</p>
Full article ">Figure 13
<p>Orientation plot of the insulated wire.</p>
Full article ">Figure 14
<p>Polarimetric images of the insulated wire synthesised in linear and circular polarisations. (<b>a</b>) Linear horizontal; (<b>b</b>) linear vertical; (<b>c</b>) linear +45°; (<b>d</b>) linear <math display="inline"><semantics> <mrow> <msup> <mrow> <mo>−</mo> <mn>45</mn> </mrow> <mo>°</mo> </msup> </mrow> </semantics></math>; (<b>e</b>) RHC; (<b>f</b>) LHC.</p>
Full article ">Figure 15
<p>Orientation plot of the curved insulated wire.</p>
Full article ">Figure 16
<p>Polarimetric images of the curved insulated wire synthesised in linear and circular polarisations. (<b>a</b>) Linear horizontal; (<b>b</b>) linear vertical; (<b>c</b>) linear +45°; (<b>d</b>) linear <math display="inline"><semantics> <mrow> <msup> <mrow> <mo>−</mo> <mn>45</mn> </mrow> <mo>°</mo> </msup> </mrow> </semantics></math>; (<b>e</b>) RHC; (<b>f</b>) LHC.</p>
Full article ">Figure 17
<p>Comparison of curved wire cross-sections for HH and RHC polarisations, shown in (<b>a</b>) linear scale and (<b>b</b>) dB scale.</p>
Full article ">Figure 18
<p>Histogram comparison for HH and RHC polarisations.</p>
Full article ">Figure 19
<p>Orientation plot of the air-filled cylinder.</p>
Full article ">Figure 20
<p>Polarimetric images of the air-filled cylinder synthesised in linear and circular polarisations. (<b>a</b>) Linear horizontal; (<b>b</b>) linear vertical; (<b>c</b>) linear +45°; (<b>d</b>) linear <math display="inline"><semantics> <mrow> <msup> <mrow> <mo>−</mo> <mn>45</mn> </mrow> <mo>°</mo> </msup> </mrow> </semantics></math>; (<b>e</b>) RHC; (<b>f</b>) LHC.</p>
Full article ">
19 pages, 4707 KiB  
Article
Chlorophyll Content Estimation of Ginkgo Seedlings Based on Deep Learning and Hyperspectral Imagery
by Zilong Yue, Qilin Zhang, Xingzhou Zhu and Kai Zhou
Forests 2024, 15(11), 2010; https://doi.org/10.3390/f15112010 - 14 Nov 2024
Viewed by 373
Abstract
Accurate estimation of chlorophyll content is essential for understanding the growth status and optimizing the cultivation practices of Ginkgo, a dominant multi-functional tree species in China. Traditional methods based on chemical analysis for determining chlorophyll content are labor-intensive and time-consuming, making them [...] Read more.
Accurate estimation of chlorophyll content is essential for understanding the growth status and optimizing the cultivation practices of Ginkgo, a dominant multi-functional tree species in China. Traditional methods based on chemical analysis for determining chlorophyll content are labor-intensive and time-consuming, making them unsuitable for large-scale dynamic monitoring and high-throughput phenotyping. To accurately quantify chlorophyll content in Ginkgo seedlings under different nitrogen levels, this study employed a hyperspectral imaging camera to capture canopy hyperspectral images of seedlings throughout their annual growth periods. Reflectance derived from pure leaf pixels of Ginkgo seedlings was extracted to construct a set of spectral parameters, including original reflectance, logarithmic reflectance, and first derivative reflectance, along with spectral index combinations. A one-dimensional convolutional neural network (1D-CNN) model was then developed to estimate chlorophyll content, and its performance was compared with four common machine learning methods, including Gaussian Process Regression (GPR), Partial Least Squares Regression (PLSR), Support Vector Regression (SVR), and Random Forest (RF). The results demonstrated that the 1D-CNN model outperformed others with the first derivative spectra, achieving higher CV-R2 and lower RMSE values (CV-R2 = 0.80, RMSE = 3.4). Furthermore, incorporating spectral index combinations enhanced the model’s performance, with the 1D-CNN model achieving the best performance (CV-R2 = 0.82, RMSE = 3.3). These findings highlight the potential of the 1D-CNN model in strengthening the chlorophyll estimations, providing strong technical support for the precise cultivation and the fertilization management of Ginkgo seedlings. Full article
(This article belongs to the Section Forest Inventory, Modeling and Remote Sensing)
Show Figures

Figure 1

Figure 1
<p>The workflow for estimating chlorophyll content in <span class="html-italic">Ginkgo</span> canopies based on hyperspectral imaging and 1D-CNN.</p>
Full article ">Figure 2
<p>Schematic representation of the experimental layout for <span class="html-italic">Ginkgo biloba</span> seedlings under five nitrogen treatments (N0–N4). Each treatment was replicated three times (R1–R3), resulting in 15 experimental units in total. Nitrogen was applied as a topdressing in three equal doses during the growing season.</p>
Full article ">Figure 3
<p>Hyperspectral images of <span class="html-italic">Ginkgo biloba</span> seedlings under different nitrogen levels (N0–N4) across growth stages (T1–T5). T1 corresponds to April (early bud development stage), T2 corresponds to May (early rapid growth stage), T3 corresponds to June (middle rapid growth stage), T4 corresponds to July (late rapid growth stage), and T5 corresponds to August (plant maturity stage).</p>
Full article ">Figure 4
<p>A suitable 1D-CNN model for spectral reflectance is proposed.</p>
Full article ">Figure 5
<p>Changes in canopy reflectance and SPAD chlorophyll content of <span class="html-italic">Ginkgo</span> seedlings across different growth stages and nitrogen fertilization levels (<b>A</b>) Spectral reflectance curves of the <span class="html-italic">Ginkgo</span> canopy at different SPAD levels. The figure includes three forms of reflectance spectra: (i) original reflectance, (ii) logarithmic reflectance, and (iii) first derivative reflectance. SPAD chlorophyll content is divided into low (SPAD_low), medium (SPAD_medium), and high (SPAD_high) levels. Low SPAD corresponds to values from 27 to 45, medium SPAD ranges between 45 and 55, and high SPAD corresponds to values from 55 to 65. Reflectance across the 400 to 900 nm range varies with SPAD levels, reflecting the sensitivity of different spectral regions to chlorophyll absorption and canopy structure. (<b>B</b>) Changes in SPAD chlorophyll content of <span class="html-italic">Ginkgo</span> seedlings at different growth stages (T1–T5). T1 represents April (early bud development stage), T2 represents May (early rapid growth stage), T3 represents June (middle rapid growth stage), T4 represents July (late rapid growth stage), and T5 represents August (plant maturity stage). SPAD content fluctuates across the different growth stages. (<b>C</b>) Changes in SPAD chlorophyll content of <span class="html-italic">Ginkgo</span> seedlings under different nitrogen fertilization treatments (N0–N4). SPAD content shows significant variation across the different nitrogen levels.</p>
Full article ">Figure 6
<p>Correlation coefficient curve between leaf SPAD-chlorophyll content in <span class="html-italic">Ginkgo</span> seedlings and original or transformed reflectance spectra.</p>
Full article ">Figure 7
<p>The correlation between leaf chlorophyll content and the three best-performing indices: SR<sub>708,775</sub>, GNDVI, and mCI<sub>Green</sub>.</p>
Full article ">Figure 8
<p>Optimal reflectance datasets for the correlation of DVI<sub>log</sub> ((<b>A</b>), logarithmic reflectance), RVI<sub>FD</sub> ((<b>B</b>), first derivative of reflectance), NDVI<sub>FD</sub> ((<b>C</b>), first derivative of reflectance), and mRVI<sub>log</sub> ((<b>D</b>), logarithmic reflectance) with chlorophyll content. The white arrow indicates the optimal band combination.</p>
Full article ">Figure 9
<p>Comparison of predicted versus measured SPAD values for <span class="html-italic">Ginkgo</span> seedling canopies using various modeling approaches. Best Spectrum-Orgi, Best Spectrum-log and Best Spectrum-FD represent the best-performing spectral data (Orgi: original spectra; log: logarithmic spectra; FD: first-derivative spectra) for each regression method. VI represents vegetation indices.</p>
Full article ">
16 pages, 39531 KiB  
Technical Note
A Geophysical Investigation in Which 3D Electrical Resistivity Tomography and Ground-Penetrating Radar Are Used to Determine Singularities in the Foundations of the Protected Historic Tower of Murcia Cathedral (Spain)
by María C. García-Nieto, Marcos A. Martínez-Segura, Manuel Navarro, Ignacio Valverde-Palacios and Pedro Martínez-Pagán
Remote Sens. 2024, 16(21), 4117; https://doi.org/10.3390/rs16214117 - 4 Nov 2024
Viewed by 550
Abstract
This study presents a procedure in which 3D electrical resistivity tomography (ERT) and ground-penetrating radar (GPR) were used to determine singularities in the foundations of protected historic towers, where space is limited due to their characteristics and location in highly populated areas. This [...] Read more.
This study presents a procedure in which 3D electrical resistivity tomography (ERT) and ground-penetrating radar (GPR) were used to determine singularities in the foundations of protected historic towers, where space is limited due to their characteristics and location in highly populated areas. This study was carried out on the Tower of the Cathedral “Santa Iglesia Catedral de Santa María” in Murcia, Spain. The novel distribution of a continuous nonlinear profile along the outer and inner perimeters of the Tower allowed us to obtain a 3D ERT model of the subsoil, even under its load-bearing walls. This nonlinear configuration of the electrodes allowed us to reach adequate investigation depths in buildings with limited interior and exterior space for data collection without disturbing the historic structure. The ERT results were compared with GPR measurements and with information from archaeological excavations conducted in 1999 and 2009. The geometry and distribution of the cavities in the entire foundation slab of the Tower were determined, verifying the proposed procedure. This methodology allows the acquisition of a detailed understanding of the singularities of the foundations of protected historic towers in urban areas with limited space, reducing time and costs and avoiding the use of destructive techniques, with the aim of implementing a more efficient and effective strategy for the protection of other tower foundations. Full article
(This article belongs to the Special Issue 3D Virtual Reconstruction for Cultural Heritage (Second Edition))
Show Figures

Figure 1

Figure 1
<p>Location map of the Tower of Murcia Cathedral: projection of the Tower in Murcia, showing the east façade and its main parts.</p>
Full article ">Figure 2
<p>Continuous ERT profile. (<b>a</b>) Measuring equipment; (<b>b</b>) ERT profile in the east façade of the Tower; (<b>c</b>) electrodes placed in the northwest corner of the Tower.</p>
Full article ">Figure 3
<p>The electrodes placed inside the Tower (sacristy) and a detailed view of the arrangement, consisting of an electrode with an aluminium plate, a steel spring, and a carbomer-based gel.</p>
Full article ">Figure 4
<p>Continuous ERT profile: (<b>a</b>) location of the electrodes inside and outside the Tower and measuring equipment between electrodes 28 and 29; (<b>b</b>) an example of a measurement sequence of the electrodes used in a 3D array.</p>
Full article ">Figure 5
<p>Position of the pole: (<b>a</b>) distance from the pole to the measuring equipment (80 m); (<b>b</b>) location of the pole outside the west façade of the Cathedral; (<b>c</b>) detail of the pole.</p>
Full article ">Figure 6
<p>(<b>a</b>) Location plan for the profiles made with the 250 and 500 MHz antennas; (<b>b</b>) measurements made with the 500 MHz antenna; (<b>c</b>) measurements made with the 250 MHz antenna.</p>
Full article ">Figure 7
<p>A 3D ERT model characterising the area underneath the Tower in terms of subsurface electrical resistivity values.</p>
Full article ">Figure 8
<p>Radargrams obtained inside the Tower with the 250 MHz antenna. The significant reflections found are highlighted by red rectangles.</p>
Full article ">Figure 9
<p>Radargrams obtained inside the Tower with the 500 MHz antenna. The significant reflections found are highlighted by red rectangles.</p>
Full article ">Figure 10
<p>A 3D model of the radargrams obtained in the east-west direction inside the Tower with the 500 MHz antenna.</p>
Full article ">Figure 11
<p>(<b>a</b>) A 3D ERT model positioned under the Tower; (<b>b</b>) a detailed view of the model with the location set according to the floor plan of the Tower of the highly resistive zones.</p>
Full article ">Figure 12
<p>(<b>a</b>) Chamber located in the northeast corner (Geocisa, 2009 [<a href="#B30-remotesensing-16-04117" class="html-bibr">30</a>]); (<b>b</b>) interior of one of the chambers to which access was gained [<a href="#B30-remotesensing-16-04117" class="html-bibr">30</a>]; (<b>c</b>) actions carried out in 1999 in the interior of the Antesacristía (photograph taken by Juan Antonio Molina Serrano [<a href="#B30-remotesensing-16-04117" class="html-bibr">30</a>]); (<b>d</b>) location according to the floor plan of the Tower’s cavities; (<b>e</b>) northeast corner in 2009, with the original plinth and archaeological remains of a rammed-earth wall [<a href="#B30-remotesensing-16-04117" class="html-bibr">30</a>]; (<b>f</b>,<b>g</b>) project carried out in 2009, involving the filling of the trenches with draining material and the creation of an aeration chamber. Photographs taken by José Antonio Sánchez Pravia. Images by Geocisa [<a href="#B30-remotesensing-16-04117" class="html-bibr">30</a>].</p>
Full article ">
20 pages, 17200 KiB  
Article
What Is Beyond Hyperbola Detection and Characterization in Ground-Penetrating Radar Data?—Implications from the Archaeological Site of Goting, Germany
by Tina Wunderlich, Bente S. Majchczack, Dennis Wilken, Martin Segschneider and Wolfgang Rabbel
Remote Sens. 2024, 16(21), 4080; https://doi.org/10.3390/rs16214080 - 31 Oct 2024
Viewed by 490
Abstract
Hyperbolae in radargrams are caused by a variety of small subsurface objects. The analysis of their curvature enables the determination of propagation velocity in the subsurface, which is important for exact time-to-depth conversion and migration and also yields information on the water content [...] Read more.
Hyperbolae in radargrams are caused by a variety of small subsurface objects. The analysis of their curvature enables the determination of propagation velocity in the subsurface, which is important for exact time-to-depth conversion and migration and also yields information on the water content of the soil. Using deep learning methods and fitting (DLF) algorithms, it is possible to automatically detect and analyze large numbers of hyperbola in 3D Ground-Penetrating Radar (GPR) datasets. As a result, a 3D velocity model can be established. Combining the hyperbola locations and the 3D velocity model with reflection depth sections and timeslices leads to improved archaeological interpretation due to (1) correct time-to-depth conversion through migration with the 3D velocity model, (2) creation of depthslices following the topography, (3) evaluation of the spatial distribution of hyperbolae, and (4) derivation of a 3D water content model of the site. In an exemplary study, we applied DLF to a 3D GPR dataset from the multi-phased (2nd to 12th century CE) archaeological site of Goting on the island of Föhr, Northern Germany. Using RetinaNet, we detected 38,490 hyperbolae in an area of 1.76 ha and created a 3D velocity model. The velocities ranged from approximately 0.12 m/ns at the surface to 0.07 m/ns at approx. 3 m depth in the vertical direction; in the lateral direction, the maximum velocity variation was ±0.048 m/ns. The 2D-migrated radargrams and subsequently created depthslices revealed the remains of a longhouse, which was not known beforehand and had not been visible in the unmigrated timeslices. We found hyperbola apex points aligned along linear strong reflections. They can be interpreted as stones contained in ditch fills. The hyperbola points help to differentiate between ditches and processing artifacts that have a similar appearance as the ditches in time-/depthslices. From the derived 3D water content model, we could identify the thickness of the archaeologically relevant layer across the whole site. The layer contains a lot of humus and has a high water retention capability, leading to a higher water content compared to the underlying glacial moraine sand, which is well-drained. Full article
(This article belongs to the Special Issue Advanced Ground-Penetrating Radar (GPR) Technologies and Applications)
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) The processed radargrams of all 16 channels have been cut into images of 256 × 256 pixels. In all images of the first channel, hyperbolae were labeled manually. These images and ground truth labels were used to train a RetinaNet. (<b>b</b>) During inference, all images of all channels were fed into the trained RetinaNet and hyperbolae were automatically detected (red, blue and yellow boxes), which is shown here for four example images. (<b>c</b>) Examples of the velocity determination workflow (see yellow and blue box in (<b>b</b>)): After detection of all hyperbolae, the boxes were extracted from the images and subsequently a thresholding and C3 clustering algorithm was applied. This results in central points of a cluster, which are used in a x<sup>2</sup>-(t/2)<sup>2</sup>-diagram for fitting of a linear equation. We can derive the velocity from the slope of this fit, and from the intercept, the apex time t<sub>0</sub>.</p>
Full article ">Figure 2
<p>(<b>a</b>) Mean velocity distribution along a west–east profile taken from the smoothed 3D data cube. (<b>b</b>) Missing velocity values at early or late times are extrapolated constantly from the existing ones and smoothed with a moving average filter over 150 samples. (<b>c</b>) Lateral extrapolation and smoothing with a box filter over 7 cells.</p>
Full article ">Figure 3
<p>(<b>a</b>) Location of the island of Föhr in northern Germany (Basemap: OpenStreetMap); (<b>b</b>) magnetic map of Goting, where the numbers refer to different zones of the settlement, which are further explained in the text; (<b>c</b>) zoomed in view of the magnetic map with interpretation and the GPR measurement area; (<b>d</b>) topographic overview based on a digital elevation model (Basemaps: DGM1, DOP20 ©GeoBasis-DE/LVermGeo SH/CC BY 4.0).</p>
Full article ">Figure 4
<p>Detected hyperbola in different time ranges on the complete area (Basemap: ©GeoBasis-DE/LVermGeo SH/CC BY 4.0).</p>
Full article ">Figure 5
<p>Timeslices of interpolated RMS velocities from the 3D velocity model (Basemap: ©GeoBasis-DE/LVermGeo SH/CC BY 4.0).</p>
Full article ">Figure 6
<p>Density maps of hyperbola points in different time intervals for grid cells of 2 m × 2 m in a radius of 5 m. The approximate depth ranges were calculated with a mean 1D velocity function, which was derived by calculating the average velocity in each time interval over the complete area (Basemap: ©GeoBasis-DE/LVermGeo SH/CC BY 4.0).</p>
Full article ">Figure 7
<p>(<b>a</b>) Magnetic map with location of smaller area (red) shown in (<b>b</b>). The purple line shows the position of the example profile in <a href="#remotesensing-16-04080-f009" class="html-fig">Figure 9</a>. (<b>b</b>) Cut-out from the magnetic map. (<b>c</b>,<b>e</b>) Timeslices. (<b>d</b>,<b>f</b>) Depthslices following the topography with a given depth below the surface. (<b>g</b>) Interpretation of magnetics (green), timeslice 10–15 ns (yellow), and timeslice 15–20 ns (blue) (main features). (<b>h</b>) Interpretation of magnetics (green), depthlice 60–80 cm (yellow), and depthslice 80–100 cm (blue) (main features). The red dotted box marks the location of the cut-out shown in <a href="#remotesensing-16-04080-f008" class="html-fig">Figure 8</a> (Basemap: ©GeoBasis-DE/LVermGeo SH/CC BY 4.0).</p>
Full article ">Figure 8
<p>Depthslice between 80 and 100 cm (for location, see red dotted box in <a href="#remotesensing-16-04080-f007" class="html-fig">Figure 7</a>h) and marked outline of the longhouse (blue dotted line).</p>
Full article ">Figure 9
<p>Comparison of slices and radargrams: (<b>a</b>) Timeslices between 15 and 20 ns with location of exemplary radargram, (<b>b</b>) processed radargram (Profile 109 of channel 1), (<b>c</b>) smoothed velocity model along the profile as an overlay on the radargram, (<b>d</b>) depthslice between 60 and 80 cm below the surface, following the topography, (<b>e</b>) migrated radargram using the velocity model in (<b>c</b>) and taking the topography into account, (<b>f</b>) water content distribution as an overlay on the migrated radargram. The time or depth ranges of the slices are shown as purple lines in the radargrams, respectively. The yellow color marks the remains of a pit house, whereas the blue anomaly marks a well. Green and orange arrows in (<b>a</b>,<b>b</b>) indicate the locations of isolated hyperbolae forming linear features in the time-/depthslices, interpreted as ditches. (<b>g</b>) Exemplary 1D velocity and water content functions. Their locations are marked by the blue/reg/green lines in (<b>c</b>,<b>f</b>). Dashed lines represent interval velocities, bold lines RMS velocity functions.</p>
Full article ">Figure 10
<p>(<b>a</b>) Depthslice between 80 and 100 cm with hyperbola locations. Red arrows indicate the alignment of points in ditches. (<b>b</b>) Photography of an excavation profile in the eastern part of the area with a simplified interpretation. A stone is located in a ditch filling (Photograph by J. Gebühr, NIhK Wilhelmshaven; full description in [<a href="#B22-remotesensing-16-04080" class="html-bibr">22</a>], pp. 550–551, Profile 17).</p>
Full article ">Figure 11
<p>Volumetric water content maps at different depths (Basemap: ©GeoBasis-DE/LVermGeo SH/CC BY 4.0).</p>
Full article ">Figure 12
<p>Comparison of depth calculations using a mean 1D velocity model, a mean 1D velocity model varied by 0.05 m/ns, or constant velocity models.</p>
Full article ">
20 pages, 6895 KiB  
Article
Three-Dimensional Reconstruction of Retaining Structure Defects from Crosshole Ground Penetrating Radar Data Using a Generative Adversarial Network
by Donghao Zhang, Zhengzheng Wang, Yu Tang, Shengshan Pan and Tianming Pan
Remote Sens. 2024, 16(21), 3995; https://doi.org/10.3390/rs16213995 - 28 Oct 2024
Viewed by 454
Abstract
Crosshole ground penetrating radar (GPR) is an efficient method for ensuring the quality of retaining structures without the need for excavation. However, interpreting crosshole GPR data is time-consuming and prone to inaccuracies. To address this challenge, we proposed a novel three-dimensional (3D) reconstruction [...] Read more.
Crosshole ground penetrating radar (GPR) is an efficient method for ensuring the quality of retaining structures without the need for excavation. However, interpreting crosshole GPR data is time-consuming and prone to inaccuracies. To address this challenge, we proposed a novel three-dimensional (3D) reconstruction method based on a generative adversarial network (GAN) to recover 3D permittivity distributions from crosshole GPR images. The established framework, named CGPR2VOX, integrates a fully connected layer, a residual network, and a specialized 3D decoder in the generator to effectively translate crosshole GPR data into 3D permittivity voxels. The discriminator was designed to enhance the generator’s performance by ensuring the physical plausibility and accuracy of the reconstructed models. This adversarial training mechanism enables the network to learn non-linear relationships between crosshole GPR data and subsurface permittivity distributions. CGPR2VOX was trained using a dataset generated through finite-difference time-domain (FDTD) simulations, achieving precision, recall and F1-score of 91.43%, 96.97% and 94.12%, respectively. Model experiments validate that the relative errors of the estimated positions of the defects were 1.67%, 1.65%, and 1.30% in the X-, Y-, and Z-direction, respectively. Meanwhile, the method exhibits noteworthy generalization capabilities under complex conditions, including condition variations, heterogeneous materials and electromagnetic noise, highlighting its reliability and effectiveness for practical quality assurance of retaining structures. Full article
Show Figures

Figure 1

Figure 1
<p>Schematic diagram of the crosshole GPR detection process: (<b>a</b>) antennas setup, (<b>b</b>) wavefield snapshot, (<b>c</b>) waveforms received at multiple depths and (<b>d</b>) final crosshole GPR data block.</p>
Full article ">Figure 2
<p>Structure of the CGPR2VOX network.</p>
Full article ">Figure 3
<p>Crosshole GPR dataset preparation for 3D reconstruction.</p>
Full article ">Figure 4
<p>Evaluation of the training effects of CGPR2VOX: (<b>a</b>) training losses and (<b>b</b>) defect recognition performance.</p>
Full article ">Figure 5
<p>Example of crosshole GPR 3D permittivity reconstruction on the testing dataset.</p>
Full article ">Figure 6
<p>Verification results for different defect condition scenarios: (<b>a</b>) FDTD model consistent with the training dataset, (<b>b</b>) FDTD model with a 50% deviation in defect permittivity, (<b>c</b>) FDTD model with cylindrical defect geometry, (<b>d</b>) reconstruction consistent with the training dataset, (<b>e</b>) reconstruction with a 50% deviation in defect permittivity, (<b>f</b>) reconstruction with cylindrical defect geometry.</p>
Full article ">Figure 7
<p>Numerical simulation for verification in complex scenarios: (<b>a</b>) FDTD simulation, (<b>b</b>) ground truth model.</p>
Full article ">Figure 8
<p>Results of the forward simulation and defect reconstruction under complex scenarios: (<b>a</b>) A-scan waveform, (<b>b</b>) B-scan image and (<b>c</b>) reconstruction results.</p>
Full article ">Figure 9
<p>The experimental system: (<b>a</b>) photograph, (<b>b</b>) exploded view diagram.</p>
Full article ">Figure 10
<p>Step-frequency crosshole GPR system.</p>
Full article ">Figure 11
<p>Photograph of the model box. (<b>a</b>) Visual effects, (<b>b</b>) measurement results of the electrical parameters of the filling materials.</p>
Full article ">Figure 12
<p>Crosshole GPR experiment data measurement: (<b>a</b>) photograph of the experiment site, (<b>b</b>) measured B-scan GPR images.</p>
Full article ">Figure 13
<p>Reconstruction result of the experimental data: (<b>a</b>) photograph of ground truth, (<b>b</b>) result of CGPR2VOX, (<b>c</b>) result of ray-based tomography, and (<b>d</b>) result of probabilistic FWI.</p>
Full article ">
25 pages, 13122 KiB  
Article
Comparative Study of GPR Acquisition Methods for Shallow Buried Object Detection
by Primož Smogavec, Blaž Pongrac, Andrej Sarjaš, Venceslav Kafedziski, Nabojša Dončov and Dušan Gleich
Remote Sens. 2024, 16(21), 3931; https://doi.org/10.3390/rs16213931 - 22 Oct 2024
Viewed by 480
Abstract
This paper investigates the use of ground-penetrating radar (GPR) technology for detecting shallow buried objects, utilizing an air-coupled stepped frequency continuous wave (SFCW) radar system that operates within a 2 GHz bandwidth starting at 500 MHz. Different GPR data acquisition methods for air-coupled [...] Read more.
This paper investigates the use of ground-penetrating radar (GPR) technology for detecting shallow buried objects, utilizing an air-coupled stepped frequency continuous wave (SFCW) radar system that operates within a 2 GHz bandwidth starting at 500 MHz. Different GPR data acquisition methods for air-coupled systems are compared, specifically down-looking, side-looking, and circular acquisition strategies, employing the back projection algorithm to provide focusing of the acquired GPR data. Experimental results showed that the GPR can penetrate up to 0.6 m below the surface in a down-looking mode. The developed radar and the back projection focusing algorithm were used to acquire data in the side-looking and circular mode, providing focused images with a resolution of 0.1 m and detecting subsurface objects up to 0.3 m below the surface. The proposed approach transforms B-scans of the GPR-based data into 2D images. The provided approach has significant potential for advancing shallow object detection capabilities by transforming hyperbola-based features into point-like features. Full article
(This article belongs to the Section Urban Remote Sensing)
Show Figures

Figure 1

Figure 1
<p>Geometry illustrating Snell’s law.</p>
Full article ">Figure 2
<p>Backprojection algorithm.</p>
Full article ">Figure 3
<p>(<b>a</b>) Modular SFCW radar system, (<b>b</b>) signal synthesis board, (<b>c</b>) block diagram of modular SFCW radar.</p>
Full article ">Figure 4
<p>Noise level for modular SFCW radar.</p>
Full article ">Figure 5
<p>S11, S22, (<b>a</b>) and S12 (<b>b</b>) parameters of TEM horn antennas.</p>
Full article ">Figure 6
<p>(<b>a</b>) Soil-box used as a test polygon of <math display="inline"><semantics> <mrow> <mn>3</mn> <mo>×</mo> <mn>3</mn> <mo>×</mo> <mn>0.5</mn> </mrow> </semantics></math> m. (<b>b</b>) TEM horn antennas [<a href="#B31-remotesensing-16-03931" class="html-bibr">31</a>].</p>
Full article ">Figure 7
<p>Different configuration of the transmitter and receiver. (<b>a</b>) Monostatic GPR, (<b>b</b>) monostatic side-looking GPR.</p>
Full article ">Figure 8
<p>Side-looking configuration. The platform is at the center position, and it is rotating around the <span class="html-italic">Z</span> axis.</p>
Full article ">Figure 9
<p>Circular SAR geometry. Platform is moving along circular trajectory.</p>
Full article ">Figure 10
<p>(<b>a</b>) Position of corners placed on the top of the surface. Corners were separated by 0.5 m. (<b>b</b>) Geometry of the data acquisition for side-looking SAR.</p>
Full article ">Figure 11
<p>The scenarios are depicted as top-down views within a 3 × 3 m polygon. (<b>a</b>) Trajectory of the down-looking scenario. (<b>b</b>) Scenario of side-looking experiment while rotating around <span class="html-italic">z</span>-axis. (<b>c</b>) Scenario of side-looking experiment with circular trajectory.</p>
Full article ">Figure 12
<p>Photograph of soil structure.</p>
Full article ">Figure 13
<p>B-Scans with removed background of targets buried at different depths. (<b>a</b>) B-scan of target buried at 0.4 m using VV polarization. (<b>b</b>) B-scan of target buried at 0.4 m using HH polarization. (<b>c</b>) B-scan of target buried at 0.6 m using VV polarization. (<b>d</b>) B-scan of target buried at 0.6 m using HH polarization. (<b>e</b>) B-scan of target buried at 0.8 m using VV polarization. (<b>f</b>) B-scan of target buried at 0.8 m using HH polarization.</p>
Full article ">Figure 14
<p>Monostatic down-looking B-scan with corner reflector on top of the surface: (<b>a</b>) HH polarization, (<b>b</b>) VV polarization; 20 cm below the surface: (<b>c</b>) HH polarization, (<b>d</b>) VV polarization. Background removal was applied to all images.</p>
Full article ">Figure 15
<p>Scene with surface targets and side-looking configuration: (<b>a</b>) B-scan of side-looking configuration in HH polarization using background subtraction. (<b>b</b>) B-scan of side-looking configuration in VV polarization using background subtraction. (<b>c</b>) Focused image using back-projection algorithm in HH polarization. (<b>d</b>) Focused image using back-projection algorithm in VV polarization.</p>
Full article ">Figure 16
<p>Scene with buried targets and side-looking configuration: (<b>a</b>) B-scan of side-looking configuration in HH polarization using background subtraction. (<b>b</b>) B-scan of side-looking configuration in VV polarization using background subtraction. (<b>c</b>) Focused image using BP algorithm in HH polarization. (<b>d</b>) Focused image using BP algorithm in VV polarization.</p>
Full article ">Figure 17
<p>(<b>a</b>) B-scan of circular GPR using HH polarization and background removal. (<b>b</b>) B-scan of circular GPR using VV polarization and background removal. (<b>c</b>) Focused image using HH polarization and BP algorithm. (<b>d</b>) Focused image using VV polarization and BP algorithm.</p>
Full article ">Figure 18
<p>Results using buried targets. (<b>a</b>) B-scan of circular GPR using HH polarization and background removal. (<b>b</b>) B-scan of circular GPR using VV polarization and background removal. (<b>c</b>) Focused image using HH polarization and BP algorithm. (<b>d</b>) Focused image using VV polarization and BP algorithm.</p>
Full article ">Figure 19
<p>(<b>a</b>) AP land-mine buried 20 cm below surface. (<b>b</b>) Focused image using back-projection algorithm.</p>
Full article ">
27 pages, 31281 KiB  
Article
Tracking Moisture Dynamics in a Karst Rock Formation Combining Multi-Frequency 3D GPR Data: A Strategy for Protecting the Polychrome Hall Paintings in Altamira Cave
by Vicente Bayarri, Alfredo Prada, Francisco García, Carmen De Las Heras and Pilar Fatás
Remote Sens. 2024, 16(20), 3905; https://doi.org/10.3390/rs16203905 - 21 Oct 2024
Viewed by 878
Abstract
This study addresses the features of the internal structure of the geological layers adjacent to the Polychrome Hall ceiling of the Cave of Altamira (Spain) and their link to the distribution of moisture and geological discontinuities mainly as fractures, joints, bedding planes and [...] Read more.
This study addresses the features of the internal structure of the geological layers adjacent to the Polychrome Hall ceiling of the Cave of Altamira (Spain) and their link to the distribution of moisture and geological discontinuities mainly as fractures, joints, bedding planes and detachments, using 3D Ground Penetrating Radar (GPR) mapping. In this research, 3D GPR data were collected with 300 MHz, 800 MHz and 1.6 GHz center frequency antennas. The data recorded with these three frequency antennas were combined to further our understanding of the layout of geological discontinuities and how they link to the moisture or water inputs that infiltrate and reach the ceiling surface where the rock art of the Polychrome Hall is located. The same 1 × 1 m2 area was adopted for 3D data acquisition with the three antennas, obtaining 3D isosurface (isoattribute-surface) images of internal distribution of moisture and structural features of the Polychrome Hall ceiling. The results derived from this study reveal significant insights into the overlying karst strata of Polychrome Hall, particularly the interface between the Polychrome Layer and the underlying Dolomitic Layer. The results show moisture patterns associated with geological features such as fractures, joints, detachments of strata and microcatchments, elucidating the mechanisms driving capillary rise and water infiltration coming from higher altitudes. The study primarily identifies areas of increased moisture content, correlating with earlier observations and enhancing our understanding of water infiltration patterns. This underscores the utility of 3D GPR as an essential tool for informing and putting conservation measures into practice. By delineating subsurface structures and moisture dynamics, this research contributes to a deeper analysis of the deterioration processes directly associated with the infiltration water both in this ceiling and in the rest of the Cave of Altamira, providing information to determine its future geological and hydrogeological evolution. Full article
(This article belongs to the Special Issue Multi-data Applied to Near-Surface Geophysics (Second Edition))
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) Location map of the Cave of Altamira (<b>b</b>) Projection of the Cave of Altamira on the orthoimage, indicating in a red frame the position of (<b>c</b>) the Polychrome Hall with the rock art orthoimage overlayed.</p>
Full article ">Figure 2
<p>Overview of previous GPR studies conducted in the Altamira Cave and Polychrome Hall ceiling: (<b>a</b>) 100 MHz exterior grid 3 × 3 m, (<b>b</b>) 400 MHz exterior grid 1 × 1 m, (<b>c</b>) 900 MHz exterior grid 1 × 1 m, (<b>d</b>) 400 and 900 MHz interior reflection profiles (indicated with white lines), and location of area around ALT1 (indicated with a red square), (<b>e</b>) 1.6 GHz interior grid 5 × 5 cm in 2022.</p>
Full article ">Figure 3
<p>(<b>a</b>) The location of the control area ALT1 on the ceiling of the Polychrome Hall. The dashed light-green line marks the path of the central fracture that runs across the ceiling from west to east. (<b>b</b>) A photograph showing a general view of the ALT1 study area, with the same dashed light-green line indicating the central crack, which has been sealed with cement.</p>
Full article ">Figure 4
<p>(<b>a</b>) Orthoimage of the Polychrome ceiling with the sectorisation grid and the location of the control area ALT1. (<b>b</b>) View of the ALT1 control area together with the sectorisation grids and the working areas for the 300/800 MHz and 1.6 GHz antennas and (<b>c</b>) overlapping of the working areas for 300/800 MHz and 1.6 GHz antennas with basins, streams and drip points.</p>
Full article ">Figure 5
<p>General workflow.</p>
Full article ">Figure 6
<p>Photographs showing the GPR data collection set-ups in control area ALT1. (<b>a</b>) with the 300–800 MHz dual-frequency antenna and (<b>b</b>) with the 1.6 GHz antenna.</p>
Full article ">Figure 7
<p>(<b>a</b>) Location of the two grids with an equidistance of 5 cm of the profiles using the 1.6 GHz antenna in the study areas ALT1-A (in white dashed square) and ALT1-B (in yellow dashed square). (<b>b</b>) Location of the two 10 cm × 10 cm grids using the 300/800 MHz dual frequency antenna in the study areas ALT1-C (in white dashed square) and ALT1-D (in yellow dashed square).</p>
Full article ">Figure 8
<p>3D model including basins (in red), streams (in light blue), drip points (in black) and the working areas (green boxes) for 1.6 GHz and 300/800 MHz antennas.</p>
Full article ">Figure 9
<p>GPR reflection profiles processed in the ALT1 control area showing the main subsurface features (discontinuities), in particular the central fracture, with their stratigraphic correlations of the overlying layers of the ceiling for each frequency antenna used: (<b>a</b>) Location of P63 and P8 profiles which spatially coincide, (<b>b</b>) reflection profile P63 antenna 1.6 GHz; (<b>c</b>) reflection profile P8 antenna 800 MHz; (<b>d</b>) reflection profile P8 antenna 300 MHz; (<b>e</b>) schematic stratigraphy of the Polychrome Hall ceiling according to [<a href="#B2-remotesensing-16-03905" class="html-bibr">2</a>,<a href="#B3-remotesensing-16-03905" class="html-bibr">3</a>], showing the depth range reached with the 1.6 GHz, 800 MHz and 300 MHz frequency antennas.</p>
Full article ">Figure 10
<p>(<b>a</b>) location of the drip points ALT1_1, ALT1_2 and ALT1_4 involved in migration and pigment loss within the control area ALT1. (<b>b</b>) Evolution of the internal moisture zones detected by the 1.6 GHz antenna directly involved in the above mentioned dropping points.</p>
Full article ">Figure 11
<p>South area of the ceiling of the Polychrome Hall affected by carbonate concretions produced by the access of seepage water through the fracture zones.</p>
Full article ">Figure 12
<p>Evolution of moisture in the ALT1 control area. The highlighted boxes indicate significant changes in moisture depths relative to depth. In February 2022, higher moisture zones were observed at lower depths, while in April 2023, lower moisture zones were recorded at lower depths. The moisture recorded between April and June, which was involved in paint drips, has already reached the surface of the Ceiling; therefore, it was not detected by the 1.6 GHz antenna on 26 April 2023. The blue arrow connects the moisture moisture zones that, year after year, are associated with pigment dragging drips in ALT1.</p>
Full article ">Figure 13
<p>Overlying of the 3D horizontal sections framing the Dolomitic Layer in the 70 cm to 100 cm depth interval in the ALT1-C zone, captured using a 1.6 GHz GPR antenna. The figure shows multiple zones of irregularly distributed moisture, projected on the ceiling surface of the Polychrome Hall with the micro-basins and with the active drip points (water/moisture zones in shades of purple; active drip points in blue) for zones (<b>a</b>) ALT1-C and (<b>b</b>) ALT1-D.</p>
Full article ">Figure 14
<p>Pigment particles and other mineral deposits carried away by the drip water and collected in the containers located on the ground perpendicular to the drip points. SEM analysis of the detached particles indicates the existence of microcorrosion processes of the supporting rock associated with the presence of infiltration water, condensation and CO<sub>2</sub> concentration.</p>
Full article ">Figure 15
<p>Perspective view of the discontinuities revealing zones of higher moisture, represented by a blue to white gradient, across the different strata of the Polychrome ceiling. Data were collected using three GPR antennas: a 1.6 GHz antenna for the upper 50 cm, an 800 MHz antenna for depths up to 1.2 m, and a 300 MHz antenna for depths up to 4.9 m. This complex network of discontinuities highlights the primary pathways for fluid exchange within the interior of the Polychrome Hall, as well as seepage from the active drip points located on the basal surface of the Polychrome Layer in the ALT1 control area.</p>
Full article ">Figure 16
<p>(<b>a</b>) Location of the drip points within the control area ALT1. (<b>b</b>) Mean dip of the Polychrome Hall calculated from the 3D model (<b>c</b>) Moistures associated with two large vertically developing fractures in the Polychrome ceiling. The red box shows the moistures associated with the large central fracture, the yellow box shows the moisture zones coming from the fracture located in the polychrome bison. The red arrows as well as the degree of dip of the ceiling show the relationship of this with the moistures associated with the dripping points with migration and loss of pigment from ALT1. (<b>d</b>) Reflection profile P8 obtained with the 300 MHz antenna. The dashed red line indicates the fracture (central fracture) associated with the moisture shown in red, and the dashed yellow line corresponds to the fracture associated with the moisture marked in yellow in (<b>c</b>).</p>
Full article ">
20 pages, 4472 KiB  
Article
Loureirin B Reduces Insulin Resistance and Chronic Inflammation in a Rat Model of Polycystic Ovary Syndrome by Upregulating GPR120 and Activating the LKB1/AMPK Signaling Pathway
by Jing Wang, Zheng Huang, Zhiyong Cao, Yehao Luo, Yueting Liu, Huilu Cao, Xiusong Tang and Gang Fang
Int. J. Mol. Sci. 2024, 25(20), 11146; https://doi.org/10.3390/ijms252011146 - 17 Oct 2024
Viewed by 799
Abstract
Polycystic ovary yndrome (PCOS) is a common metabolic disorder in women, which is usually associated with insulin resistance (IR) and chronic inflammation. Loureirin B (LrB) can effectively improve insulin resistance and alleviate chronic inflammation, and in order to investigate the therapeutic effect of [...] Read more.
Polycystic ovary yndrome (PCOS) is a common metabolic disorder in women, which is usually associated with insulin resistance (IR) and chronic inflammation. Loureirin B (LrB) can effectively improve insulin resistance and alleviate chronic inflammation, and in order to investigate the therapeutic effect of LrB on polycystic ovary syndrome with insulin resistance (PCOS-IR), we conducted animal experiments. A PCOS-IR rat model was established by feeding a high-fat diet combined with letrozole (1 mg/kg·d for 21 days). The rats were treated with the GPR120 agonists TUG-891 and LrB for 4 weeks. Biochemical parameters (fasting blood glucose, total cholesterol, triglycerides, high- and low-density lipoprotein), hormone levels (serum insulin, E2, T, LH, and FSH), and inflammatory cytokines (TNF-α, IL-1β, IL-6, and IL-18) were analyzed. Histopathological analyses of ovaries were performed using hematoxylin/eosin (H&E) staining. Real-time PCR and western blotting were used to assess GPR120, NLRP3, and caspase-1 expression in ovaries, and immunohistochemistry was used to evaluate LKB1 and AMPK protein expression. LrB reduced body weight, Lee’s index, ovarian index, ovarian area, and volume in PCOS-IR rats. It lowered fasting blood glucose, serum insulin, and HOMA-IR. LrB decreased total serum cholesterol, triglyceride, and LDL levels and increased HDL levels. It reduced serum T, LH, and LH/FSH and raised serum E2 and FSH levels. LrB downregulated the mRNA and protein expression levels of NLRP3 and Caspase-1, increased the protein and mRNA expression levels of GPR120 in rat ovaries, and increased LKB1 and AMPK protein expression in ovaries, ameliorating ovarian histopathological changes in PCOS-IR rats. Taken together, LrB upregulated GPR120, LKB1, and AMPK protein expression, downregulated NLRP3 and Caspase-1 protein expression, reduced insulin resistance and chronic inflammation, and ameliorated histopathological changes in ovarian tissues in PCOS rats, suggesting its potential as a treatment for PCOS. Full article
(This article belongs to the Section Molecular Immunology)
Show Figures

Figure 1

Figure 1
<p>Therapeutic effect of LrB on PCOS-IR rats.</p>
Full article ">Figure 2
<p>Effect of LrB on body weight and Lee’s index in PCOS-IR rats. The data are mean ± SDs, n = 10. Body weights of rats at (<b>A</b>) week 0 of treatment, (<b>B</b>) week 1, (<b>C</b>) week 2, (<b>D</b>) week 3, (<b>E</b>) week 4. (<b>F</b>) Lee’s indices of the rats in each group after four weeks of treatment. In comparison with the normal group, <sup>##</sup> <span class="html-italic">p &lt;</span> 0.01; in comparison with the PCOS-IR group, * <span class="html-italic">p &lt;</span> 0.05; ** <span class="html-italic">p &lt;</span> 0.01.</p>
Full article ">Figure 3
<p>Effect of LrB on the ovarian area, ovarian volume, and ovarian index in PCOS-IR rats. The values are means ± SDs; n = 10. (<b>A</b>) Ovarian area at four weeks of treatment. (<b>B</b>) Ovarian volume at four weeks. (<b>C</b>) Ovarian indices at four weeks. In comparison with the normal group, <sup>##</sup> <span class="html-italic">p &lt;</span> 0.01; in comparison with the PCOS-IR group, * <span class="html-italic">p &lt;</span> 0.05; ** <span class="html-italic">p &lt;</span> 0.01.</p>
Full article ">Figure 4
<p>Vaginal cytology images (magnification 10×, 40×). LrB improved the estrous cycle in PCOS-IR rats.</p>
Full article ">Figure 5
<p>Effect of LrB on lipid metabolic activity in PCOS-IR rats. The values are means ± SDs; n = 10. (<b>A</b>) TC levels at four weeks of treatment. (<b>B</b>) TG levels at four weeks. (<b>C</b>) LDL levels at four weeks. (<b>D</b>) HDL levels at four weeks. Relative to the normal group, <sup>##</sup> <span class="html-italic">p &lt;</span> 0.01; relative to the PCOS-IR group, ** <span class="html-italic">p &lt;</span> 0.01.</p>
Full article ">Figure 6
<p>Effect of LrB on fasting blood glucose and serum insulin in PCOS-IR rats; values are presented as the means ± SDs, n = 10. (<b>A</b>) INS levels at four weeks of treatment. (<b>B</b>) FBG levels at four weeks. (<b>C</b>) HOMA-IR at four weeks. Relative to the normal group, <sup>##</sup> <span class="html-italic">p &lt;</span> 0.01; relative to the PCOS-IR group, * <span class="html-italic">p &lt;</span> 0.05; ** <span class="html-italic">p &lt;</span> 0.01.</p>
Full article ">Figure 7
<p>Effects of LrB on hormone levels in PCOS-IR rats; values are means ± SDs; n = 10. (<b>A</b>) LH levels at four weeks of treatment. (<b>B</b>) T levels at four weeks. (<b>C</b>) LH/FSH ratio at four weeks. (<b>D</b>) FSH levels at four weeks. (<b>E</b>) E2 levels at four weeks. Relative to the normal group, <sup>##</sup> <span class="html-italic">p &lt;</span> 0.01; relative to the PCOS-IR group, * <span class="html-italic">p &lt;</span> 0.05; ** <span class="html-italic">p &lt;</span> 0.01.</p>
Full article ">Figure 8
<p>LrB alleviates the serum inflammatory response in PCOS-IR rats. Values are means ± SDs; n = 10. (<b>A</b>) Serum TNF-α content at four weeks of treatment. (<b>B</b>) Serum IL-1β at four weeks. (<b>C</b>) Serum IL-6 at four weeks. (<b>D</b>) Serum IL-18 at four weeks. Relative to the normal group, <sup>##</sup> <span class="html-italic">p &lt;</span> 0.01; relative to the PCOS-IR group, * <span class="html-italic">p &lt;</span> 0.05, ** <span class="html-italic">p &lt;</span> 0.01.</p>
Full article ">Figure 9
<p>Effects of LrB treatment on morphological changes in ovarian tissue (magnification 4×, 10×). Representative H&amp;E-stained ovarian tissue sections after four weeks of treatment. C, cystic follicles; CL, corpus luteum; ANF, antral follicles; ATF, atretic follicles.</p>
Full article ">Figure 10
<p>Effects of LrB on GPR120 expression in PCOS-IR rat ovaries. Data are means ± SDs; n = 3. (<b>A</b>) Western blot analysis of GPR120 and GAPDH in ovarian tissues at four weeks of treatment. (<b>B</b>) Quantitative analysis of GPR120 levels in ovarian tissues at four weeks of treatment. (<b>C</b>) Expression levels of GPR120 mRNA in ovarian tissues at four weeks of treatment. Relative to the normal group, <sup>##</sup> <span class="html-italic">p &lt;</span> 0.01; relative to the PCOS-IR group, * <span class="html-italic">p &lt;</span> 0.05; ** <span class="html-italic">p &lt;</span> 0.01.</p>
Full article ">Figure 11
<p>Effects of LrB on NLRP3 and Caspase-1 levels in the ovaries of PCOS-IR rats. The values are means ± SDs; n = 3. (<b>A</b>) Western blot analysis of NLRP3 and GAPDH in ovarian tissues at four weeks of treatment. (<b>B</b>) Quantitative analysis of NLRP3 levels in ovarian tissues at four weeks of treatment. (<b>C</b>) Expression levels of NLRP3 mRNA in ovarian tissues at four weeks of treatment. (<b>D</b>) Western blot analysis of Caspase-1 and GAPDH in ovarian tissues at four weeks of treatment. (<b>E</b>) Quantitative analysis of Caspase-1 levels in ovarian tissues at four weeks of treatment. (<b>F</b>) Expression levels of Caspase-1 mRNA in ovarian tissues at four weeks of treatment. Relative to the normal group, <sup>##</sup> <span class="html-italic">p &lt;</span> 0.01; relative to the PCOS-IR group, * <span class="html-italic">p &lt;</span> 0.05; ** <span class="html-italic">p &lt;</span> 0.01.</p>
Full article ">Figure 12
<p>The influence of LrB on expression of LKB1 and AMPK in the ovaries of PCOS-IR rats. The values are means ± SDs; n = 5. (<b>A</b>) Expression of AMPK in the ovarian tissues of each group at four weeks of treatment (magnification 10×, 40×). (<b>B</b>) Expression of LKB1 in the ovarian tissues of each group at four weeks of treatment (magnification 10×, 40×). (<b>C</b>) Positive expression of AMPK was analyzed. (<b>D</b>) Positive expression of LKB1 was analyzed. Relative to the normal group, <sup>##</sup> <span class="html-italic">p &lt;</span> 0.01; relative to the PCOS-IR group, * <span class="html-italic">p &lt;</span> 0.05; ** <span class="html-italic">p &lt;</span> 0.01.</p>
Full article ">Figure 12 Cont.
<p>The influence of LrB on expression of LKB1 and AMPK in the ovaries of PCOS-IR rats. The values are means ± SDs; n = 5. (<b>A</b>) Expression of AMPK in the ovarian tissues of each group at four weeks of treatment (magnification 10×, 40×). (<b>B</b>) Expression of LKB1 in the ovarian tissues of each group at four weeks of treatment (magnification 10×, 40×). (<b>C</b>) Positive expression of AMPK was analyzed. (<b>D</b>) Positive expression of LKB1 was analyzed. Relative to the normal group, <sup>##</sup> <span class="html-italic">p &lt;</span> 0.01; relative to the PCOS-IR group, * <span class="html-italic">p &lt;</span> 0.05; ** <span class="html-italic">p &lt;</span> 0.01.</p>
Full article ">Figure 13
<p>Animal grouping and handling procedures.</p>
Full article ">
13 pages, 3778 KiB  
Article
Preliminary Insights on Moisture Content Measurement in Square Timbers Using GPR Signals and 1D-CNN Models
by Jiaxing Guo, Huadong Xu, Yan Zhong and Kuanjie Yu
Forests 2024, 15(10), 1800; https://doi.org/10.3390/f15101800 - 14 Oct 2024
Viewed by 711
Abstract
Accurately measuring the moisture content (MC) of square timber is crucial for ensuring the quality and performance of wood products in wood processing. Traditional MC detection methods have certain limitations. Therefore, this study developed a one-dimensional convolutional neural network (1D-CNN) model based on [...] Read more.
Accurately measuring the moisture content (MC) of square timber is crucial for ensuring the quality and performance of wood products in wood processing. Traditional MC detection methods have certain limitations. Therefore, this study developed a one-dimensional convolutional neural network (1D-CNN) model based on the first 8 nanoseconds of ground-penetrating radar (GPR) signals to predict the MC of square timber. The study found that the mixed-species model exhibited effective predictive performance (R2 = 0.9864, RMSE = 0.0393) across the tree species red spruce, Dahurian larch, European white birch, and Manchurian ash (MC range 0%–133.1%), while single-species models showed even higher accuracy (R2 ≥ 0.9876, RMSE ≤ 0.0358). Additionally, the 1D-CNN model outperformed other algorithms in automatically capturing complex patterns in GPR full-waveform amplitude data. Moreover, the algorithms based on full-waveform amplitude data demonstrated significant advantages in detecting wood MC compared to those based on a traditional time–frequency feature parameter. These results indicate that the 1D-CNN model can be used to optimize the drying process and detect the MC of load-bearing timber in construction and bridge engineering. Future work will focus on expanding the dataset, further optimizing the algorithm, and validating the models in industrial applications to enhance their reliability and applicability. Full article
Show Figures

Figure 1

Figure 1
<p>Flowchart for detecting moisture content in square timber using tree radar.</p>
Full article ">Figure 2
<p>Electromagnetic wave time domain variations at different moisture content levels (10%–90%) for Manchurian ash square timbers, with an error range of ±0.5%.</p>
Full article ">Figure 3
<p>Flowchart of the 1D-CNN algorithm based on full-waveform amplitude data from GPR signals.</p>
Full article ">Figure 4
<p>1D-CNN model loss curves (<b>a</b>) and residuals (<b>b</b>) for moisture content in mixed-species square timber.</p>
Full article ">Figure 5
<p>1D-CNN model loss curves (<b>a</b>,<b>c</b>,<b>e</b>,<b>g</b>) and residuals (<b>b</b>,<b>d</b>,<b>f</b>,<b>h</b>) for moisture content in red spruce, Dahurian larch, European white birch, and Manchurian ash square timbers.</p>
Full article ">
15 pages, 2753 KiB  
Article
Assessing Soil Physical Quality in a Layered Agricultural Soil: A Comprehensive Approach Using Infiltration Experiments and Time-Lapse Ground-Penetrating Radar Surveys
by Simone Di Prima, Gersende Fernandes, Maria Burguet, Ludmila Ribeiro Roder, Vittoria Giannini, Filippo Giadrossich, Laurent Lassabatere and Alessandro Comegna
Appl. Sci. 2024, 14(20), 9268; https://doi.org/10.3390/app14209268 - 11 Oct 2024
Viewed by 854
Abstract
Time-lapse ground-penetrating radar (GPR) surveys, combined with automated infiltration experiments, provide a non-invasive approach for investigating the distribution of infiltrated water within the soil medium and creating three-dimensional images of the wetting bulb. This study developed and validated an experimental protocol aimed at [...] Read more.
Time-lapse ground-penetrating radar (GPR) surveys, combined with automated infiltration experiments, provide a non-invasive approach for investigating the distribution of infiltrated water within the soil medium and creating three-dimensional images of the wetting bulb. This study developed and validated an experimental protocol aimed at quantifying and visualizing water distribution fluxes in layered soils under both unsaturated and saturated conditions. The 3D images of the wetting bulb significantly enhanced the interpretation of infiltration data, enabling a detailed analysis of water movement through the layered system. We used the infiltrometer data and the Beerkan Estimation of Soil Transfer parameters (BEST) method to determine soil capacitive indicators and evaluate the physical quality of the upper soil layer. The field survey involved conducting time-lapse GPR surveys alongside infiltration experiments between GPR repetitions. These experiments included both tension and ponding tests, designed to sequentially activate the soil matrix and the full pore network. The results showed that the soil under study exhibited significant soil aeration and macroporosity (represented by AC and pMAC), while indicators related to microporosity (such as PAWC and RFC) were notably low. The RFC value of 0.55 m3 m−3 indicated the soil’s limited capacity to retain water relative to its total pore volume. The PAWC value of 0.10 m3 m−3 indicated a scarcity of micropores ranging from 0.2 to 30 μm in diameter, which typically hold water accessible to plant roots within the total porosity. The saturated soil hydraulic conductivity, Ks, values ranged from 192.2 to 1031.0 mm h−1, with a mean of 424.4 mm h−1, which was 7.9 times higher than the corresponding unsaturated hydraulic conductivity measured at a pressure head of h = −30 mm (K−30). The results indicated that the upper soil layer supports root proliferation and effectively drains excess water to the underlying limestone layer. However, this layer has limited capacity to store and supply water to plant roots and acts as a restrictive barrier, promoting non-uniform downward water movement, as revealed by the 3D GPR images. The observed difference in hydraulic conductivity between the two layers suggests that surface ponding and overland flow are generated through a saturation excess mechanism. Water percolating through the soil can accumulate above the limestone layer, creating a shallow perched water table. During extreme rainfall events, this water table may rise, leading to the complete saturation of the soil profile. Full article
Show Figures

Figure 1

Figure 1
<p>Flowchart outlining the process to generate a 3D image of the wetting bulb. The arrow indicates the funneling flow path through the limestone layer.</p>
Full article ">Figure 2
<p>Three-dimensional representations of the wetting zones obtained from ground-penetrating radar surveys conducted before and after wetting, during (<b>a</b>) tension and (<b>e</b>) ponding infiltrometer experiments at the Ottava site. Panels (<b>b</b>,<b>f</b>) illustrate horizontal cross-sections taken from the 3D models at a depth of −0.1m from the soil surface. Panels (<b>c</b>,<b>g</b>) present vertical cross-sections oriented north–south with a view to the east, while panels (<b>d</b>,<b>h</b>) show vertical cross-sections oriented west–east within a view to the north. The red arrows highlight the detected flow channeling through the limestone layer (see <a href="#applsci-14-09268-f001" class="html-fig">Figure 1</a> for reference).</p>
Full article ">Figure 3
<p>Example of the procedure adopted for detecting flow impedance owing to the hydraulic resistance exerted by the underlying limestone layer. (<b>a</b>): Entire cumulative infiltration curve [<span class="html-italic">I</span>(<span class="html-italic">t</span>) vs. <span class="html-italic">t</span>]. (<b>b</b>): Data linearized according to the cumulative linearization (CL, Smiles and Knight, 1976) method (<span class="html-italic">I</span>√<span class="html-italic">t</span> vs. √<span class="html-italic">t</span>). The abscissa (√<span class="html-italic">t</span>) of the intersection point of the two straight lines splits the infiltration data into two subsets. (<b>c</b>): Cumulative infiltration data representative of the first stage when water infiltrates into the upper layer.</p>
Full article ">Figure 4
<p>θ<span class="html-italic"><sub>PWP</sub></span> [m<sup>3</sup> m<sup>−3</sup>] is the permanent wilting point soil water content, corresponding to <span class="html-italic">h</span> = −150 m. θ<span class="html-italic"><sub>FC</sub></span> [m<sup>3</sup> m<sup>−3</sup>] is the field capacity (gravity drained) soil water content, corresponding to <span class="html-italic">h</span> = −1 m. θ<span class="html-italic"><sub>m</sub></span> [m<sup>3</sup> m<sup>−3</sup>] is the saturated volumetric water content of the soil matrix, corresponding to <span class="html-italic">h</span> = −0.1 m. θ<span class="html-italic"><sub>TI</sub></span> [m<sup>3</sup> m<sup>−3</sup>] is the final volumetric water content at the end of the TI test (corresponding to <span class="html-italic">h</span> = −0.03 m), θ<span class="html-italic"><sub>s</sub></span> [m<sup>3</sup> m<sup>−3</sup>] is the saturated volumetric water content. <span class="html-italic">AC</span> [m<sup>3</sup> m<sup>−3</sup>] is the air capacity. <span class="html-italic">PAWC</span> [m<sup>3</sup> m<sup>−3</sup>] is the plant-available water capacity. <span class="html-italic">RFC</span> [−] is the relative field capacity. <span class="html-italic">p<sub>MAC</sub></span> [m<sup>3</sup> m<sup>−3</sup>] is the soil macroporosity. <sup>†</sup> Water content values determined from wet soil samples collected after the tension (θ<span class="html-italic"><sub>TI</sub></span>) and Beerkan (θ<span class="html-italic"><sub>s</sub></span>) infiltration tests.</p>
Full article ">
19 pages, 9337 KiB  
Article
Investigating the Internal Deterioration of the Auriga Statue of Mozia Island, Sicily, through Ultrasonic and Ground-Penetrating Radar Studies
by Patrizia Capizzi, Raffaele Martorana and Alessandra Carollo
Sensors 2024, 24(19), 6450; https://doi.org/10.3390/s24196450 - 5 Oct 2024
Viewed by 1158
Abstract
The Greek marble statue of the Auriga of Mozia Island, in Sicily, is the most important artwork displayed at the Whitaker Foundation Archaeological Museum. It underwent geophysical investigations twice, in 2012 and 2021, to assess the marble’s degradation. The 2012 investigation prepared the [...] Read more.
The Greek marble statue of the Auriga of Mozia Island, in Sicily, is the most important artwork displayed at the Whitaker Foundation Archaeological Museum. It underwent geophysical investigations twice, in 2012 and 2021, to assess the marble’s degradation. The 2012 investigation prepared the statue for transfer to the Paul Getty Museum in New York and repositioning on an anti-seismic pedestal. The 2021 investigation evaluated potential new damage before another transfer. Both investigations utilized 3D ultrasonic tomography (UST) to detect degraded marble areas and ground-penetrating radar (GPR) to identify internal discontinuities, such as fractures or lesions, and locate metal pins that were previously inserted to reassemble the statue and its pedestal. Results from the UST indicate an average marble velocity of approximately 4700 m/s, suggesting good mechanical strength, with some areas showing lower velocities (~3000 m/s) within the material’s variability range. The GPR profiles demonstrated internal signal homogeneity, excluding internal fracture surfaces or lesions, and confirmed the presence of metallic pins. This study highlights the effectiveness of integrating UST and GPR for non-invasive diagnostics of marble sculptures, providing detailed insights into the marble’s condition and identifying hidden defects or damage. Full article
Show Figures

Figure 1

Figure 1
<p>Geographical map of Sicily (<b>a</b>), where the black rectangle indicates the coastal lagoon of Stagnone (<b>b</b>), in the center of which is the island of Mozia (<b>c</b>).</p>
Full article ">Figure 2
<p>(<b>a</b>) Frontal view of the statue of the “Auriga” (charioteer), housed at the Whitaker Museum of Mozia; (<b>b</b>) a 3D digital reconstruction of the statue’s surface, used for the correct positioning of sensors and for the graphical rendering of the tomographies.</p>
Full article ">Figure 3
<p>Various stages of UST data acquisition on the statue of the Auriga of Mozia covered with transparent film: (<b>a</b>) view of adhesive paper circles during the 2012 measurements; (<b>b</b>) US measurements in 2021; (<b>c</b>) some measurement points on the head of the statue in 2021.</p>
Full article ">Figure 4
<p>GPR data acquisition on the statue of the Auriga (charioteer) of Mozia: (<b>a</b>) measurements taken in 2012 along horizontal closed paths marked with white adhesive tape; (<b>b</b>) a moment of the acquisition in 2021 along the vertical direction.</p>
Full article ">Figure 5
<p>Measured ultrasonic travel times versus source–receiver distance. The 2012 data (in magenta), acquired only in the lower part of the statue, are compared with those acquired in 2021 on the whole surface of the statue (white diamonds) and with those, among these, acquired in the lower part of the statue (yellow diamonds). The solid line represents the linear regression of the 2012 data, while the dotted and dashed lines represent, respectively, the linear regression of the entire 2021 dataset and of only the data acquired in the lower part of the statue.</p>
Full article ">Figure 6
<p>Results of the UST performed in 2012 on the statue of the Auriga (charioteer) of Mozia: (<b>a</b>) a 3D visualization of the tomographic model and (<b>b</b>) horizontal sections of the tomographic model.</p>
Full article ">Figure 7
<p>Results of the UST performed in 2021 on the statue of the Auriga (charioteer) of Mozia. Seventeen horizontal sections, spaced ten centimeters apart, are shown.</p>
Full article ">Figure 8
<p>Three-dimensional rendering of the UST model from data acquired in 2021 on the Auriga (charioteer) of Mozia.</p>
Full article ">Figure 9
<p>GPR sections related to the data acquired in 2012, in the horizontal direction, on approximately elliptical circuits with an interdistance of 5 cm. The height of each section from the base (in centimeters) is written in each section. The top part of the sections refers to the rear side of the statue.</p>
Full article ">Figure 10
<p>GPR profile collected in 2021 on the Auriga (charioteer) of Mozia. The reflection amplitudes are represented in a color scale from blue (maximum negative amplitude) to magenta (maximum positive amplitude). The red arrows indicate the traces of each profile on the surface of the statue. The dashed red lines highlight, in each GPR section, the reflection caused by the marble/air contact, opposite to the acquisition surface.</p>
Full article ">Figure 11
<p>Joint interpretation of GPR and UST results. The red dashed arrows highlight the correspondences between US anomalies and GPR reflections. In the GPR profiles, he reflection amplitudes are represented in a color scale from blue (maximum negative amplitude) to magenta (maximum positive amplitude); the red continuous lines indicate the traces of GPR profiles on the surface of the statue; the continuous red line in P29 highlight the reflection of a metal pin.</p>
Full article ">
19 pages, 7626 KiB  
Article
Measuring Biophysical Parameters of Wheat Canopy with MHz- and GHz-Frequency Range Impulses Employing Contactless GPR
by Konstantin Muzalevskiy, Sergey Fomin, Andrey Karavayskiy, Julia Leskova, Alexey Lipshin and Vasily Romanov
Remote Sens. 2024, 16(19), 3547; https://doi.org/10.3390/rs16193547 - 24 Sep 2024
Viewed by 505
Abstract
In this paper, the advantages of the joint use of MHz- and GHz-frequency band impulses when employing contactless ground penetration radar (GPR) for the remote sensing of biomass, the height of the wheat canopy, and underlying soil moisture were experimentally investigated. A MHz-frequency [...] Read more.
In this paper, the advantages of the joint use of MHz- and GHz-frequency band impulses when employing contactless ground penetration radar (GPR) for the remote sensing of biomass, the height of the wheat canopy, and underlying soil moisture were experimentally investigated. A MHz-frequency band nanosecond impulse with a duration of 1.2 ns (average frequency of 750 MHz and spectrum bandwidth of 580 MHz, at a level of –6 dB) was emitted and received by a GPR OKO-3 equipped with an AB-900 M3 antenna unit. A GHz-frequency band sub-nanosecond impulse with a duration of 0.5 ns (average frequency of 3.2 GHz and spectral bandwidth of 1.36 GHz, at a level of −6 dB) was generated using a horn antenna and a Keysight FieldFox N9917B 18 GHz vector network analyzer. It has been shown that changes in the relative amplitudes and time delays of nanosecond impulses, reflected from a soil surface covered with wheat at a height from 0 to 87 cm and fresh above-ground biomass (AGB) from 0 to 1.5 kg/m2, do not exceed 6% and 0.09 ns, respectively. GPR nanosecond impulses reflected/scattered by the wheat canopy have not been detected. In this research, sub-nanosecond impulses reflected/scattered by the wheat canopy have been confidently identified and make it possible to measure the wheat height (fresh AGB up to 2.3 kg/m2 and height up to 104 cm) with a determination coefficient (R2) of ~0.99 and a bias of ~−7 cm, as well as fresh AGB where R2 = 0.97, with a bias = −0.09 kg/m2, and a root-mean-square error of 0.1 kg/m2. The joint use of impulses in two different MHz- and GHz-frequency bands will, in the future, make it possible to create UAV-based reflectometers for simultaneously mapping the soil moisture, height, and biomass of vegetation for precision farming systems. Full article
Show Figures

Figure 1

Figure 1
<p>Geographic location of the test field (red rectangle).</p>
Full article ">Figure 2
<p>(<b>a</b>–<b>d</b>) The process used in the experiment on the remote sensing of the wheat canopy in the MHz-frequency range with an OKO-3 GPR (23 August 2023); (<b>d</b>) measurement over a metal screen. Free space calibration is not shown.</p>
Full article ">Figure 3
<p>(<b>a</b>–<b>d</b>) The process used in the experiment on the remote sensing of the wheat canopy in the GHz-frequency range with a horn antenna (3 September 2023); (<b>d</b>) measurement over a metal screen. Free space calibration is not shown.</p>
Full article ">Figure 4
<p>Lambda dependency of the <span class="html-italic">b</span>-factors: 1—corn [<a href="#B64-remotesensing-16-03547" class="html-bibr">64</a>,<a href="#B65-remotesensing-16-03547" class="html-bibr">65</a>], 2—soybean [<a href="#B64-remotesensing-16-03547" class="html-bibr">64</a>,<a href="#B65-remotesensing-16-03547" class="html-bibr">65</a>,<a href="#B66-remotesensing-16-03547" class="html-bibr">66</a>,<a href="#B67-remotesensing-16-03547" class="html-bibr">67</a>], 3—wheat [<a href="#B64-remotesensing-16-03547" class="html-bibr">64</a>,<a href="#B65-remotesensing-16-03547" class="html-bibr">65</a>], 4—alfalfa [<a href="#B64-remotesensing-16-03547" class="html-bibr">64</a>,<a href="#B65-remotesensing-16-03547" class="html-bibr">65</a>], 5—wheat grains [<a href="#B60-remotesensing-16-03547" class="html-bibr">60</a>], 6—cereals, and sorghum [<a href="#B66-remotesensing-16-03547" class="html-bibr">66</a>] are <span class="html-italic">b<sub>rad</sub></span>-factors, estimated based on radiometric measurements, and the 7–<span class="html-italic">b<sub>refr</sub></span>-factor, calculated based on the refractive mixing model (4).</p>
Full article ">Figure 5
<p>Flowchart of algorithms used for soil geophysical and canopy biometric parameters, based on the measured OKO-3 GPR data in the time domain (TD) and VNA data in the frequency domain (FD).</p>
Full article ">Figure 6
<p>(<b>a</b>) Time shapes of the GPR impulse MHz-frequency range, reflected from the vegetation–soil cover <span class="html-italic">s<sub>sv</sub></span>(<span class="html-italic">t</span>) (color lines) and metal reflector <span class="html-italic">s<sub>ref</sub></span>(<span class="html-italic">t</span>) (gray solid line), with the normalized envelope of impulse, reflected from the metal screen |<span class="html-italic">R<sub>ref</sub></span>(<span class="html-italic">t</span>)|=<math display="inline"><semantics> <mrow> <mtext> </mtext> <mfenced open="|" close="|" separators="|"> <mrow> <msub> <mrow> <mover accent="true"> <mrow> <mi>s</mi> </mrow> <mo>˙</mo> </mover> </mrow> <mrow> <mi>r</mi> <mi>e</mi> <mi>f</mi> </mrow> </msub> <mfenced separators="|"> <mrow> <mi>t</mi> </mrow> </mfenced> </mrow> </mfenced> <mo>/</mo> <msub> <mrow> <mfenced open="|" close="|" separators="|"> <mrow> <msub> <mrow> <mover accent="true"> <mrow> <mi>s</mi> </mrow> <mo>˙</mo> </mover> </mrow> <mrow> <mi>r</mi> <mi>e</mi> <mi>f</mi> </mrow> </msub> <mfenced separators="|"> <mrow> <mi>t</mi> </mrow> </mfenced> </mrow> </mfenced> </mrow> <mrow> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> </msub> </mrow> </semantics></math> (gray dashed line); (<b>b</b>) the corresponding normalized module of the impulse spectrum |<span class="html-italic">S<sub>sv,ref</sub></span>(<span class="html-italic">f</span>)|; (<b>c</b>) the module of the reflection coefficient |<span class="html-italic">R<sub>sv</sub></span>(<span class="html-italic">f</span><sub>0</sub>)| and the delay time Δ<span class="html-italic">t</span> (calculated from the maximum envelopes) of GPR impulses (see <a href="#remotesensing-16-03547-f006" class="html-fig">Figure 6</a>a), depending on the measured fresh AGB value calculated by the thermostat-weight method.</p>
Full article ">Figure 7
<p>(<b>a</b>) The measured (solid lines) and retrieved (dash lines) modules <math display="inline"><semantics> <mrow> <mo>|</mo> <msub> <mrow> <mi>R</mi> </mrow> <mrow> <mi>s</mi> <mi>v</mi> </mrow> </msub> <mfenced separators="|"> <mrow> <mi>f</mi> <mo>,</mo> <mi>B</mi> </mrow> </mfenced> <mo>|</mo> <mtext> </mtext> </mrow> </semantics></math> and arguments <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>a</mi> <mi>r</mi> <mi>g</mi> <mtext> </mtext> <mi>R</mi> </mrow> <mrow> <mi>s</mi> <mi>v</mi> </mrow> </msub> <mfenced separators="|"> <mrow> <mi>f</mi> <mo>,</mo> <mi>B</mi> </mrow> </mfenced> </mrow> </semantics></math> of the spectrum of reflection coefficients; (<b>b</b>) the correlation between measured and retrieved fresh AGB values. The optimally found parameters while solving the inverse problem (6) were <span class="html-italic">W<sub>retr</sub> </span>= 23.4 ± 2.6% and σ<span class="html-italic"><sub>r,retr</sub></span>= 1.4 ± 0.3 cm. The color scheme of solid and dashed color lines is the same (various colors of lines correspond to the different values for fresh AGB in [kg/m<sup>2</sup>]).</p>
Full article ">Figure 8
<p>(<b>a</b>) Time shapes of the GPR impulses GHz-frequency range, reflected from the vegetation–soil cover <span class="html-italic">s<sub>sv</sub></span>(<span class="html-italic">t</span>) (color lines) and metal reflector <span class="html-italic">s<sub>ref</sub></span>(<span class="html-italic">t</span>) (gray solid line), normalized envelope of impulse, reflected from the metal screen |<span class="html-italic">R<sub>ref</sub></span>(<span class="html-italic">t</span>)| (gray dashed line); (<b>b</b>) the corresponding normalized module of the impulse spectrum |<span class="html-italic">S<sub>sv,ref</sub></span>(<span class="html-italic">f</span>)|. Various colors of lines correspond to the different values of fresh AGB in [kg/m<sup>2</sup>]. The color scheme is the same as in the pictures.</p>
Full article ">Figure 9
<p>(<b>a</b>,<b>b</b>) Normalized time shapes of impulse <span class="html-italic">s<sub>sv</sub></span>(<span class="html-italic">t</span>), reflected from the vegetation–soil cover and (<b>c</b>,<b>d</b>) their envelopes <math display="inline"><semantics> <mrow> <mfenced open="|" close="|" separators="|"> <mrow> <msub> <mrow> <mi>R</mi> </mrow> <mrow> <mi>s</mi> <mi>v</mi> </mrow> </msub> <mfenced separators="|"> <mrow> <mi>t</mi> </mrow> </mfenced> </mrow> </mfenced> <mo>=</mo> <mfenced open="|" close="|" separators="|"> <mrow> <msub> <mrow> <mover accent="true"> <mrow> <mi>s</mi> </mrow> <mo>˙</mo> </mover> </mrow> <mrow> <mi>s</mi> <mi>v</mi> </mrow> </msub> <mfenced separators="|"> <mrow> <mi>t</mi> </mrow> </mfenced> </mrow> </mfenced> <mo>/</mo> <msub> <mrow> <mfenced open="|" close="|" separators="|"> <mrow> <msub> <mrow> <mover accent="true"> <mrow> <mi>s</mi> </mrow> <mo>˙</mo> </mover> </mrow> <mrow> <mi>r</mi> <mi>e</mi> <mi>f</mi> </mrow> </msub> <mfenced separators="|"> <mrow> <mi>t</mi> </mrow> </mfenced> </mrow> </mfenced> </mrow> <mrow> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> </msub> </mrow> </semantics></math> as wheat is cut down from the top to the soil surface. In the legend, the numbers indicate fresh AGB in [kg/m<sup>2</sup>], measured in situ by the thermostat–weight method. The key for the black circles will be made clear further on.</p>
Full article ">Figure 10
<p>Retrieved local values of the volumetric contents of vegetation elements Δ<span class="html-italic">v<sub>v,retr</sub></span>(<span class="html-italic">z</span>) and the refractive index <span class="html-italic">n<sub>can,retr</sub></span>(<span class="html-italic">z</span>) in the canopy. The numbers in the legend indicate the values of fresh AGB in [kg/m<sup>2</sup>], measured by the thermostat-weight method.</p>
Full article ">Figure 11
<p>Dependence of the retrieved data on the measured canopy heights.</p>
Full article ">Figure 12
<p>(<b>a</b>) The retrieved <span class="html-italic">B<sub>retr</sub></span> (color lines) and measured <span class="html-italic">B<sub>meas</sub></span> (line with black circles) for fresh AGB, depending on the retrieved <span class="html-italic">h<sub>can,retr</sub></span> and measured <span class="html-italic">h<sub>can,meas</sub></span> canopy heights, respectively; (<b>b</b>) correlation between the retrieved and measured fresh AGB values.</p>
Full article ">Figure 13
<p>Relief of the norm of the difference between measured and retrieved fresh AGB values in total for all wheat-height cuts (see <a href="#remotesensing-16-03547-f012" class="html-fig">Figure 12</a>a), depending on the half-bandwidth Δ<span class="html-italic">f</span> and the central frequency <span class="html-italic">f</span><sub>0</sub> of the Gaussian window function (see <a href="#sec2dot2-remotesensing-16-03547" class="html-sec">Section 2.2</a> and the text for Formulas (1)–(2)).</p>
Full article ">
18 pages, 7717 KiB  
Article
Development of a Background Filtering Algorithm to Improve the Accuracy of Determining Underground Cavities Using Multi-Channel Ground-Penetrating Radar and Deep Learning
by Dae Wook Park, Han Eung Kim, Kicheol Lee and Jeongjun Park
Remote Sens. 2024, 16(18), 3454; https://doi.org/10.3390/rs16183454 - 18 Sep 2024
Viewed by 526
Abstract
In the process of using multi-channel ground-penetrating radar (GPR) for underground cavity exploration, the acquired 3D data include reflection data from underground cavities or various underground objects (structures). Reflection data from unspecified structures can interfere with the identification process of underground cavities. This [...] Read more.
In the process of using multi-channel ground-penetrating radar (GPR) for underground cavity exploration, the acquired 3D data include reflection data from underground cavities or various underground objects (structures). Reflection data from unspecified structures can interfere with the identification process of underground cavities. This study aims to identify underground cavities using a C-GAN model with an applied ResBlock technique. This deep learning model demonstrates excellent performance in the image domain and can automatically classify the presence of cavities by analyzing 3D GPR data, including reflection waveforms (A-scan), cross-sectional views (B-scan), and plan views (C-scan) measured from the ground under roads. To maximize the performance of the C-GAN model, a background filtering algorithm (BFA) was developed and applied to enhance the visibility and clarity of underground cavities. To verify the performance of the developed BFA, 3D data collected from roads in Seoul, Republic of Korea, using 3D GPR equipment were transformed, and the C-GAN model was applied. As a result, it was confirmed that the recall, an indicator of cavity prediction, improved by approximately 1.15 times compared to when the BFA was not applied. This signifies the verification of the effectiveness of the BFA. This study developed a special algorithm to distinguish underground cavities. This means that in the future, not only the advancement of separate equipment and systems but also the development of specific algorithms can contribute to the cavity exploration process. Full article
(This article belongs to the Special Issue Advanced Ground-Penetrating Radar (GPR) Technologies and Applications)
Show Figures

Figure 1

Figure 1
<p>Schematic of data acquisition.</p>
Full article ">Figure 2
<p>A-scan images: (<b>a</b>) underground cavity −1; (<b>b</b>) underground cavity −2; (<b>c</b>) loose gravel layer; (<b>d</b>) buried pipes.</p>
Full article ">Figure 3
<p>B-scan images: (<b>a</b>) longitudinal section; (<b>b</b>) cross-section.</p>
Full article ">Figure 4
<p>C-scan images with depth from the ground surface: (<b>a</b>) underground cavity at a depth of 32 cm; (<b>b</b>) underground cavity at a depth of 42 cm; (<b>c</b>) buried pipes at a depth of 26 cm; (<b>d</b>) buried pipes at a depth of 36 cm.</p>
Full article ">Figure 5
<p>Workflow of the BFA.</p>
Full article ">Figure 6
<p>Example of reconstructing data in the form of a 3-dimensional matrix.</p>
Full article ">Figure 7
<p>Process of ground alignment.</p>
Full article ">Figure 8
<p>Data amplification function: (<b>a</b>) error function; (<b>b</b>) improved error correction function.</p>
Full article ">Figure 9
<p>Data images: (<b>a</b>) non-application of the improved error correction function; (<b>b</b>) application of the improved error correction function.</p>
Full article ">Figure 10
<p>Noise removal using local average subtraction: (<b>a</b>) non-application; (<b>b</b>) application.</p>
Full article ">Figure 11
<p>Data processing images of the deep learning model: (<b>a</b>) input data image; (<b>b</b>) data processing; (<b>c</b>) output of the cavity section.</p>
Full article ">Figure 12
<p>C-GAN architecture tailored for 3D GPR data classification.</p>
Full article ">Figure 13
<p>Field validation tests using a van with a mounted 3D GPR.</p>
Full article ">Figure 14
<p>GPR radargrams: (<b>a</b>) cross-section in the longitudinal direction; (<b>b</b>) planar at a depth of 0.67 m.</p>
Full article ">Figure 15
<p>Comparison of B-Scan data with the BFA applied: (<b>a</b>) longitudinal section; (<b>b</b>) cross-section.</p>
Full article ">Figure 16
<p>Comparison of plan view (C-Scan) data with the BFA applied.</p>
Full article ">Figure 17
<p>Comparison of the plan view of cavities at different depths without the BFA: (<b>a</b>) 15 cm; (<b>b</b>) 18 cm; (<b>c</b>) 21 cm; (<b>d</b>) 24 cm; (<b>e</b>) 27 cm.</p>
Full article ">Figure 18
<p>Comparison of the plan view of cavities at different depths with the BFA: (<b>a</b>) 15 cm; (<b>b</b>) 18 cm; (<b>c</b>) 21 cm; (<b>d</b>) 24 cm; (<b>e</b>) 27 cm.</p>
Full article ">Figure 19
<p>Comparison of the longitudinal view of cavities at cChannels without the BFA: (<b>a</b>) Channel 5; (<b>b</b>) Channel 6; (<b>c</b>) Channel 7; (<b>d</b>) Channel 8; (<b>e</b>) Channel 9.</p>
Full article ">Figure 20
<p>Comparison of the longitudinal view of cavities at cChannels with the BFA: (<b>a</b>) Channel 5; (<b>b</b>) Channel 6; (<b>c</b>) Channel 7; (<b>d</b>) Channel 8; (<b>e</b>) Channel 9.</p>
Full article ">Figure 21
<p>Results of the underground object classification: (<b>a</b>) aApplication of the BFA; (<b>b</b>) nNon-application of the BFA.</p>
Full article ">Figure 22
<p>Evaluation index with and without BFA aApplication.</p>
Full article ">
21 pages, 3642 KiB  
Article
Transcriptome Profiling of Phenylalanine-Treated Human Neuronal Model: Spotlight on Neurite Impairment and Synaptic Connectivity
by Sara Stankovic, Andrijana Lazic, Marina Parezanovic, Milena Stevanovic, Sonja Pavlovic, Maja Stojiljkovic and Kristel Klaassen
Int. J. Mol. Sci. 2024, 25(18), 10019; https://doi.org/10.3390/ijms251810019 - 18 Sep 2024
Viewed by 3646
Abstract
Phenylketonuria (PKU) is the most common inherited disorder of amino acid metabolism, characterized by high levels of phenylalanine (Phe) in the blood and brain, leading to cognitive impairment without treatment. Nevertheless, Phe-mediated brain dysfunction is not fully understood. The objective of this study [...] Read more.
Phenylketonuria (PKU) is the most common inherited disorder of amino acid metabolism, characterized by high levels of phenylalanine (Phe) in the blood and brain, leading to cognitive impairment without treatment. Nevertheless, Phe-mediated brain dysfunction is not fully understood. The objective of this study was to address gene expression alterations due to excessive Phe exposure in the human neuronal model and provide molecular advances in PKU pathophysiology. Hence, we performed NT2/D1 differentiation in culture, and, for the first time, we used Phe-treated NT2-derived neurons (NT2/N) as a novel model for Phe-mediated neuronal impairment. NT2/N were treated with 1.25 mM, 2.5 mM, 5 mM, 10 mM, and 30 mM Phe and subjected to whole-mRNA short-read sequencing. Differentially expressed genes (DEGs) were analyzed and enrichment analysis was performed. Under three different Phe concentrations (2.5 mM, 5 mM, and 10 mM), DEGs pointed to the PREX1, LRP4, CDC42BPG, GPR50, PRMT8, RASGRF2, and CDH6 genes, placing them in the context of PKU for the first time. Enriched processes included dendrite and axon impairment, synaptic transmission, and membrane assembly. In contrast to these groups, the 30 mM Phe treatment group clearly represented the neurotoxicity of Phe, exhibiting enrichment in apoptotic pathways. In conclusion, we established NT2/N as a novel model for Phe-mediated neuronal dysfunction and outlined the Phe-induced gene expression changes resulting in neurite impairment and altered synaptic connectivity. Full article
Show Figures

Figure 1

Figure 1
<p>Immunostaining of NT2/N with the anti-MAP2 antibody (red). Mature neuronal phenotype of NT2-derived neurons was confirmed by observing MAP2, a neuron-specific cytoskeletal protein. Magnification used 20×, scale bar 100 μm.</p>
Full article ">Figure 2
<p>Gene co-expression between treatment groups in comparison to control cells (untreated NT2/N). Group treated with 1.25 mM Phe shows the lowest number of uniquely expressed genes, while group treated with 30 mM Phe has far more uniquely expressed genes compared to other groups (<b>A</b>). There are 5 genes co-expressed between all treatments (<b>A</b>). Groups treated with 2.5 mM, 5 mM, and 10 mM Phe show increased number of shared genes (<b>B</b>).</p>
Full article ">Figure 3
<p>Principal component analysis biplot. PCA shows clear clustering of samples treated with 30 mM Phe. Samples treated with lower concentrations are clustered together, along with untreated cells. PCA points out the difference between highest Phe concentration group and the rest of the samples, explaining around 85% of the variability. Top 5 loadings contributing to PC1 are <span class="html-italic">CCND1</span>, <span class="html-italic">CPA4</span>, <span class="html-italic">RGS5</span>, <span class="html-italic">ADAMTS1</span>, and <span class="html-italic">CCND2</span>, while, for PC2, top 5 loadings are <span class="html-italic">FOSB</span>, <span class="html-italic">EGR1</span>, <span class="html-italic">FOS</span>, <span class="html-italic">COL4A2</span>, and <span class="html-italic">NES</span> (<a href="#app1-ijms-25-10019" class="html-app">Figure S3</a>).</p>
Full article ">Figure 4
<p>Significantly up- and down-regulated DEGs observed for each treatment group in comparison to untreated NT2/N. Genes with significantly altered expression upon exposure to 2.5 mM (<b>A</b>), 5 mM (<b>B</b>), 10 mM (<b>C</b>), and 30 mM (<b>D</b>) Phe visualized on volcano plots. Horizontal dashed lines represent adjusted <span class="html-italic">p</span>-value cutoff (0.05) on a logarithmic scale for base 10, while vertical dashed lines represent absolute log<sub>2</sub>FoldChange cutoff (0.5). Dots visualized in lighter colors indicate down-regulated genes, while darker-colored dots indicate up-regulated genes. If possible, all DEGs were shown on volcano plots (<b>A</b>,<b>B</b>), and, if not, at least top 20 DEGs were represented on the plot (<b>C</b>,<b>D</b>).</p>
Full article ">Figure 5
<p>Significant GO terms obtained by GSEA of 2.5 mM, 5 mM, 10 mM, and 30 mM Phe treatments in comparison to control group. Upon GSEA for 2.5 mM treatment, GO terms that emerged include G protein-coupled receptor signaling pathway, actin cytoskeleton organization, and neuron projection morphogenesis, along with terms connected to mitochondrial respirasome and ATP synthesis (<b>A</b>). In 5 mM Phe-treated group, among enriched GO terms were actin filament polymerization, synaptic vesicle recycling, and synaptic signaling (<b>B</b>). In NT2/N treated with 10 mM Phe, overrepresented GO terms include neurotransmitter loading into synaptic vesicle, neuron projection guidance, and axonogenesis (<b>C</b>). In the group treated with highest Phe concentration (30 mM), enriched terms were cell migration, cell adhesion, and response to external biotic stimulus, but also actin cytoskeleton and neuron to neuron synapse (<b>D</b>). Dotplots represent the top 20 most overrepresented GO terms for each treatment group.</p>
Full article ">Figure 6
<p>DEGs observed upon treatment with 2.5 mM, 5 mM, and 10 mM Phe and their known roles in human neuron dysfunction. Two distinct processes emerged as important in Phe-mediated effect on neurons: neurite impairment and synaptic connectivity. Neurite impairment due to Phe exposure is observed in gene expression alterations contributing to actin cytoskeleton regulation, dendritic arborization and dendritic spines mobility, and neurite outgrowth and guidance. Phe-mediated effect on synaptic connectivity was observed as gene expression changes connected to calcium-mediated signaling, synaptic membrane and vesicle assembly, and synaptic potentiation and neurotransmitter loading. Created with <a href="http://BioRender.com" target="_blank">BioRender.com</a> (accessed on 13 August 2024).</p>
Full article ">
5 pages, 2272 KiB  
Proceeding Paper
Three-Dimensional Reconstruction of Water Leaks in Water Distribution Networks from Ground-Penetrating Radar Images by Exploring New Influencing Factors with Multi-Agent and Intelligent Data Analysis
by Samira Islam and David Ayala-Cabrera
Eng. Proc. 2024, 69(1), 121; https://doi.org/10.3390/engproc2024069121 - 10 Sep 2024
Viewed by 201
Abstract
This paper promotes water distribution networks’ (WDNs) sustainability and efficiency by integrating intelligent data analysis with ground-penetrating radar (GPR) to better interpret GPR images for detecting water leaks, favouring their asset assessment. This work uses GPR data from a laboratory setting to investigates [...] Read more.
This paper promotes water distribution networks’ (WDNs) sustainability and efficiency by integrating intelligent data analysis with ground-penetrating radar (GPR) to better interpret GPR images for detecting water leaks, favouring their asset assessment. This work uses GPR data from a laboratory setting to investigates the effects of various parameters on image interpretability across pipes. This methodology aims to advance the automation of leak and pipe identification, improving data interpretation and reducing dependency on human experts for leakage detection purposes. The findings suggest the possibility of uncovering new features enhancing GPR image interpretability, presented in 3D models. Full article
Show Figures

Figure 1

Figure 1
<p>Outlines the methodology.</p>
Full article ">Figure 2
<p>Case 1. (<b>a</b>–<b>c</b>) 3D reconstructions. (<b>d</b>–<b>f</b>) 2D cross-sections.</p>
Full article ">Figure 3
<p>Case 2. (<b>a</b>–<b>c</b>) 3D reconstructions.</p>
Full article ">
Back to TopTop