# Camera Calibration for Long-Distance Photogrammetry Using Unmanned Aerial Vehicles

Authors: Yang Zhang; Jun Yang; Guoliang Li; Tianqing Zhao; Xiaokai Song; Suoqi Zhang; Ang Li; Hui Bian; Jin Li; Min Zhang
Journal: Journal of Sensors
 (2022)

---

## Abstract

The traditional target-dependent camera calibration method has been widely used in close-distance and small field of view scenes. However, in view of the field coordinate measurement in the large-scale monitoring area under the complex field environment, the standard target can hardly meet the requirements of covering most of the camera’s field of view. In view of the above problem, a stereo camera calibration method is studied, using the unmanned aerial vehicles (UAV) as feature points, combined with the high-precision position information measured by the real-time kinematic (RTK) positioning system it carries. The measured UAV coordinates are unified in World Geodetic System 1984 (WGS-84). Therefore, through several preset points, the measurement reference coordinate system which is the new world coordinate system we need can be established in any monitoring area, which greatly improves the flexibility of measurement. The experimental results show that the measurement accuracy of the proposed method can reach 0.5% in the monitoring area with a diameter of 100 m. The calibration method has a wide range of application and does not need the traditional standard target, and the measurement reference coordinate system can be established according to the actual needs. It is suitable for field spatial coordinate measurement in long-distance and complex terrain environment.


---

## Body

## 1. Introduction

Binocular stereo vision simulates human eye vision to realize the mapping from two-dimensional (2D) images to three-dimensional (3D) space and realizes the use of 3D information. At present, this method has been widely used in autonomous driving, robot navigation, virtual reality, and industrial production [1–6]. The process of solving this mapping relationship is called camera calibration, which involves some parameters, including both intrinsic and extrinsic parameters. Intrinsic parameters consist of principal points, focal lengths, and lens distortion. Extrinsic parameters include a rotation matrix and a translation vector between the two cameras.Various effective calibration methods have been proposed, including traditional calibration methods, self-calibration methods, and camera calibration based on active vision. In the traditional calibration method, the intrinsic and extrinsic parameters of the camera are obtained by mathematical transformation of 3D coordinates and 2D image coordinates by presetting some targets. Faig [7] proposed an imaging model based on the optimization algorithm, which has a complex solution process and high initial value requirements. Abdel-Aziz and Karara [8] proposed the direct linear transformation (DLT) method, which ignored the effects of distortion and obtained unknown parameters of the equation by solving the linear equations. Tsai [9] proposed a two-step method based on radial constraint by combining optimization algorithm and direct linear transformation method on the basis of only considering radial distortion. Zhang [10] is best known for his flexible calibration method, in which he provided a good method for estimating the initial parameters of the camera by using the constraints of the homography between planes. The premise of the application of the above methods is to manufacture specific targets, such as checkerboard or circular targets, which is difficult to achieve in the field with a large field of view due to the size limitation. In response to the above problems, Faugeras et al. [11] and Maybank and Faugeras [12, 13] proposed a camera self-calibration method, which calibrated the camera by taking multiple images with distinct features and relative motion. However, this method had great limitations in sky, desert, sea, and other environments, and its robustness was poor and data reliability was insufficient. Similarly, Ma and Zhang [14–17] proposed a calibration method based on active vision, which required the camera to make specific movement and was not suitable for the occasion when the camera was fixed in the field with a large field of view. Besides, many scholars have proposed camera calibration methods in large field of view environment. For example, Kong et al. [18] proposed a method of camera calibration based on the Global Positioning System (GPS), which directly took the GPS instrument as the feature points, which limited the flexibility of the method in practical use. Xiao et al. [19] proposed a binocular 3D measurement system that uses a cross target with ring coded points. Shang et al. [20] proposed a large field of view calibration method in which the optical center and control point of the camera are close to the coplanar, which has many limitations. Sun et al. [21] proposed a baseline-based camera calibration method in which the calibration target must be randomly placed in the field of view several times. Wang et al. [22–24] proposed a stereo calibration method for out-of-focus cameras when acquiring images for long- and short-distance photogrammetry, which has high robustness and high accuracy. None of these methods enable precise and fast camera calibration at large field of view.In this paper, a calibration method using the UAV with RTK as a high-precision mobile calibration target is proposed. This method does not need to manufacture large-scale calibration target, which reduces the requirement of calibration conditions, and is suitable for large scene field environment. In addition, by using the WGS-84 earth coordinate system as the intermediary, the measurement reference coordinate system can be flexibly converted to any desired position through several preset coordinate points, even if the position cannot be observed by the binocular cameras simultaneously, which is very suitable for some complex field scenes where the view is partially obscured by trees or hills. Experimental results show that the proposed method performs well in the monitoring area with a diameter of 50-100 m at the distance of 500-1000 m from the cameras.The subsequent compositions of this article are as follows: Section2 introduces the basic principles, Section 3 introduces the calibration process and experimental results, and Section 4 summarizes this article.


## 2. Calibration Theory

### 2.1. Camera Imaging Model

This paper focuses on where the camera is 500 m-1000 m away from the center of the monitoring area; therefore, the telephoto lens is used. Considering that the telephoto lens of the camera has very little distortion, the ideal pinhole imaging mode [25] is chosen to describe the mapping relationship between the object space and the image space, as is shown in Figure 1.According to the pinhole imaging model, the coordinates from the world coordinate system, camera coordinate system, camera physical coordinate system and image pixel coordinate system have undergone three parts of rigid body transformation, projection transformation, and rigid body transformation, respectively. One point in the camera coordinate system is expressed asXc,Yc,Zc and in the world coordinate system is expressed as Xw,Yw,Zw, which are named x,y and u,v, respectively, in the camera physical and image pixel coordinate systems. According to the relationship between each coordinate system, the linear transformation relationship between the world coordinate system Xw,Yw,Zw and the pixel coordinate system u,v is established by the following equation.
(1)Zcuv1=1dx0u001dyv0001f0000f000010RT01XwYwZw1,where R and T are the rotation matrix and translation vector between the world coordinate system and the camera coordinate system. f is the focal length of the lens. dx and dy are the physical size of the pixel. u0 and v0 are the camera principal points. Camera calibration requires the solution of these parameters.


### 2.2. Coordinate System Conversion

As can be seen from Section2.1, obtaining the correspondence between the pixel coordinates and the world coordinates of the feature point is the key to estimating the camera parameters. With the help of the UAV working in RTK mode, we can obtain the UAV’s current GPS navigation coordinates PGB,L,H, which can be converted into earth rectangular coordinates PEXE,YE,ZE [18]:
(2)XE=N+HcosBcosL,YE=N+HcosBsinL,ZE=N1−E2+HsinB,whereN is the radius of curvature of the ellipsoid and E is the first eccentricity of the ellipsoid. Let a,b be the long and short semiaxes of the Earth, respectively, and χ be the ellipsoidal flattening rate of the Earth. Without losing generality [26],
(3)a=6378137,χ=1298.257223563,b=a×1−χ.Also, we know(4)E=a2+b2a,N=a1−E2sin2B.By combining equations (2)–(4), the representation of PG in the WGS-84 earth rectangular coordinate system, PE, can be obtained. However, the coordinates obtained by the above steps are based on the Earth’s center of mass, which has two disadvantages: first, the scale of the obtained coordinates is too large to estimate the camera parameters; second, the origin and direction of the current coordinate axes have been fixed, which is not conducive to further measurement. Therefore, we need to set the origin and direction of the world coordinate system according to our own needs and complete the camera calibration in this coordinate system which is also called the preset coordinate system. Since coordinates of all measuring points are in the WGS-84 earth rectangular coordinate system, this transformation is not difficult. According to different application scenarios, the preset coordinate systems can be established by the following two ways.
(1)
Establishment of the preset coordinate system for a rectangular region of interestAs is shown in Figure2, the latitude and longitude of A, B, C, and D four points are measured at the four corners of a rectangle and converted to the earth rectangular coordinate system by the above steps. Take the coordinates of the intersection of lines AC and BD as the origin of the preset coordinate system Ow, and the vector between Ow and the midpoint of CD is the direction vector X of the X-axis. Then, the direction vector of the Z-axis can be expressed as OwD⟶×OwC⟶, and the direction vector of the Y-axis can be represented as Z×X.
(2)
Establishment of the preset coordinate system for a region of interest with a center pointMeasure the latitude and longitude of a point as the originOw of the preset coordinates and convert the coordinate to the earth rectangular coordinate system under WGS-84. Without loss of generality, in the wild, due north is usually used as the Y-direction vector, due east is the X-direction vector, and the Z direction is perpendicular to them.The preset coordinate system can be established by the above rules, as is shown in Figure3. And the transformation vector between the preset coordinate system and the earth rectangular coordinate system is obtained as TE=Ow.Normalize the three orthogonal direction vectors:(5)I=XX,YY,ZZ,where  represents the modular arithmetic. Then, we can obtain the rotation matrix RE between the preset coordinate system and the earth rectangular coordinate system:
(6)RE=I−1.Thus, all space coordinatesPEUAV of UAV in the earth rectangular coordinate system can be converted to the preset coordinate system (new world coordinate system):
(7)P=RE·PEUAV−TE.


### 2.3. Single-Camera Calibration

The basic condition of parameter estimation is to find the matching relationship between image coordinates and 3D coordinates. In this paper, the centroid of the UAV is designated as the feature points in the left and right cameras, as shown in Figure4.Figure 1
Camera pinhole imaging mode.First, the initial values of intrinsic parameters are given based on the theoretical values:(8)fx=fdx,fy=fdy,u0=umax2,v0=vmax2,where f is the theoretical focal length and fx and fy are the focal lengths (in pixels). dx and dy represent the physical size of pixels in the x and y directions, respectively, and umax and vmax represent the resolution of the image in the x and y directions, respectively. Subsequently, initial solutions for other parameters (such as extrinsic parameters) can be obtained by DLT [8]. Finally, the constrained adjustment method is used to minimize the reprojection errors:
(9)argmin∑n=1Nmn−m∧fx,fy,u0,v0,R,T,Mn2,where mn represents the image coordinate of the n-th point, Mn denotes its corresponding spatial coordinate, m^ is the projection of point Mn in image n according to equation (1), and R and T represent the rotation matrix and translation vector, respectively. It is worth noting that the principal points u0 and v0 are set to constant values and do not participate in the iterative process because their values are an order of magnitude smaller than the other parameters. Otherwise, although small reprojection errors can be obtained, these values have no physical significance and cause instability of other parameters. It has been experimentally proven in Reference [27] that the fixation of the principal points has little effect on the final reconstruction accuracy.


### 2.4. Binocular Calibration

According to Section2.2, we can convert the GPS navigation coordinates of UAV to the required preset coordinate system. Suppose P is the coordinates of a point in the preset coordinate system and PL and PR are its corresponding coordinates in the left and right camera coordinate system, as is shown in Figure 5, their relationship can be described as
(10)PL=RLP+TL,PR=RRP+TR,where RL,TL and RR,TR describe the extrinsic parameters of the left and right cameras, respectively. Obviously, it is easy to obtain the extrinsic parameters between the left and right cameras as
(11)Rw=RRRL−1,Tw=TR−RRRL−1TL.In this way, the binocular camera can be calibrated. The advantage of establishing the preset world coordinate system is that we can quickly convert the coordinatesPL (generally in the left camera coordinate system) reconstructed by the binocular camera to the preset coordinate system:
(12)P=RL−1PL−TL.


## 2.1. Camera Imaging Model

This paper focuses on where the camera is 500 m-1000 m away from the center of the monitoring area; therefore, the telephoto lens is used. Considering that the telephoto lens of the camera has very little distortion, the ideal pinhole imaging mode [25] is chosen to describe the mapping relationship between the object space and the image space, as is shown in Figure 1.According to the pinhole imaging model, the coordinates from the world coordinate system, camera coordinate system, camera physical coordinate system and image pixel coordinate system have undergone three parts of rigid body transformation, projection transformation, and rigid body transformation, respectively. One point in the camera coordinate system is expressed asXc,Yc,Zc and in the world coordinate system is expressed as Xw,Yw,Zw, which are named x,y and u,v, respectively, in the camera physical and image pixel coordinate systems. According to the relationship between each coordinate system, the linear transformation relationship between the world coordinate system Xw,Yw,Zw and the pixel coordinate system u,v is established by the following equation.
(1)Zcuv1=1dx0u001dyv0001f0000f000010RT01XwYwZw1,where R and T are the rotation matrix and translation vector between the world coordinate system and the camera coordinate system. f is the focal length of the lens. dx and dy are the physical size of the pixel. u0 and v0 are the camera principal points. Camera calibration requires the solution of these parameters.


## 2.2. Coordinate System Conversion

As can be seen from Section2.1, obtaining the correspondence between the pixel coordinates and the world coordinates of the feature point is the key to estimating the camera parameters. With the help of the UAV working in RTK mode, we can obtain the UAV’s current GPS navigation coordinates PGB,L,H, which can be converted into earth rectangular coordinates PEXE,YE,ZE [18]:
(2)XE=N+HcosBcosL,YE=N+HcosBsinL,ZE=N1−E2+HsinB,whereN is the radius of curvature of the ellipsoid and E is the first eccentricity of the ellipsoid. Let a,b be the long and short semiaxes of the Earth, respectively, and χ be the ellipsoidal flattening rate of the Earth. Without losing generality [26],
(3)a=6378137,χ=1298.257223563,b=a×1−χ.Also, we know(4)E=a2+b2a,N=a1−E2sin2B.By combining equations (2)–(4), the representation of PG in the WGS-84 earth rectangular coordinate system, PE, can be obtained. However, the coordinates obtained by the above steps are based on the Earth’s center of mass, which has two disadvantages: first, the scale of the obtained coordinates is too large to estimate the camera parameters; second, the origin and direction of the current coordinate axes have been fixed, which is not conducive to further measurement. Therefore, we need to set the origin and direction of the world coordinate system according to our own needs and complete the camera calibration in this coordinate system which is also called the preset coordinate system. Since coordinates of all measuring points are in the WGS-84 earth rectangular coordinate system, this transformation is not difficult. According to different application scenarios, the preset coordinate systems can be established by the following two ways.
(1)
Establishment of the preset coordinate system for a rectangular region of interestAs is shown in Figure2, the latitude and longitude of A, B, C, and D four points are measured at the four corners of a rectangle and converted to the earth rectangular coordinate system by the above steps. Take the coordinates of the intersection of lines AC and BD as the origin of the preset coordinate system Ow, and the vector between Ow and the midpoint of CD is the direction vector X of the X-axis. Then, the direction vector of the Z-axis can be expressed as OwD⟶×OwC⟶, and the direction vector of the Y-axis can be represented as Z×X.
(2)
Establishment of the preset coordinate system for a region of interest with a center pointMeasure the latitude and longitude of a point as the originOw of the preset coordinates and convert the coordinate to the earth rectangular coordinate system under WGS-84. Without loss of generality, in the wild, due north is usually used as the Y-direction vector, due east is the X-direction vector, and the Z direction is perpendicular to them.The preset coordinate system can be established by the above rules, as is shown in Figure3. And the transformation vector between the preset coordinate system and the earth rectangular coordinate system is obtained as TE=Ow.Normalize the three orthogonal direction vectors:(5)I=XX,YY,ZZ,where  represents the modular arithmetic. Then, we can obtain the rotation matrix RE between the preset coordinate system and the earth rectangular coordinate system:
(6)RE=I−1.Thus, all space coordinatesPEUAV of UAV in the earth rectangular coordinate system can be converted to the preset coordinate system (new world coordinate system):
(7)P=RE·PEUAV−TE.


## 2.3. Single-Camera Calibration

The basic condition of parameter estimation is to find the matching relationship between image coordinates and 3D coordinates. In this paper, the centroid of the UAV is designated as the feature points in the left and right cameras, as shown in Figure4.Figure 1
Camera pinhole imaging mode.First, the initial values of intrinsic parameters are given based on the theoretical values:(8)fx=fdx,fy=fdy,u0=umax2,v0=vmax2,where f is the theoretical focal length and fx and fy are the focal lengths (in pixels). dx and dy represent the physical size of pixels in the x and y directions, respectively, and umax and vmax represent the resolution of the image in the x and y directions, respectively. Subsequently, initial solutions for other parameters (such as extrinsic parameters) can be obtained by DLT [8]. Finally, the constrained adjustment method is used to minimize the reprojection errors:
(9)argmin∑n=1Nmn−m∧fx,fy,u0,v0,R,T,Mn2,where mn represents the image coordinate of the n-th point, Mn denotes its corresponding spatial coordinate, m^ is the projection of point Mn in image n according to equation (1), and R and T represent the rotation matrix and translation vector, respectively. It is worth noting that the principal points u0 and v0 are set to constant values and do not participate in the iterative process because their values are an order of magnitude smaller than the other parameters. Otherwise, although small reprojection errors can be obtained, these values have no physical significance and cause instability of other parameters. It has been experimentally proven in Reference [27] that the fixation of the principal points has little effect on the final reconstruction accuracy.


## 2.4. Binocular Calibration

According to Section2.2, we can convert the GPS navigation coordinates of UAV to the required preset coordinate system. Suppose P is the coordinates of a point in the preset coordinate system and PL and PR are its corresponding coordinates in the left and right camera coordinate system, as is shown in Figure 5, their relationship can be described as
(10)PL=RLP+TL,PR=RRP+TR,where RL,TL and RR,TR describe the extrinsic parameters of the left and right cameras, respectively. Obviously, it is easy to obtain the extrinsic parameters between the left and right cameras as
(11)Rw=RRRL−1,Tw=TR−RRRL−1TL.In this way, the binocular camera can be calibrated. The advantage of establishing the preset world coordinate system is that we can quickly convert the coordinatesPL (generally in the left camera coordinate system) reconstructed by the binocular camera to the preset coordinate system:
(12)P=RL−1PL−TL.


## 3. Experiments and Analysis

To verify the effectiveness of the proposed method, we set up a series of experiments. Five groups of camera-lens pairs were calibrated independently. Details of the camera-lens pairs are shown in Table1.Figure 2
Establishment of the preset coordinate system for a rectangular region of interest.Table 1
Details of the camera-lens pairs.
GroupCamera resolutionCamera modelLens modelPixel size (μm)Focal length (mm)11920×1080Phantom V341Nikon 70-200 mm1010021024×1024Photron Nova s12Nikon 70-200 mm2013031280×800Phantom VEO 310Nikon 200-500 mm2035041280×800Phantom VEO 310Nikon 200-500 mm2035051920×1080Phantom VEO 440Nikon 70-200 mm10170In each group, identical camera-lens pairs were used to form a stereo camera, with the two cameras placed vertically, while monitoring an area 500-1000 meters away. The area covered by the cameras varies in diameter from 50 m to 100 m, depending on the focal lengths.In the experiment, the UAV (DJI M300) with RTK (DJI RTK-2) was used as the high-precision mobile calibration target. The RTK master station was arranged on the ground, and the fuselage was equipped with the RTK slave station. In the range of 10 km, the measuring accuracy of the slave station can reach the order of centimeters [28], which is a satisfactory accuracy compared with the camera monitoring diameter of tens of meters.Control the UAV navigate over the monitoring area, and confirm that the UAV is in the field of view of the cameras. At 8 m, 16 m, 24 m, 32 m, and 40 m above the planeX‐O‐Y in the preset coordinate system, 10 points were suspended to record the GPS navigation coordinates and corresponding image coordinates of the UAV. Figure 6 illustrates the UAV images taken by two cameras. Convert the GPS coordinates to the preset coordinate system, and the position distribution of the UAV is shown in Figure 7.Figure 3
Establishment of the preset coordinate system for a region of interest with a center point.Figure 4
UAV in a single camera.Figure 5
Binocular vision model [27].Figure 6
UAV in camera images.Figure 7
Position distribution of UAV in the preset coordinate system.Figure 8
Relationship between errors and the number of feature points.


### 3.1. Influence of the Feature Point Number on Calibration Results

As we know, the camera parameters can be correctly estimated only if there are at least six sets of 2D and 3D coordinates corresponding to each other. Adding a feature point means that the UAV needs to fly one more time, which will undoubtedly increase our workload. Therefore, it is meaningful to explore the appropriate number of feature points to reduce the work. Five independent experiments were carried out for the five camera-lens pairs described in Table1.In each experiment, 6, 10, and 40 UAV images (one image corresponds to a feature point position) were used to calibrate the stereo cameras. Then, the calibration results were used to reconstruct the space positions of another 10 UAVs. It is worth noting that the navigation coordinates measured by the GPS on the fuselage were used as the real space position of the UAV positions. Table2 reveals the influence of different numbers of feature points on calibration results, in which the mean Euclidean distances of the reconstructed space positions and ideal ones of UAV are used to evaluate the accuracy of the results.Table 2
Errors between reconstructed positions and ideal positions (unit: m).
Point numberGroup 1Group 2Group 3Group 4Group 5612.690.398.368.683.2101.440.313.320.350.21120.160.320.110.360.19150.130.300.090.280.21200.120.280.090.220.16250.120.280.090.170.17300.110.260.080.150.15350.110.250.080.140.16400.100.220.080.130.15As is shown in Figure8, the results of five experiments show that when the number of feature points is less than 12, the reconstruction errors decrease rapidly with the increase in the number of feature points. However, when the number of feature points is greater than 12, the impact of the number of feature points on the accuracy becomes smaller and the reconstruction accuracy only improves slightly. Therefore, 15~30 points are a good choice to balance efficiency and accuracy in practical applications.


### 3.2. Reconstruction Accuracy

The actual measurement accuracy is an important criterion to evaluate the calibration accuracy. Two markers were placed in the monitoring area of the cameras, and the actual distance between them can be measured by RTK. The same steps were used to calibrate the two cameras, and the coordinates of the two markers were reconstructed according to the calibration results, and then, the distance between them was calculated. Experiments were carried out on the five groups of camera-lens configurations, and the reconstruction errors are shown in Table3.Table 3
Distance errors of the two marks (unit: m).
Group numberMeasured lengthReconstruction lengthAbsolute errorsRelative errors124.2524.360.110.45%232.8332.770.060.18%346.6046.480.120.26%446.6046.500.100.21%530.0029.950.050.17%It can be seen that the reconstruction results are stable in accuracy, the maximum absolute error is less than 0.12 m, and the relative error is less than 0.5%. This is satisfactory when the monitoring diameter ranges from 50 m to 100 m. The results show that the proposed method is accurate and flexible in calibrating cameras with large field of view in the wild.


## 3.1. Influence of the Feature Point Number on Calibration Results

As we know, the camera parameters can be correctly estimated only if there are at least six sets of 2D and 3D coordinates corresponding to each other. Adding a feature point means that the UAV needs to fly one more time, which will undoubtedly increase our workload. Therefore, it is meaningful to explore the appropriate number of feature points to reduce the work. Five independent experiments were carried out for the five camera-lens pairs described in Table1.In each experiment, 6, 10, and 40 UAV images (one image corresponds to a feature point position) were used to calibrate the stereo cameras. Then, the calibration results were used to reconstruct the space positions of another 10 UAVs. It is worth noting that the navigation coordinates measured by the GPS on the fuselage were used as the real space position of the UAV positions. Table2 reveals the influence of different numbers of feature points on calibration results, in which the mean Euclidean distances of the reconstructed space positions and ideal ones of UAV are used to evaluate the accuracy of the results.Table 2
Errors between reconstructed positions and ideal positions (unit: m).
Point numberGroup 1Group 2Group 3Group 4Group 5612.690.398.368.683.2101.440.313.320.350.21120.160.320.110.360.19150.130.300.090.280.21200.120.280.090.220.16250.120.280.090.170.17300.110.260.080.150.15350.110.250.080.140.16400.100.220.080.130.15As is shown in Figure8, the results of five experiments show that when the number of feature points is less than 12, the reconstruction errors decrease rapidly with the increase in the number of feature points. However, when the number of feature points is greater than 12, the impact of the number of feature points on the accuracy becomes smaller and the reconstruction accuracy only improves slightly. Therefore, 15~30 points are a good choice to balance efficiency and accuracy in practical applications.


## 3.2. Reconstruction Accuracy

The actual measurement accuracy is an important criterion to evaluate the calibration accuracy. Two markers were placed in the monitoring area of the cameras, and the actual distance between them can be measured by RTK. The same steps were used to calibrate the two cameras, and the coordinates of the two markers were reconstructed according to the calibration results, and then, the distance between them was calculated. Experiments were carried out on the five groups of camera-lens configurations, and the reconstruction errors are shown in Table3.Table 3
Distance errors of the two marks (unit: m).
Group numberMeasured lengthReconstruction lengthAbsolute errorsRelative errors124.2524.360.110.45%232.8332.770.060.18%346.6046.480.120.26%446.6046.500.100.21%530.0029.950.050.17%It can be seen that the reconstruction results are stable in accuracy, the maximum absolute error is less than 0.12 m, and the relative error is less than 0.5%. This is satisfactory when the monitoring diameter ranges from 50 m to 100 m. The results show that the proposed method is accurate and flexible in calibrating cameras with large field of view in the wild.


## 4. Conclusion

In this paper, a camera calibration method for long-distance photogrammetry using unmanned aerial vehicles is studied. Instead of traditional targets, the GPS carried by UAV is used to obtain the spatial coordinate information, so as to complete camera calibration. This method overcomes the problem that standard target cannot cover most of the camera’s field of view and enhances the environmental adaptability. In addition, by using the WGS-84 coordinate system as the intermediary, the preset coordinate system can be established in any area of interest, improving the flexibility of measurement. Experimental results show that the absolute measurement error of the proposed method is less than 0.5% in the monitoring area with a diameter of 50-100 m and at the distance of 500-1000 m from the cameras.



---

*Source: 8573315-2022-05-06.xml*