High Accuracy Reconstruction of Airborne Streak Tube Imaging LiDAR Using Particle Swarm Optimization

: Airborne streak tube imaging LiDAR (STIL) consists of several different data-generating subsystems and introduces system errors each time it is installed on an aircraft. These errors change with each installation, which makes the parametric calibration of the LiDAR meaningless. In this study, we propose a high-precision reconstruction method for point clouds that can be used without calibrating the system parameters. In essence, after each remote sensing measurement, a self-checking process is performed with experimental data to replace the fixed system parameters. In this process, the splicing error of the same region measured under different conditions is used as a criterion to optimize the reconstruction parameters via a particle swarm optimization (PSO) algorithm. For a detection distance of 3000 m, the elevation error of the point cloud reconstruction reaches more than 1 m if the placement parameters are not optimized; after optimization, the elevation error can be controlled within 0.3 m.


Introduction
Airborne LiDAR is an active remote sensing method that has the benefits of 24 h operation, high penetration, precise range finding, and a quick turnaround.These characteristics have made it very popular for surveying a variety of regions, including urban areas [1,2], powerlines [3,4], forestry resources [5,6], glaciers and snow [7,8], and other areas [9][10][11].
We have developed an STIL as an airborne laser system (ALS) to assess the local landscape and structures at a working distance of 3000 m in Zhaodong, China.According to the mapping mission requirements, the altitude accuracy should be better than 0.3 m.However, the error is more sensitive as the operating distance increases.For instance, a 0.02-degree deviation can cause a one-meter offset.
In the reconstruction process, there are three categories of error: computational, random, and systematic [21].During our data processing, the computational errors were consistently smaller than 10 −6 m.The random error for an ALS was also less than 0.2 m within a range of 30 • scan angles and at the operational height of 3000 m [22].Both error effects are compounded accuracy requirements.The primary element affecting the mapping mission goal is systematic error.
If installation errors were not taken into account, the elevation accuracy of the point cloud would be lower than one meter during data processing.This falls far short of the predicted outcome.
There are two essential steps to modifying the effect of installations.First, the whole point cloud generated by a single experiment is too large to be processed within a limited time, so we have to select some representative points for optimization.There are generally two options for picking points: one is from the overlapping bands [23], and the other requires layout homonymy points on the ground [24].
On the one hand, it is impossible to find a one-to-one correspondence on the ground from the rough point cloud that has not been checked and calibrated; on the other hand, measuring and averaging all points in an area is too time-consuming.Given these two factors, we decided to calibrate the system using data from overlapping bands.
However, the whole set of data from a flight test is too large to be processed in a given period of time.Thus, we need to match the same regions for different point clouds using target features and average them as control points to check parameters.
The second step is solving all installation errors using the points we selected.The most classic algorithm is the least square method (LSM) for this problem [25].Unfortunately, this method did not work for our airborne STIL.We have to use nonlinear optimization with global search capabilities.
Particle swarm optimization (PSO) is one of the best-known global search algorithms [26]; it is widely used in many fields, such as path planning [27] and neural network training [28].
This algorithm, based on particle competition, is insensitive to initial values and monotonicity and thus has a strong global solution search capability.These features make it ideal for parametric calibration [29][30][31][32][33].
Based on PSO, we finally obtained multiple sets of checking parameters, reducing the elevation error from more than 1 m to 0.3 m.

Streak Tube Imaging LiDAR
The STIL workflow is depicted in Figure 1a.The laser is activated by the sequential pulse and is received by a streak tube image sensor.In this process, four subsystems generate data for reconstruction.The streak tube image detector records the time when echo photons come back, the scanning system provides the scanning angle using a grating ruler, the positioning and orientation system (POS) records the position and posture information at the same time, and the pulse sequential system triggers the three subsystems to start with a delay after the laser beam is emitted.The STIL works as an ALS, which is operated as shown in Figure 1b.It is located by the POS, which integrates a GPS and an inertial navigation system (INS).The POS calculates the GPS signals from satellites to its own position and obtains posture information from the INS.
In addition, the scanning mirror swings back along the perpendicular direction of flight, and the laser footprint is formed on the ground, as shown in Figure 1c, where the green lines represent the laser streaks, and red indicates the scan trajectory.We designed different scanning speeds for forward and reverse directions and only used the forward scanning data.There are some overlapping areas between valid scans, which play an important role when checking the parameters.
The POS, delay, and angle data can be calculated for reconstruction, but the streak image should be converted to distance and field angle.The working principle of the streak tube detector is shown in Figure 2. Currently, the impulse delay, incident time, incidence angle, scanning angle, and POS data may all be used to reconstruct the point cloud.
However, there are inclinations among the scanner, the streak tube detector, and the POS during installation.Their installation mistakes make the point cloud reconstruction more challenging and less accurate.
For low-altitude ALS, installation errors have a low impact on the inversion accuracy, such that a rough measurement is sufficient to ensure that the effect of systematic errors is smaller than the effect of random errors.Unfortunately, as the working height increases, the installation errors are amplified during the point cloud-inversion process.The placement parameters obtained from ground measurements cannot meet the accuracy requirements of point cloud reconstruction.
Therefore, we had to check the system parameters with experimental flight data.These parameters are included in the coordinate transformation matrixes, under which the raw data are reconstructed into a point cloud.This is clearly a global nonlinear problem.

Particle Swarm Optimization
The PSO algorithm is a computational optimization technique inspired by the social behavior of birds flocking or fish schooling.It was first introduced by James Kennedy and Russell Eberhart in 1995 [26].
Figure 3 shows an overview of how the PSO works.Firstly, it initializes a swarm of particles with random positions and velocities within the search space.Secondly, it evaluates each particle's position and determines its fitness value based on the objective function.Then, it starts a cycle that updates the personal best and global best, as well as the velocity if the best solution is worse than the cut-off point.

Reconstruction Process
To convert from a laser signal to a real point in the WGS-84 system, the target point is determined using several coordinate transitions.The process of point cloud reconstruction can be abstracted as coordinate transformations of a vector.According to Figure 1, the streak signal undergoes a reflection as it travels from the target to the streak tube detector.When data such as stripe signals and POS are synthesized, the rotation of coordinates is inevitably involved due to placement errors.
A vector should be multiplied by a rotation matrix while being transformed between different Cartesian coordinates.The three-dimensional rotation matrix can be expressed by Euler angles as R(A, B, C), whose matrix form is shown as A vector should be multiplied by a reflection matrix while reflected by a mirror.The three-dimensional reflection matrix can be calculated by the mirror norm vector u.We define M as the mirror reflection matrix, which can be expressed by We define the normalized vector of incident light as L S , the photon travel distance as d, and the whole transformation from this vector to the point in world frame as This equation contains one reflection and three rotations: R S L is the rotation matrix between the sensor and laser footprint, R L P is the rotation matrix from the POS to the footprint, and the matrix R P W expresses the rotation from the POS frame to the real world.
The parameter P p is the POS position calculated from the GPS and IMU.The specific form of the position is P = x p y p z p T .(4)

The Sensor Coordinate
The first-hand data are located on the pixel coordinate of the streak images shown in Figure 4a.We need to extract the location information from it using the center-of-mass method.The centroid in row i of the image can be calculated using the following formula: where I ij is the grayscale value of the position (i, j), N is the column count, and the centroid array {x i } is shown in Figure 4b.We rewrite x i in the coordinate form as (x, y) where the variable y is just the row number i.After calculating the gray centroid horizontally, the image is converted to points with pixel units.It is the first coordinate transition that maps pixel points from an image to a distance-angle coordinate using calibration matrixes [34].Then, the y-axis data are mapped to the incident angle θ, and the x-axis is converted to the time of flight (TOF) t s .
The calibration matrixes whose pseudo-color figure is shown in Figure 4c,d are datamapping tables between the (x, y) pixel frame and the (t s , θ) system.We denote T and Θ as the matrixes of the time delay and incident angle, respectively.
The subscripts of the elements in the matrix are integers, while the centroid coordinates in the x-direction are often fractional, and their corresponding delay parameters need to be obtained via interpolation.We define x L = ⌊x⌋ as the floor of x and x R = ⌈x⌉ as the ceiling of x.The correspondence from the image coordinate system to the sensor coordinate system can be expressed as Its accuracy meets the inversion requirements, where one pixel corresponds to a maximum of 0.004 degrees in the spatially calibrated matrix and 1.1 nanoseconds in the temporal calibration matrix.After using the center-of-mass method and interpolation of the calibrated matrix, the accuracy can reach 0.1 pixel.Thus, the (t s , θ) data calibrated from the two matrixes are 0.0004 • and 0.11 nanoseconds.Both cause an error of less than 0.02 m, which is negligible for a plotting accuracy of 0.3 m.
In the sensor frame (t s , θ), t denotes the time between when the camera starts working and when it receives the light signal, and θ is the incident angle.Then, the normalized vector of incident light can be expressed by θ as follows: t s is just one part of TOF, and another part of the TOF t d is the delay between the time when the laser and camera were triggered.Thus, the whole TOF t = t S + t d and the distance measurement value d can be shown as where c is the speed of light.At this point, the first coordinate transition was completed, and the next goal was the position in the scanning coordinate system.The errors introduced by this transition from the streak image and time delay to the incident vector L s and ranging value d are less than the accuracy requirements for mapping, so in the subsequent calculations, we just used the vector L s and the distance d as a definite value.

The Scanning Coordinate
Laser streaks on the target object were reflected using a scanner mirror and incident to the sensor shown in Figure 5, so the distance and angle inversed using the calibration matrixes represent the image point P i instead of the target P t .Ideally, there would be a simple relationship between P i and P t .However, the scanning motor and image sensor are not coaxial, which leads to rotation between the two axes systems.In addition, the normal oscillating mirror and the spin axis are not at a 45 • angle, leading to a more complex reflection between image and object systems.
We define the Euler angles between the sensor and mirror as r, s, and t and the normal space angles as Φ and Ω.The parameter Φ is the angle between the normal mirror and the motor's rotational axis, and Ω is the placement error in the direction of motor rotation.
A high-precision circular grating ruler was fixed on the motor, from which we could read the mirror rotation angle.Not surprisingly, there was also an error between the grating ruler and scanning system.Thus, we should not only define the raster reading φ but also another angle Ψ to modify it.Now, the matrix R(r, s, t) can be used to describe the first rotation transformation R S L , and R(0, 0, Ψ + φ) can indicate the second.Then, the norm vector can be deduced as This step converts the coordinate from the sensor to the scanning system and contains two rotations and one reflection transformation.

The POS Coordinate
At this point, we should migrate the coordinate from the scanning system to the POS via rotation and translation.
First, we define u, v, and w as the installation angles of the POS, and the rotation matrix ) is used to transform the scanning system to a navigational system.It is so difficult to measure these three angles that we should use them in optimization.
Second, we define o, p, and q as the sensor's x, y, and z axis position in the POS coordinate system.We use this position to shift values from the sensor to the navigational system, which connects the real world with longitude, latitude, and altitude.The shift vector is simpler than those matrixes and can be written as the vector P p = [o, p, q] T .
As the last process in the local relative coordinate system, we should solve the position P with those transition matrixes and vectors.The distance d plays a key role as the length component.

The World Coordinate
The POS outputs position and posture data, both located in the real world, or rather, in the WGS-84 coordinate system.
We define α, β, and γ as the POS angles, which expresses the real POS posture to the ground plane of reference.These angles are formulated as the rotation matrix R(α, β, γ).
We use L 0 , B 0 , and H 0 to indicate the longitude, latitude, and altitude of the POS, respectively.L 0 and B 0 continue to rotate the POS axes using the following matrix: To combine with the POS coordinate, longitude and latitude should be converted from angle to length variables using the following conversion: The function f T converts the positions (L, B, H) from the WGS-84 system to an ECEF (Earth-centered, Earth-fixed point (x, y, z), where a is the Earth radius, and e is the first eccentricity in the Earth model.
With the same coordinate format, the point x p , y p , z p is converted to the position on Earth by Finally, we convert the ECEF coordinate back to the WGS-84 system using the iterative process f W (x, y, z) shown in (13)  [35], which is the standard conversion relationship between ECEF and WGS84 systems.
The end condition is that the difference between the old and new value is less than 1 × 10 −9• .The error in length due to a difference of 1 • in latitude and longitude is not greater than 111 km, so this process results in an error of less than 1 mm.

Variable Types
There are many variables in the point cloud-reconstruction process, some of which are experimental data, some of which are reconfiguration parameters that can be measured, and others of which are parameters that cannot be measured well enough.The parameters that cannot be measured need to be searched using the PSO algorithm.
In Table 1, these parameters are classified into three types: experimental data (E), measurable parameters (M), and parameters that need to be checked and calibrated (C).The roll, pitch, and head angle of the POS E L 0 , B 0 , H 0 The longitude, latitude, and altitude of the POS E In addition to the calibration matrixes, the distance between the two subsystems is a are measurable parameter, because even with the coarsest measuring methods, a measuring accuracy higher than 1 mm can be obtained.
However, the angular placement error behaves quite differently, and it varies with the working distance.The mounting angle error causes distance errors in two directions with the form of where d is the dummy ranging value in the current coordinate system, ϕ is the measured mounting angle, ϵ is the angular placement error, and ∆ is the distance error due to the angle placement.Since ∆ 1 = ∆ 2 = 0.3 m, and D = 3000 m, ϵ has the minimum value of 0.3 3000 = 1 × 10 −4 rad, while ϕ equals 0 or π 2 .This means that an angular accuracy of 10 −4 rad must be achieved; otherwise, a distance accuracy of 0.3 m cannot be reached.
Thus, the main parameters to be checked are from the placement errors between the image, scanning, and POS coordinate systems.

Point Selection
The whole point cloud is too big for an intelligence algorithm to optimize.As a result, we must choose a few checking sites for iterative calculation.First, the point cloud is reconstructed using the uncalibrated parameters, and then the composite required points are selected from it.After target point selection, the selected points need to be denoised via radius filtering.Finally, the raw data are determined one-to-one with the point cloud using subscripts and selected for self-checking.
Different areas play different roles in the self-checking process, and these points need to be selected using different methods.Target buildings have definite edges that can be boxed to an area with determinant latitudes and longitudes.Further, the flat ground's altitude slowly changes, so we can use it to adjust the vertical accuracy.
Target buildings are real buildings in a city, which means they appear randomly in the point cloud.It is difficult to find a complete one-to-one mapping building in adjacent strips, but it is easy in two air lines.Therefore, we have three sets of control points: P b : target building points between different air routes.P f : flat area points among different air routes.P a : flat area points between adjacent scanning bands.

Target Buildings
To ensure that the altitude accuracy is meaningful, it is important to have the control points aligned in the latitude and longitude directions.The role of the target buildings is to align the two sets of point clouds in the horizontal direction.
In the inspection area, the buildings suitable for a parametric check have a regular geometry and are at a certain height above the ground.A typical building suitable for checking schools is shown in Figure 6b.In order to choose the target building points, a rectangle should be boxed around the building.Then, the points higher than the altitude threshold in the rectangular area should be selected.It must be ensured that the rectangular selection box is flat except for the target building.
Target buildings are used to match the points from different air lines, so the strip splicing error should be eliminated through point-selection work.We have to select buildings that stand on one strip completely, which are always small, as shown in Figure 6.
There are three sets of targets: the first contains points generated from east to west, the second from north to south, and the third from west to east.The three airlines constitute two control groups: white numbers on the first picture are compared against the second, and black numbers are subtracted by the third.
We checked the accuracy of the latitude and longitude directions using two groups, one of which is reverse-checking and the other is crosschecking.Both of them need four vertexes selected from imprecise point cloud manually.
Depending on the point of the target building, it is possible to make a judgment about the error in the horizontal direction as where the subscript i denotes the average value of the ith set of point clouds, c1, c2, d, e correspond to the four different labeled builds in the three figures; x is the relative distance in longitude; and y is the latitude.

Flat Areas among Different Air Routes
There is a further characteristic of the flat mark point that distinguishes it from its building counterpoint.The flat object has no clear edge to locate its horizontal position, and we cannot even make sure that it is an object that contains the same area as another.However, this characteristic results from the slowly changing height, which means horizontal misalignment does not affect the altitude accuracy.
Thus, we selected a square area on streets, especially on crossroads, and made four points in every cloud.The square is more convenient than a target building because it is small enough to stay within the same scanning band.We can find the same area among the three clouds generated from different heading airlines, as shown in Figure 7.There are also three sets in Figure 7, but we painted them in different pseudo colors to show the street contour.Every subgraph has eight marked areas, which contain approximately thirty points and are averaged to just one height.We subtracted the second and third height from the first and gained two groups of altitude difference values as follows: where h a i , h b i , and h c i are the altitude averages of the points in the ith target area in Figure 7a-c.

Flat Area for Adjacent Scanning Bands
Error-checking parameters not only lead to misalignment between buildings generated from different airlines but also cause altitude differences between adjacent scanning bands.The points across two strips have a greater field angle, so their photons are incident on the sensor more marginally.Sometimes, half of them are wasted in the reconstruction process, but the continuously high trend from the middle to sides is still useful for checking the installation error.
The area is not wide enough to accommodate complete buildings on overlapped areas, so we have to select flat targets as in the second selection method.Figure 8 shows the selection positions across scanning bands beneath the east-west airline.We select a point that spans both strips and draw a circle horizontally with that point as the center and then average the heights of all the points within the circle.For the ith circle, we obtain an average value on each of the front and back strips, denoted as h f i and h b j ; the error criterion can be expressed as

Parameter Optimization Using the PSO Algorithm
There are nine parameters that should be optimized, as listed in Table 1, and we defined them as a checking vector p = (r, s, t, Ψ, Φ, Ω, u, v, w).
It is obvious that the reconstruction of STIL is a nonlinear process whose inverse problem cannot be solved using the least square method (LSM) based on the matrix form.Furthermore, it is a global optimization process, where the nonlinear LSM falls into the local optimal solution through the direction of the gradient descent.
One of the best ways to avoid local optimums is to introduce randomness, and most intelligent heuristic algorithms work with it.If we set a random step size for the nonlinear LSM, the program evolves into an intelligent algorithm.
The PSO not only uses random steps but also introduces a competition mechanism, a population-based stochastic optimization algorithm inspired by the social behavior of birds flocking or fish schooling.The PSO is commonly used to solve optimization problems where the objective is to find the optimal solution in a search space.
During the optimization of placement parameters, the checking vector p is abstracted into a particle.Since the PSO starts running, 50 checking vectors are generated randomly and comprise the particle swarm.Three types of checking points are reconstructed with every particle, and the weighted average is determined to obtain a matching error shown as e = w b e b + w f e f + w a e a . ( The average value e i is the distance between two groups of check points in point set P i , and w i is its weight.
The particle p i records its own minimum matching error as p i h , and there is a best particle q j in the jth iteration.Every particle has a v in every iteration and moves along a straight line without a matching error.However, it exists in a swarm and is burdened with the task of finding an optimal solution, so it has to move to the historical optimal solution p i h or the ethnic optimal solution q j .Thus, the particle p i j , which is the ith particle in the swarm, in its jth iteration, would have the new speed v i j+1 .
The vectors p i j , p i h , and q j are the positions of their corresponding particles; c 0 and c 1 are the global and local constants; and r 1 , r 2 are randoms.The pseudocode is shown in Algorithm 1. function pso(points, max_iter, target_precision) 2.
if test_one_parameter(points, p i h ) > e i 10.
end for 13.g = argmin i e i 14. q j = p i g 15. if e g < target_precision 16.
return q j 17. v j+1 = v j + c 0 r j q j − p j + c 1 r 1 p g − p j Fifteen optimizations were performed, and the nine checking parameters in each optimization are shown in Figure 9.The maximum angle error is approximately 5.5 • , which is impossible.These large angular errors imply coupling between the nine param-eters.Coupling but not simple merging is the biggest difficulty in checking the error of airborne STILs.
Their linear fitting diagram is shown as Figure 10.The parameter coupling phenomenon indicates a correlation before different placement angles and verifies that the particle swarm algorithm can solve the multi-parameter optimization problem.
Based on the definitions in Table 1, the coupling parameters can be found to be related in orientation.The nine parameters come from three devices, the imaging system, the scanning system and the POS.Although the orientation of each parameter was not specified during modelling, based on the optimization results, it can be seen that ϕ, s, and v represent the same orientation in all three coordinate systems.The same is true for other parameters.In addition, the scanning system was modelled with redundancy, resulting in Ψ and Ω playing the same role.
The coefficients between the parameters w and r in these three sets of couplings converge to 1, suggesting that the two may be condensed.The linear coupling coefficients for the other two sets of relationships are not 1, indicating that these parameters cannot be optimized by fitting the relationships.However, according to the model's definition, the parameters Ψ and Ω are both used to describe the error of the motor rotation angle.The grating ruler has been installed co-axially to the rotating shaft so that its zero position is virtually the same physical quantity as the motor placement error in the tangential direction of the rotating shift.
When one or more of these parameters are fixed, the parameters associated with them will also converge to fixed values during the optimization process.In traditional calibration methods, it is necessary to obtain the exact placement error.The phenomenon of parameter coupling, on the other hand, makes the number of parameters that have to be obtained by measurements much lower, which is very beneficial for the traditional checking methods.
Considering that the purpose of checking is to improve the accuracy of the point cloud, the non-uniqueness of the optimum caused by parameter coupling is not a major issue.Even the exact physical meaning of the nine parameters seems less important, as long as they are substituted into the conformal equations to be able to invert a highly accurate point cloud.

Reconstruction Accuracy
An ALS topographic mapping task is essentially a measurement of the topographic relief of an area.In its reconstructed point clouds, latitude and longitude often act as the coordinate frame, while elevation is the measured value.Therefore, we locate the proximity points in different point clouds by giving the latitude and longitude and use the elevation difference between the two as the reconstruction error.
Though three types of points are used in global optimization, the criteria for determining the error are used for elevation.
Finally, the vertical errors are around 0.3 m, different every time when solved through the PSO, and the minimum value is as low as 0.271 m, as shown in Figure 11.In contrast, whether assuming placement errors of 0 or optimizing using nonlinear LSM starting with random numbers, the elevation errors of the generated point cloud are all greater than 1 m.On the one hand, the optimization results of the placement errors have huge jitter and are coupled with each other; on the other hand, the point cloud reconstruction results are relatively stable, and their reconstruction errors are controlled within a small range.In fact, although we have been referring to the parameters to be optimized as placement errors, they are clearly not true installation angles.In other words, we have achieved highly accurate point cloud reconstruction without checking parameters.
The placement error obtained by the PSO algorithm produces a significant effect on point cloud reconstruction.As far as the horizontal direction is concerned, it solves the misalignment problem between the strips and reflects the validity of the building target criterion shown in Figure 12.The difference before and after PSO is more intuitively shown in Figure 13, where a clear misalignment before the nonlinear optimization is shown in Figure 13b.The PSO has smoothed this issue out, as shown in Figure 13c.

Conclusions
In this study, we propose a high-accuracy reconstruction method for airborne STIL.Three types of points are selected from different airlines or scanning bands, and iterated as weighted variables using the PSO algorithm.The optimization results show that the elevation accuracy is always maintained in a small range despite the complex coupling relationship between different parameters.This means that high-precision reconstruction of 3D point clouds can be achieved without the need to measure placement errors of the LiDAR system during mapping missions.Among several optimization parameters, the best optimization setting increases the elevation accuracy to within 0.3 m.

Figure 1 .
Figure 1.STIL working mechanism; (a) system compositions and workflow of STILs; (b) work status of an ALS; (c) laser line footprints.

Figure 2 .
Figure 2. Schematic of the STIL system.(a) Schematic diagram of the data collection process.(b) The working principle of the streak array detector.(c) Streak image on CCD.

Figure 3 .
Figure 3. Diagram of the PSO algorithm.

Figure 4 .
Figure 4. Factors of conversion from a pixel coordinate to a distance-angle coordinate: (a) is the original streak image, (b) is the centroid of the streak, (c) is the calibration matrix T mapping the time delay, and (d) is the calibration matrix Θ mapping the incident angle.

Figure 6 .
Figure 6.Target buildings: (a) top view of target buildings; (b) side view of the building in red in (a); (c-e) selected buildings, with (c) showing the airline from east to west, (d) from north to south, and (e) from west to east.

Figure 7 .
Figure 7. Target flat areas from different air routes; (a) the east-west airline; (b) the south-north airline; (c) the west-east airline.

Figure 8 .
Figure 8.(a) Marked positions across strips; (b) diagram of the strip overlap.

Figure 9 .
Figure 9.The distribution of each parameter after multiple PSO optimizations.In this boxplot, the boxes represent the range of data from the first to the third quartile.There are three directions in the rotation operation, corresponding to the three sets of coupling relationships of the checking parameters, shown as      w= 0.9979r − 0.0355 Φ= 0.5632(s + v) + 45.228 Ψ + Ω= −0.4965(t − u) + 0.5391 .(20)

Figure 10 .
Figure 10.The coupling diagram of placement parameters.

Figure 11 .
Figure 11.The altitude errors in fifteen PSO solutions.

Figure 12 .
Figure 12.A comparison of misalignment with and without PSO parameters; (a) before PSO; (b) after PSO.

Figure 13 .
Figure 13.Coincidence of point clouds under different air routes; (a) the test zone; (b) before PSO; (c) after PSO.

Table 1 .
Variables in the reconstruction process.