[go: up one dir, main page]

 
 
Sign in to use this feature.

Years

Between: -

Subjects

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Journals

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
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

Article Types

Countries / Regions

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Search Results (1,758)

Search Parameters:
Keywords = GPR

Order results
Result details
Results per page
Select all
Export citation of selected articles as:
36 pages, 5088 KiB  
Article
Eco-Friendly Synthesis of ZnO Nanoparticles for Quinoline Dye Photodegradation and Antibacterial Applications Using Advanced Machine Learning Models
by Hayet Chelghoum, Noureddine Nasrallah, Hichem Tahraoui, Mahmoud F. Seleiman, Mustapha Mounir Bouhenna, Hayet Belmeskine, Meriem Zamouche, Souhila Djema, Jie Zhang, Amina Mendil, Fayçal Dergal, Mohammed Kebir and Abdeltif Amrane
Catalysts 2024, 14(11), 831; https://doi.org/10.3390/catal14110831 - 19 Nov 2024
Viewed by 115
Abstract
Community drinking water sources are increasingly contaminated by various point and non-point sources, with emerging organic contaminants and microbial strains posing health risks and disrupting ecosystems. This study explores the use of zinc oxide nanoparticles (ZnO-NPs) as a non-specific agent to address groundwater [...] Read more.
Community drinking water sources are increasingly contaminated by various point and non-point sources, with emerging organic contaminants and microbial strains posing health risks and disrupting ecosystems. This study explores the use of zinc oxide nanoparticles (ZnO-NPs) as a non-specific agent to address groundwater contamination and combat microbial resistance effectively. The ZnO-NPs were synthesized via a green chemistry approach, employing a sol-gel method with lemon peel aqueous extract. The catalyst was characterized using techniques including XRD, ATR-FTIR, SEM-EDAX, UV-DRS, BET, and Raman spectroscopy. ZnO-NPs were then tested for photodegradation of quinoline yellow dye (QY) under sunlight irradiation, as well as for their antibacterial and antioxidant properties. The ZnO-NP photocatalyst showed significant photoactivity, attributed to effective separation of photogenerated charge carriers. The efficiency of sunlight dye photodegradation was influenced by catalyst dosage (0.1–0.6 mg L−1), pH (3–11), and initial QY concentration (10–50 mg L−1). The study developed a first-order kinetic model for ZnO-NPs using the Langmuir–Hinshelwood equation, yielding kinetic constants of equilibrium adsorption and photodegradation of Kc = 6.632 × 10−2 L mg−1 and kH = 7.104 × 10−2 mg L−1 min−1, respectively. The results showed that ZnO-NPs were effective against Gram-positive bacterial strains and showed moderate antioxidant activity, suggesting their potential in wastewater disinfection to achieve sustainable development goals. A potential antibacterial mechanism of ZnO-NPs involving interactions with microbial cells is proposed. Additionally, Gaussian Process Regression (GPR) combined with an improved Lévy flight distribution (FDB-LFD) algorithm was used to model QY photodegradation by ZnO-NPs. The ARD-Exponential kernel function provided high accuracy, validated through residue analysis. Finally, an innovative MATLAB-based application was developed to integrate the GPR_FDB-LFD model and FDB-LFD algorithm, streamlining optimization for precise photodegradation rate predictions. The results obtained in this study show that the GPR and FDB-LFD approaches offer efficient and cost-effective methods for predicting dye photodegradation, saving both time and resources. Full article
(This article belongs to the Special Issue Cutting-Edge Photocatalysis)
Show Figures

Figure 1

Figure 1
<p>TGA and DTA plot of the reagent’s mixture used for ZnO under N<sub>2</sub> environment.</p>
Full article ">Figure 2
<p>XRD patterns of ZnO-NPs synthesized by sol-gel mediated green chemistry.</p>
Full article ">Figure 3
<p>Room temperature Raman spectra of ZnO (bottom) powders.</p>
Full article ">Figure 4
<p>SEM analysis (<b>a</b>) and particle size distribution (<b>b</b>) of the prepared ZnO-NPs.</p>
Full article ">Figure 5
<p>EDS spectra of ZnO-NPs.</p>
Full article ">Figure 6
<p>ATR-FTIR spectra of ZnO-NPs.</p>
Full article ">Figure 7
<p>UV-visible diffuse reflectance spectra (<b>a</b>); the Tauc-plot variation of (αhν)<sup>2</sup> versus photon energy of ZnO-NPs (<b>b</b>).</p>
Full article ">Figure 8
<p>Progression of QY removal efficiency via photodegradation, as a function of different (<b>a</b>) ZnO-NP masses, (<b>b</b>) pH, (<b>c</b>) initial dye concentrations. (Experimental conditions: catalyst dose = 0.3 g L<sup>−1</sup>, pH = 9, 25 °C.)</p>
Full article ">Figure 8 Cont.
<p>Progression of QY removal efficiency via photodegradation, as a function of different (<b>a</b>) ZnO-NP masses, (<b>b</b>) pH, (<b>c</b>) initial dye concentrations. (Experimental conditions: catalyst dose = 0.3 g L<sup>−1</sup>, pH = 9, 25 °C.)</p>
Full article ">Figure 9
<p>Kinetics of photodegradation of different initial concentrations of QY dye.</p>
Full article ">Figure 10
<p>DPPH radical scavenging activity of (<b>a</b>) prepared ZnO-NPs and (<b>b</b>) ascorbic acid.</p>
Full article ">Figure 11
<p>Schematic representation of bacterial death mechanism caused by ZnO-NPs disrupting bacterial cell membrane integrity and protein damages.</p>
Full article ">Figure 12
<p>Residuals relating to the models established by the different techniques according to the estimated values: (<b>a</b>) relationship between experimental data and the predicted data of samples, and (<b>b</b>) histogram of model errors.</p>
Full article ">Figure 12 Cont.
<p>Residuals relating to the models established by the different techniques according to the estimated values: (<b>a</b>) relationship between experimental data and the predicted data of samples, and (<b>b</b>) histogram of model errors.</p>
Full article ">Figure 13
<p>Application for prediction and optimization of QY photodegradation rate using the GPR_FDB-LFD and FDB-LFD.</p>
Full article ">Figure 14
<p>Organization chart for the development and optimization of the GPR_FDB-LFD.</p>
Full article ">
18 pages, 6212 KiB  
Article
A Method to Detect Concealed Damage in Concrete Tunnels Using a Radar Feature Vector and Bayesian Analysis of Ground-Penetrating Radar Data
by Junfang Wang, Heng Chen, Jianfu Lin and Xiangxiong Li
Buildings 2024, 14(11), 3662; https://doi.org/10.3390/buildings14113662 - 18 Nov 2024
Viewed by 266
Abstract
Many machine learning (ML)-based detection methods for interpreting ground-penetrating radar (GPR) data of concrete tunnels require extensive labeled damage-state data for model training, limiting their practical use in concealed damage detection of in-service tunnels. This study presents a probabilistic, data-driven method for GPR-based [...] Read more.
Many machine learning (ML)-based detection methods for interpreting ground-penetrating radar (GPR) data of concrete tunnels require extensive labeled damage-state data for model training, limiting their practical use in concealed damage detection of in-service tunnels. This study presents a probabilistic, data-driven method for GPR-based damage detection, which exempts the requirement in the training process of supervised ML models. The approach involves extracting a radar feature vector (RFV), building a Bayesian baseline model with healthy data, and quantifying damage severity with the Bayes factor. The RFV is a complex vector obtained by radargram data fusion. Bayesian regression is applied to build a model for the relationship between real and imaginary parts of the RFV. The Bayes factor is employed for defect identification and severity assessment, by quantifying the difference between the RFV built with new observations and the baseline RFV predicted by the baseline model with new input. The probability of damage is calculated to reflect the influence of uncertainties on the detection result. The effectiveness of the proposed method is validated through simulated data with random noise and physical model tests. This method facilitates GPR-based hidden damage detection of in-service tunnels when lacking labeled damage-state data in the model training process. Full article
(This article belongs to the Special Issue Structural Health Monitoring and Vibration Control)
Show Figures

Figure 1

Figure 1
<p>A flowchart of the proposed methodology.</p>
Full article ">Figure 2
<p>Illustration of radargram simulation by DEM-FDTD.</p>
Full article ">Figure 3
<p>Simulation of a tunnel structure without voids.</p>
Full article ">Figure 4
<p>The comparison between observed and predicted historical patterns with simulated data: (<b>a</b>) pattern comparison based on the predicted real part as a function of the imaginary part; (<b>b</b>) pattern comparison based on the predicted imaginary part as a function of the real part.</p>
Full article ">Figure 5
<p>The comparison between observed and predicted new patterns in simulated scenario I: (<b>a</b>) pattern comparison based on the predicted real part as a function of the imaginary part; (<b>b</b>) pattern comparison based on the predicted imaginary part as a function of the real part.</p>
Full article ">Figure 6
<p>Simulation of the tunnel structure with a 0.2 m void.</p>
Full article ">Figure 7
<p>The comparison between observed and predicted new patterns in simulated scenario II: (<b>a</b>) pattern comparison based on the predicted real part as a function of the imaginary part; (<b>b</b>) pattern comparison based on the predicted imaginary part as a function of the real part.</p>
Full article ">Figure 8
<p>Simulation of the tunnel structure with a 0.5 m void.</p>
Full article ">Figure 9
<p>The comparison between observed and predicted new patterns in simulated scenario III: (<b>a</b>) pattern comparison based on the predicted real part as a function of the imaginary part; (<b>b</b>) pattern comparison based on the predicted imaginary part as a function of the real part.</p>
Full article ">Figure 10
<p>The relationship of the Bayes factor and cavity size.</p>
Full article ">Figure 11
<p>Physical model of a metro tunnel: (<b>a</b>) side view of the tunnel model; (<b>b</b>) top view of the tunnel model with the approximate location of a cavity.</p>
Full article ">Figure 12
<p>The comparison between observed and predicted historical patterns with test data: (<b>a</b>) pattern comparison based on the predicted real part as a function of the imaginary part; (<b>b</b>) pattern comparison based on the predicted imaginary part as a function of the real part.</p>
Full article ">Figure 13
<p>The comparison between observed and predicted new patterns in test scenario I: (<b>a</b>) pattern comparison based on the predicted real part as a function of the imaginary part; (<b>b</b>) pattern comparison based on the predicted imaginary part as a function of the real part.</p>
Full article ">Figure 14
<p>Linear regression results for the baseline models: (<b>a</b>) the model with the real part as a function of the imaginary part and (<b>b</b>) the model with the imaginary part as a function of the real part.</p>
Full article ">Figure 15
<p>The comparison between observed and predicted new patterns in test scenario II: (<b>a</b>) pattern comparison based on the predicted real part as a function of the imaginary part; (<b>b</b>) pattern comparison based on the predicted imaginary part as a function of the real part.</p>
Full article ">
24 pages, 14942 KiB  
Article
The Ground-Penetrating Radar Image Matching Method Based on Central Dense Structure Context Features
by Jie Xu, Qifeng Lai, Dongyan Wei, Xinchun Ji, Ge Shen and Hong Yuan
Remote Sens. 2024, 16(22), 4291; https://doi.org/10.3390/rs16224291 - 18 Nov 2024
Viewed by 212
Abstract
Subsurface structural distribution can be detected using Ground-Penetrating Radar (GPR). The distribution can be considered as road fingerprints for vehicle positioning. Similar to the principle of visual image matching for localization, the position coordinates of the vehicle can be calculated by matching real-time [...] Read more.
Subsurface structural distribution can be detected using Ground-Penetrating Radar (GPR). The distribution can be considered as road fingerprints for vehicle positioning. Similar to the principle of visual image matching for localization, the position coordinates of the vehicle can be calculated by matching real-time GPR images with pre-constructed reference GPR images. However, GPR images, due to their low resolution, cannot extract well-defined geometric features such as corners and lines. Thus, traditional visual image processing algorithms perform inadequately when applied to GPR image matching. To address this issue, this paper innovatively proposes a GPR image matching and localization method based on a novel feature descriptor, termed as central dense structure context (CDSC) features. The algorithm utilizes the strip-like elements in GPR images to improve the accuracy of GPR image matching. First, a CDSC feature descriptor is designed. By applying threshold segmentation and extremum point extraction to the GPR image, stratified strip-like elements and pseudo-corner points are obtained. The pseudo-corner points are treated as the centers, and the surrounding strip-like elements are described in context to form the GPR feature descriptors. Then, based on the feature description method, feature descriptors for both the real-time image and the reference image are calculated separately. By searching for the nearest matching point pairs and removing erroneous pairs, GPR image matching and localization are achieved. The proposed algorithm was evaluated on datasets collected from urban roads and railway tracks, achieving localization errors of 0.06 m (RMSE) and 1.22 m (RMSE), respectively. Compared to the traditional Speeded Up Robust Features (SURF) visual image matching algorithm, localization errors were reduced by 86.6% and 95.7% in urban road and railway track scenarios, respectively. Full article
(This article belongs to the Special Issue Advanced Ground-Penetrating Radar (GPR) Technologies and Applications)
Show Figures

Figure 1

Figure 1
<p>Vehicle-borne GPR detects underground information.</p>
Full article ">Figure 2
<p>Three data types of ground-penetrating radar. (<b>a</b>) A-scan: the single-channel data signal acquired by GPR. (<b>b</b>) B-scan: data acquired by the GPR antenna through continuous scanning in the direction of movement. (<b>c</b>) C-scan: data composed of multiple B-scans.</p>
Full article ">Figure 3
<p>Stripe and spot features in GPR images.</p>
Full article ">Figure 4
<p>Data were collected from the same road segment at different times and under different weather conditions using GPR. Subfigure (<b>a</b>,<b>b</b>,<b>e</b>,<b>f</b>) are collected in sunny days, and subfigure (<b>c</b>,<b>d</b>) are collected in rainy days.</p>
Full article ">Figure 5
<p>Algorithm flow.</p>
Full article ">Figure 6
<p>Comparison of feature point extraction results. (<b>a</b>) The SITF algorithm was used to extract feature points in GPR images. (<b>b</b>) The ORB algorithm was used to extract feature points in GPR images. (<b>c</b>) The SURF algorithm was used to extract feature points in GPR images.</p>
Full article ">Figure 7
<p>Shape context algorithm feature extraction.</p>
Full article ">Figure 8
<p>The steps of CDSC feature extraction. Descriptors are obtained by leveraging the stripes in the binarized image surrounding the pseudo-corner points. Subgraph (<b>a</b>) is the filtered image, subgraph (<b>b</b>) is the Binarized image, and subgraph (<b>c</b>) is the Center Dense Struct Context Feature.</p>
Full article ">Figure 9
<p>CDSC features were extracted from GPR images collected twice at the same location and compared. The nine feature points on each image on the left correspond to the nine features on the right.</p>
Full article ">Figure 10
<p>Image matching schematic diagram.</p>
Full article ">Figure 11
<p>Comparison of GPR images before and after preprocessing; (<b>a</b>) displays the original output image, (<b>b</b>) shows the filtered image.</p>
Full article ">Figure 12
<p>Urban test trajectory and equipment setup. (<b>a</b>) Urban test trajectory, the yellow line represents the test trajectory. (<b>b</b>) The equipment setup for the road experiment.</p>
Full article ">Figure 13
<p>Train test trajectory and equipment setup. (<b>a</b>) Train test trajectory total length is approximately 6.7 km, the yellow line represents the test trajectory. (<b>b</b>) The equipment setup for the railway experiment.</p>
Full article ">Figure 14
<p>Railway test trajectory positioning error.</p>
Full article ">Figure 15
<p>Comparison of matching results of different methods in the railway test trajectory. (<b>a</b>) Matching GPR image with strong interference; (<b>b</b>) matching GPR image without interference.</p>
Full article ">Figure 16
<p>Railway test trajectory positioning error CDF.</p>
Full article ">Figure 17
<p>Urban test trajectory positioning error.</p>
Full article ">Figure 18
<p>Comparison of matching results of different methods in the urban test trajectory. (<b>a</b>) Matching GPR image with strong interference; (<b>b</b>) matching GPR image without interference.</p>
Full article ">Figure 19
<p>Urban test trajectory positioning error CDF.</p>
Full article ">
23 pages, 19347 KiB  
Article
Georadar Survey and Simulation for Subsurface Investigation at Historical Mosque of Sorghatmesh, Cairo, Egypt
by Mohamed Elkarmoty, Hussien E. Allam, Khalid Helal, Fathy Ahmed, Stefano Bonduà and Sherif A. Mourad
Buildings 2024, 14(11), 3653; https://doi.org/10.3390/buildings14113653 - 17 Nov 2024
Viewed by 309
Abstract
Sorghatmesh mosque is a historical structure that was constructed in Cairo, Egypt, by Prince Saif El-Din Sorghatmesh in 1356. A dual-frequency ground-penetrating radar (GPR) with 250–700 MHz was used to investigate the subsurface of the Sorghatmesh mosque for restoration purposes. A total of [...] Read more.
Sorghatmesh mosque is a historical structure that was constructed in Cairo, Egypt, by Prince Saif El-Din Sorghatmesh in 1356. A dual-frequency ground-penetrating radar (GPR) with 250–700 MHz was used to investigate the subsurface of the Sorghatmesh mosque for restoration purposes. A total of 37 lines were surveyed on the ground floor of the mosque. The subsurface utilities were detected, and the status of the concrete base and the medium of the ground floor were assessed. A set of subsurface anomalies were detected and interpreted within the ground floor area of the mosque. In order to validate the interpretation, a trial pit was drilled on the ground floor, allowing for the visual inspection of the subsurface, and a Georadar numerical simulation was carried out to study the responses of the subsurface materials and conditions. For a better comprehension of the results, the ground floor area was categorized into five zones where the GPR interpretations between survey lines are almost similar. This work not only demonstrates the effectiveness of GPR as a non-invasive investigation tool but also highlights the potential of integrating advanced technologies into cultural heritage preservation by offering refined methodologies and insights for future research and restoration efforts. Full article
(This article belongs to the Special Issue Advanced Research on Cultural Heritage)
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) The location of Sorghatmesh mosque; (<b>b</b>) the outside architecture of the mosque with the presence of a unique spiral minaret architect; and (<b>c</b>) the inside structure of the mosque (open hall with a dome in the middle and 4 iwans facing each other).</p>
Full article ">Figure 2
<p>The ground floor structural plan of the mosque.</p>
Full article ">Figure 3
<p>(<b>a</b>) A diagram of the ground-penetrating radar and (<b>b</b>) GPR measurements in Sorghatmesh mosque.</p>
Full article ">Figure 4
<p>The ground floor of the Sorghatmesh mosque with the relative location of the 37 GPR lines (blue lines).</p>
Full article ">Figure 5
<p>Results of radargram 1 with interpretation.</p>
Full article ">Figure 6
<p>The location of the trial pit (red color) relative to the lines of GPR (blue color).</p>
Full article ">Figure 7
<p>(<b>a</b>) The content of the trial pit and (<b>b</b>) the nearest radargram to the trial pit (radargram 14).</p>
Full article ">Figure 8
<p>(<b>a</b>) The design of the Georadar simulation model; (<b>b</b>) the results of the simulation with interpretation; (<b>c</b>) the actual measured Georadar profile (radargram 16) with similar responses found by simulation.</p>
Full article ">Figure 8 Cont.
<p>(<b>a</b>) The design of the Georadar simulation model; (<b>b</b>) the results of the simulation with interpretation; (<b>c</b>) the actual measured Georadar profile (radargram 16) with similar responses found by simulation.</p>
Full article ">Figure 9
<p>The five main zones of GPR measurements (GPR profiles are in blue) are zone A (red polygon), zone B (orange polygon), zone C (purple polygon), zone D (dark blue polygon), and zone E (yellow polygon).</p>
Full article ">Figure 10
<p>The discontinuity closed shape form (red polygon) with respect to GPR profiles (in blue).</p>
Full article ">Figure A1
<p>(<b>a</b>) Radargram 2 shows presence of corroded rebar (at X = 7.2–8.6 m, Z = 0.35–0.55 m), water/saturated zone (at X = 7.0–8.8 m, Z = 0.2–0.8 m), reinforced concrete (at Z = 0.25–0.6 m), and rocks/concrete (at Z = 0.02–0.25 m) (700 MHz). (<b>b</b>) Radargram 3 shows presence of water/saturated zones (at X = 3.6–4.7 m, Z = 0.2–0.7 m), cable/tube (at X = 2.2 m, Z = 0.25 m), and rocks/concrete (at Z = 0.02–0.2 m) (700 MHz). (<b>c</b>) Radargram 4 shows presence of water/saturated zone (at X = 0–1.6 m, Z = 0.3–0.85 m), rocks/concrete (at Z = 0.02–0.3 m) (700 MHz), sink drain (at X = 1.9 m), and anomaly (at X = 3.2–3.8 m, Z = 7.0–7.3 m) (250 MHz). (<b>d</b>) Radargram 5 shows presence of cable/tube (at X = 1.3 m, Z = 0.4 m), rocks/concrete (at Z = 0.02–0.25 m) (700 MHz), and sink drain (at X = 1.1 m, Z = 0.2 m) (250 MHz).</p>
Full article ">Figure A2
<p>(<b>a</b>) Radargram 6 shows presence of void (at X = 3.0 m, Z = 0.3–1.3 m), anomaly (at X = 6.8–7.3 m, Z = 0.25–1.0 m), rocks/concrete (at Z = 0.02–0.25 m) (700 MHz), and anomaly (at X = 0.6–7.8 m, Z = 7.0–7.3 m) (250 MHz). (<b>b</b>) Radargram 7 shows presence of anomaly (at X = 0.3–0.7 m, Z = 0.2–0.7 m), rocks/concrete (at Z = 0.02–0.2 m) (700 MHz), and another anomaly (at X = 2.1–4.4 m, Z = 7.0–7.25 m) (250 MHz). (<b>c</b>) Radargram 8 shows presence of rocks/concrete (at Z = 0.02–0.2 m) (700 MHz) and anomaly (at X = 5.2–6.6 m, Z = 7.0–7.3 m) (250 MHz). (<b>d</b>) Radargram 9 shows presence of anomaly (at X = 6.7–7.1 m, Z = 0.1–0.6 m), rocks/concrete (at X = 0–4.2, Z = 0.02–0.2 m and at X = 5.7–7.6 m Z = 0.02–0.2 m) (700 MHz), and another anomaly (at X = 6.7–7.6 m, Z = 6.9–7.2 m) (250 MHz).</p>
Full article ">Figure A3
<p>(<b>a</b>) Radargram 10 shows presence of rocks/concrete (at Z = 0.02–0.2 m) (700 MHz) and anomaly (at X = 6.5–6.8 m, Z = 6.9–7.2 m) (250 MHz). (<b>b</b>) Radargram 11 shows presence of inclined blocks of concrete intersected with void presence between them (at X = 0.8–3.1 m, Z = 0.6–2.1 m) (700–250 MHz), rocks/concrete (at Z = 0.02–0.2 m) (700 MHz), and anomaly (at X = 5.0–7.0 m, Z = 7.0–7.2 m) (250 MHz). (<b>c</b>) Radargram 12 shows presence of part of hyperbolic reflection (at X = 1.8–2.1 m, Z = 0.55 m) and rocks/concrete (at Z = 0.02–0.2 m) (700 MHz). (<b>d</b>) Radargram 13 shows presence of water/saturated zones (at X = 1.6–3.6 m, Z = 0.2–0.7 m, and X = 5.4–7.9 m, Z = 0.2–0.7 m), reinforced concrete (at Z = 0.2–0.6 m), rocks/concrete (at Z = 0.02–0.2 m) (700 MHz), and anomalies (at X = 1.2–1.6 m, Z = 7.0–7.2 m, at X = 3.6–4.0 m, Z = 7.0–7.2 m, at X = 4.5–5.0 m, Z = 7.0–7.2 m and at X = 6.5–6.9 m, Z = 7.0–7.2 m) (250 MHz).</p>
Full article ">Figure A4
<p>(<b>a</b>) Radargram 14 shows presence of corroded rebar (at X = 0–1.8 m, Z = 0.3–0.5 m, at X = 9.5–10.1 m, Z = 0.3–0.5 m, at X = 11.2–12.0 m, Z = 0.3–0.5 m, and at X = 14–15.8 m, Z = 0.3–0.5 m); cable/tube (at X = 2.2 m, Z = 0.45 m and at X = 7.8 m, Z = 0.45 m); reinforced concrete (at Z = 0.3–0.75 m); rocks/concrete (at X = 0–1.8, Z = 0.02–0.3 m, and at X = 3.0–15.8 m, Z = 0.02–0.3 m) (700 MHz); water/saturated zones (at X = 0.2–1.6 m, Z = 0.2–1.4 m, at X = 9.6–10.4 m, Z = 0.2–1.4 m, at X = 11.5–12.4 m, Z = 0.2–1.4 m, and at X = 14.4–15.8 m, Z = 0.1–1.4 m); and anomaly (at X = 0.4–1.6 m, Z = 7.0–7.2) (250 MHz). (<b>b</b>) Radargram 15 shows presence of corroded rebar (at X = 1.3–2.4 m, Z = 0.2–0.5 m); cable/tube (at X = 0.8 m, Z = 0.75 m, at X = 2.7 m, Z = 0.3 m, and at X = 6.15 m, Z = 0.3 m); reinforced concrete (at Z = 0.2–0.6 m); rocks/concrete (at Z = 0.02–0.2 m) (700 MHz); water/saturated zone (at X = 1.1–2.2 m, Z = 0.4–1.4 m, and at X = 4.5–5.6 m, Z = 0.4–1.4 m); and anomaly (at X = 12.7–16.0 m, Z = 7.0–7.4 m) (250 MHz). (<b>c</b>) Radargram 16 shows presence of water/saturated zone (at X = 4.0–8.2 m, Z = 0.3–0.8 m); reinforced concrete (at X = 0–4.0 m, Z = 0.3–0.65 m); rocks/concrete (at Z = 0.02–0.3 m) (700 MHz); and anomaly (at Z = 7.0–7.2 m) (250 MHz). (<b>d</b>) Radargram 17 shows presence of voids (at X = 2.4 m, Z = 0.35–1.5 m and at X = 6.9 m, Z = 0.5–1.35 m); water/saturated zones (at X = 0.6–1.8 m, Z = 0.5–1.0 m and at X = 4.2–5.9 m, Z = 0.35–1.0 m); reinforced concrete (at Z = 0.2–0.55 m); rocks/concrete (at Z = 0.02–0.2 m) (700 MHz); and anomalies (at X = 0–0.6 m, Z = 7.0–7.2 m) (250 MHz).</p>
Full article ">Figure A5
<p>(<b>a</b>) Radargram 18 shows presence of water/saturated zones (at X = 2.4–4.4 m, Z = 0.3–0.8 m, and X = 7.2–10.0 m, Z = 0.3–0.65 m), reinforced concrete (at X = 0–2.4 m, Z = 0.3–0.65 m and at X = 4.4–7.2 m, Z = 0.3–0.65 m ), and rocks/concrete (at Z = 0.02–0.3 m) (700 MHz). (<b>b</b>) Radargram 19 shows presence of cable/tube (at X = 0.9 m, Z = 0.3 m), water/saturated zone (at X = 1.7–4.8 m, Z = 0.25–0.75 m), and rocks/concrete (at Z = 0.02–0.25 m) (700 MHz). (<b>c</b>) Radargram 20 shows presence of cable/tube (at X = 1.0 m, Z = 0.45 m), corroded rebar (at X = 10.0–11.4 m, Z = 0.3–0.55 m), reinforced concrete (at Z = 0.2–0.6 m), and rocks/concrete (at Z = 0.02–0.2 m) (700 MHz). (<b>d</b>) Radargram 21 shows presence of cable/tube (at X = 1.7 m, Z = 0.5 m and at X = 2.7 m, Z = 0.35 m) and rocks/concrete (at Z = 0.02–0.2 m) (700 MHz).</p>
Full article ">Figure A6
<p>(<b>a</b>) Radargram 22 shows presence of rocks/concrete (at Z = 0.02–0.25 m) (700 MHz). (<b>b</b>) Radargram 23 shows presence of metallic object (at X = 0.8 m, Z = 0.7 m); cable/tube (at X = 2.4 m, Z = 0.6 m, X = 2.8 m, Z = 0.7 m, X = 3.2 m, Z = 0.6 m, and X = 3.7 m, Z = 0.3 m); and rocks/concrete (at Z = 0.02–0.25 m) (700 MHz). (<b>c</b>) Radargram 24 shows presence of cable/tube (at X = 0.1 m, Z = 0.55 m), discontinuity (at Z = 0.4), and rocks/concrete (at Z = 0.02–0.2 m) (700 MHz). (<b>d</b>) Radargram 25 shows presence of metallic object (at X = 0.75 m, Z = 0.8 m), cable/tube (at X = 0.1 m, Z = 0.5 m), water/saturated zones (at X = 1.0–3.2 m, Z = 0.3–0.75 m), discontinuity (at Z = 0.6 m), and rocks/concrete (at Z = 0.02–0.25 m) (700 MHz).</p>
Full article ">Figure A7
<p>(<b>a</b>) Radargram 26 shows presence of water/saturated zone (at Z = 0.25–0.7 m) and rocks/concrete (at Z = 0.02–0.25 m) (700 MHz). (<b>b</b>) Radargram 27 shows presence of water/saturated zone (at Z = 0.2–1.0 m) and rocks/concrete (at X = 0–2.4 m, Z = 0.02–0.2 m and at X = 2.85–4.5 m, Z = 0.02–0.2 m) (700 MHz). (<b>c</b>) Radargram 28 shows presence of cable/tube (at X = 0.45 m, Z = 0.35 m), water/saturated zone (at X = 1.0–2.75 m, Z = 0.25–0.75 m), and rocks/concrete (at Z = 0.02–0.2 m) (700 MHz). (<b>d</b>) Radargram 29 shows presence of anomaly (at X = 2.8–7.8 m, Z = 1.2–1.5 m), discontinuity (at Z = 1.0 m), and rocks/concrete (at Z = 0.02–0.6 m) (700 MHz).</p>
Full article ">Figure A8
<p>(<b>a</b>) Radargram 30 shows presence of sink drains (at X = 0.8 m and X = 6.4 m) and rocks/concrete (at Z = 0.02–0.2 m) (700 MHz). (<b>b</b>) Radargram 31 shows presence of sink drain (at X = 0.55 m), water/saturated zone (at X = 0.8–1.6 m, Z = 0.2–0.8 m), anomaly (at X = 1.8–2.8 m, Z = 1.5–1.75 m), and rocks/concrete (at Z = 0.02–0.2 m) (700 MHz). (<b>c</b>) Radargram 32 shows presence of strong reflection from solid material (at Z = 0.3–0.9 m) and rocks/concrete (at Z = 0.08–0.3 m) (700 MHz). (<b>d</b>) Radargram 33 shows presence of strong reflection from solid material (at Z = 0.4–0.8 m) and rocks/concrete (at Z = 0.05–0.3 m) (700 MHz).</p>
Full article ">Figure A9
<p>(<b>a</b>) Radargram 34 shows presence of rocks/concrete (at Z = 0.02–0.2 m) (700 MHz). (<b>b</b>) Radargram 35 shows presence of and rocks/concrete (at Z = 0.02–0.2 m) (700 MHz). (<b>c</b>) Radargram 36 shows presence of rocks/concrete (at Z = 0.02–0.2 m) (700 MHz). (<b>d</b>) Radargram 37 shows presence of anomaly (at X = 3.0 m, Z = 0.3 m) and rocks/concrete (at Z = 0.02–0.2 m) (700 MHz).</p>
Full article ">
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 - 17 Nov 2024
Viewed by 446
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 ">
13 pages, 3882 KiB  
Article
Machine Learning-Based Software for Predicting Pseudomonas spp. Growth Dynamics in Culture Media
by Fatih Tarlak
Life 2024, 14(11), 1490; https://doi.org/10.3390/life14111490 - 15 Nov 2024
Viewed by 340
Abstract
In predictive microbiology, both primary and secondary models are widely used to estimate microbial growth, often applied through two-step or one-step modelling approaches. This study focused on developing a tool to predict the growth of Pseudomonas spp., a prominent bacterial genus in food [...] Read more.
In predictive microbiology, both primary and secondary models are widely used to estimate microbial growth, often applied through two-step or one-step modelling approaches. This study focused on developing a tool to predict the growth of Pseudomonas spp., a prominent bacterial genus in food spoilage, by applying machine learning regression models, including Support Vector Regression (SVR), Random Forest Regression (RFR) and Gaussian Process Regression (GPR). The key environmental factors—temperature, water activity, and pH—served as predictor variables to model the growth of Pseudomonas spp. in culture media. To assess model performance, these machine learning approaches were compared with traditional models, namely the Gompertz, Logistic, Baranyi, and Huang models, using statistical indicators such as the adjusted coefficient of determination (R2adj) and root mean square error (RMSE). Machine learning models provided superior accuracy over traditional approaches, with R2adj values from 0.834 to 0.959 and RMSE values between 0.005 and 0.010, showcasing their ability to handle complex growth patterns more effectively. GPR emerged as the most accurate model for both training and testing datasets. In external validation, additional statistical indices (bias factor, Bf: 0.998 to 1.047; accuracy factor, Af: 1.100 to 1.167) further supported GPR as a reliable alternative for microbial growth prediction. This machine learning-driven approach bypasses the need for the secondary modelling step required in traditional methods, highlighting its potential as a robust tool in predictive microbiology. Full article
(This article belongs to the Collection Feature Papers in Microbiology)
Show Figures

Figure 1

Figure 1
<p>A flow chart outlining the main steps followed in the present study.</p>
Full article ">Figure 2
<p>Histograms of the variables for (<b>a</b>) temperature (°C), (<b>b</b>) water activity and (<b>c</b>) pH.</p>
Full article ">Figure 3
<p>The observed and predicted <span class="html-italic">Pseudomonas</span> spp. in culture medium using traditional models: (<b>a</b>) modified Gompertz, (<b>b</b>) Logistic, (<b>c</b>) Baranyi and (<b>d</b>) Huang for training process.</p>
Full article ">Figure 4
<p>The observed and predicted <span class="html-italic">Pseudomonas</span> spp. in culture medium using machine learning models (<b>a</b>) Support Vector Regression, (<b>b</b>) Random Forest Regression and (<b>c</b>) Gaussian Process Regression for training process.</p>
Full article ">Figure 5
<p>The relative importance of predictor variables to microorganism populations in culture medium using Gaussian Process Regression.</p>
Full article ">Figure 6
<p>The observed and predicted <span class="html-italic">Pseudomonas</span> spp. in culture medium using traditional models: (<b>a</b>) modified Gompertz, (<b>b</b>) Logistic, (<b>c</b>) Baranyi and (<b>d</b>) Huang for testing process.</p>
Full article ">Figure 7
<p>The observed and predicted <span class="html-italic">Pseudomonas</span> spp. in culture medium using machine learning models: (<b>a</b>) Support Vector Regression, (<b>b</b>) Random Forest Regression and (<b>c</b>) Gaussian Process Regression for testing process.</p>
Full article ">Figure 8
<p>Illustration of developed software and its parts.</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 429
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 ">
39 pages, 1899 KiB  
Review
Cannabis, Endocannabinoids and Brain Development: From Embryogenesis to Adolescence
by Ricardo J. Rodrigues, Joana M. Marques and Attila Köfalvi
Cells 2024, 13(22), 1875; https://doi.org/10.3390/cells13221875 - 13 Nov 2024
Viewed by 923
Abstract
The endocannabinoid signalling system (ECS) plays a critical role from the very beginning of embryogenesis. Accordingly, the ECS is engaged early on in nervous system development, starting from neurulation, supported by the identification of ECS components—both receptors and enzymes controlling endocannabinoid metabolism—at these [...] Read more.
The endocannabinoid signalling system (ECS) plays a critical role from the very beginning of embryogenesis. Accordingly, the ECS is engaged early on in nervous system development, starting from neurulation, supported by the identification of ECS components—both receptors and enzymes controlling endocannabinoid metabolism—at these early stages. In particular, regarding the brain, the ECS is involved in the tightly regulated sequence of events that comprise brain development, from neurogenesis to neuronal migration, morphological guidance for neuronal connectivity, and synaptic circuitry refinement. The importance of this broad role of the ECS across various brain development processes is further underscored by the growing understanding of the consequences of cannabis exposure at different developmental stages. Despite the considerable knowledge we have on the role of the ECS in brain development, significant gaps in our understanding remain, particularly regarding the long-term impact and underlying mechanisms of cannabis exposure at different developmental stages. This review provides an overview of the current state of knowledge on the role of the ECS throughout brain development, from embryogenesis to adulthood, and discusses the impact of cannabis exposure, especially during adolescence—a critical period of circuitry maturation and refinement coinciding with an increased risk of cannabis use. Full article
Show Figures

Figure 1

Figure 1
<p>Overview of the endocannabinoid system in the brain. The endocannabinoid system (ECS) was uncovered through research investigating the molecular targets of key phytocannabinoids found in Cannabis sativa, particularly Δ<sup>9</sup>-tetrahydrocannabinol (Δ<sup>9</sup>-THC), the psychoactive component, and cannabidiol (CBD), a non-psychoactive compound. Both Δ<sup>9</sup>-THC and CBD interact with numerous targets within the brain, and here we focus on four key receptors: the cannabinoid receptors CB<sub>1</sub> and CB<sub>2</sub> (CB<sub>1</sub>R and CB<sub>2</sub>R), GPR55, and the transient receptor potential vanilloid type 1 (TRPV<sub>1</sub>) receptor. CB<sub>1</sub>R, CB<sub>2</sub>R, and GPR55 are G protein-coupled receptors (GPCRs) with seven transmembrane-spanning domains. These four receptors are expressed across various brain cell types, including astrocytes, microglia, oligodendrocytes [<a href="#B22-cells-13-01875" class="html-bibr">22</a>], glutamatergic neurons, GABAergic interneurons, and projection neurons (GABAergic, monoaminergic, and cholinergic), at highly variable densities, depending on cell types and factors like brain region, age, and neuropsychiatric conditions. For instance, CB<sub>1</sub>Rs are present at high levels in CCK<sup>+</sup> cortical and hippocampal GABAergic interneurons and at moderate levels in VGLUT1<sup>+</sup> pyramidal cells but are virtually absent in parvalbumin<sup>+</sup> interneurons and VGLUT2<sup>+</sup> pyramidal cells. Although synaptic pruning is assisted by resting microglia, which expresses low levels of CB<sub>1</sub>Rs and CB<sub>2</sub>Rs, when activated, these cells express much higher amounts of both receptors. All four receptors are typically found in the cytoplasmic membrane—primarily in nerve terminals, dendrites, and cell bodies— albeit there is substantial evidence for their intracellular localisation, too. While Δ<sup>9</sup>-THC acts as a partial agonist at these GPCRs, CBD’s pharmacological actions are more complex, often resembling negative allosteric modulation at the CB<sub>1</sub>R and the CB<sub>2</sub>R and weak partial agonism at the CB<sub>1</sub>R-CB<sub>2</sub>R heterodimer [<a href="#B10-cells-13-01875" class="html-bibr">10</a>,<a href="#B23-cells-13-01875" class="html-bibr">23</a>]. CBD is a functionally selective antagonist at the GPR55 because it inhibits agonist-induced G-protein signalling, while it does not affect β-arrestin-mediated pathways [<a href="#B24-cells-13-01875" class="html-bibr">24</a>]. Whereas CBD at pharmacologically high concentrations (≥10 µM) can activate and desensitise the ionotropic TRPV<sub>1</sub>R, at therapeutically relevant concentrations (≤1 µM), CBD only deactivates the TRPV<sub>1</sub>R [<a href="#B25-cells-13-01875" class="html-bibr">25</a>]. In addition to receptors, the ECS includes enzymes responsible for synthesising lipid ligands that activate these receptors. One of the most well-studied eCBs, anandamide (N-arachidonoyl-ethanolamine or AEA), is synthesised from N-acylphosphatidylethanolamine (NAPE) via NAPE-specific phospholipase D (PLD). Several alternative pathways also contribute to AEA production. Diacylglycerol lipase α (DAGLα) is the primary enzyme that synthesises 2-arachidonoyl-glycerol (2-AG), another major eCB. Both AEA and 2-AG activate all four receptors, though other ligands exhibit more receptor-selective actions. For example, N-arachidonoyl-dopamine (NADA), likely produced by fatty acid amide hydrolase (FAAH) in dopaminergic cells, acts as a hybrid agonist for CB<sub>1</sub>R and TRPV<sub>1</sub>R [<a href="#B26-cells-13-01875" class="html-bibr">26</a>]. Similarly, L-α-lysophosphatidyl-inositol (LPI) and its congeners resemble classical eCBs but selectively activate GPR55. The activation of these receptors can influence virtually all functions of the brain cells expressing them, but their actions are highly context-dependent. The effects depend on factors such as receptor splice variants, heteromeric interactions with other receptors (e.g., TrkB, insulin receptor, or EGF receptor), the cell’s metabolic state and age, and the ontogenetic stage of the organism. For instance, homo(di)meric CB<sub>1</sub>Rs in neurons are mostly coupled to G<sub>i/o</sub> proteins and inhibit cAMP production, but in astrocytes, CB<sub>1</sub>R activation stimulates G<sub>q</sub> and, consequently, [Ca<sup>2+</sup>]<sub>i</sub> levels. Heteromeric CB<sub>1</sub>Rs, such as the adenosine A<sub>2A</sub> receptor-CB<sub>1</sub>R heterotetramer, may also couple to G<sub>s</sub> and stimulate cAMP synthesis [<a href="#B10-cells-13-01875" class="html-bibr">10</a>]. Many receptor-mediated effects are tied to brain cell processes, such as differentiation, maturation, migration, circuit formation, and plasticity, which are key topics in this review. Finally, after eCBs activate their receptors, they are primarily metabolised intracellularly by a variety of enzymes. The key enzymes for this review are FAAH and cyclooxygenase-2 (COX-2), which degrade anandamide, and monoacylglycerol lipase (MAGL), which metabolises 2-AG. Cytochrome P450 (P450) enzymes may also contribute to eCB metabolism. LPI is broken down by various lysophospholipases (A, C, and D).</p>
Full article ">Figure 2
<p>Schematic representation of the involvement of endocannabinoid signalling system (ECS) in corticogenesis. ECS through CB<sub>1</sub>R may be involved in cortical cell proliferation [<a href="#B92-cells-13-01875" class="html-bibr">92</a>,<a href="#B94-cells-13-01875" class="html-bibr">94</a>,<a href="#B133-cells-13-01875" class="html-bibr">133</a>] and intermediate precursor cell generation [<a href="#B101-cells-13-01875" class="html-bibr">101</a>]. CB<sub>1</sub>R is expressed in Cajal-Retzius cells and may control early-born cortical projection neuron positioning [<a href="#B93-cells-13-01875" class="html-bibr">93</a>,<a href="#B100-cells-13-01875" class="html-bibr">100</a>,<a href="#B104-cells-13-01875" class="html-bibr">104</a>]. CB<sub>1</sub>R is involved in radial migration [<a href="#B92-cells-13-01875" class="html-bibr">92</a>] in the transition from the intermediate zone (IZ) towards the cortical plate (CP) [<a href="#B106-cells-13-01875" class="html-bibr">106</a>] by controlling the neuronal polarisation [<a href="#B102-cells-13-01875" class="html-bibr">102</a>] and eventually through the control of cell movement [<a href="#B102-cells-13-01875" class="html-bibr">102</a>,<a href="#B139-cells-13-01875" class="html-bibr">139</a>], controlling the distribution of neurons across the different cortical layers [<a href="#B92-cells-13-01875" class="html-bibr">92</a>,<a href="#B106-cells-13-01875" class="html-bibr">106</a>]. CB<sub>1</sub>R also regulates tangential migration of cortical interneurons [<a href="#B37-cells-13-01875" class="html-bibr">37</a>,<a href="#B104-cells-13-01875" class="html-bibr">104</a>]. ECS, via CB<sub>1</sub>R, is also engaged in the development of cortical excitatory cytoarchitecture by controlling the differentiation of cortical projection neurons of layer 5 (Ctip2<sup>+</sup>) [<a href="#B94-cells-13-01875" class="html-bibr">94</a>,<a href="#B109-cells-13-01875" class="html-bibr">109</a>,<a href="#B128-cells-13-01875" class="html-bibr">128</a>]. Next, ECS controls the guidance of corticofugal axons [<a href="#B92-cells-13-01875" class="html-bibr">92</a>,<a href="#B94-cells-13-01875" class="html-bibr">94</a>,<a href="#B103-cells-13-01875" class="html-bibr">103</a>,<a href="#B110-cells-13-01875" class="html-bibr">110</a>,<a href="#B124-cells-13-01875" class="html-bibr">124</a>] by regulating growth cone steering through autocrine signalling by 2-AG via CB<sub>1</sub>R at the growth cones, through the regulation of Robo1 receptor and the concomitant CB<sub>2</sub>R-induced release of Slit2 by oligodendrocytes [<a href="#B124-cells-13-01875" class="html-bibr">124</a>], whose differentiation was shown to entail CB<sub>1</sub>R and CB<sub>2</sub>R [<a href="#B124-cells-13-01875" class="html-bibr">124</a>,<a href="#B140-cells-13-01875" class="html-bibr">140</a>,<a href="#B141-cells-13-01875" class="html-bibr">141</a>,<a href="#B142-cells-13-01875" class="html-bibr">142</a>,<a href="#B143-cells-13-01875" class="html-bibr">143</a>]. ECS may also control axon pathfinding through the regulation of the trafficking of deleted in colorectal cancer (DCC) receptors, which tethers the action of the guidance cue netrin [<a href="#B144-cells-13-01875" class="html-bibr">144</a>]. 2-AG signalling through CB<sub>1</sub>R also controls growth cone steering of cortical interneurons via RhoA activation [<a href="#B108-cells-13-01875" class="html-bibr">108</a>]. ECS is later involved in cortical synaptic refinement by controlling synaptic weakening/pruning through CB<sub>1</sub>R-mediated long-term depression (LTD), observed in afferent inputs at layer 2/3 and layer 4-layer 2/3 synapses [<a href="#B145-cells-13-01875" class="html-bibr">145</a>,<a href="#B146-cells-13-01875" class="html-bibr">146</a>,<a href="#B147-cells-13-01875" class="html-bibr">147</a>,<a href="#B148-cells-13-01875" class="html-bibr">148</a>,<a href="#B149-cells-13-01875" class="html-bibr">149</a>,<a href="#B150-cells-13-01875" class="html-bibr">150</a>,<a href="#B151-cells-13-01875" class="html-bibr">151</a>]. VZ, Ventricular Zone; SVZ, Subventricular zone; PP, Preplate.</p>
Full article ">
19 pages, 10723 KiB  
Article
Potential Release of Phosphorus by Runoff Loss and Stabilization of Arsenic and Cadmium in Mining-Contaminated Soils with Exogenous Phosphate Fertilizers
by Meng Zhang, Chaoyang Wei, Fen Yang, Yujian Lai, Xuemei Wang, Menglu Wang, Wei Han, Xinlian Zhong, Jian Wang, Hongbing Ji and Zhiling Guo
Sustainability 2024, 16(22), 9783; https://doi.org/10.3390/su16229783 - 9 Nov 2024
Viewed by 449
Abstract
Phosphate has been proven to be effective in remediating soils contaminated with potentially toxic elements (PTEs); however, the potential release of phosphorus (P) through runoff and the impact on PTEs’ transport in this process have never been assessed. A rainfall simulation study was [...] Read more.
Phosphate has been proven to be effective in remediating soils contaminated with potentially toxic elements (PTEs); however, the potential release of phosphorus (P) through runoff and the impact on PTEs’ transport in this process have never been assessed. A rainfall simulation study was conducted to investigate P runoff loss and its impact on the stability of arsenic (As) and cadmium (Cd) after applying potassium dihydrogen phosphate (PDP), superphosphate (SSP), and ground phosphate rock (GPR) in soil trays packed with As–Cd-contaminated soil. The phosphorus loss through runoff and sedimentary phases followed the order of SSP > PDP > GPR > control. Phosphate fertilizers’ application reduced the mobility of As and Cd. In the first rainfall, the enrichment ratios (ERs) of As and Cd in the sedimentary phase after PDP, SSP, and GPR treatment were 0.12, 0.04, and 0.08 and 0.24, 0.16, and 0.07 units lower than the control, respectively. The <53 μm fraction in the sedimentary phase accounted for 53.06–75.95%, and phosphate fertilizers significantly enhanced the As and Cd stability in this fraction. The XPS analysis showed that the conversion of As(III) to As(V) and the generation of Cd–phosphate compounds were important reasons for enhancing As and Cd stability. This study demonstrated that PDP might be capable of the remediation of As–Cd contamination with the least release of P to watersheds. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Variation curves of P (<b>a</b>–<b>c</b>), As (<b>d</b>–<b>f</b>), and Cd (<b>g</b>–<b>i</b>) concentrations in the runoff phase with rainfall duration. (PDP—potassium dihydrogen phosphate, SSP—superphosphate, GPR—ground phosphate rock). The points in the figure are average values, and the error bar represents the standard deviation (SD) (<span class="html-italic">n</span> = 3).</p>
Full article ">Figure 2
<p>Variation curves of P (<b>a</b>–<b>c</b>), As (<b>d</b>–<b>f</b>), and Cd (<b>g</b>–<b>i</b>) concentrations in the sedimentary phase with rainfall time (PDP—potassium dihydrogen phosphate, SSP—superphosphate, GPR—ground phosphate rock). The points in the figure are average values, and the error bar represents the standard deviation (SD) (<span class="html-italic">n</span> = 3).</p>
Full article ">Figure 3
<p>Variation curves of As (<b>a</b>–<b>c</b>) and Cd (<b>d</b>–<b>f</b>) enrichment ratio in sedimentary phase with rainfall time. (PDP—potassium dihydrogen phosphate, SSP—superphosphate, GPR—ground phosphate rock) (The points in the figure are average values, and the error bar represents the standard deviation (SD) (<span class="html-italic">n</span> = 3)).</p>
Full article ">Figure 4
<p>P (<b>a</b>–<b>d</b>), As (<b>e</b>–<b>h</b>), and Cd (<b>i</b>–<b>l</b>) content of each particle size at different periods in the third runoff experiment (I = initial: 0–24 min; M = middle: 25–40 min; F = final: 41–60 min, PDP—potassium dihydrogen phosphate, SSP—superphosphate, GPR—ground phosphate rock). The points in the figure are average values, and the error bar represents the standard deviation (SD) (<span class="html-italic">n</span> = 3). Different letters indicate significant differences among the treatments at <span class="html-italic">p</span> &lt; 0.05 according to Duncan’s test.</p>
Full article ">Figure 5
<p>Speciation of P (<b>a</b>–<b>d</b>), As (<b>e</b>–<b>h</b>), and Cd (<b>i</b>–<b>l</b>) in different particle sizes during different rainfall periods in the third experiment (I = initial: 0–24 min; M = middle: 25–40 min; F = final: 41–60 min).</p>
Full article ">Figure 6
<p>XPS spectra of P with particle sizes of &lt;53 μm, 53–250 μm, and &gt;250 μm in the sedimentary phase at the initial stage (0–24 min) of rainfall under different phosphate fertilizers treatments: (<b>a</b>,<b>e</b>) Control, (<b>b</b>,<b>f</b>) PDP, (<b>c</b>,<b>g</b>) SSP, (<b>d</b>,<b>h</b>) GPR (PDP—potassium dihydrogen phosphate, SSP—superphosphate, GPR—ground phosphate rock).</p>
Full article ">Figure 7
<p>XPS spectra of As and Cd with particle sizes of &lt;53 μm, 53–250 μm, and &gt;250 μm in the sedimentary phase at the initial stage (0–24 min) of rainfall under different phosphate fertilizer treatments ((<b>a</b>,<b>e</b>) Control, (<b>b</b>,<b>f</b>) PDP, (<b>c</b>,<b>g</b>) SSP, (<b>d</b>,<b>h</b>) GPR) (PDP—potassium dihydrogen phosphate, SSP—superphosphate, GPR—ground phosphate rock).</p>
Full article ">
21 pages, 10448 KiB  
Article
Research on Airborne Ground-Penetrating Radar Imaging Technology in Complex Terrain
by Yuelong Chi, Su Pang, Lifeng Mao, Qiang Zhou and Yuehua Chi
Remote Sens. 2024, 16(22), 4174; https://doi.org/10.3390/rs16224174 - 8 Nov 2024
Viewed by 430
Abstract
The integration of ground-penetrating radar (GPR) with unmanned aerial vehicles (UAVs) enables efficient non-contact detection, performing exceptionally well in complex terrains and extreme environments. However, challenges in data processing and interpretation remain significant obstacles to fully utilizing this technology. To mitigate the effects [...] Read more.
The integration of ground-penetrating radar (GPR) with unmanned aerial vehicles (UAVs) enables efficient non-contact detection, performing exceptionally well in complex terrains and extreme environments. However, challenges in data processing and interpretation remain significant obstacles to fully utilizing this technology. To mitigate the effects of numerical dispersion, this paper develops a high-order finite-difference time-domain (FDTD(2,4)) three-dimensional code suitable for airborne GPR numerical simulations. The simulation results are compared with traditional FDTD methods, validating the accuracy of the proposed approach. Additionally, a Kirchhoff migration algorithm that considers the influence of the air layer is developed for airborne GPR. Different processing strategies are applied to flat and undulating terrain models, significantly improving the identification of shallowly buried targets. Particularly under undulating terrain conditions, the energy ratio method is introduced, effectively suppressing the interference of surface reflections caused by terrain variations. This innovative approach offers a new technical pathway for efficient GPR data processing in complex terrains. The study provides new insights and methods for the practical application of airborne GPR. Full article
Show Figures

Figure 1

Figure 1
<p>Schematic diagram of using a UAV-based airborne ground-penetrating radar system.</p>
Full article ">Figure 2
<p>Three-dimensional validation model for airborne ground-penetrating radar (GPR).</p>
Full article ">Figure 3
<p>Source waveform in time domain and frequency domain. (<b>a</b>) Time-domain waveform; (<b>b</b>) frequency-domain waveform.</p>
Full article ">Figure 4
<p>Comparison of FDTD(2,4) and FDTD(2,2) simulation results for the two-layer model. (<b>a</b>) Direct wave and surface reflection; (<b>b</b>) reflection from the stratigraphic interface.</p>
Full article ">Figure 5
<p>Relative error of FDTD(2,2) and FDTD(2,4).</p>
Full article ">Figure 6
<p>Three-dimensional plots of wave-field snapshots at four different moments. (<b>a</b>) Time: 10 ns; (<b>b</b>) time: 24 ns; (<b>c</b>) time: 36 ns; (<b>d</b>) time: 60 ns.</p>
Full article ">Figure 6 Cont.
<p>Three-dimensional plots of wave-field snapshots at four different moments. (<b>a</b>) Time: 10 ns; (<b>b</b>) time: 24 ns; (<b>c</b>) time: 36 ns; (<b>d</b>) time: 60 ns.</p>
Full article ">Figure 7
<p>Three-dimensional slice plots of wave-field snapshots at four different times. (<b>a</b>) Time: 10 ns; (<b>b</b>) time: 24 ns; (<b>c</b>) time: 36 ns; (<b>d</b>) time: 60 ns.</p>
Full article ">Figure 8
<p>Schematic diagram of Kirchhoff migration considering the air layer.</p>
Full article ">Figure 9
<p>Schematic diagram of the time-window rolling energy-ratio method.</p>
Full article ">Figure 10
<p>Schematic diagram of the energy-ratio method with a sliding time window.</p>
Full article ">Figure 11
<p>Schematic diagram of the 3D model with three shallowly buried anomalous bodies under flat terrain conditions.</p>
Full article ">Figure 12
<p>Three-dimensional forward modeling results and 3D slice diagrams. (<b>a</b>) Three-dimensional forward modeling results; (<b>b</b>) slice diagram in the x-direction, y = 5 m; (<b>c</b>) x–o–y plane diagram, time = 38 ns; (<b>d</b>) x–o–y plane diagram, yime = 50 ns.</p>
Full article ">Figure 13
<p>Three-dimensional forward modeling results and slice diagrams after removing surface reflections. (<b>a</b>) Three-dimensional forward modeling results; (<b>b</b>) slice diagram in the x-direction, y = 5 m, and x–o–y plane diagram, time = 38 ns; (<b>c</b>) slice diagram in the x-direction, y = 5 m, and x–o–y plane diagram, time = 50 ns.</p>
Full article ">Figure 14
<p>Migration results for the three shallowly buried anomalies in horizontal terrain. (<b>a</b>) Migration results for the original data; (<b>b</b>) migration results after removing surface reflections from the original data.</p>
Full article ">Figure 15
<p>Three-dimensional model of a random undulating surface with a shallowly buried cubic anomaly.</p>
Full article ">Figure 16
<p>Three-dimensional forward modeling results and slice diagrams of the random undulating surface model. (<b>a</b>) Original 3D forward modeling results; (<b>b</b>) 3D slice diagram.</p>
Full article ">Figure 17
<p>Three-dimensional forward modeling results and slice diagrams of the random undulating surface model after removing surface reflections. (<b>a</b>) Three-dimensional forward modeling results; (<b>b</b>) 3D slice diagram.</p>
Full article ">Figure 18
<p>Migration results of the 3D random undulating surface model. (<b>a</b>) Migration results of the original data; (<b>b</b>) migration results after removing surface reflections.</p>
Full article ">Figure 19
<p>Three-dimensional model of a random undulating surface with a shallowly buried cubic anomaly.</p>
Full article ">Figure 20
<p>Numerical and physical simulation results (y = 1.5 m). (<b>a</b>) Numerical simulation results; (<b>b</b>) physical simulation results.</p>
Full article ">Figure 21
<p>Migration results after removing surface reflections (y = 1.5 m). (<b>a</b>) Numerical simulation results; (<b>b</b>) physical simulation results.</p>
Full article ">
20 pages, 10237 KiB  
Article
A Leaf Chlorophyll Content Estimation Method for Populus deltoides (Populus deltoides Marshall) Using Ensembled Feature Selection Framework and Unmanned Aerial Vehicle Hyperspectral Data
by Zhulin Chen, Xuefeng Wang, Shijiao Qiao, Hao Liu, Mengmeng Shi, Xingjing Chen, Haiying Jiang and Huimin Zou
Forests 2024, 15(11), 1971; https://doi.org/10.3390/f15111971 - 8 Nov 2024
Viewed by 384
Abstract
Leaf chlorophyll content (LCC) is a key indicator in representing the photosynthetic capacity of Populus deltoides (Populus deltoides Marshall). Unmanned aerial vehicle (UAV) hyperspectral imagery provides an effective approach for LCC estimation, but the issue of band redundancy significantly impacts model accuracy [...] Read more.
Leaf chlorophyll content (LCC) is a key indicator in representing the photosynthetic capacity of Populus deltoides (Populus deltoides Marshall). Unmanned aerial vehicle (UAV) hyperspectral imagery provides an effective approach for LCC estimation, but the issue of band redundancy significantly impacts model accuracy and computational efficiency. Commonly used single feature selection algorithms not only fail to balance computational efficiency with optimal set search but also struggle to combine different regression algorithms under dynamic set conditions. This study proposes an ensemble feature selection framework to enhance LCC estimation accuracy using UAV hyperspectral data. Firstly, the embedded algorithm was improved by introducing the SHapley Additive exPlanations (SHAP) algorithm into the ranking system. A dynamic ranking strategy was then employed to remove bands in steps of 10, with LCC models developed at each step to identify the initial band subset based on estimation accuracy. Finally, the wrapper algorithm was applied using the initial band subset to search for the optimal band subset and develop the corresponding model. Three regression algorithms including gradient boosting regression trees (GBRT), support vector regression (SVR), and gaussian process regression (GPR) were combined with this framework for LCC estimation. The results indicated that the GBRT-Optimal model developed using 28 bands achieved the best performance with R2 of 0.848, RMSE of 1.454 μg/cm2 and MAE of 1.121 μg/cm2. Compared with a model performance that used all bands as inputs, this optimal model reduced the RMSE value by 24.37%. In addition to estimating biophysical and biochemical parameters, this method is also applicable to other hyperspectral imaging tasks. Full article
(This article belongs to the Special Issue Panoptic Segmentation of Tree Scenes from Mobile LiDAR Data)
Show Figures

Figure 1

Figure 1
<p>The location of study area.</p>
Full article ">Figure 2
<p>The DJI 350M UAV with 300TC hyperspectral camera.</p>
Full article ">Figure 3
<p><span class="html-italic">Populus deltoides</span> leaves collection, storage and LCC extraction.</p>
Full article ">Figure 4
<p>Framework of the ensembled feature selection method.</p>
Full article ">Figure 5
<p>(<b>a</b>) RMSE value of the GBRT model during the stepwise dimensionality reduction process; (<b>b</b>) RMSE value of the SVR model during the stepwise dimensionality reduction process; (<b>c</b>) RMSE value of the GPR model during the stepwise dimensionality reduction process.</p>
Full article ">Figure 5 Cont.
<p>(<b>a</b>) RMSE value of the GBRT model during the stepwise dimensionality reduction process; (<b>b</b>) RMSE value of the SVR model during the stepwise dimensionality reduction process; (<b>c</b>) RMSE value of the GPR model during the stepwise dimensionality reduction process.</p>
Full article ">Figure 6
<p>Performances of the GBRT-50V model, SVR-50V model and GRP-50V model.</p>
Full article ">Figure 7
<p>Residual distribution in different LCC range of the GBRT-50V model, SVR-50V model and GPR-50V model.</p>
Full article ">Figure 8
<p>Performances of the GBRT-Optimal model, SVR-Optimal model and GRP-Optimal model.</p>
Full article ">Figure 9
<p>Residual distribution in different LCC range of the GBRT-Optimal model, SVR-Optimal model, and GPR-Optimal model.</p>
Full article ">Figure 10
<p>(<b>a</b>) Hyperspectral bands used in the GBRT-50V and GBRT-Optimal models; (<b>b</b>) Hyperspectral bands used in the SVR-50V and SVR-Optimal models; (<b>c</b>) Hyperspectral bands used in the GPR-50V and GPR-Optimal models.</p>
Full article ">Figure 11
<p>Bands selected according to APCC value.</p>
Full article ">Figure 12
<p>Performances of the C-GBRT-49V model, C-SVR-49V model, and C-GRP-49V model.</p>
Full article ">Figure 13
<p>(<b>a</b>) RMSE value of GBRT model during the stepwise dimensionality reduction process based on GBRT ranking; (<b>b</b>) RMSE value of SVR model during the stepwise dimensionality reduction process based on GBRT ranking; (<b>c</b>) RMSE value of GPR model during the stepwise dimensionality reduction process based on GBRT ranking.</p>
Full article ">Figure 14
<p>Performances of the C-GBRT-Optimal model, C-SVR-Optimal model and C-GRP-Optimal model.</p>
Full article ">
19 pages, 9917 KiB  
Article
Full-Waveform Inversion of Two-Parameter Ground-Penetrating Radar Based on Quadratic Wasserstein Distance
by Kai Lu, Yibo Wang, Heting Han, Shichao Zhong and Yikang Zheng
Remote Sens. 2024, 16(22), 4146; https://doi.org/10.3390/rs16224146 - 7 Nov 2024
Viewed by 679
Abstract
Full-waveform inversion (FWI) is one of the most promising techniques in current ground-penetrating radar (GPR) inversion methods. The least-squares method is usually used, minimizing the mismatch between the observed signal and the simulated signal. However, the cycle-skipping problem has become an urgent focus [...] Read more.
Full-waveform inversion (FWI) is one of the most promising techniques in current ground-penetrating radar (GPR) inversion methods. The least-squares method is usually used, minimizing the mismatch between the observed signal and the simulated signal. However, the cycle-skipping problem has become an urgent focus of this method because of the nonlinearity of the inversion problem. To mitigate the issue of local minima, the optimal transport problem has been introduced into full-waveform inversion in this study. The Wasserstein distance derived from the optimal transport problem is defined as the mismatch function in the FWI objective function, replacing the L2 norm. In this study, the Wasserstein distance is computed by using entropy regularization and the Sinkhorn algorithm to reduce computational complexity and improve efficiency. Additionally, this study presents the objective function for dual-parameter full-waveform inversion of ground-penetrating radar, with the Wasserstein distance as the mismatch function. By normalizing with the Softplus function, the electromagnetic wave signals are adjusted to meet the non-negativity and mass conservation assumptions of the Wasserstein distance, and the convexity of the method has been proven. A multi-scale frequency-domain Wasserstein distance full-waveform inversion method based on the Softplus normalization approach is proposed, enabling the simultaneous inversion of relative permittivity and conductivity from ground-penetrating radar data. Numerical simulation cases demonstrate that this method has low initial model dependency and low noise sensitivity, allowing for high-precision inversion of relative permittivity and conductivity. The inversion results show that it, in particular, significantly improves the accuracy of conductivity inversion. Full article
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) The <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>L</mi> </mrow> <mrow> <mn>2</mn> </mrow> </msub> </mrow> </semantics></math> norm measures the mismatch between two signals by calculating the difference point by point. (<b>b</b>) The Wasserstein distance computes the minimum cost required to transfer the mass of one distribution to another.</p>
Full article ">Figure 2
<p>(<b>a</b>) Ricker wavelets m and n. (<b>b</b>) The objective function with the <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>L</mi> </mrow> <mrow> <mn>2</mn> </mrow> </msub> </mrow> </semantics></math> norm. (<b>c</b>) The Softplus normalized objective function with the <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>W</mi> </mrow> <mrow> <mn>2</mn> </mrow> </msub> </mrow> </semantics></math> distance.</p>
Full article ">Figure 3
<p>(<b>a</b>) The true relative permittivity model, where the red asterisks represent the transmitter antennas. (<b>b</b>) The true conductivity model, where red triangles represent the receiver antennas.</p>
Full article ">Figure 4
<p>Initial model I. Initial relative permittivity model (<b>a</b>) and initial conductivity model (<b>b</b>) with true model background medium.</p>
Full article ">Figure 5
<p>Initial model II. Initial relative permittivity model (<b>a</b>) and initial conductivity model (<b>b</b>) with true model background medium.</p>
Full article ">Figure 6
<p>The inverted images of relative permittivity and conductivity obtained from the <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>L</mi> </mrow> <mrow> <mn>2</mn> </mrow> </msub> </mrow> </semantics></math> norm method based on initial model I. (<b>a</b>–<b>d</b>) are the reconstructed relative permittivity images from frequencies in batch 1, batch 5, batch 9, and batch 14, respectively. (<b>e</b>–<b>h</b>) are the reconstructed conductivity images from frequencies for batch 1, batch 2, batch 3, and batch 4, respectively.</p>
Full article ">Figure 7
<p>The inverted images of relative permittivity and conductivity obtained from the <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>W</mi> </mrow> <mrow> <mn>2</mn> </mrow> </msub> </mrow> </semantics></math> distance method based on initial model I. (<b>a</b>–<b>d</b>) are the reconstructed relative permittivity images from frequencies in batch 1, batch 5, batch 9, and batch 14, respectively. (<b>e</b>–<b>h</b>) are the reconstructed conductivity images from frequencies for batch 1, batch 2, batch 3, and batch 4, respectively.</p>
Full article ">Figure 8
<p>The inverted images of relative permittivity and conductivity obtained from the <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>L</mi> </mrow> <mrow> <mn>2</mn> </mrow> </msub> </mrow> </semantics></math> norm method based on initial model II. (<b>a</b>–<b>d</b>) are the reconstructed relative permittivity images from frequencies in batch 1, batch 5, batch 9, and batch 14, respectively. (<b>e</b>–<b>h</b>) are the reconstructed conductivity images from frequencies for batch 1, batch 2, batch 3, and batch 4, respectively.</p>
Full article ">Figure 9
<p>The inverted images of relative permittivity and conductivity obtained from the <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>W</mi> </mrow> <mrow> <mn>2</mn> </mrow> </msub> </mrow> </semantics></math> distance method based on initial model II. (<b>a</b>–<b>d</b>) are the reconstructed relative permittivity images from frequencies in batch 1, batch 5, batch 9, and batch 14, respectively. (<b>e</b>–<b>h</b>) are the reconstructed conductivity images from frequencies for batch 1, batch 2, batch 3, and batch 4, respectively.</p>
Full article ">Figure 10
<p>(<b>a</b>) The original multi-offset forward simulation data with the transmitter antenna located at 7 m. (<b>b</b>) The image after adding random noise to the original data (<b>a</b>).</p>
Full article ">Figure 11
<p>(<b>a</b>) The relative permittivity inversion result based on the <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>L</mi> </mrow> <mrow> <mn>2</mn> </mrow> </msub> </mrow> </semantics></math> norm for noisy data. (<b>b</b>) The relative permittivity inversion result based on the <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>W</mi> </mrow> <mrow> <mn>2</mn> </mrow> </msub> </mrow> </semantics></math> distance for noisy data.</p>
Full article ">Figure 12
<p>(<b>a</b>) The conductivity inversion result based on the <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>L</mi> </mrow> <mrow> <mn>2</mn> </mrow> </msub> </mrow> </semantics></math> norm for noisy data. (<b>b</b>) The conductivity inversion result based on the <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>W</mi> </mrow> <mrow> <mn>2</mn> </mrow> </msub> </mrow> </semantics></math> distance for noisy data.</p>
Full article ">Figure 13
<p>(<b>a</b>) True relative permittivity model. (<b>b</b>) True conductivity model. In (<b>a</b>), red stars represent transmitter locations. In (<b>b</b>), red triangles represent receiver locations. (<b>c</b>) Initial relative permittivity model. (<b>d</b>) Initial conductivity model.</p>
Full article ">Figure 14
<p>The relative permittivity inversion images obtained by using the <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>L</mi> </mrow> <mrow> <mn>2</mn> </mrow> </msub> </mrow> </semantics></math> norm as the mismatch function. (<b>a</b>–<b>d</b>) represent the relative permittivity images reconstructed from frequencies in batch 1, batch 4, batch 7, and batch 11, respectively.</p>
Full article ">Figure 15
<p>The relative permittivity inversion images obtained by using the <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>W</mi> </mrow> <mrow> <mn>2</mn> </mrow> </msub> </mrow> </semantics></math> distance as the mismatch function. (<b>a</b>–<b>d</b>) represent the relative permittivity images reconstructed from frequencies in batch 1, batch 4, batch 7, and batch 11, respectively.</p>
Full article ">Figure 16
<p>The conductivity inversion images obtained by using the <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>L</mi> </mrow> <mrow> <mn>2</mn> </mrow> </msub> </mrow> </semantics></math> norm as the mismatch function. (<b>a</b>–<b>d</b>) represent the conductivity images reconstructed from frequencies in batch 1, batch 2, batch 3, and batch 5, respectively.</p>
Full article ">Figure 17
<p>The conductivity inversion images obtained by using the <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>W</mi> </mrow> <mrow> <mn>2</mn> </mrow> </msub> </mrow> </semantics></math> distance as the mismatch function. (<b>a</b>–<b>d</b>) represent the conductivity images reconstructed from frequencies in batch 1, batch 2, batch 3, and batch 5, respectively.</p>
Full article ">
20 pages, 29993 KiB  
Article
Research on the Forward Simulation and Intelligent Detection of Defects in Highways Using Ground-Penetrating Radar
by Pengxiang Li, Mingzhou Bai, Xin Li and Chenyang Liu
Appl. Sci. 2024, 14(22), 10183; https://doi.org/10.3390/app142210183 - 6 Nov 2024
Viewed by 390
Abstract
The increasing variety and frequency of subgrade defects in operational highways have led to a rise in road safety incidents. This study employed ground-penetrating radar (GPR) detection and forward simulation to analyze the characteristic patterns of common subgrade defects, such as looseness, voids, [...] Read more.
The increasing variety and frequency of subgrade defects in operational highways have led to a rise in road safety incidents. This study employed ground-penetrating radar (GPR) detection and forward simulation to analyze the characteristic patterns of common subgrade defects, such as looseness, voids, and cavities. Through the integration of instantaneous feature information from different defect patterns with complex signal techniques, the boundary judgment of structural layers and anomalies in GPR images of various subgrade defects was improved. An intelligent recognition platform was established, and a radar image dataset was created and trained to evaluate the recognition performance of the You Only Look Once (YOLO) v3 and Single-Shot Multi-Box Detector (SSD) algorithms. Evaluation metrics such as precision, recall, F1-score, average precision (AP), and mean average precision (mAP) were used to assess the detection efficiency and accuracy for subgrade defect images. The results showed that YOLO v3 achieved an average detection accuracy of 76.69%, while the SSD achieved 75.07%. This study demonstrates that the reliability of the intelligent recognition and classification of highway subgrade defects can be enhanced by using GPR for non-destructive testing. Full article
(This article belongs to the Section Civil Engineering)
Show Figures

Figure 1

Figure 1
<p>Principle of geological radar detection.</p>
Full article ">Figure 2
<p>SSD target detection algorithm network architecture.</p>
Full article ">Figure 3
<p>YOLO target detection network architecture.</p>
Full article ">Figure 4
<p>Subgrade defect model.</p>
Full article ">Figure 5
<p>Forward simulation results of non-defect highway base model.</p>
Full article ">Figure 6
<p>Loose defect instantaneous amplitude features.</p>
Full article ">Figure 7
<p>Loose defect instantaneous phase features.</p>
Full article ">Figure 8
<p>Loose defect instantaneous frequency features.</p>
Full article ">Figure 9
<p>Cavity defect instantaneous amplitude features.</p>
Full article ">Figure 10
<p>Cavity defect instantaneous phase features.</p>
Full article ">Figure 11
<p>Cavity defect instantaneous frequency features.</p>
Full article ">Figure 12
<p>Void defect instantaneous amplitude features.</p>
Full article ">Figure 13
<p>Void defect instantaneous phase features.</p>
Full article ">Figure 14
<p>Void defect instantaneous frequency features.</p>
Full article ">Figure 15
<p>Image processing of void defect.</p>
Full article ">Figure 16
<p>Partial samples of defect data.</p>
Full article ">Figure 17
<p>SSD target detection algorithm: training parameter index changes.</p>
Full article ">Figure 18
<p>SSD target detection algorithm: roadbed defect target prediction results in GPR files.</p>
Full article ">Figure 19
<p>YOLO v3 target detection algorithm: roadbed defect target prediction results in GPR files.</p>
Full article ">Figure 20
<p>AP trend charts for different roadbed defects.</p>
Full article ">Figure 21
<p>Field roadbed defect test.</p>
Full article ">Figure 22
<p>Geological radar image.</p>
Full article ">Figure 23
<p>Intelligent detection results for on-site roadbed defects.</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 561
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 509
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 ">
Back to TopTop