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

You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Sign in to use this feature.

Years

Between: -

Subjects

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Journals

Article Types

Countries / Regions

Search Results (49)

Search Parameters:
Keywords = local track geometry

Order results
Result details
Results per page
Select all
Export citation of selected articles as:
27 pages, 3382 KiB  
Article
DOT-SLAM: A Stereo Visual Simultaneous Localization and Mapping (SLAM) System with Dynamic Object Tracking Based on Graph Optimization
by Yuan Zhu, Hao An, Huaide Wang, Ruidong Xu, Zhipeng Sun and Ke Lu
Sensors 2024, 24(14), 4676; https://doi.org/10.3390/s24144676 - 18 Jul 2024
Cited by 1 | Viewed by 561
Abstract
Most visual simultaneous localization and mapping (SLAM) systems are based on the assumption of a static environment in autonomous vehicles. However, when dynamic objects, particularly vehicles, occupy a large portion of the image, the localization accuracy of the system decreases significantly. To mitigate [...] Read more.
Most visual simultaneous localization and mapping (SLAM) systems are based on the assumption of a static environment in autonomous vehicles. However, when dynamic objects, particularly vehicles, occupy a large portion of the image, the localization accuracy of the system decreases significantly. To mitigate this challenge, this paper unveils DOT-SLAM, a novel stereo visual SLAM system that integrates dynamic object tracking through graph optimization. By integrating dynamic object pose estimation into the SLAM system, the system can effectively utilize both foreground and background points for ego vehicle localization and obtain a static feature points map. To rectify the inaccuracies in depth estimation from stereo disparity directly on the foreground points of dynamic objects due to their self-similarity characteristics, a coarse-to-fine depth estimation method based on camera–road plane geometry is presented. This method uses rough depth to guide fine stereo matching, thereby obtaining the 3 dimensions (3D)spatial positions of feature points on dynamic objects. Subsequently, by establishing constraints on the dynamic object’s pose using the road plane and non-holonomic constraints (NHCs) of the vehicle, reducing the initial pose uncertainty of dynamic objects leads to more accurate dynamic object initialization. Finally, by considering foreground points, background points, the local road plane, the ego vehicle pose, and dynamic object poses as optimization nodes, through the establishment and joint optimization of a nonlinear model based on graph optimization, accurate six degrees of freedom (DoFs) pose estimations are obtained for both the ego vehicle and dynamic objects. Experimental validation on the KITTI-360 dataset demonstrates that DOT-SLAM effectively utilizes features from the background and dynamic objects in the environment, resulting in more accurate vehicle trajectory estimation and a static environment map. Results obtained from a real-world dataset test reinforce the effectiveness. Full article
Show Figures

Figure 1

Figure 1
<p>The pose representation of a dynamic object. Cubes represent the same dynamic object in different frames, solid lines are the pose transformations in the world frame, dashed lines are transformations between camera frames, and dotted lines originating from the camera optical center.</p>
Full article ">Figure 2
<p>The pipeline of the proposed system. The inputs are stereo images and instance segmentation from the left camera’s images with instances of dynamic objects and masks of road. The outputs are the vehicle poses and the global map. The system consists of three parts, namely front-end, dynamic object management, and back-end.</p>
Full article ">Figure 3
<p>The pipeline of dynamic object management module. The red points are the foreground features and the green points are the background features in the upper left image.</p>
Full article ">Figure 4
<p>The results of SGBM with different initial values of the disparity search range.</p>
Full article ">Figure 5
<p>Depth estimation using camera–road plane geometry.</p>
Full article ">Figure 6
<p>Tracking-by-detection method for object-level tracking. Dynamic objects from the <math display="inline"><semantics> <msup> <mi>i</mi> <mi>th</mi> </msup> </semantics></math> camera frame, represented by green triangles, are predicted in the <math display="inline"><semantics> <msup> <mrow> <mo>(</mo> <mi>i</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> <mi>th</mi> </msup> </semantics></math> camera frame as blue triangles, using a constant velocity model. Association gates are then established around these predicted dynamic objects, within which instances detected in the <math display="inline"><semantics> <msup> <mrow> <mo>(</mo> <mi>i</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> <mi>th</mi> </msup> </semantics></math> camera frame are potentially associated with the predicted dynamic objects.</p>
Full article ">Figure 7
<p>Feature-level tracking by optical flow. Incorrect tracking is marked in red boxes.</p>
Full article ">Figure 8
<p>The initial position and orientation of dynamic objects. The blue points are the features in dynamic objects. The cubes represent the positions of the vehicle in consecutive frames, with the vehicle moving closely along the local road plane.</p>
Full article ">Figure 9
<p>Factor graph of the nonlinear optimization in local bundle adjustment.</p>
Full article ">Figure 10
<p>ORB features were detected in frame 4218 of KITTI-360 sequence 05. Features of dynamic objects are marked in red, while other features are marked in green.</p>
Full article ">Figure 11
<p>Estimated trajectories along with the ground truth for KITTI-360 sequence 05.</p>
Full article ">Figure 12
<p>Estimated trajectories along with the ground truth for KITTI-360 sequence 10.</p>
Full article ">Figure 13
<p>Estimated positions along with ground truth in three directions for KITTI-360 sequence 10.</p>
Full article ">Figure 14
<p>The data collection vehicle and its equipped sensors. The sub-image shows coordinate systems of different sensors: the red represents the LiDAR coordinate system, the green represents the stereo camera coordinate system, and the blue represents the IMU coordinate system.</p>
Full article ">Figure 15
<p>Estimated trajectories along with the ground truth for real-world sequence 01.</p>
Full article ">Figure 16
<p>Estimated trajectories along with the ground truth for real-world sequence 00.</p>
Full article ">Figure 17
<p>Estimated positions along with ground truth in three directions for real-world sequence 00.</p>
Full article ">
20 pages, 1567 KiB  
Article
Dynamic SNR, Spectral Efficiency, and Rate Characterization in 5G/6G mmWave/sub-THz Systems with Macro- and Micro-Mobilities
by Darya Ostrikova, Elizaveta Golos, Vitalii Beschastnyi, Egor Machnev, Yuliya Gaidamaka and Konstantin Samouylov
Future Internet 2024, 16(7), 240; https://doi.org/10.3390/fi16070240 - 6 Jul 2024
Viewed by 3175
Abstract
The performance of 5G/6G cellular systems operating in millimeter wave (mmWave, 30–100 GHz) and sub-terahertz (sub-THz, 100–300 GHz) bands is conventionally assessed by utilizing the static distributions of user locations. The rationale is that the use of the beam tracking procedure allows for [...] Read more.
The performance of 5G/6G cellular systems operating in millimeter wave (mmWave, 30–100 GHz) and sub-terahertz (sub-THz, 100–300 GHz) bands is conventionally assessed by utilizing the static distributions of user locations. The rationale is that the use of the beam tracking procedure allows for keeping the beams of a base station (BS) and user equipment (UE) aligned at all times. However, by introducing 3GPP Reduced Capability (RedCap) UEs utilizing the Radio Resource Management (RRM) Relaxation procedure, this may no longer be the case, as UEs are allowed to skip synchronization signal blocks (SSB) to improve energy efficiency. Thus, to characterize the performance of such UEs, methods explicitly accounting for UE mobility are needed. In this paper, we will utilize the tools of the stochastic geometry and random walk theory to derive signal-to-noise ratio (SNR), spectral efficiency, and rate as an explicit function of time by accounting for mmWave/sub-THZ specifics, including realistic directional antenna radiation patterns and micro- and macro-mobilities causing dynamic antenna misalignment. Different from other studies in the field that consider time-averaged performance measures, these metrics are obtained as an explicit function of time. Our numerical results illustrate that the macro-mobility specifies the overall trend of the time-dependent spectral efficiency, while local dynamics at 1–3 s scales are mainly governed by micro-mobility. The difference between spectral efficiency corresponding to perfectly synchronized UE and BS antennas and time-dependent spectral efficiency in a completely desynchronized system is rather negligible for realistic cell coverages and stays within approximately 5–10% for a wide range of system parameters. These conclusions are not affected by the utilized antenna array at the BS side. However, accounting for realistic radiation patterns is critical for a time-dependent performance analysis of 5G/6G mmWave/sub-THz systems. Full article
(This article belongs to the Section Internet of Things)
Show Figures

Figure 1

Figure 1
<p>The considered deployment for a time-dependent performance model with macro- and micro-mobility impairments.</p>
Full article ">Figure 2
<p>Antenna radiation pattern example [<a href="#B46-futureinternet-16-00240" class="html-bibr">46</a>].</p>
Full article ">Figure 3
<p>Illustration of the beam misalignment due to macro-mobility.</p>
Full article ">Figure 4
<p>Side view of the considered scenario.</p>
Full article ">Figure 5
<p>Density <math display="inline"><semantics> <mrow> <mi>p</mi> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>t</mi> <mo>|</mo> <msub> <mi>x</mi> <mn>0</mn> </msub> <mo>,</mo> <msub> <mi>t</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> </semantics></math> for different <math display="inline"><semantics> <msub> <mi>x</mi> <mn>0</mn> </msub> </semantics></math>, <span class="html-italic">t</span>, <math display="inline"><semantics> <mrow> <mi>L</mi> <mo>=</mo> <mn>10</mn> </mrow> </semantics></math>. (<b>a</b>) <math display="inline"><semantics> <mrow> <msub> <mi>x</mi> <mn>0</mn> </msub> <mo>=</mo> <mn>5</mn> </mrow> </semantics></math>. (<b>b</b>) <math display="inline"><semantics> <mrow> <msub> <mi>x</mi> <mn>0</mn> </msub> <mo>=</mo> <mn>1</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 6
<p>Top view of the considered scenario.</p>
Full article ">Figure 7
<p>Spectral efficiency as a function of the user location with respect to the BS. (<b>a</b>) Within 120 s. (<b>b</b>) Within 20 s. (<b>c</b>) Within 5 s. (<b>d</b>) Within 1 s.</p>
Full article ">Figure 8
<p>Beam misalignment as a function of mobility.</p>
Full article ">Figure 9
<p>Spectral efficiency as a function of deployment dimensions (compartment size).</p>
Full article ">Figure 10
<p>Spectral efficiency as a function of micro-mobility.</p>
Full article ">Figure 11
<p>Spectral efficiency as a function of macro-mobility. (<b>a</b>) Macro-mobility. (<b>b</b>) Mean distance to the BS.</p>
Full article ">Figure 12
<p>Spectral efficiency as a function of antenna pattern at the BS. (<b>a</b>) Within 60 s. (<b>b</b>) Within 20 s.</p>
Full article ">Figure 13
<p>Time-dependent gain as a function of the antenna type.</p>
Full article ">
27 pages, 10638 KiB  
Article
Influence of the Road Model on the Optimal Maneuver of a Racing Motorcycle
by Jan Biniewicz and Mariusz Pyrz
Appl. Sci. 2024, 14(10), 4006; https://doi.org/10.3390/app14104006 - 8 May 2024
Viewed by 637
Abstract
Motorcycle motion is largely influenced by the road geometry, which alters the allowable accelerations in longitudinal and lateral directions and influences the vertical wheel loads. Recently, a method for three-dimensional road reconstruction and its incorporation into transient and quasi-steady-state (QSS) minimum lap time [...] Read more.
Motorcycle motion is largely influenced by the road geometry, which alters the allowable accelerations in longitudinal and lateral directions and influences the vertical wheel loads. Recently, a method for three-dimensional road reconstruction and its incorporation into transient and quasi-steady-state (QSS) minimum lap time simulations (MLTSs) has been proposed. The main purpose of this work is to demonstrate how significantly different results from a minimum lap time optimal control problem can be obtained when using inappropriate elevation data sources in the track reconstruction problem, and how the road model reconstructed using poor input data can lead to misleading conclusions when analyzing real vehicle and driver performances. Two road models derived from high- and low-resolution digital elevation models (DEMs) are compared and their impact on the optimal maneuver of a racing motorcycle is examined. The essentials of track identification are presented, as well as vehicle positioning on the 3D road and the generalized QSS motorcycle model. Obtained 3D and 2D road models are analyzed in detail, on a case example of the Road Atlanta racetrack, and used in minimum lap time simulations, which are validated by the experimental data recorded on the Supersport motorcycle. The comparative analysis showed that great care should be taken when selecting the elevation dataset in the track reconstruction process, and that the 1 m resolution local DEMs seem to be sufficient to obtain MLTS results close to the measured ones. The example of using the 3D free-trajectory QSS minimum lap time problem to localize the track segments where real driver actions can be improved is also presented. The differences between simulation results on different road models of the same racetrack can be large and influence the interpretation of optimal maneuver. Full article
(This article belongs to the Section Transportation and Future Mobility)
Show Figures

Figure 1

Figure 1
<p>The ribbon frame <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>O</mi> </mrow> <mrow> <mi>t</mi> <mi>n</mi> <mi>m</mi> </mrow> </msub> </mrow> </semantics></math> and the vehicle frame <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>A</mi> </mrow> <mrow> <mover accent="true"> <mrow> <mi>x</mi> </mrow> <mo stretchy="false">^</mo> </mover> <mover accent="true"> <mrow> <mi>y</mi> </mrow> <mo stretchy="false">^</mo> </mover> <mover accent="true"> <mrow> <mi>z</mi> </mrow> <mo stretchy="false">^</mo> </mover> </mrow> </msub> </mrow> </semantics></math>.</p>
Full article ">Figure 2
<p>Vehicle positioning on the three-dimensional road.</p>
Full article ">Figure 3
<p>Generalized rigid-body motorcycle model. <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>F</mi> </mrow> <mrow> <mi>x</mi> </mrow> </msub> <mo>,</mo> </mrow> </semantics></math> <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>F</mi> </mrow> <mrow> <mi>y</mi> </mrow> </msub> <mo>,</mo> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mi>N</mi> </mrow> </semantics></math> are the longitudinal, lateral, and vertical forces acting on the front (subscript <math display="inline"><semantics> <mrow> <mi>f</mi> </mrow> </semantics></math>) and the rear (subscript <math display="inline"><semantics> <mrow> <mi>r</mi> </mrow> </semantics></math>) tire. <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>F</mi> </mrow> <mrow> <mi>d</mi> </mrow> </msub> </mrow> </semantics></math> is the drag force. The center of gravity and center of pressure are indicated by the acronyms CoG and CoP, respectively.</p>
Full article ">Figure 4
<p>Vehicle performance envelope: (<b>a</b>) exemplary g-g diagram of motorcycle at selected speed <math display="inline"><semantics> <mrow> <mi>V</mi> </mrow> </semantics></math> modified by road camber <math display="inline"><semantics> <mrow> <mi>ϕ</mi> </mrow> </semantics></math>, road inclination <math display="inline"><semantics> <mrow> <mi>μ</mi> </mrow> </semantics></math> and road normal curvature variations represented by the term <math display="inline"><semantics> <mrow> <msub> <mrow> <mover accent="true"> <mrow> <mi>ω</mi> </mrow> <mo stretchy="false">^</mo> </mover> </mrow> <mrow> <mi>y</mi> </mrow> </msub> <mi>V</mi> </mrow> </semantics></math>; (<b>b</b>) hypersurfaces of the adhesion radius generated for various values of apparent gravity <math display="inline"><semantics> <mrow> <mover accent="true"> <mrow> <mi mathvariant="normal">g</mi> </mrow> <mo stretchy="false">~</mo> </mover> </mrow> </semantics></math>.</p>
Full article ">Figure 5
<p>Enlarged view of turn 12. Measured points on track edges (circles), with the reconstructed track boundaries using penalty factor <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>w</mi> </mrow> <mrow> <mi>θ</mi> </mrow> </msub> </mrow> </semantics></math> equal to <math display="inline"><semantics> <mrow> <mn>4</mn> <mo>×</mo> <msup> <mrow> <mn>10</mn> </mrow> <mrow> <mn>4</mn> </mrow> </msup> </mrow> </semantics></math> (thick line in light blue) and <math display="inline"><semantics> <mrow> <mn>4</mn> <mo>×</mo> <msup> <mrow> <mn>10</mn> </mrow> <mrow> <mn>6</mn> </mrow> </msup> </mrow> </semantics></math> (thin line in purple). The dashed lines symbolize the calculated spine curves.</p>
Full article ">Figure 6
<p>Track reconstruction results: (<b>a</b>) plain view of Road Atlanta circuit, (<b>b</b>) track curvatures.</p>
Full article ">Figure 7
<p>Track boundaries altitude: (<b>a</b>) left edge, (<b>b</b>) right edge.</p>
Full article ">Figure 8
<p>Height difference between track boundaries.</p>
Full article ">Figure 9
<p>Comparison of the Euler angles: (<b>a</b>) road camber <math display="inline"><semantics> <mrow> <mi>ϕ</mi> </mrow> </semantics></math> (<b>b</b>) road inclination <math display="inline"><semantics> <mrow> <mi>μ</mi> </mrow> </semantics></math>.</p>
Full article ">Figure 10
<p>Discrepancies in the road geometry between 3D-local and 3D-global cases: (<b>a</b>) track sector between turns 5 and 8, (<b>b</b>) turn 11 and chicane composed of corners 10a and 10b. The z-axis is exaggerated by a factor of three to emphasize the differences between the models.</p>
Full article ">Figure 11
<p>Track width in particular road models.</p>
Full article ">Figure 12
<p>Plots related to the vehicle powertrain: (<b>a</b>) engine torque and power measured at the rear wheel, (<b>b</b>) driving force for gears 2–6 and total resistance force.</p>
Full article ">Figure 13
<p>Velocity profiles (<b>top</b>) and time difference (<b>bottom</b>) as a function of elapsed distance: (<b>a</b>) comparison between experimental data and 3D-local case, (<b>b</b>) comparison between GPS speed and vehicle speed computed in 3D-global and 2D cases.</p>
Full article ">Figure 14
<p>Longitudinal acceleration (<b>top</b>) and lateral acceleration (<b>bottom</b>). Experimental data compared with (<b>a</b>) 3D-local simulation, (<b>b</b>) 3D-global and 2D cases.</p>
Full article ">Figure 15
<p>g-g diagrams. Experimentally measured accelerations compared with (<b>a</b>) 3D-local simulation, (<b>b</b>) 3D-global simulation, (<b>c</b>) 2D road case.</p>
Full article ">Figure 16
<p>Comparison of (<b>a</b>) vehicle lateral position, (<b>b</b>) radius of curvature modulus.</p>
Full article ">Figure 17
<p>Comparison of the optimal trajectories; enlarged view of (<b>a</b>) turn 11, (<b>b</b>) chicane 10a–10b.</p>
Full article ">Figure 18
<p>From top to bottom: vertical force acting on the front wheel, front suspension deflection and (in blue) wheel speed, throttle position.</p>
Full article ">Figure 19
<p>Total vertical force depending on the adopted road model.</p>
Full article ">Figure 20
<p>Enlarged view of the track section including turn 5 and the straight between turns 5 and 6 (<b>left</b>). The right side of the figure shows (from top to bottom) vehicle speed, longitudinal acceleration, lateral acceleration, radius of curvature modulus, throttle position, front brake pressure, engine speed, and gear selected.</p>
Full article ">
25 pages, 3433 KiB  
Article
Analysis of Local Track Discontinuities and Defects in Railway Switches Based on Track-Side Accelerations
by Susanne Reetz, Taoufik Najeh, Jan Lundberg and Jörn Groos
Sensors 2024, 24(2), 477; https://doi.org/10.3390/s24020477 - 12 Jan 2024
Viewed by 1199
Abstract
Switches are an essential, safety-critical part of the railway infrastructure. Compared to open tracks, their complex geometry leads to increased dynamic loading on the track superstructure from passing trains, resulting in high maintenance costs. To increase efficiency, condition monitoring methods specific to railway [...] Read more.
Switches are an essential, safety-critical part of the railway infrastructure. Compared to open tracks, their complex geometry leads to increased dynamic loading on the track superstructure from passing trains, resulting in high maintenance costs. To increase efficiency, condition monitoring methods specific to railway switches are required. A common approach to track superstructure monitoring is to measure the acceleration caused by vehicle track interaction. Local interruptions in the wheel–rail contact, caused for example by local defects or track discontinuities, appear in the data as transient impact events. In this paper, such transient events are investigated in an experimental setup of a railway switch with track-side acceleration sensors, using frequency and waveform analysis. The aim is to understand if and how the origins of these impact events can be distinguished in the data of this experiment, and what the implications for condition monitoring of local track discontinuities and defects with wayside acceleration sensors are in practice. For the same experimental configuration, individual impact events are shown to be reproducible in waveform and frequency content. Nevertheless, with this track-side sensor setup, the different types of track discontinuities and defects (squats, joints, crossing) could not be clearly distinguished using characteristic frequencies or waveforms. Other factors, such as the location of impact event origin relative to the sensor, are shown to have a much stronger influence. The experimental data suggest that filtering the data to narrow frequency bands around certain natural track frequencies could be beneficial for impact event detection in practice, but differentiating between individual impact event origins requires broadband signals. A multi-sensor setup with time-synchronized acceleration sensors distributed over the switch is recommended. Full article
(This article belongs to the Special Issue Real-Time Monitoring Technology for Built Infrastructure Systems)
Show Figures

Figure 1

Figure 1
<p>Schematic representation of the experimental setup (not to scale). The squats and joints on the straight track are labeled as A to K and X to Z respectively.</p>
Full article ">Figure 2
<p>Experimental setup with track discontinuities and defects (cf. <a href="#sensors-24-00477-f001" class="html-fig">Figure 1</a>). Point machine with bogie (<b>upper left</b>), joint X (<b>upper right</b>), squat D (<b>lower left</b>) and crossing (<b>lower right</b>).</p>
Full article ">Figure 3
<p>Pre-processed measurements obtained from stop block to stop block on the straight track.</p>
Full article ">Figure 4
<p>Power spectral density of whole acceleration measurement series, colored by average speed (calculated over the whole speed measurement series). The number of samples per window in the calculation of Welch’s method is 4096 and the windows overlap by 2048 samples. The lower two plots zoom in on the frequency range up to 1100 Hz, with densities displayed in log and linear scale, respectively.</p>
Full article ">Figure 5
<p>A single measurement on the straight track in facing direction. All known impact events according to <a href="#sensors-24-00477-f001" class="html-fig">Figure 1</a> are labeled (“impact event origin: axle”) at their start times.</p>
Full article ">Figure 6
<p>Hypothetical impact event signature of an 8-axle passenger train (Bombardier Regina X52) over the experimental switch with all its track discontinuities and defects at 100 km/h. Squats and joints are assumed to have a signal duration of 0.04 s (value adapted from [<a href="#B8-sensors-24-00477" class="html-bibr">8</a>]), and the crossing is 0.08 s.</p>
Full article ">Figure 7
<p>Mean correlation coefficient between impact event time series of different directions and axle for each impact event origin. The data are band-pass filtered to 10–1000 Hz and the maximal allowed timelag is set to 0.01 s, to compensate for the manually labeled start times.</p>
Full article ">Figure 8
<p>Measurement data of squat G, colored by direction and axle. The data are band-pass filtered to 10–1000 Hz after standard pre-processing and aligned using the timelag derived from the correlation coefficient (maximal allowed timelag set to 0.01 s).</p>
Full article ">Figure 9
<p>Energy spectral densities of squats in the area of the intermediate rails on a logarithmic scale, colored by axle. Axle A2 is in between axle A1 and the sensor.</p>
Full article ">Figure 10
<p>Squats observed in facing measurements, with axle A2. Each column corresponds to one measurement (i.e., one bogie passage over the switch), each row to a specific squat. The data are filtered to 1–8000 Hz during pre-processing.</p>
Full article ">Figure 11
<p>Energy spectral densities of squats observed in facing measurements, with axle A2. Each column corresponds to one measurement (i.e., one bogie passage over the switch), each row to a specific squat, as in <a href="#sensors-24-00477-f010" class="html-fig">Figure 10</a>.</p>
Full article ">Figure 12
<p>Joint Y and Z (blue) and adjacent squats (gray), observed in facing measurements, caused by axle A2. Each column corresponds to one measurement (i.e., one bogie passage over the switch), each row to a specific joint or squat, as in <a href="#sensors-24-00477-f010" class="html-fig">Figure 10</a>. The data are filtered to 1–8000 Hz during pre-processing.</p>
Full article ">Figure 13
<p>Energy spectral densities of joint Y and Z (blue) and adjacent squats (gray), observed in facing measurements, caused by axle A2. Each column corresponds to one measurement (i.e., one bogie passage over the switch), each row to a specific joint or squat, compared to <a href="#sensors-24-00477-f011" class="html-fig">Figure 11</a>.</p>
Full article ">Figure 14
<p>Impact events caused by the crossing. The data are filtered to 1–8000 Hz during pre-processing.</p>
Full article ">Figure 15
<p>Energy spectral densities of impact events caused by the crossing, compared to <a href="#sensors-24-00477-f014" class="html-fig">Figure 14</a>.</p>
Full article ">Figure 16
<p>Mean correlation coefficient between impact event time series of the track discontinuities and defects in the experiment. For the top plots, the data are band-pass filtered at 10–1000 Hz after standard pre-processing, for the bottom plots at 200–400 Hz. The plots on the left contain all combinations of all driving directions and axle, the plots on the right only in the facing direction and axle A2. For all plots, the maximum allowed timelag is set to 0.01 s, to compensate for the manually labeled start times.</p>
Full article ">
26 pages, 12200 KiB  
Article
TSG-SLAM: SLAM Employing Tight Coupling of Instance Segmentation and Geometric Constraints in Complex Dynamic Environments
by Yongchao Zhang, Yuanming Li and Pengzhan Chen
Sensors 2023, 23(24), 9807; https://doi.org/10.3390/s23249807 - 13 Dec 2023
Cited by 1 | Viewed by 1089
Abstract
Although numerous effective Simultaneous Localization and Mapping (SLAM) systems have been developed, complex dynamic environments continue to present challenges, such as managing moving objects and enabling robots to comprehend environments. This paper focuses on a visual SLAM method specifically designed for complex dynamic [...] Read more.
Although numerous effective Simultaneous Localization and Mapping (SLAM) systems have been developed, complex dynamic environments continue to present challenges, such as managing moving objects and enabling robots to comprehend environments. This paper focuses on a visual SLAM method specifically designed for complex dynamic environments. Our approach proposes a dynamic feature removal module based on the tight coupling of instance segmentation and multi-view geometric constraints (TSG). This method seamlessly integrates semantic information with geometric constraint data, using the fundamental matrix as a connecting element. In particular, instance segmentation is performed on frames to eliminate all dynamic and potentially dynamic features, retaining only reliable static features for sequential feature matching and acquiring a dependable fundamental matrix. Subsequently, based on this matrix, true dynamic features are identified and removed by capitalizing on multi-view geometry constraints while preserving reliable static features for further tracking and mapping. An instance-level semantic map of the global scenario is constructed to enhance the perception and understanding of complex dynamic environments. The proposed method is assessed on TUM datasets and in real-world scenarios, demonstrating that TSG-SLAM exhibits superior performance in detecting and eliminating dynamic feature points and obtains good localization accuracy in dynamic environments. Full article
(This article belongs to the Special Issue Advanced Sensing and Control Technologies for Autonomous Robots)
Show Figures

Figure 1

Figure 1
<p>Overall SLAM system framework.</p>
Full article ">Figure 2
<p>Framework of dynamic feature removal algorithm.</p>
Full article ">Figure 3
<p>Framework of SOLOv2 network.</p>
Full article ">Figure 4
<p>Multi-view epipolar geometry constraints.</p>
Full article ">Figure 5
<p>Framework of scenario semantic map construction algorithm.</p>
Full article ">Figure 6
<p>Comparison of dynamic feature removal method. (<b>a</b>) Original grayscale image. (<b>b</b>) ORB feature extraction. (<b>c</b>) Dynamic feature removal based on instance segmentation. (<b>d</b>) Dynamic feature removal based on tightly coupled method.</p>
Full article ">Figure 6 Cont.
<p>Comparison of dynamic feature removal method. (<b>a</b>) Original grayscale image. (<b>b</b>) ORB feature extraction. (<b>c</b>) Dynamic feature removal based on instance segmentation. (<b>d</b>) Dynamic feature removal based on tightly coupled method.</p>
Full article ">Figure 7
<p>Selected image frame sequence diagram.</p>
Full article ">Figure 8
<p>Comparison of point cloud map before and after filtering.</p>
Full article ">Figure 9
<p>Octree map before integrating target semantic color information.</p>
Full article ">Figure 10
<p>Octree map after integrating target semantic color information.</p>
Full article ">Figure 11
<p>Experimental platform.</p>
Full article ">Figure 12
<p>Estimated trajectory vs. true trajectory for fr1/desk.</p>
Full article ">Figure 13
<p>Estimated trajectory vs. true trajectory for fr1/room.</p>
Full article ">Figure 14
<p>Estimated trajectory vs. true trajectory for fr3/sitting_static.</p>
Full article ">Figure 15
<p>Estimated trajectory vs. true trajectory for fr3/sitting_xyz.</p>
Full article ">Figure 16
<p>Estimated trajectory vs. true trajectory for fr3/walking_static.</p>
Full article ">Figure 17
<p>Estimated trajectory vs. true trajectory for fr3/walking_xyz.</p>
Full article ">Figure 18
<p>Octree semantic map construction results.</p>
Full article ">Figure 19
<p>Partial image sequence of real scenario.</p>
Full article ">Figure 20
<p>Comparison of estimated trajectory for real static scenarios.</p>
Full article ">Figure 21
<p>Comparison of estimated trajectory for real dynamic scenarios.</p>
Full article ">Figure 22
<p>Semantic map of real scenario.</p>
Full article ">
19 pages, 10784 KiB  
Article
Numerical and Analytical Determination of Rockburst Characteristics: Case Study from Polish Deep Copper Mine
by Witold Pytel, Krzysztof Fuławka, Bogumiła Pałac-Walko and Piotr Mertuszka
Appl. Sci. 2023, 13(21), 11881; https://doi.org/10.3390/app132111881 - 30 Oct 2023
Cited by 3 | Viewed by 843
Abstract
A simplified analytical method useful for ductile ground support design in underground mine workings is presented. This approach allows for maintaining the stability of sidewalls in rectangular openings extracted in competent and homogeneous rocks, especially in high-pressure conditions, favoring rockburst event occurrence. The [...] Read more.
A simplified analytical method useful for ductile ground support design in underground mine workings is presented. This approach allows for maintaining the stability of sidewalls in rectangular openings extracted in competent and homogeneous rocks, especially in high-pressure conditions, favoring rockburst event occurrence. The proposed design procedure involves the typical assumptions governing the limit equilibrium method (LEM) with respect to a triangular rock block expelled from a sidewall of a long mine excavation subjected to normal stresses of the values determined based on the Maugis’s analytical solution concerned with stress distribution around the elliptical opening extracted within the homogeneous infinite elastic space. This stage of the local assessment of rock susceptibility to ejection from the walls of the excavation allowed for determining the geometry of the block whose ejection is most likely in a given geological and mining situation. Having extensive information about the geometry of the excavations and the properties of the surrounding rocks, it was possible to make an exemplary map of the risk from rockburst hazard, developed as the 2D contours of safety indexes’ values, for special-purpose excavations such as heavy machinery chambers, main excavations, etc. in conditions of selected mining panel of the deep copper mine at Legnica-Głogów Copper Basin, Poland. Another important element of the obtained results is the calculated values of the horizontal forces potentially pushing out the predetermined rock blocks. These forces are the surplus over the potential of frictional resistance and cohesion on the surfaces of previously identified discontinuities or on new cracks appearing as a result of overloading of the sidewalls. Finally, the presented algorithm allows us to perform quantitative tracking of rockburst phenomena as a function of time by determination of acceleration, velocity, and displacement of expelled rocks. Such information may be useful at the stage of designing the support for underground workings. Full article
Show Figures

Figure 1

Figure 1
<p>Scheme of utilized room and pillar-mining system.</p>
Full article ">Figure 2
<p>Example of room-and-pillar mining with roof deflection in Polish copper mines—the site where the proposed approach was validated.</p>
Full article ">Figure 3
<p>General flowchart diagram explaining the proposed methodology.</p>
Full article ">Figure 4
<p>Physical model of the sidewall loading for the case of a deep opening excavated within a homogeneous or joined rock mass.</p>
Full article ">Figure 5
<p>Views of damaged mine workings with inclined sidewall surfaces: (<b>left</b>)—long pillar at an intersection in Rudna mine, (<b>right</b>)—sidewall of a rectangular gallery with inclusion of competent rocks—Polkowice-Sieroszowice mine after blasting works in the close neighborhood, February 2016.</p>
Full article ">Figure 6
<p>Rock mass ejection from the excavation sidewall.</p>
Full article ">Figure 7
<p>Values of unbalanced forces <span class="html-italic">P<sub>k</sub></span> which should be compensated by the appropriate ground support system (<b>top</b>), angle α<sub>2</sub> (<b>middle</b>), and angle α<sub>1</sub> (<b>bottom</b>), assessed for a specific case of excavation geometry, load, and strength characteristics of surrounding rocks.</p>
Full article ">Figure 8
<p>Modulus of deformation of surrounding rock mass, scaled down using the Hoek–Diederichs approach [<a href="#B41-applsci-13-11881" class="html-bibr">41</a>].</p>
Full article ">Figure 9
<p>Maximum shear strain distribution obtained numerically (model 1, top; model 4, bottom).</p>
Full article ">Figure 10
<p>Quarter of the entire intersection; contour of the values of safety margin within sidewalls of the gallery (<b>left</b>); contours of slip planes of the failure hexahedron at the underground galleries’ intersection (<b>right</b>).</p>
Full article ">Figure 11
<p>Predicted changes in the acceleration (<b>left</b>) and velocity (<b>right</b>) of the detached rock wedge movement (t<sub>0</sub> = 0.0131 s).</p>
Full article ">Figure 12
<p>Critical depth for different types of rocks within a drift’s sidewall in selected mining-geologic conditions (<span class="html-italic">B</span> = 6 m, <span class="html-italic">H</span> = 3 m, <span class="html-italic">p<sub>x</sub></span> = 10 MPa).</p>
Full article ">Figure 13
<p>Rockburst safety factor distribution over one of the analyzed regions in the Rudna mine.</p>
Full article ">
13 pages, 2892 KiB  
Article
Kinematic Precise Point Positioning Performance-Based Cost-Effective Robot Localization System
by Ashraf Farah and Mehdi Tlija
Appl. Sci. 2023, 13(18), 10408; https://doi.org/10.3390/app131810408 - 18 Sep 2023
Cited by 1 | Viewed by 1472
Abstract
The use of high-precision positioning systems in modern navigation applications is crucial since location data is one of the most important pieces of information in Industry 4.0, especially for robots operating outdoors. In the modernization process of global navigation satellite system (GNSS) positioning, [...] Read more.
The use of high-precision positioning systems in modern navigation applications is crucial since location data is one of the most important pieces of information in Industry 4.0, especially for robots operating outdoors. In the modernization process of global navigation satellite system (GNSS) positioning, precise point positioning (PPP) has demonstrated its effectiveness in comparison to traditional differential positioning methods over the past thirty years. However, various challenges hinder the integration of PPP techniques into Internet of Things (IoT) systems for robot localization, with accuracy being a primary concern. This accuracy is impacted by factors such as satellite availability and signal disruptions in outdoor environments, resulting in less precise determination of satellite observations. Effectively addressing various GNSS errors is crucial when collecting PPP observations. The paper investigates the trade-off between kinematic PPP accuracy and cost effectiveness, through the examination of various influencing factors, including the choice of GNSS system (single or mixed), observation type (single or dual frequency), and satellite geometry. This research investigates kinematic PPP accuracy variation on a 10.4 km observed track based on different factors, using the GNSS system (single or mixed), and observation type (single or dual frequency). It can be concluded that mixed (GPS/GLONASS) dual frequency offers a 3D position accuracy of 9 cm, while mixed single frequency offers a 3D position accuracy of 13 cm. In industry, the results enable manufacturers to select suitable robot localization solutions according to the outdoor working environment (number of available satellites), economical constraint (single or dual frequency), and 3D position accuracy. Full article
(This article belongs to the Special Issue Advanced Robotics and Mechatronics)
Show Figures

Figure 1

Figure 1
<p>The architecture of CPS-based mobile robots for localization function.</p>
Full article ">Figure 2
<p>Study scope.</p>
Full article ">Figure 3
<p>The study’s observed kinematic track (KSU campus), Riyadh, KSA (8 October 2022).</p>
Full article ">Figure 4
<p>The rover setup for the study’s observed track (KSU campus), Riyadh, KSA (8 October 2022).</p>
Full article ">Figure 5
<p>Kinematic PPP coordinate differences using (<b>a</b>) GPS single-frequency observations, (<b>b</b>) GLONASS single-frequency observations, (<b>c</b>) mixed GPS/GLONASS single-frequency observations, (<b>d</b>) GPS dual-frequency observations, (<b>e</b>) GLONASS dual-frequency observations, and (<b>f</b>) mixed GPS/GLONASS dual-frequency observations.</p>
Full article ">Figure 5 Cont.
<p>Kinematic PPP coordinate differences using (<b>a</b>) GPS single-frequency observations, (<b>b</b>) GLONASS single-frequency observations, (<b>c</b>) mixed GPS/GLONASS single-frequency observations, (<b>d</b>) GPS dual-frequency observations, (<b>e</b>) GLONASS dual-frequency observations, and (<b>f</b>) mixed GPS/GLONASS dual-frequency observations.</p>
Full article ">Figure 6
<p>Kinematic PPP RMSE from GPS, GLONASS, and mixed GPS/GLONASS using (<b>a</b>) single-frequency observations, and (<b>b</b>) dual-frequency observations.</p>
Full article ">
29 pages, 6132 KiB  
Article
Investigation of the Causes of Railway Track Gauge Narrowing
by Péter Bocz, Nándor Liegner, Ákos Vinkó and Szabolcs Fischer
Vehicles 2023, 5(3), 949-977; https://doi.org/10.3390/vehicles5030052 - 10 Aug 2023
Cited by 1 | Viewed by 2138
Abstract
On behalf of MÁV Hungarian State Railways Ltd., the authors carried out a research and development (R&D) project on behalf of the Budapest University of Technology and Economics, Department of Highway and Railway Engineering, on the subject of “Research and investigation of the [...] Read more.
On behalf of MÁV Hungarian State Railways Ltd., the authors carried out a research and development (R&D) project on behalf of the Budapest University of Technology and Economics, Department of Highway and Railway Engineering, on the subject of “Research and investigation of the causes of gauge narrowing by finite-element modeling in running track and turnout, and under operational and laboratory conditions”. The main objective of the research was to investigate the causes of localized defects of gauge narrowing in railway tracks based on machine and manual track measurements, laboratory measurements, and theoretical considerations. The measures proposed as a consequence of identifying the causes could significantly contribute to reducing the number and extent of local defects in the future. Furthermore, the research aims to develop new theories in less scientifically mature areas and provide procedures and instructions that professional engineers and practitioners can easily apply. The main areas of research, which are not exhaustive, are as follows: (i) the evaluation of the measurement results provided by track geometry measuring and recording cars; (ii) on-site investigations in the railway track in terms of gauge and rail profile measurements; and, based on these, (iii) the selection of concrete sleepers, which were removed from the track and subjected to more detailed geometrical investigations in the laboratory, together with the components of the rail reinforcement; (iv) the track–vehicle connection, tight running in straight and curved track sections under track confinement; (v) modeling of the stability and deflection of the rail when the rail fastenings lose part of their supporting function; and (vi) finite element modeling of the concrete sleepers under operating conditions such as slow deformation of the concrete, temperature variation effects, and lateral support on the ballast. In the already-narrowed track section, the tight vehicle running is not the cause of the track gauge narrowing but a consequence, so it is not investigated in this paper. Full article
(This article belongs to the Special Issue Railway Vehicles and Infrastructure)
Show Figures

Figure 1

Figure 1
<p>Example of diagram indicating local gauge-narrowing errors (“mh.” means station stop, “ipvk.” means siding connection).</p>
Full article ">Figure 2
<p>Example of dividing the track into sections by the cumulative deviation from the mean method (the red vertical lines illustrate the considered section borders for the calculation of cumulative deviation between Budapest-Nyugati and Szajol railway stations).</p>
Full article ">Figure 3
<p>Change of gauge in straight sections of line 100a in the function of time.</p>
Full article ">Figure 4
<p>Change of gauge in straight sections of line 100a in the function of time (classified according to rail system).</p>
Full article ">Figure 5
<p>Schematic figures on the modeled <span class="html-italic">L4 SV-type</span> pre-tensioned reinforced concrete railway sleeper (above part: a 3D model with the steel tendons; below part: the applied mesh in the FE model).</p>
Full article ">Figure 6
<p>Central loads and supports in normal operation (the red–white tetrahedron symbolizes the load distribution from the rail head to rail foot to be able to consider realistic loading onto the modeled sleeper). Below the sleeper, a surface was considered to support the sleeper; see the meshed “plate”.</p>
Full article ">Figure 7
<p>Eccentric loads and supports (the red–white tetrahedron symbolizes the load distribution from the rail head to rail foot to be able to consider realistic loading onto the modeled sleeper). Below the sleeper, a surface was considered to support the sleeper; see the meshed “plate”.</p>
Full article ">Figure 8
<p>Calculated shrinkage of the concrete until the first 10,000 days after manufacture.</p>
Full article ">Figure 9
<p>Calculated shrinkage of the concrete until the first ten days after manufacture.</p>
Full article ">Figure 10
<p>The <span class="html-italic">x-direction</span> (<b>above part</b>) and total displacement (<b>below part</b>) plots, in the case of (a), are based on <a href="#sec4dot3-vehicles-05-00052" class="html-sec">Section 4.3</a>. The legend shows values with 10<sup>−3</sup> mm units.</p>
Full article ">Figure 11
<p>Result of gauge narrowing (all units are in mm) in the different load conditions (the legend from (a) to (h) are in accordance with <a href="#sec4dot3-vehicles-05-00052" class="html-sec">Section 4.3</a>).</p>
Full article ">Figure 12
<p>Forces acting on the rail with initial distortions (<span class="html-italic">l<sub>dist</sub></span> = 9 m, 54E1; the horizontal arrows symbolize the longitudinal forces in the track; hence, the vertical arrows mean the transversal forces).</p>
Full article ">Figure 13
<p>Bending moment in the rail with initial distortions (<span class="html-italic">l</span> = 9 m, 54E1).</p>
Full article ">Figure 14
<p>Lateral displacement of the examined section of rail.</p>
Full article ">
18 pages, 7112 KiB  
Article
Mobile Collaborative Heatmapping to Infer Self-Guided Walking Tourists’ Preferences for Geomedia
by Iori Sasaki, Masatoshi Arikawa, Min Lu and Ryo Sato
ISPRS Int. J. Geo-Inf. 2023, 12(7), 283; https://doi.org/10.3390/ijgi12070283 - 15 Jul 2023
Cited by 2 | Viewed by 1429
Abstract
This paper proposes a model-less feedback system driven by tourist tracking data that are automatically collected through mobile applications to visualize the gap between geomedia recommendations and the actual routes selected by tourists. High-frequency GPS data essentially make it difficult to interpret the [...] Read more.
This paper proposes a model-less feedback system driven by tourist tracking data that are automatically collected through mobile applications to visualize the gap between geomedia recommendations and the actual routes selected by tourists. High-frequency GPS data essentially make it difficult to interpret the semantic importance of hot spots and the presence of street-level features on a density map. Our mobile collaborative framework reorganizes tourist trajectories. This processing comprises (1) extracting the location of the user-generated content (UGC) recording, (2) abstracting the locations where tourists stay, (3) discarding locations where users remain stationary, and (4) simplifying the remaining points of location. Then, our heatmapping system visualizes heatmaps for hot streets, UGC-oriented hot spots, and indoor-oriented hot spots. According to our experimental study, this method can generate a trajectory that is more adaptable for hot street visualization than the raw trajectory and a simplified trajectory according to its geometry. This paper extends our previous work at the 2022 IEEE International Conference on Big Data, providing deeper discussions on application for local tourism. The framework allows us to derive insights for the development of guide content from mobile sensor data. Full article
(This article belongs to the Special Issue Urban Geospatial Analytics Based on Crowdsourced Data)
Show Figures

Figure 1

Figure 1
<p>Example of a heatmap with high-frequency GPS trajectories. There are too many factors that cause locally dense areas to properly judge their semantic importance. As the research subject area is Akita City in Japan, all background maps are in Japanese in this paper.</p>
Full article ">Figure 2
<p>Density maps using raw trajectories based on three values of <math display="inline"><semantics><mrow><msub><mrow><mi>T</mi><mi>h</mi></mrow><mrow><mi>c</mi><mo>.</mo><mi>r</mi><mo>.</mo></mrow></msub></mrow></semantics></math>. These maps are not compatible with hot street visualizations, as the topology of streets is not visible even after adjusting the color range.</p>
Full article ">Figure 3
<p>Structure realizing the feedback system on the basis of current mobile environments for walking tourism businesses. Our proposal for a novel heatmapping framework focuses on two sub-systems: (1) semi-ready data construction on the user side and (2) thematic heatmap generation to visualize hot spots and hot streets on the analyst side.</p>
Full article ">Figure 4
<p>A walking route in the experiments. A walker traced the blue line at a constant speed and stopped at each red point A, B, C, and D for one or two minutes. Gray rectangles depict indoor areas.</p>
Full article ">Figure 5
<p>Diagram of resampling process for calculating synchronous Euclidean distances between the ground truth and a target trajectory. A point <math display="inline"><semantics><mrow><msubsup><mrow><mi>p</mi></mrow><mrow><mi>s</mi></mrow><mrow><mo>′</mo></mrow></msubsup></mrow></semantics></math> is added to maintain time ratio.</p>
Full article ">Figure 6
<p>Total SED of the target trajectory data (red line: <math display="inline"><semantics><mrow><msub><mrow><mi>T</mi></mrow><mrow><mi>r</mi><mi>a</mi><mi>w</mi></mrow></msub></mrow></semantics></math>; brown dashed line: <math display="inline"><semantics><mrow><msub><mrow><mi>T</mi></mrow><mrow><mi>D</mi><mi>P</mi></mrow></msub></mrow></semantics></math>; blue dashed line <math display="inline"><semantics><mrow><msub><mrow><mi>T</mi></mrow><mrow><mi>S</mi><mi>R</mi></mrow></msub></mrow></semantics></math> ). This implies that the proposed method can decrease total SED with a small tolerance parameter.</p>
Full article ">Figure 7
<p>Trajectory shape (<b>left</b>) and time series changes in the SED (<b>right</b>) of <math display="inline"><semantics><mrow><msub><mrow><mi>T</mi></mrow><mrow><mi>r</mi><mi>a</mi><mi>w</mi></mrow></msub></mrow></semantics></math>. Orange areas in the graph of time series changes represent the periods when the user is stationary outdoor and indoor, as indicated by the red points in <a href="#ijgi-12-00283-f004" class="html-fig">Figure 4</a> (A, B, C, and D, in order).</p>
Full article ">Figure 8
<p>Trajectory shape (<b>left</b>) and time series changes in the SED (<b>right</b>) of <math display="inline"><semantics><mrow><msub><mrow><mi>T</mi></mrow><mrow><mi>D</mi><mi>P</mi></mrow></msub></mrow></semantics></math>. The tolerance parameter <math display="inline"><semantics><mrow><mi>ε</mi></mrow></semantics></math> is set to 12.0 m. Orange areas in the graph of time series changes represent the periods when the user is stationary outdoor and indoor, as indicated by the red points in <a href="#ijgi-12-00283-f004" class="html-fig">Figure 4</a> (A, B, C, and D, in order).</p>
Full article ">Figure 9
<p>Trajectory shape (<b>left</b>) and time series changes in the SED (<b>right</b>) of <math display="inline"><semantics><mrow><msub><mrow><mi>T</mi></mrow><mrow><mi>S</mi><mi>R</mi></mrow></msub></mrow></semantics></math>. The tolerance parameter <math display="inline"><semantics><mrow><mi>ε</mi></mrow></semantics></math> is set to 1.0 m. Orange areas in the graph of time series changes represent the periods when the user is stationary outdoor and indoor, as indicated by the red points in <a href="#ijgi-12-00283-f004" class="html-fig">Figure 4</a> (A, B, C, and D, in order).</p>
Full article ">Figure 10
<p>Recommended spots with IDs from one to nine and walking routes in the walking guidebook that is available on [<a href="#B44-ijgi-12-00283" class="html-bibr">44</a>] for Japanese tourists. Red pins are facilities where tourists can stay, and green pins are monuments or viewpoints they can look at.</p>
Full article ">Figure 11
<p>Location-based services: (<b>a</b>) positioning the current location on the illustrated maps which is provided in a Japanese tourist guidebook published by Akita City; (<b>b</b>) location-based push services that automatically display geomedia, such as Japanese guide scripts and pictures, on the screen when the user gets close to the registered spots.</p>
Full article ">Figure 12
<p>Example of the distribution of horizontal GPS accuracy values, obtained by monitoring twelve subjects within the dataset The device used was iPhone 11, manufactured by Apple Inc., based in Cupertino, California, USA. The kCLLocationAccuracyBest setting was applied, which is specified when very high accuracy is required in Core Location framework. The left-side graph represents an outdoor condition, i.e., street between spots 7 and 9 in <a href="#ijgi-12-00283-f010" class="html-fig">Figure 10</a>, and the right-side graph represents an indoor condition, i.e., spot 7 in <a href="#ijgi-12-00283-f010" class="html-fig">Figure 10</a>.</p>
Full article ">Figure 13
<p>An example of a hot street heatmap. Equalizing the density per area enables visualization of the presence of polyline-shaped features, such as walking routes and streets.</p>
Full article ">Figure 14
<p>An example of a UGC-oriented hot spot heatmap that considers point data drawn only from <span class="html-italic">ugc</span> tags. Dense areas represent attractive photo spots and places that are worth sharing.</p>
Full article ">Figure 15
<p>An example of an indoor-oriented hot spot heatmap that considers point data drawn only from <span class="html-italic">indoor</span> tags. Dense areas represent attractive buildings and facilities visited by many tourists.</p>
Full article ">Figure 16
<p>Heatmaps that were used for a user experiment. The experiment involved the generation of heatmaps from raw data and semi-ready data using different values for <math display="inline"><semantics><mrow><mi>T</mi><msub><mrow><mi>h</mi></mrow><mrow><mi>c</mi><mo>.</mo><mi>r</mi><mo>.</mo></mrow></msub></mrow></semantics></math>.</p>
Full article ">Figure 17
<p>Stacked bar chart of the selection distribution of heatmaps ranked as the top three.</p>
Full article ">
18 pages, 164375 KiB  
Article
Evaluating Urban Building Damage of 2023 Kahramanmaras, Turkey Earthquake Sequence Using SAR Change Detection
by Xiuhua Wang, Guangcai Feng, Lijia He, Qi An, Zhiqiang Xiong, Hao Lu, Wenxin Wang, Ning Li, Yinggang Zhao, Yuedong Wang and Yuexin Wang
Sensors 2023, 23(14), 6342; https://doi.org/10.3390/s23146342 - 12 Jul 2023
Cited by 13 | Viewed by 3778
Abstract
On February 6, 2023 (local time), two earthquakes (Mw7.8 and Mw7.7) struck central and southern Turkey, causing extensive damage to several cities and claiming a toll of 40,000 lives. In this study, we propose a method for seismic building damage assessment and analysis [...] Read more.
On February 6, 2023 (local time), two earthquakes (Mw7.8 and Mw7.7) struck central and southern Turkey, causing extensive damage to several cities and claiming a toll of 40,000 lives. In this study, we propose a method for seismic building damage assessment and analysis by combining SAR amplitude and phase coherence change detection. We determined building damage in five severely impacted urban areas and calculated the damage ratio by measuring the urban area and the damaged area. The largest damage ratio of 18.93% is observed in Nurdagi, and the smallest ratio of 7.59% is found in Islahiye. We verified the results by comparing them with high-resolution optical images and AI recognition results from the Microsoft team. We also used pixel offset tracking (POT) technology and D-InSAR technology to obtain surface deformation using Sentinel-1A images and analyzed the relationship between surface deformation and post-earthquake urban building damage. The results show that Nurdagi has the largest urban average surface deformation of 0.48 m and Antakya has the smallest deformation of 0.09 m. We found that buildings in the areas with steeper slopes or closer to earthquake faults have higher risk of collapse. We also discussed the influence of SAR image parameters on building change recognition. Image resolution and observation geometry have a great influence on the change detection results, and the resolution can be improved by various means to raise the recognition accuracy. Our research findings can guide earthquake disaster assessment and analysis and identify influential factors of earthquake damage. Full article
(This article belongs to the Special Issue Sensing and Signal Analysis in Synthetic Aperture Radar Systems)
Show Figures

Figure 1

Figure 1
<p>Flowchart used in this study.</p>
Full article ">Figure 2
<p>Surface deformation of the five urban areas caused by the 2023 Mw7.7 and Mw7.8 Turkey earthquakes obtained by POT. (<b>a</b>) Surface deformation of central southern Turkey. The blue box indicates the Sentinel-1 image coverage; the dark gray lines are active faults; the magenta lines are seismogenic faults; and the red pentagrams indicate the epicenter of the two main earthquakes and the Mw6.8 aftershock. The epicenter and focal mechanism are cited from GCMT. The red dots are the location of the affected cities. The inset shows the regional seismotectonic background: The red box is the study area shown in (<b>a</b>); the black line is the fault zone; the north one is North Anatolian fault zone (NAFZ); and the south are East Anatolian fault zone (EAFZ) and the Death Sea fault zone (DSFZ). The red pentagonal stars are the epicenters, (<b>b</b>–<b>f</b>) zoom-in of the deformation map of the five cities. The red polygon delineates the urban boundary.</p>
Full article ">Figure 3
<p>Urban building damage identification results with descending data. The red polygon delineates the urban boundary. Yellow-filled polygons are the identified damaged buildings.</p>
Full article ">Figure 4
<p>Comparison of the building damage identification results with optical images in urban areas. The first column is the pre-earthquake optical images; the second column is the post-earthquake optical images superimposed damage area contour (cyan polygon). The third column is the post-earthquake optical image superimposed damage area contours and recognition results (yellow filled polygon), and the forth column is the post-earthquake optical image superimposed damage area contours and coherence results. (<b>a</b>–<b>d</b>) is in Turkoglu; (<b>e</b>–<b>h</b>) is in Islahiye; (<b>i</b>–<b>p</b>) is in Marash; and (<b>q</b>–<b>t</b>) is in Nurdagi.</p>
Full article ">Figure 5
<p>Comparison of urban building damage identification results with Microsoft artificial intelligence identification results. The first column is the post-earthquake optical images; the second column is the post-earthquake optical images superimposed on the artificial intelligence identification results of the Microsoft team (red block is the damaged building; blue block is the undamaged building). The third column is the post-earthquake optical images superimposed on the identification results in this study. (<b>a</b>) is in Marash; (<b>b</b>) is in Turkoglu; (<b>c</b>) is in Nurdagi; (<b>d</b>) is in Islahiye.</p>
Full article ">Figure 6
<p>Building damage identification results. The number of damaged buildings comes from the artificial intelligence identification results of the Microsoft team.</p>
Full article ">Figure 7
<p>Radar chart of influencing factors of building damage.</p>
Full article ">Figure 8
<p>Comprehensive analysis chart of influencing factors. The red box above is the range of the damage ratio of the five cities. The following red box is the range of the superposition results of the average slope and the distance between the urban center and the seismogenic fault of the five cities.</p>
Full article ">Figure 9
<p>Building echo examples. Gray-filled rectangles are the buildings that receive radar microwave signal, and the white-filled rectangle is the building that does not receive radar microwave signal. Black arrows are radar signals. (<b>a</b>) The microwave signal direction is on the left of the building. (<b>b</b>) The microwave signal direction is on the right of the building.</p>
Full article ">Figure 10
<p>Urban building damage identification results with ascending data. The red polygon delineates the urban boundary. Yellow-filled polygons are the identified damaged buildings.</p>
Full article ">Figure A1
<p>Surface deformation of the five urban areas caused by the 2023 Mw7.7 and Mw7.8 Turkey earthquakes which were obtained by D-InSAR. (<b>a</b>) Surface deformation of central southern Turkey. The blue box indicates the Sentinel-1 image cover; the dark gray lines are active faults; the yellow lines are seismogenic faults; and the magenta pentagrams indicate the epicenter of the two main earthquakes and the Mw6.8 aftershock. The epicenter and focal mechanism are cited from GCMT. The magenta circles with white outline are the location of the affected cities. The inset shows surface deformation results with the unwrapping of central southern Turkey: the yellow lines are seismogenic faults, and the magenta pentagrams indicate the epicenter of the two main earthquakes and the Mw6.8 aftershock. The magenta dots are the location of the affected cities. (<b>b</b>–<b>f</b>) zoom-in of the deformation map of the five cities. The red polygon delineates the urban boundary.</p>
Full article ">
18 pages, 8822 KiB  
Article
Data Comparison and Cross-Calibration between Level 1 Products of DPC and POSP Onboard the Chinese GaoFen-5(02) Satellite
by Xuefeng Lei, Zhenhai Liu, Fei Tao, Hao Dong, Weizhen Hou, Guangfeng Xiang, Lili Qie, Binghuan Meng, Congfei Li, Feinan Chen, Yanqing Xie, Miaomiao Zhang, Lanlan Fan, Liangxiao Cheng and Jin Hong
Remote Sens. 2023, 15(7), 1933; https://doi.org/10.3390/rs15071933 - 4 Apr 2023
Cited by 2 | Viewed by 1956
Abstract
The Polarization CrossFire (PCF) suite onboard the Chinese GaoFen-5(02) satellite has been sophisticatedly composed by the Particulate Observing Scanning Polarimeter (POSP) and the Directional Polarimetric Camera (DPC). Among them, DPC is a multi-angle sequential measurement polarization imager, while POSP is a cross-track scanning [...] Read more.
The Polarization CrossFire (PCF) suite onboard the Chinese GaoFen-5(02) satellite has been sophisticatedly composed by the Particulate Observing Scanning Polarimeter (POSP) and the Directional Polarimetric Camera (DPC). Among them, DPC is a multi-angle sequential measurement polarization imager, while POSP is a cross-track scanning simultaneous polarimeter with corresponding radiometric and polarimetric calibrators, which can theoretically be used for cross comparison and calibration with DPC. After the data preprocessing of these two sensors, we first select local homogeneous cluster scenes by calculating the local variance-to-mean ratio in DPC’s Level 1 product projection grids to reduce the influence of scale differences and geometry misalignment between DPC and POSP. Then, taking the observation results after POSP data quality assurance as the abscissa and taking the DPC observation results under the same wavelength band and geometric conditions as the same ordinate, a two-dimensional radiation/polarization feature space is established. Results show that the normalized top of the atmosphere (TOA) radiances of DPC and POSP processed data at the nadir are linearly correlated. The normalized TOA radiance root mean square errors (RMSEs) look reasonable in all common bands. The DPC and POSP normalized radiance ratios in different viewing zenith angle ranges at different times reveal the temporal drift of the DPC relative radiation response. The RMSEs, mean absolute errors (MAEs), relative errors (REs), and scatter percentage of DPC degree of linear polarization (DoLP) falling within the expected error (EE = ±0.02) of POSP measured DoLP are better than 0.012, 0.009, 0.066, and 91%, respectively. Full article
Show Figures

Figure 1

Figure 1
<p>The sensor assembly of the PCF suite installed on the GaoFen-5(02) satellite.</p>
Full article ">Figure 2
<p>The layout diagram of the POSP onboard calibrators.</p>
Full article ">Figure 3
<p>The sampling schematic diagram of PCF [<a href="#B27-remotesensing-15-01933" class="html-bibr">27</a>].</p>
Full article ">Figure 4
<p>Normalized relative spectral responses of DPC and POSP common bands.</p>
Full article ">Figure 5
<p>The linear fitting coefficients of DPC and POSP normalized radiance in common spectral bands at nadir. The data pairs are a collection selected from May 2022 of (<b>a</b>) 443 nm; (<b>b</b>) 490 nm; (<b>c</b>) 670 nm; and (<b>d</b>) 865 nm bands in DPC and POSP Level 1 products. Yellow solid and red dashed lines are the 1:1 lines and fit lines, respectively. The reason for the fewer matching data points is that the amount of raw data points at nadir is small.</p>
Full article ">Figure 5 Cont.
<p>The linear fitting coefficients of DPC and POSP normalized radiance in common spectral bands at nadir. The data pairs are a collection selected from May 2022 of (<b>a</b>) 443 nm; (<b>b</b>) 490 nm; (<b>c</b>) 670 nm; and (<b>d</b>) 865 nm bands in DPC and POSP Level 1 products. Yellow solid and red dashed lines are the 1:1 lines and fit lines, respectively. The reason for the fewer matching data points is that the amount of raw data points at nadir is small.</p>
Full article ">Figure 6
<p>The time-varying characteristics of DPC and POSP normalized radiance linear fitting results (<b>a</b>) <math display="inline"><semantics> <mrow> <msup> <mrow> <mi>A</mi> </mrow> <mrow> <mi>k</mi> </mrow> </msup> </mrow> </semantics></math>; (<b>b</b>) R<sup>2</sup>; and (<b>c</b>) RMSE in common spectral bands at nadir. These data pairs are selected from collections from November 2021 to May 2022 in DPC and POSP Level 1 products.</p>
Full article ">Figure 7
<p>The linear fitting slope and corresponding quadratic fitting curve of DPC and POSP normalized radiance in common spectral bands with different post-launch days and different VZA ranges. The data pairs are collected in each day selected from October 2021 to July 2022 of (<b>a</b>) 443 nm; (<b>b</b>) 490 nm; (<b>c</b>) 670 nm; and (<b>d</b>) 865 nm bands in DPC and POSP Level 1 products.</p>
Full article ">Figure 8
<p>The ratio of DPC and POSP normalized radiance in the 443 nm band with different VZAs. These data pairs are selected for collection from (<b>a</b>) November 2021; (<b>b</b>) December 2021; (<b>c</b>) January 2022; (<b>d</b>) February 2022; (<b>e</b>) April 2022; and (<b>f</b>) May 2022 in DPC and POSP Level 1 products.</p>
Full article ">Figure 9
<p>The ratio of DPC to POSP normalized radiance in the 865 nm band with different VZAs. These data pairs are selected for collection from (<b>a</b>) November 2021; (<b>b</b>) December 2021; (<b>c</b>) January 2022; (<b>d</b>) February 2022; (<b>e</b>) April 2022; and (<b>f</b>) May 2022 in DPC and POSP Level 1 products.</p>
Full article ">Figure 10
<p>The ratio of DPC to POSP normalized radiance in the 490 nm band with different VZAs. These data pairs are selected for collection from (<b>a</b>) November 2021; (<b>b</b>) December 2021; (<b>c</b>) January 2022; (<b>d</b>) February 2022; (<b>e</b>) April 2022; and (<b>f</b>) May 2022 in DPC and POSP Level 1 products.</p>
Full article ">Figure 11
<p>The ratio of DPC to POSP normalized radiance in the 670 nm band with different VZAs. These data pairs are selected for collection from (<b>a</b>) November 2021; (<b>b</b>) December 2021; (<b>c</b>) January 2022; (<b>d</b>) February 2022; (<b>e</b>) April 2022; and (<b>f</b>) May 2022 in DPC and POSP Level 1 products.</p>
Full article ">Figure 11 Cont.
<p>The ratio of DPC to POSP normalized radiance in the 670 nm band with different VZAs. These data pairs are selected for collection from (<b>a</b>) November 2021; (<b>b</b>) December 2021; (<b>c</b>) January 2022; (<b>d</b>) February 2022; (<b>e</b>) April 2022; and (<b>f</b>) May 2022 in DPC and POSP Level 1 products.</p>
Full article ">Figure 12
<p>The linear fitting and statistical distribution results of DPC and POSP DoLP in common polarimetric bands. The data pairs are selected for collection from October 2021 to May 2022 for (<b>a</b>) and (<b>d</b>) 490 nm; (<b>b</b>) and (<b>e</b>) 670 nm; and (<b>c</b>,<b>f</b>) 865 nm bands in DPC and POSP Level 1 products. Yellow solid, green, and red dashed lines are the 1:1 lines, EE envelope lines, and fit lines, respectively.</p>
Full article ">Figure 12 Cont.
<p>The linear fitting and statistical distribution results of DPC and POSP DoLP in common polarimetric bands. The data pairs are selected for collection from October 2021 to May 2022 for (<b>a</b>) and (<b>d</b>) 490 nm; (<b>b</b>) and (<b>e</b>) 670 nm; and (<b>c</b>,<b>f</b>) 865 nm bands in DPC and POSP Level 1 products. Yellow solid, green, and red dashed lines are the 1:1 lines, EE envelope lines, and fit lines, respectively.</p>
Full article ">
16 pages, 4654 KiB  
Article
Residual Stress Formation Mechanisms in Laser Powder Bed Fusion—A Numerical Evaluation
by Moritz Kaess, Martin Werz and Stefan Weihe
Materials 2023, 16(6), 2321; https://doi.org/10.3390/ma16062321 - 14 Mar 2023
Cited by 6 | Viewed by 2214
Abstract
Additive manufacturing methods, such as the laser powder bed fusion, do not need any special tool or casting mold. This enables the fast realization of complex and individual geometries with integrated functions. However, the local heat input during the manufacturing process often leads [...] Read more.
Additive manufacturing methods, such as the laser powder bed fusion, do not need any special tool or casting mold. This enables the fast realization of complex and individual geometries with integrated functions. However, the local heat input during the manufacturing process often leads to residual stresses and distortion. This in turn causes poor quality, scrap parts or can even terminate a job prematurely if the powder recoating mechanism collides with a distorted part during the process. This study investigates the generation mechanisms of residual stresses and distortion during laser powder bed fusion (LPBF) of stainless steel 316L in order to reduce these effects and thus contribute to improved process safety and efficiency. Therefore, numerical investigations with a finite element model on the scale of a few melt tracks and layers serve to develop a detailed understanding of the mechanisms during production. The work includes an investigation of the build plate temperature, the laser power and speed and the layer thickness. The results show a strong dependency on the build plate preheating and energy per unit length. A higher build plate temperature and a reduction of the energy per unit length both lead to lower residual stresses. Full article
(This article belongs to the Section Manufacturing Processes and Systems)
Show Figures

Figure 1

Figure 1
<p>Distortions and delamination during manufacturing process leading to process failure.</p>
Full article ">Figure 2
<p>Approach used: weakly-coupled thermomechanical simulation with boundary conditions.</p>
Full article ">Figure 3
<p>Finite element model—total dimensions and model cut with influence area of heat source.</p>
Full article ">Figure 4
<p>Temperature dependent material properties of 316L: (<b>a</b>) Young’s Modulus (value at room temperature [<a href="#B33-materials-16-02321" class="html-bibr">33</a>] vertical; slope from 100 °C to 500 °C [<a href="#B36-materials-16-02321" class="html-bibr">36</a>,<a href="#B37-materials-16-02321" class="html-bibr">37</a>]; values at melting temperature: assumption) (<b>b</b>) Flow curve (values at room temperature [<a href="#B33-materials-16-02321" class="html-bibr">33</a>] vertical—rounded; values at higher temperatures: assumption) (<b>c</b>) Thermal conductivity (derived from [<a href="#B32-materials-16-02321" class="html-bibr">32</a>,<a href="#B37-materials-16-02321" class="html-bibr">37</a>,<a href="#B38-materials-16-02321" class="html-bibr">38</a>,<a href="#B39-materials-16-02321" class="html-bibr">39</a>,<a href="#B40-materials-16-02321" class="html-bibr">40</a>]) (<b>d</b>) Specific heat (derived from [<a href="#B32-materials-16-02321" class="html-bibr">32</a>,<a href="#B37-materials-16-02321" class="html-bibr">37</a>,<a href="#B38-materials-16-02321" class="html-bibr">38</a>,<a href="#B39-materials-16-02321" class="html-bibr">39</a>,<a href="#B40-materials-16-02321" class="html-bibr">40</a>]) (<b>e</b>) Thermal expansion (values up to 500 °C [<a href="#B36-materials-16-02321" class="html-bibr">36</a>,<a href="#B37-materials-16-02321" class="html-bibr">37</a>]; values from 500 °C: slope assumed according to [<a href="#B39-materials-16-02321" class="html-bibr">39</a>]).</p>
Full article ">Figure 4 Cont.
<p>Temperature dependent material properties of 316L: (<b>a</b>) Young’s Modulus (value at room temperature [<a href="#B33-materials-16-02321" class="html-bibr">33</a>] vertical; slope from 100 °C to 500 °C [<a href="#B36-materials-16-02321" class="html-bibr">36</a>,<a href="#B37-materials-16-02321" class="html-bibr">37</a>]; values at melting temperature: assumption) (<b>b</b>) Flow curve (values at room temperature [<a href="#B33-materials-16-02321" class="html-bibr">33</a>] vertical—rounded; values at higher temperatures: assumption) (<b>c</b>) Thermal conductivity (derived from [<a href="#B32-materials-16-02321" class="html-bibr">32</a>,<a href="#B37-materials-16-02321" class="html-bibr">37</a>,<a href="#B38-materials-16-02321" class="html-bibr">38</a>,<a href="#B39-materials-16-02321" class="html-bibr">39</a>,<a href="#B40-materials-16-02321" class="html-bibr">40</a>]) (<b>d</b>) Specific heat (derived from [<a href="#B32-materials-16-02321" class="html-bibr">32</a>,<a href="#B37-materials-16-02321" class="html-bibr">37</a>,<a href="#B38-materials-16-02321" class="html-bibr">38</a>,<a href="#B39-materials-16-02321" class="html-bibr">39</a>,<a href="#B40-materials-16-02321" class="html-bibr">40</a>]) (<b>e</b>) Thermal expansion (values up to 500 °C [<a href="#B36-materials-16-02321" class="html-bibr">36</a>,<a href="#B37-materials-16-02321" class="html-bibr">37</a>]; values from 500 °C: slope assumed according to [<a href="#B39-materials-16-02321" class="html-bibr">39</a>]).</p>
Full article ">Figure 5
<p>Temperature distribution [°C] around hemispherical heat source—only elements in solid/liquid state are displayed.</p>
Full article ">Figure 6
<p>Scanning strategy—67° rotation.</p>
Full article ">Figure 7
<p>Model cut: cantilever after partial separation—path for distortion evaluation marked in red and path for stress evaluation marked in black.</p>
Full article ">Figure 8
<p>Results: Bending (<b>a</b>) and residual stress (<b>b</b>) for a variation of build plate temperature; bending (<b>c</b>) and residual stress (<b>d</b>) for a variation of layer thickness; bending (<b>e</b>) and residual stress (<b>f</b>) for a variation of laser power; bending (<b>g</b>) and residual stress (<b>h</b>) for a variation of laser power &amp; speed.</p>
Full article ">Figure 9
<p>Maximum bending at cantilever front position dependent on build plate temperature (<b>a</b>) and laser power (<b>b</b>)—fixed speed brings variation in LED—varying speed keeps LED constant.</p>
Full article ">Figure 10
<p>Cantilever bending (<b>a</b>,<b>b</b>) and residual stress (<b>c</b>,<b>d</b>) for laser power variation from layer to layer.</p>
Full article ">
14 pages, 4652 KiB  
Article
A Methodology Linking Tamping Processes and Railway Track Behaviour
by Stefan Offenbacher, Christian Koczwara, Matthias Landgraf and Stefan Marschnig
Appl. Sci. 2023, 13(4), 2137; https://doi.org/10.3390/app13042137 - 7 Feb 2023
Cited by 3 | Viewed by 2731
Abstract
Today’s railway transport is built upon high-performance infrastructure. Cost-effective yet sustainable infrastructure presumes tracks with a precise and durable geometry. At ballasted tracks, the geometry is created and restored through tamping machines, which position the track panel and compact the ballast beneath the [...] Read more.
Today’s railway transport is built upon high-performance infrastructure. Cost-effective yet sustainable infrastructure presumes tracks with a precise and durable geometry. At ballasted tracks, the geometry is created and restored through tamping machines, which position the track panel and compact the ballast beneath the sleepers. It is commonly agreed that the ballast compaction plays an important role in the long-term stability of the track. Yet, there is no method available which allows a direct correlation between the compactness of the ballast and the stability of the track geometry. Available studies either model track behaviour without considering the bedding, or analyse ballast compactness locally while disregarding its influence on the track geometry. This paper presents a new methodology which establishes a relation between these two topics—ballast compaction during tamping and subsequent track behaviour. A state-of-the-art tamping machine has been equipped with an experimental measurement setup, constantly recording relevant data during every tamping process. These data can be used to derive an indication for the achieved compaction under every sleeper. Utilising the tamping machine’s internal measuring system for track geometry documentation, every tamping process (every sleeper) is assigned to the precise position along the track. The data set is merged and synchronised with regular track geometry measurements of the infrastructure manager. The result is a comprehensive data set which allows precise analyses between tamping machine measurements and track behaviour. This data set provides the foundation for future research, aiming towards a better understanding of the tamping process and its influence on the quality and durability of the established track geometry. Full article
(This article belongs to the Special Issue Railway Infrastructures Engineering: Latest Advances and Prospects)
Show Figures

Figure 1

Figure 1
<p>Concept and scope of this paper.</p>
Full article ">Figure 2
<p>Principal track geometry parameters according to EN 13848-1; red lines illustrate the deviation from the ideal geometry.</p>
Full article ">Figure 3
<p>Theoretical track behaviour by means of the vertical track geometry quality.</p>
Full article ">Figure 4
<p>Main phases of a tamping process.</p>
Full article ">Figure 5
<p>(<b>a</b>) Tamping machine Unimat 09-4x4/4S E<sup>3</sup> and (<b>b</b>,<b>c</b>) measurement equipment installed on its four tamping unit segments.</p>
Full article ">Figure 6
<p>(<b>a</b>) Raw measurements and low-pass filtered data of the squeezing distance and (<b>b</b>) squeezing force and calculated parameters; the data originate from two different squeezing processes.</p>
Full article ">Figure 7
<p>Illustration of the outlier detection algorithm.</p>
Full article ">Figure 8
<p>(<b>a</b>) Three metre sine wave and (<b>b</b>) 25 m sine wave and respective standard deviations with influence lengths of 100, 25, and 10 m.</p>
Full article ">Figure 9
<p>Longitudinal level (LL) D1 signals and standard deviations (SD) calculated with a moving window over 25 m; (<b>a</b>) before, (<b>b</b>) after the synchronisation process. Each colour represents an individual measurement run.</p>
Full article ">Figure 10
<p>Merging process of sleeper-specific tamping data with the DRP control measurement.</p>
Full article ">Figure 11
<p>Merging tamping data and track geometry data by means of longitudinal level measurements; (<b>a</b>) before, (<b>b</b>) after the synchronisation process.</p>
Full article ">Figure 12
<p>Application of the tamping machine measurements: Ballast condition assessment (penetration force per tamping tine; modified from [<a href="#B48-applsci-13-02137" class="html-bibr">48</a>]).</p>
Full article ">
30 pages, 14431 KiB  
Article
AdaSplats: Adaptive Splatting of Point Clouds for Accurate 3D Modeling and Real-Time High-Fidelity LiDAR Simulation
by Jean Pierre Richa, Jean-Emmanuel Deschaud, François Goulette and Nicolas Dalmasso
Remote Sens. 2022, 14(24), 6262; https://doi.org/10.3390/rs14246262 - 10 Dec 2022
Cited by 7 | Viewed by 5140
Abstract
LiDAR sensors provide rich 3D information about their surroundings and are becoming increasingly important for autonomous vehicles tasks such as localization, semantic segmentation, object detection, and tracking. Simulation accelerates the testing, validation, and deployment of autonomous vehicles while also reducing cost and eliminating [...] Read more.
LiDAR sensors provide rich 3D information about their surroundings and are becoming increasingly important for autonomous vehicles tasks such as localization, semantic segmentation, object detection, and tracking. Simulation accelerates the testing, validation, and deployment of autonomous vehicles while also reducing cost and eliminating the risks of testing in real-world scenarios. We address the problem of high-fidelity LiDAR simulation and present a pipeline that leverages real-world point clouds acquired by mobile mapping systems. Point-based geometry representations, more specifically splats (2D oriented disks with normals), have proven their ability to accurately model the underlying surface in large point clouds, mainly with uniform density. We introduce an adaptive splat generation method that accurately models the underlying 3D geometry to handle real-world point clouds with variable densities, especially for thin structures. Moreover, we introduce a fast LiDAR sensor simulator, working in the splatted model, that leverages the GPU parallel architecture with an acceleration structure while focusing on efficiently handling large point clouds. We test our LiDAR simulation in real-world conditions, showing qualitative and quantitative results compared to basic splatting and meshing techniques, demonstrating the interest of our modeling technique. Full article
(This article belongs to the Special Issue Point Cloud Processing in Remote Sensing Technology)
Show Figures

Figure 1

Figure 1
<p>Starting with a point cloud acquired using a mobile mapping system (MMS), we obtain point-wise semantic labels by performing semantic segmentation. Using the semantic labels, we remove dynamic objects in the scene and perform our splats generation method. The splatted scene can then be used to simulate the different sensors. Dynamic objects can be added to the splatted scene either in the form of splatted point clouds or using a bank of CAD meshed models.</p>
Full article ">Figure 2
<p>Splats generation starts by including points in the neighborhood, until the error bounds are exceeded, then the center of the splat is moved along the normal vector to minimize the distance from the splat to the neighboring points.</p>
Full article ">Figure 3
<p>Illustrating the stopping cases to ensure the preservation of sharp features and avoid classes interference in the splats generation.</p>
Full article ">Figure 4
<p>Parallel ray casting and accelerated ray–splat intersection are achieved using OptiX. The BVH is created from the splats primitives; then, the rays are cast in parallel on the device. Each ray traverses the BVH, and an intersection is reported back if a hit is found. Otherwise, the intersection for the specific ray is ignored.</p>
Full article ">Figure 5
<p>Our pipeline is split into different modules. The first generates accurate 3D modeling of the static environment using our adaptive splatting method. The second module takes the sensor model as input and simulates the sensor, including but not limited to camera and LiDAR, and generates the corresponding rays. In the third module, the rays are cast in parallel using the GPU architecture, then a bounding volume hierarchy (BVH) structure containing the generated splats is traversed, and ray–splat intersection point or color information is reported to generate the LiDAR or camera output, respectively.</p>
Full article ">Figure 6
<p>Point clouds used in the experiments (left to right:) PC3D-Paris, SemanticKITTI, and M-City. In red, we show the trajectory used for simulation.</p>
Full article ">Figure 7
<p>Rendering results on different choices of <math display="inline"><semantics> <mi mathvariant="script">K</mi> </semantics></math>-nn on PC3D-Paris dataset. A small <math display="inline"><semantics> <mi mathvariant="script">K</mi> </semantics></math> (e.g., 10 or 20) results in holes on the surface and ground groups while also resulting in a better approximation on the non-surface and linear groups. A large <math display="inline"><semantics> <mi mathvariant="script">K</mi> </semantics></math> (40 to 120) results in a hole-free approximation of the surface and ground semantic groups while creating artifacts on small structures belonging to the linear and surface groups.</p>
Full article ">Figure 8
<p>Rendering the different surface representations on PC3D-Paris. The top row shows the meshed scene using IMLS (left) and Poisson (right). The middle row shows the splatted scene using basic splats (left) and AdaSplats using KPConv semantics (right). The bottom row shows the splatted scene using AdaSplats-GT, which contains the ground truth point-wise semantic information (left) and the original point cloud (right). We show the ability of AdaSplats to recover a better geometry, especially on fine structures (in green, red and yellow boxes).</p>
Full article ">Figure 9
<p>Comparison of simulated LiDAR data using different reconstruction and modeling methods on PC3D-Paris. The top row shows the simulation in meshed IMLS (left) and Poisson (right). The middle row shows the simulation with Basic Splats (left) and AdaSplats-KPConv (right). The bottom row shows the simulation with AdaSplats-GT (left) and original point cloud (right).</p>
Full article ">Figure 10
<p>Rendering the different surface representations on SemanticKITTI. The top row shows the meshed scene using IMLS (<b>left</b>) and Poisson (<b>right</b>). The middle row shows the splatted scene using basic splats (<b>left</b>) and AdaSplats using KPConv semantics (<b>right</b>). The bottom row shows the splatted scene using AdaSplats-GT, which contains the ground truth point-wise semantic information (<b>left</b>) and the original point cloud (<b>right</b>).</p>
Full article ">Figure 11
<p>Comparison of simulated LiDAR data using different reconstruction and modeling methods on SemanticKITTI. The top row shows the simulation in meshed IMLS (<b>left</b>) and Poisson (<b>right</b>). The middle row shows the simulation with Basic Splats (<b>left</b>) and AdaSplats-KPConv (<b>right</b>). The bottom row shows the simulation with AdaSplats-GT (<b>left</b>) and original point cloud (<b>right</b>).</p>
Full article ">Figure 12
<p>Rendering the different surface representations on M-City. The top row shows the manually meshed scene (<b>left</b>) and basic splats (<b>right</b>). The bottom row shows the results of rendering AdaSplats using GT semantics (<b>left</b>) and the original point cloud (<b>right</b>).</p>
Full article ">Figure 13
<p>Comparison of simulated LiDAR data using different reconstruction and modeling methods on M-City. The top row shows the simulation in the manually meshed scene (<b>left</b>) and the modeled scene with Basic Splats (<b>right</b>). The bottom row shows the simulation with AdaSplats-GT (<b>left</b>) and original point cloud (<b>right</b>). Modeling vegetation is not an easy task and usually requires different ray–primitive intersection methods.</p>
Full article ">Figure 14
<p>Showing an original frame from the SemanticKITTI [<a href="#B18-remotesensing-14-06262" class="html-bibr">18</a>] sequence 08 dataset with dynamic objects (<b>top</b>). The simulated HDL-64 LiDAR at the same position with dynamic objects (<b>middle</b>). The simulated HDL-32 LiDAR translated by −0.5 m on the <span class="html-italic">z</span>-axis (<b>bottom</b>).</p>
Full article ">
15 pages, 40805 KiB  
Article
Experimental Study on the Target–Receiver Formation Problem with the Exploitation of Coherent and Non-Coherent Bearing Information
by Lu Wang, Shiliang Fang, Yixin Yang and Xionghou Liu
J. Mar. Sci. Eng. 2022, 10(12), 1922; https://doi.org/10.3390/jmse10121922 - 6 Dec 2022
Viewed by 1212
Abstract
Localization of emitting sources is a fundamental task in sonar applications. One of the most important factors that affect the localization performance is the sensor–target geometry. The sensor formation problem is usually addressed in related work assuming that the target is static and [...] Read more.
Localization of emitting sources is a fundamental task in sonar applications. One of the most important factors that affect the localization performance is the sensor–target geometry. The sensor formation problem is usually addressed in related work assuming that the target is static and the location is known to a certain degree, but this is not the case for many underwater surveillance problems. In this paper, we deal with the target–receiver formation problem from a different perspective, and propose to investigate the effect of target–receiver geometry on localization performance by exploiting the spatial spectrum of the direct position determination (DPD) methods. For a given multi-array system, the transformation of geometrical patterns can be explicitly demonstrated as the target moves along the track. Meaningful characteristics of the DPD methods are obtained from the experimental results, where coherent and non-coherent bearing information is used and compared. The feasibility of the DPD approaches in the ocean environments is also investigated by comparing with a matched filter processing (MFP)-based multi-array processor in order to validate the credibility of the results in this paper. Full article
(This article belongs to the Special Issue Application of Sensing and Machine Learning to Underwater Acoustic)
Show Figures

Figure 1

Figure 1
<p>(Color online) Geometric plot of the 27-element subsets of the two horizontal line arrays and the source track for Event S5. Filled green and yellow circles indicate the north large array (HLA North) and the south large array (HLA South), respectively. The four small-aperture arrays 1 to 4 are indicated by element indices. Dashed black line connects the center of the two HLAs. Solid black line indicates the source track for the first 30 min. Filled blue circles indicate true source positions for the 19 min data processed, and filled red circles mark several key points along the track. Dashed box indicates the 2D search area.</p>
Full article ">Figure 2
<p>(Color online) Localization results of GMA and N-MUSIC with 198 Hz for the 20 segments processed. (<b>a</b>) Position estimates of GMA (red cross symbols) and N-MUSIC (filled blue circles). Solid black line indicates true source track. (<b>b</b>) Range estimates of GMA to HLA North (red cross symbols) and HLA South (dotted red line with filled circles), and range estimates of N-MUSIC to HLA North (blue cross symbols) and HLA South (dotted blue line with filled circles). Solid green line and yellow line indicate true range to HLA North and HLA South, respectively. (<b>c</b>) Localization errors of GMA (solid red line with red circles) and N-MUSIC (solid blue line with plus symbols).</p>
Full article ">Figure 3
<p>(Color online) Localization results of GMA and N-MUSIC with nine tone frequencies at segment 2. (<b>a</b>,<b>b</b>) Spatial spectra of GMA with 127 Hz and 145 Hz. (<b>d</b>,<b>e</b>) Spatial spectra of N-MUSIC with 127 Hz and 145 Hz. True source position is indicated by red cross symbol inside the square. Estimated source position is indicated by white circle. (<b>c</b>,<b>f</b>) Position estimates of GMA (red symbols) and N-MUSIC (blue symbols) with nine shallow tone frequencies. (<b>g</b>) Localization errors of GMA (solid red line with red circles) and N-MUSIC (solid blue line with plus symbols).</p>
Full article ">Figure 4
<p>(Color online) 3D spatial spectrum of GMA (<b>a</b>) and N-MUSIC (<b>b</b>) with 127 Hz at segment 2.</p>
Full article ">Figure 5
<p>(Color online) (<b>a</b>–<b>c</b>) Spatial spectra of GMA at segments 3, 4, and 5, respectively. (<b>d</b>–<b>f</b>) Spatial spectra of N-MUSIC at segments 3, 4, and 5, respectively. (<b>g</b>,<b>h</b>) Spatial spectra of GMA and N-MUSIC at 17.5 s past segment 4. Results are with 232 Hz. True source position is indicated by red cross symbol inside the square. Estimated source position is indicated by white circle.</p>
Full article ">Figure 6
<p>(Color online) Localization errors with nine tone frequencies at (<b>a</b>) segment 3, (<b>b</b>) segment 4, (<b>c</b>) segment 5, and (<b>d</b>) 17.5 s past segment 4 of GMA (solid red line with red circles) and N-MUSIC (solid blue line with plus symbols). Note that the estimation error of GMA at segment 3 with 109 Hz is 1681 m and is indicated by the filled circle in (<b>a</b>).</p>
Full article ">Figure 7
<p>(Color online) Demonstration of typical behavior of GMA and N-MUSIC in section 2. Results are with 198 Hz at segment 8 (first row) and segment 14 (second row), respectively. (<b>a</b>,<b>e</b>) Spatial spectra of GMA. (<b>b</b>,<b>f</b>) Spatial spectra of N-MUSIC. (<b>c</b>,<b>g</b>) Position estimates of GMA (red symbols) and N-MUSIC (blue symbols). True source position is indicated by red cross symbol inside the square. (<b>d</b>,<b>h</b>) Estimation results of GMA (solid red line with red circles) and N-MUSIC (solid blue line with plus symbols).</p>
Full article ">Figure 8
<p>(Color online) 3D spatial spectra of GMA (<b>a</b>) and N-MUSIC (<b>b</b>) with 198 Hz at segment 8.</p>
Full article ">Figure 9
<p>(Color online) Localization results of GMA and N-MUSIC at segment 18 (first row) and segment 20 (second row). (<b>a</b>,<b>d</b>) Localization errors with nine tone frequencies of GMA (solid red line with red circles) and N-MUSIC (solid blue line with plus symbols). (<b>b</b>,<b>e</b>) Position estimates of GMA (red symbols) and N-MUSIC (blue symbols). True source position is indicated by red cross symbol inside the square. (<b>c</b>,<b>f</b>) Spatial spectra of N-MUSIC with 280 Hz. Note that the estimation error of GMA and N-MUSIC at segment 20 with 109 Hz are 1728 m and 1750 m. They are indicated by filled red circle and filled blue circle in (<b>d</b>) respectively.</p>
Full article ">Figure 10
<p>(Color online) (<b>a</b>–<b>c</b>) Position estimates of GMA (red symbols), N-MUSIC (blue symbols) and the relative-amplitude MFP processor (green symbols) for <a href="#sec1-jmse-10-01922" class="html-sec">Section 1</a>, <a href="#sec2-jmse-10-01922" class="html-sec">Section 2</a> and <a href="#sec3-jmse-10-01922" class="html-sec">Section 3</a>, respectively. Solid black line indicates true source track in corresponding section. (<b>d</b>) Localization errors of GMA (red symbols), N-MUSIC (blue symbols) and the relative-amplitude MFP processor (green symbols). Results are with 198 Hz. Filled green symbols in (<b>c</b>) indicate the estimate with <span class="html-italic">x</span> value beyond this area. Filled green symbols in (<b>d</b>) on the dotted lines correspond to the estimates with error larger than 1 km.</p>
Full article ">
Back to TopTop