S1, obtaining aerial photographs and ground control point data from a survey area, with a number of ground control points of the ground control point data denoted as n, and performing point-placing and aerial triangulation densification on the aerial photographs;

S2, randomly selecting 40% to 60% of the number of the ground control points as control points to participate in bundle adjustment optimization, with a number of the control points denoted as m, and using remaining (n−m) ground control points as check points;

S3, keeping the number of the control points unchanged at m, selecting a 1^st ground control point of the ground control points as a fixed control point, and randomly selecting remaining (m−1) control points from the ground control points more than 2n times, performing the bundle adjustment optimization to the control points after each time of selection, calculating and recording errors of the check points and the control points after each bundle adjustment optimization;

wherein a formula for calculating the errors of the check points and the control points after each bundle adjustment optimization is as follows:

where Δx_i represents a difference between an observed value of each ground control point after the aerial triangulation densification and a coordinate value of each ground control point in an x direction, Δy_i represents a difference between an observed value of each ground control point after the aerial triangulation densification and a coordinate value of each ground control point in a y direction, and Δz_i represents a difference between an observed value of each ground control point after the aerial triangulation densification and a coordinate value of each ground control point in a z direction;

S4, sequentially selecting a 2^nd ground control point to a n-th ground control point of the ground control points as the fixed control point, repeating the step S3 until iteration for the n-th ground control point is completed; and

S5, calculating and evaluating a quality coefficient Q of each ground control point to obtain a quality evaluation result of each ground control point; wherein the step S5 specifically comprises following sub-steps:

S5.1, calculating average errors of the ground control points serving as the control points and the average errors of the ground control points serving as the check points to obtain an average error graph of the ground control points as the control points and the check points;

S5.2, calculating average values u of the average errors of the ground control points as the control points and the check points, calculating standard deviations σ of the ground control points as the control points and the check points, and calculating the quality coefficient Q of each ground control point; and

S5.3, evaluating the quality coefficient Q of each ground control point to obtain the quality evaluation result of each ground control point.