Figure 1.
Schematic diagram of STMF imaging algorithm flow.
Figure 1.
Schematic diagram of STMF imaging algorithm flow.
Figure 2.
Schematic diagram of the data preprocessing steps. The information derived by the system is a two-dimensional array, and the pixels are not arranged by the target pixel position. We first reconstruct the order of the echo pixel data according to the system scanning law, and then superimpose the third-dimension data of all pixels containing the echo information of the current pixel. Since the reflection of the target is the most intense, the moment with the largest number of echo photons is the approximate distance of the target determined by all pixels “voting”. This is similar to the clustering rule.
Figure 2.
Schematic diagram of the data preprocessing steps. The information derived by the system is a two-dimensional array, and the pixels are not arranged by the target pixel position. We first reconstruct the order of the echo pixel data according to the system scanning law, and then superimpose the third-dimension data of all pixels containing the echo information of the current pixel. Since the reflection of the target is the most intense, the moment with the largest number of echo photons is the approximate distance of the target determined by all pixels “voting”. This is similar to the clustering rule.
Figure 3.
Pixel value selection diagram. Select the group with the largest number of photons on the third-dimension time histogram of each pixel, and find the peak corresponding to the bin moment and the total number of photons in the group. The range of the group is based on the accumulation law of the reflected photons of the target and the empirical value after many experiments. Select 5 bin ranges to meet the measurement of various targets.
Figure 3.
Pixel value selection diagram. Select the group with the largest number of photons on the third-dimension time histogram of each pixel, and find the peak corresponding to the bin moment and the total number of photons in the group. The range of the group is based on the accumulation law of the reflected photons of the target and the empirical value after many experiments. Select 5 bin ranges to meet the measurement of various targets.
Figure 4.
Pixel replacement strategy in STMF imaging algorithm.
Figure 4.
Pixel replacement strategy in STMF imaging algorithm.
Figure 5.
Schematic diagram of pixel value filtering and replacement. To scan the edge pixels as well, a layer of 0 values is extended without introducing added value. The scanning window scans horizontally from left to right in step 1 and pixels are replaced according to Equation (6), then scanned for the next row until the scan is completed.
Figure 5.
Schematic diagram of pixel value filtering and replacement. To scan the edge pixels as well, a layer of 0 values is extended without introducing added value. The scanning window scans horizontally from left to right in step 1 and pixels are replaced according to Equation (6), then scanned for the next row until the scan is completed.
Figure 6.
Schematic diagram of the total variational smoothing constraint flow. The total variation regularized objective function plays a key role in the optimization of depth maps and reflection intensity maps, enhancing the block-ness of the image by minimizing the total variation while maintaining edge sharpness.
Figure 6.
Schematic diagram of the total variational smoothing constraint flow. The total variation regularized objective function plays a key role in the optimization of depth maps and reflection intensity maps, enhancing the block-ness of the image by minimizing the total variation while maintaining edge sharpness.
Figure 7.
Single-photon lidar working schematic. 1: Upper computer, 2: Heat exchanger, 3: Beam expander, 4: Open-aperture reflector, 5: Galvanometer, 6: Beam and target, 7: Radar window glass, 8: Lens, 9: Filter.
Figure 7.
Single-photon lidar working schematic. 1: Upper computer, 2: Heat exchanger, 3: Beam expander, 4: Open-aperture reflector, 5: Galvanometer, 6: Beam and target, 7: Radar window glass, 8: Lens, 9: Filter.
Figure 8.
Schematic diagram of regular target plate and dimensions (Unit/mm).
Figure 8.
Schematic diagram of regular target plate and dimensions (Unit/mm).
Figure 9.
Schematic diagram of complex geometry target and dimensions (Unit/mm).
Figure 9.
Schematic diagram of complex geometry target and dimensions (Unit/mm).
Figure 10.
Schematic diagram of multi-depth parameters target and dimensions (Unit/mm).
Figure 10.
Schematic diagram of multi-depth parameters target and dimensions (Unit/mm).
Figure 11.
Depth (top) and reflection intensity (bottom) maps of the regular target plate at pulse number 5, first-photon group imaging algorithm (left column, μ = 2), STMF imaging algorithm (right column).
Figure 11.
Depth (top) and reflection intensity (bottom) maps of the regular target plate at pulse number 5, first-photon group imaging algorithm (left column, μ = 2), STMF imaging algorithm (right column).
Figure 12.
Depth (top) and reflection intensity (bottom) maps of complex geometry target at pulse number 5, first-photon group imaging algorithm (left column, μ = 2), STMF imaging algorithm (right column).
Figure 12.
Depth (top) and reflection intensity (bottom) maps of complex geometry target at pulse number 5, first-photon group imaging algorithm (left column, μ = 2), STMF imaging algorithm (right column).
Figure 13.
Depth (top) maps and reflection intensity (bottom) maps of multi-depth parametric targets at pulse number 5, first-photon group imaging algorithm (left column, μ = 2), STMF imaging algorithm (right column).
Figure 13.
Depth (top) maps and reflection intensity (bottom) maps of multi-depth parametric targets at pulse number 5, first-photon group imaging algorithm (left column, μ = 2), STMF imaging algorithm (right column).
Figure 14.
High-resolution detailed display of the imaged 5000-pulse map as a reference map. The three targets are imaged separately. The first row shows the depth maps and the second row shows the reflection intensity maps.
Figure 14.
High-resolution detailed display of the imaged 5000-pulse map as a reference map. The three targets are imaged separately. The first row shows the depth maps and the second row shows the reflection intensity maps.
Figure 15.
Depth maps of the three targets after adjusting the pixel display range at a single pixel pulse count of 5 through the STMF imaging algorithm.
Figure 15.
Depth maps of the three targets after adjusting the pixel display range at a single pixel pulse count of 5 through the STMF imaging algorithm.
Figure 16.
Reference Depth Image (pulse number 5000, top left). First-photon group imaging algorithm with parameters μ = 1 (top center), μ = 2 (top right), μ = 3 (bottom left), μ = 4 (bottom center), and μ = 5 (bottom right) for depth imaging results; the tests (MSE, RMSE, SSIM) are (2233.168382, 47.256411, 0.720168) (1944.474829, 44.096200, 0.719074) (3149.886326, 56.123848, 0.706244) (3636.324166, 60.301942, 0.723175) (3733.789040, 61.104738, 0.726569).
Figure 16.
Reference Depth Image (pulse number 5000, top left). First-photon group imaging algorithm with parameters μ = 1 (top center), μ = 2 (top right), μ = 3 (bottom left), μ = 4 (bottom center), and μ = 5 (bottom right) for depth imaging results; the tests (MSE, RMSE, SSIM) are (2233.168382, 47.256411, 0.720168) (1944.474829, 44.096200, 0.719074) (3149.886326, 56.123848, 0.706244) (3636.324166, 60.301942, 0.723175) (3733.789040, 61.104738, 0.726569).
Figure 17.
Reference reflection intensity image (pulse number 5000, top left). First-photon group imaging algorithm with parameters μ = 1 (top center), μ = 2 (top right), μ = 3 (bottom left), μ = 4 (bottom center), μ = 5 (bottom right) for reflected intensity imaging results; the tests (MSE, RMSE, SSIM) are (2940.371267, 54.225190, 0.636347) (4291.454902, 65.509197, 0.550790) (5879.315111, 76.676692, 0.568072) (6842.660833, 82.720377, 0.577880) (7041.941686, 83.916278, 0.578380).
Figure 17.
Reference reflection intensity image (pulse number 5000, top left). First-photon group imaging algorithm with parameters μ = 1 (top center), μ = 2 (top right), μ = 3 (bottom left), μ = 4 (bottom center), μ = 5 (bottom right) for reflected intensity imaging results; the tests (MSE, RMSE, SSIM) are (2940.371267, 54.225190, 0.636347) (4291.454902, 65.509197, 0.550790) (5879.315111, 76.676692, 0.568072) (6842.660833, 82.720377, 0.577880) (7041.941686, 83.916278, 0.578380).
Figure 18.
Trend plot of MSE, RMSE, and SSIM evaluation values of depth map with parameter μ values by using the first-photon imaging algorithm.
Figure 18.
Trend plot of MSE, RMSE, and SSIM evaluation values of depth map with parameter μ values by using the first-photon imaging algorithm.
Figure 19.
Trend plot of MSE, RMSE, and SSIM evaluation values of reflection intensity map with parameter μ values by using the first-photon imaging algorithm.
Figure 19.
Trend plot of MSE, RMSE, and SSIM evaluation values of reflection intensity map with parameter μ values by using the first-photon imaging algorithm.
Figure 20.
Regular target maximum group estimation depth imaging map (left), depth optimization map after pixel value screening replacement + total variance smoothing constraints (middle), depth map of 5000-pulse STMF algorithm (right).
Figure 20.
Regular target maximum group estimation depth imaging map (left), depth optimization map after pixel value screening replacement + total variance smoothing constraints (middle), depth map of 5000-pulse STMF algorithm (right).
Figure 21.
Complex geometry target maximum group estimation depth imaging map (left), optimization map after pixel value screening replacement + total variance smoothing constraints (middle), depth map of 5000-pulse STMF algorithm (right).
Figure 21.
Complex geometry target maximum group estimation depth imaging map (left), optimization map after pixel value screening replacement + total variance smoothing constraints (middle), depth map of 5000-pulse STMF algorithm (right).
Figure 22.
Multi-depth parameters target maximum group estimation depth imaging map (left), optimization map after pixel value screening replacement + total variance smoothing constraints (middle), depth map of 5000-pulse STMF algorithm (right).
Figure 22.
Multi-depth parameters target maximum group estimation depth imaging map (left), optimization map after pixel value screening replacement + total variance smoothing constraints (middle), depth map of 5000-pulse STMF algorithm (right).
Figure 23.
Curves of MSE, RMSE, and SSIM values of the imaging results varying with external light intensity.
Figure 23.
Curves of MSE, RMSE, and SSIM values of the imaging results varying with external light intensity.
Figure 24.
Curves of MSE, RMSE, and SSIM values of the imaging results varying with scene distance.
Figure 24.
Curves of MSE, RMSE, and SSIM values of the imaging results varying with scene distance.
Figure 25.
Curves of MSE, RMSE, and SSIM values of the imaging results varying with the number of pulses.
Figure 25.
Curves of MSE, RMSE, and SSIM values of the imaging results varying with the number of pulses.
Figure 26.
Underwater environment and target placement. (a) Underwater photography, (b) Target placement in water, (c) Small block target.
Figure 26.
Underwater environment and target placement. (a) Underwater photography, (b) Target placement in water, (c) Small block target.
Figure 27.
Depth imaging map of underwater multi-depth parameters target with pulse number 50 (Left picture is before noise reduction, right picture is after noise reduction).
Figure 27.
Depth imaging map of underwater multi-depth parameters target with pulse number 50 (Left picture is before noise reduction, right picture is after noise reduction).
Figure 28.
Reflection intensity imaging map of underwater multi-depth parameters target with pulse number 50 (Left picture is before noise reduction, right picture is after noise reduction).
Figure 28.
Reflection intensity imaging map of underwater multi-depth parameters target with pulse number 50 (Left picture is before noise reduction, right picture is after noise reduction).
Figure 29.
Depth imaging map of underwater multi-depth parameters target with pulse number 500 (Left picture is before noise reduction, right picture is after noise reduction).
Figure 29.
Depth imaging map of underwater multi-depth parameters target with pulse number 500 (Left picture is before noise reduction, right picture is after noise reduction).
Figure 30.
Reflection intensity imaging map of underwater multi-depth parameters target with pulse number 500 (Left picture is before noise reduction, right picture is after noise reduction).
Figure 30.
Reflection intensity imaging map of underwater multi-depth parameters target with pulse number 500 (Left picture is before noise reduction, right picture is after noise reduction).
Figure 31.
Depth imaging map of underwater multi-depth parameters target with pulse number 50 (Left picture is before noise reduction, right picture is after noise reduction, μ = 4).
Figure 31.
Depth imaging map of underwater multi-depth parameters target with pulse number 50 (Left picture is before noise reduction, right picture is after noise reduction, μ = 4).
Figure 32.
Reflection intensity imaging map of underwater multi-depth parameters target with pulse number 50 (Left picture is before noise reduction, right picture is after noise reduction, μ = 4).
Figure 32.
Reflection intensity imaging map of underwater multi-depth parameters target with pulse number 50 (Left picture is before noise reduction, right picture is after noise reduction, μ = 4).
Figure 33.
Depth imaging map of underwater multi-depth parameters target with pulse number 500 (Left picture is before noise reduction, right picture is after noise reduction, μ = 29).
Figure 33.
Depth imaging map of underwater multi-depth parameters target with pulse number 500 (Left picture is before noise reduction, right picture is after noise reduction, μ = 29).
Figure 34.
Reflection intensity imaging map of underwater multi-depth parameters target with pulse number 500 (Left picture is before noise reduction, right picture is after noise reduction, μ = 29).
Figure 34.
Reflection intensity imaging map of underwater multi-depth parameters target with pulse number 500 (Left picture is before noise reduction, right picture is after noise reduction, μ = 29).
Table 1.
Comparison of MSE, RMSE, and SSIM evaluation values of depth maps of first-photon group and STMF imaging algorithms.
Table 1.
Comparison of MSE, RMSE, and SSIM evaluation values of depth maps of first-photon group and STMF imaging algorithms.
Index | MSE | RMSE | SSIM |
---|
Methods |
---|
First-photon group | I | 1944.474829 | 44.096200 | 0.719074 |
II | 2597.687134 | 50.967511 | 0.677240 |
III | 4475.418699 | 66.898570 | 0.672388 |
STMF | I | 1877.190681 | 43.326559 | 0.799389 |
II | 3072.724840 | 55.432164 | 0.760533 |
III | 1508.434584 | 38.838571 | 0.786565 |
Table 2.
Comparison of MSE, RMSE, and SSIM evaluation values of reflected intensity maps for first-photon group and STMF imaging algorithms.
Table 2.
Comparison of MSE, RMSE, and SSIM evaluation values of reflected intensity maps for first-photon group and STMF imaging algorithms.
Index | MSE | RMSE | SSIM |
---|
Methods |
---|
First-photon group | I | 4291.454902 | 65.509197 | 0.550790 |
II | 3478.640814 | 58.980004 | 0.677963 |
III | 4152.372382 | 64.438904 | 0.627105 |
STMF | I | 3064.493265 | 55.357865 | 0.577002 |
II | 2620.730375 | 51.193070 | 0.625744 |
III | 4482.540446 | 66.951777 | 0.567506 |
Table 3.
Evaluation index values of three types of targets before and after the STMF noise reduction step.
Table 3.
Evaluation index values of three types of targets before and after the STMF noise reduction step.
Index | MSE | RMSE | SSIM |
---|
Methods |
---|
STMF | I (before noise reduction) | 1236.794840 | 35.168094 | 0.789318 |
I (after noise reduction) | 937.125737 | 30.612509 | 0.912257 |
II (before noise reduction) | 878.321448 | 29.636488 | 0.749681 |
II (after noise reduction) | 601.322533 | 24.521879 | 0.897178 |
III (before noise reduction) | 1448.678334 | 38.061507 | 0.743624 |
III (after noise reduction) | 402.982322 | 20.074420 | 0.910434 |
Table 4.
Depth imaging MSE, RMSE, and SSIM values with external light intensity (some example points are shown).
Table 4.
Depth imaging MSE, RMSE, and SSIM values with external light intensity (some example points are shown).
Index | MSE | RMSE | SSIM |
---|
Intensity/Lux |
---|
0.1 | 1805.239107 | 42.488105 | 0.862387 |
16146 | 2094.045478 | 45.760742 | 0.865294 |
50828 | 3542.234712 | 59.516676 | 0.747120 |
Table 5.
Depth imaging MSE, RMSE, and SSIM values a with scene distance (some example points are shown).
Table 5.
Depth imaging MSE, RMSE, and SSIM values a with scene distance (some example points are shown).
Index | MSE | RMSE | SSIM |
---|
Distance/m |
---|
25 | 1895.817939 | 43.540991 | 0.869116 |
40 | 2094.045478 | 45.760742 | 0.865294 |
55 | 2873.154976 | 53.601819 | 0.837796 |
Table 6.
Depth imaging MSE, RMSE, and SSIM values with number of pulses (some example points are shown).
Table 6.
Depth imaging MSE, RMSE, and SSIM values with number of pulses (some example points are shown).
Index | MSE | RMSE | SSIM |
---|
Pulse Number |
---|
5 | 1888.402562 | 43.455754 | 0.874237 |
20 | 1350.848840 | 36.753896 | 0.903139 |
50 | 1059.433603 | 32.548942 | 0.926321 |
Table 7.
MSE, RMSE, and SSIM evaluation metrics of underwater depth map (STMF).
Table 7.
MSE, RMSE, and SSIM evaluation metrics of underwater depth map (STMF).
Index | MSE | RMSE | SSIM |
---|
Pulse Number |
---|
50 | before noise reduction | 4967.051578 | 70.477313 | 0.656152 |
after noise reduction | 5162.811476 | 71.852707 | 0.689871 |
500 | before noise reduction | 3974.910685 | 63.046893 | 0.685066 |
after noise reduction | 3978.879936 | 63.078363 | 0.716883 |
Table 8.
MSE, RMSE, and SSIM evaluation metrics of underwater depth map (first-photon group).
Table 8.
MSE, RMSE, and SSIM evaluation metrics of underwater depth map (first-photon group).
Index | MSE | RMSE | SSIM |
---|
Pulse Number |
---|
50 | before noise reduction | 9343.743232 | 96.663040 | 0.611259 |
after noise reduction | 10,459.830375 | 102.273312 | 0.615554 |
500 | before noise reduction | 10,699.728375 | 103.439491 | 0.631418 |
after noise reduction | 10,820.452469 | 104.021404 | 0.617580 |