Vertical Accuracy Simulation of Stereo Mapping Using a Small Matrix Charge-Coupled Device

At present, without ground control points (GCPs), the positioning accuracy of remote sensing images often fails to meet the growing requirements for mapping accuracy. Multi-load synergy to improve accuracy without GCPs by eliminating the impact of stereo accuracy, which is caused by on-orbit measurement error, is urgently needed to improve large-scale mapping. In this study, we analyzed error sources in stereo imaging mode and found that vertical accuracy depends on the relative accuracy of attitude during symmetric stereoscopic mapping. With the assistance of small matrix charge-coupled device (CCD) images and the block adjustment method, relative accuracy of attitude was improved, allowing for the improvement in vertical accuracy without GCPs. The simulation results show that vertical accuracy in symmetric stereo mode is not affected by attitude system error. After the restoration of imaging attitude processed by a sequence of matrix CCD images, the relative accuracy of the attitude increased, and the accuracy of the elevation without GCPs improved significantly. The results demonstrate the feasibility of small matrix CCD-assisted stereo mapping.


Introduction
Optical satellites with stereo mapping observation capability are among the main means of data acquisition for global mapping.A number of high-resolution optical satellites with mapping capability have been launched worldwide.As in the case of the French SPOT-5 satellite, stereo observations can be achieved using a high-resolution stereo (HRS) camera with a planimetric accuracy of less than 50 m, and a vertical accuracy of less than 15 m, without the use of ground control points (GCPs) [1,2].By using a three-linear-array of stereo cameras, the Japanese advanced land observation satellite (ALOS) satellite can implement global mapping with planimetric and vertical accuracies, reaching 8 and 10 m, respectively, without GCPs [3][4][5].The Indian Cartosat stereo satellite can achieve 7.3 and 4.7 m planimetric and vertical accuracy, respectively, without GCPs, which satisfies the requirements of 1:50,000 mapping [6].The satellites of the IKONOS and Worldview series, launched by the United States, can achieve stereo mapping with an accuracy of around 10 m, again without the use of GCPs.In recent years, the Worldview-3 satellite has achieved a planimetric accuracy of 2.16 m and a vertical accuracy of −1.61 m without GCPs, which satisfies the needs of 1:5000 mapping [7][8][9].Optical satellites used for surveying and mapping in China belong to the ZiYuan-3 (ZY3) and TianHui-1 (TH1) series.The planimetric and vertical accuracies of ZY3-01 are 10 and 5 m, respectively, without GCPs [10][11][12][13], and those for TH1-01 are 10.3 and 5.7 m, respectively [14,15].Thus, current mapping satellites can only meet the requirements of 1:50,000 stereo mapping.
Large-scale surveying and mapping places greater requirements on the hardware measurement accuracy and the stability of the platform.The accuracy by which attitude is measured by star sensors is of key importance.Some remote sensing satellites, which are equipped with advanced star sensors, can perform at a higher level.For example, Worldview-3 carries Ball Aerospace's High Accuracy Star Tracker, which has a 0.2" total error performance (1σ) [16,17]; hence, it can satisfy the accuracy need for 1:5000 or even larger-scale surveying and mapping.However, owing to technical restrictions, most satellites are equipped with relatively lower-accuracy star sensors.For example, the star sensor of Yaogan-24 in China is equipped with only 5" error in the optical axis (3σ).Although the attitude accuracy can be improved a little by post-processing, it is not sufficient to meet the accuracy requirements for large-scale surveying and mapping [18,19].
The motivation of many researchers arises from the fact that overcoming the measurement limitations of the hardware is difficult in the short term.Instead, large-scale mapping technologies have been developed using multi-load synergy, a comprehensive approach to surveying and mapping in which various types of sensors are carried on a satellite platform to achieve high measurement accuracy.For example, in the TH1 satellite series, Wang et al. proposed matching line CCD images with matrix CCD images using the equivalent frame photo (EFP) bundle adjustment, which can solve the systematic distortion of the strip mode [20,21].The vertical accuracy can reach the shuttle radar topography mission (SRTM) level without GCPs [15].The ZY3-02 satellite carries a laser altimeter with a three-linear-array camera.After calibration through multiple areas, it has a laser spot vertical accuracy of 2-3 m in areas where the terrain slope is less than 2 • .The absolute accuracy of height measurements for flat areas is better than 1 m [22,23].However, how to combine laser data with image data to improve the accuracy of three-linear-array stereo mapping requires further study.
Given that accuracy requirements for mapping are increasing, and current measurements do not always live up to those expectations, in this study we developed a method to restore satellite imaging attitude and to improve vertical accuracy without GCPs by using a combination of small matrix CCD cameras with linear array CCD cameras.First, our formula derivation proves that, under the symmetrical stereo mode, the vertical accuracy of satellite mapping is determined by the relative precision of the attitude.High-precision image attitude can be restored by using sequences of small matrix CCD images, after which vertical accuracy can be improved without the use of GCPs.We collected sequences of matrix CCD images and two-linear-array stereo images by simulation for joint processing.The results show that, when using small matrix CCD images to recover imaging attitude, the accuracy of elevation without GCPs is significantly improved, which confirms the feasibility of this method.

Accuracy Analysis of Symmetrical Stereo Mode
The error sources that affect the accuracy of surveying and mapping satellites include exterior orientation elements (e.g., position or attitude) and inner orientation elements.The influence of orbit error can be ignored, because orbit post-determination is accurate to the centimeter level [24].Inner orientation elements can be eliminated by high-precision on-orbit calibration.As a result, attitude errors (e.g., pitch, roll, and yaw) are the main error sources affecting stereo mapping without GCPs [25].Global surveying and mapping rely on the use of stereo cameras to photograph the same area from different pitch perspectives.Figure 1 shows a schematic diagram of stereo observations for both satellite forward-view and backward-view cameras.Research has shown that roll and yaw error will lead to forward intersections not being at the same point; however, this can be effectively eliminated using the free net adjustment [21,26].Therefore, the vertical accuracy of stereo mapping without GCPs depends on the accuracy of pitch, and hence our discussion will be in the XOZ plane (Figure 1).
As shown in Figure 1, the forward-view and backward-view cameras take images of the same area at times T 0 and T 1 , respectively.Assuming that the orbit coordinate system at time T 0 is a benchmark, the Z-axis of the orbit coordinate system is defined as pointing to the center of the Earth, while X-axis is defined as pointing to the satellite flight direction (along-track).Under the XOZ coordinate system, satellite images at point P have a pitch angle of ϕ 0 , ϕ 1 and coordinates (Px, Pz).Satellite coordinates are (Bx, Bz) at time T 1 .Defining counterclockwise rotation as positive, according to the relationship shown in Figure 1, we can obtain the coordinates of point P from the equation: Assuming that the altitude of the satellite orbit is H, and the average radius of Earth is R, it is obvious that: By differentiating Pz in Equation ( 1), we obtain: For the vertical accuracy of the forward intersection, the best condition is dPz = 0, such that: The base:height ratio will be controlled around 1 in the normal stereo mode, where ϕ 0 ≈ ϕ 1 ≈ 26.5 • .As R ≈ 6, 378, 140 m and H ≈ 500, 000 m, we obtain ω 0 ≈ ω 1 ≈ 2.3 • based on Equation (2).Therefore, the angle of the Z-axis between times T 0 and T 1 is about 4.6 • .Since this is so small, we can regard the Z-axis as being in the same direction at times T 0 and T 1 , with Bz = 0.When combined with Equation (3), we get: If there exist systematic errors, then dϕ 1 = dϕ 0 .Fulfilling Equation (5) in the stereo mode of the same orbit requires ϕ 0 = ϕ 1 or ϕ 0 = −ϕ 1 , based on a data range of ϕ 0 , ϕ 1 .As the perspective of the forward view and backward view cannot be the same, the situation of ϕ 0 = ϕ 1 cannot exist; therefore, ϕ 0 = −ϕ 1 is the symmetric stereo mode.When ϕ 0 = −ϕ 1 , B Z = 0, Equation (3) can be calculated as: When the attitude errors of the two stereo times are the same size and direction, the elevation errors of the forward intersection are 0 (Figure 2).As shown in Figure 1, the forward-view and backward-view cameras take images of the same area at times T0 and T1, respectively.Assuming that the orbit coordinate system at time T0 is a benchmark, the Z-axis of the orbit coordinate system is defined as pointing to the center of the Earth, while X-axis is defined as pointing to the satellite flight direction (along-track).Under the XOZ coordinate system, satellite images at point P have a pitch angle of , and coordinates (Px, Pz).Satellite coordinates are (Bx, Bz) at time T1.Defining counterclockwise rotation as positive, according  Where the stereo mode has different sides but the same track, ≥ ≥ 0 and − cot ≥ 0.Then, ( ) and ( ) have opposite signs, while and have the same sign.In this case, errors will be canceled out, and vertical accuracy is high.Where the stereo mode has the same side and track, if ≥ ≥ 0 (Figure 3) and − cot ≤ 0, then ( ) and ( ) have opposite signs, while and have the same sign.Again, errors will cancel out and the vertical accuracy will be high.In summary, for stereo modes with the same track, the influence of attitude errors will be canceled out in the forward intersection, the resulting symmetric stereo imaging will have a high elevation accuracy, and the vertical error of the forward intersection will only be affected by random noise in the attitude.Similarly, the planimetric accuracy of forward intersection is derived as follows:

Δφ Δφ
Under the conditions of stereo symmetry, if = − , = 0, we have: And the planimetric accuracy is high only when = − .This requires that attitude errors be in the opposite direction for stereo imaging, which is generally difficult to achieve.Therefore, under stereo conditions without GCPs, the possibility that planimetric accuracy is influenced by systematic attitude errors cannot be eliminated.
Based on the method above, the vertical accuracy of the forward intersection depends on the relative accuracy of attitude.The main elements that reduce the relative accuracy of attitude are (1) random error, such as the noise resulting from attitude measurement equipment (e.g., star sensor noise and gyro noise); (2) attitude high-frequency error, such as high-frequency jitter of the platform (e.g., that caused by momentum wheel activities) [27]; (3) attitude low-frequency drift, such as angle change between the CCD camera and star sensor caused by drift of the star sensor's optic axis or by Where the stereo mode has different sides but the same track, Bx ≥ Px ≥ 0 and −Px cot ϕ 1 ≥ 0.Then, f (ϕ 0 ) and f (ϕ 1 ) have opposite signs, while dϕ 0 and dϕ 1 have the same sign.In this case, errors will be canceled out, and vertical accuracy is high.Where the stereo mode has the same side and track, if Px ≥ Bx ≥ 0 (Figure 3) and −Px cot ϕ 1 ≤ 0, then f (ϕ 0 ) and f (ϕ 1 ) have opposite signs, while dϕ 0 and dϕ 1 have the same sign.Again, errors will cancel out and the vertical accuracy will be high.In summary, for stereo modes with the same track, the influence of attitude errors will be canceled out in the forward intersection, the resulting symmetric stereo imaging will have a high elevation accuracy, and the vertical error of the forward intersection will only be affected by random noise in the attitude.Where the stereo mode has different sides but the same track, ≥ ≥ 0 and − cot ≥ 0.Then, ( ) and ( ) have opposite signs, while and have the same sign.In this case, errors will be canceled out, and vertical accuracy is high.Where the stereo mode has the same side and track, if ≥ ≥ 0 (Figure 3) and − cot ≤ 0, then ( ) and ( ) have opposite signs, while and have the same sign.Again, errors will cancel out and the vertical accuracy will be high.In summary, for stereo modes with the same track, the influence of attitude errors will be canceled out in the forward intersection, the resulting symmetric stereo imaging will have a high elevation accuracy, and the vertical error of the forward intersection will only be affected by random noise in the attitude.
Similarly, the planimetric accuracy of forward intersection is derived as follows: Under the conditions of stereo symmetry, if = − , = 0, we have: And the planimetric accuracy is high only when = − .This requires that attitude errors be in the opposite direction for stereo imaging, which is generally difficult to achieve.Therefore, under stereo conditions without GCPs, the possibility that planimetric accuracy is influenced by systematic attitude errors cannot be eliminated.
Based on the method above, the vertical accuracy of the forward intersection depends on the relative accuracy of attitude.The main elements that reduce the relative accuracy of attitude are (1) random error, such as the noise resulting from attitude measurement equipment (e.g., star sensor noise and gyro noise); (2) attitude high-frequency error, such as high-frequency jitter of the platform (e.g., that caused by momentum wheel activities) [27]; (3) attitude low-frequency drift, such as angle change between the CCD camera and star sensor caused by drift of the star sensor's optic axis or by Similarly, the planimetric accuracy of forward intersection is derived as follows: Under the conditions of stereo symmetry, if ϕ 0 = −ϕ 1 , B Z = 0, we have: And the planimetric accuracy is high only when dϕ 0 = −dϕ 1 .This requires that attitude errors be in the opposite direction for stereo imaging, which is generally difficult to achieve.Therefore, under stereo conditions without GCPs, the possibility that planimetric accuracy is influenced by systematic attitude errors cannot be eliminated.
Based on the method above, the vertical accuracy of the forward intersection depends on the relative accuracy of attitude.The main elements that reduce the relative accuracy of attitude are (1) random error, such as the noise resulting from attitude measurement equipment (e.g., star sensor noise and gyro noise); (2) attitude high-frequency error, such as high-frequency jitter of the platform (e.g., that caused by momentum wheel activities) [27]; (3) attitude low-frequency drift, such as angle change between the CCD camera and star sensor caused by drift of the star sensor's optic axis or by on-orbit thermal environment change.The third factor has the greatest influence when GCPs are not used [28][29][30][31].
As shown in Figure 4, the period of on-orbit low-frequency drift error is difficult to eliminate by means of a small number of calibration fields.During the period of stereo imaging (≤100 s), attitude error was mainly systematic and was about −1.7"; however, there were also minor changes of about 0.1".Based on the above derivation, the systematic part of attitude low-frequency drift can be offset by symmetric stereo mode during the stereo time, but changeable parts of the attitude errors affect stereo vertical accuracy.Through processing a sequence of small matrix CCD images, the vertical accuracy of forward intersection caused by random errors of attitude measurement and low-frequency drift could be reduced, thereby improving the accuracy of elevation without GCPs.
As shown in Figure 4, the period of on-orbit low-frequency drift error is difficult to eliminate by means of a small number of calibration fields.During the period of stereo imaging (≤100 s), attitude error was mainly systematic and was about −1.7″; however, there were also minor changes of about 0.1″.Based on the above derivation, the systematic part of attitude low-frequency drift can be offset by symmetric stereo mode during the stereo time, but changeable parts of the attitude errors affect stereo vertical accuracy.Through processing a sequence of small matrix CCD images, the vertical accuracy of forward intersection caused by random errors of attitude measurement and low-frequency drift could be reduced, thereby improving the accuracy of elevation without GCPs.

Vertical Accuracy Improvement Using Small Matrix CCD
Figure 5 shows synchronous observations of a small matrix CCD camera and a linear array CCD camera.The satellite design ensures that the adjacent images of a small matrix CCD have an overlap.Generally, the small matrix CCD in orbit needs to be calibrated first because of the potential for lens distortion or optical perturbations.The method currently used for calibration of matrix CCDs has been extensively developed using images covering calibration fields, and the interior orientation is better than 0.3 pixels [1][2][3][4][5]32].Therefore, we assume that the matrix CCDs are perfectly planar after calibration.Under on-orbit conditions, attitude error is not considered, as shown in Figure 6.The homonymous points (p and q) of the two adjacent images should intersect at the same ground position (S); however, because of different attitude errors for moments and , the two points cannot intersect in the same position [33].Therefore, the attitude of moment relative to moment can be restored by positioning the consistency constraints of homonymous points at moments and .

Vertical Accuracy Improvement Using Small Matrix CCD
Figure 5 shows synchronous observations of a small matrix CCD camera and a linear array CCD camera.The satellite design ensures that the adjacent images of a small matrix CCD have an overlap.Generally, the small matrix CCD in orbit needs to be calibrated first because of the potential for lens distortion or optical perturbations.The method currently used for calibration of matrix CCDs has been extensively developed using images covering calibration fields, and the interior orientation is better than 0.3 pixels [1][2][3][4][5]32].Therefore, we assume that the matrix CCDs are perfectly planar after calibration.Under on-orbit conditions, attitude error is not considered, as shown in Figure 6.The homonymous points (p and q) of the two adjacent images should intersect at the same ground position (S); however, because of different attitude errors for moments t 0 and t 1 , the two points cannot intersect in the same position [33].Therefore, the attitude of moment t 1 relative to moment t 0 can be restored by positioning the consistency constraints of homonymous points at moments t 0 and t 1 .
As shown in Figure 4, the period of on-orbit low-frequency drift error is difficult to eliminate by means of a small number of calibration fields.During the period of stereo imaging (≤100 s), attitude error was mainly systematic and was about −1.7″; however, there were also minor changes of about 0.1″.Based on the above derivation, the systematic part of attitude low-frequency drift can be offset by symmetric stereo mode during the stereo time, but changeable parts of the attitude errors affect stereo vertical accuracy.Through processing a sequence of small matrix CCD images, the vertical accuracy of forward intersection caused by random errors of attitude measurement and low-frequency drift could be reduced, thereby improving the accuracy of elevation without GCPs.

Vertical Accuracy Improvement Using Small Matrix CCD
Figure 5 shows synchronous observations of a small matrix CCD camera and a linear array CCD camera.The satellite design ensures that the adjacent images of a small matrix CCD have an overlap.Generally, the small matrix CCD in orbit needs to be calibrated first because of the potential for lens distortion or optical perturbations.The method currently used for calibration of matrix CCDs has been extensively developed using images covering calibration fields, and the interior orientation is better than 0.3 pixels [1][2][3][4][5]32].Therefore, we assume that the matrix CCDs are perfectly planar after calibration.Under on-orbit conditions, attitude error is not considered, as shown in Figure 6.The homonymous points (p and q) of the two adjacent images should intersect at the same ground position (S); however, because of different attitude errors for moments and , the two points cannot intersect in the same position [33].Therefore, the attitude of moment relative to moment can be restored by positioning the consistency constraints of homonymous points at moments and .For matrix CCD images at moments and , we extracted the homonymous points ( , ) and ( , ); [ , , ] represents the position vector of ground object S corresponding to and ; [ , , ] represents the satellite's position in the Cartesian coordinates of WGS84 at time ; and represents the scaling factor at time .Assuming that ( , , ) represents the principal point and focal length of the matrix CCD, the consistency constraints of positioning mode [32,[34][35][36] can be built as: where represents the rotation matrix from camera coordinate to WGS84 coordinate at time (t stands for or in Equation ( 9)); and represents the attitude compensation of the small matrix CCD image at moment , which constitutes ∆ , ∆ , and ∆ by: where = ∆ , = ∆ , and so on.For any pair of homonymous points and , the unknown parameters in Equation ( 9) are the attitude compensation of two time points ( , ) and the ground position [ , , ] of the homonymous points.
For any time in Equation ( 9), for example, , let: From which we get: After linearizing the above equation, the error equation can be calculated as [10,26,36,37]: where , , … , are the Jacobi matrices of each unknown parameter, and ΔX, ΔY, ΔZ are the correction numbers of the ground positions of the homonymous points.By transforming the equation to the matrix form: For matrix CCD images at moments t 0 and t 1 , we extracted the homonymous points p(x t 0 , y t 0 ) and q(x t 1 , y t 1 ); [X, Y, Z] T represents the position vector of ground object S corresponding to p and q; [Xs, Ys, Zs] T t represents the satellite's position in the Cartesian coordinates of WGS84 at time t; and m t represents the scaling factor at time t.Assuming that (x 0 , y 0 , f ) represents the principal point and focal length of the matrix CCD, the consistency constraints of positioning mode [32,[34][35][36] can be built as: where R wgs84 cam t represents the rotation matrix from camera coordinate to WGS84 coordinate at time t (t stands for t 0 or t 1 in Equation ( 9)); and Ru t represents the attitude compensation of the small matrix CCD image at moment t, which constitutes ∆ϕ t , ∆ω t , and ∆κ t by: where cω = cos∆ω t , sω = sin∆ω t , and so on.For any pair of homonymous points p and q, the unknown parameters in Equation ( 9) are the attitude compensation Ru t of two time points (t 0 , t 1 ) and the ground position [X, Y, Z] T of the homonymous points.
For any time t in Equation ( 9), for example, t 0 , let: From which we get: processing, the attitude error is random.When using A2, the attitude error is systematic (dA = dA = … = dA ).As derived from Equation ( 5), when the satellite uses the symmetric stereo mode, the attitude systematic error does not affect vertical accuracy without GCPs.Therefore, restoring relative attitude through sequences of small matrix CCD images actually improves vertical accuracy by raising the relative accuracy of attitude data.Finally, by applying the attitude compensation matrix to the original attitude measurement parameters, we obtain the corrected attitude, which can be used in the production of the forward intersection between the forward-and backward-view linear array camera.

Experimental Data Simulation
The experimental parameters were simulated based on existing satellite design.Orbit, attitude, and linear array CCD parameters were from the ZY3 platform [38].The digital orthophoto map (DOM) and digital elevation model (DEM) data of the Songshan (Henan province, China) calibration field were used as the base image.Specific parameters settings are shown in Tables 1 and 2. Figure 8 shows the simulation results for the forward and backward view of the linear array CCD images and small matrix CCD images.

Symmetry and Vertical Accuracy
First, we verified the relationship between symmetry and vertical accuracy.During the initial simulation conditions, the attitude systematic error was set to omega =5", while the attitude stability and random error were set to 0. Under different symmetry conditions, the vertical accuracies of stereo image production according to the forward-and backward-view cameras are shown in Table 3 and Figure 9.At symmetry condition value 1, at which the pitch angles of the forward and backward cameras were symmetrical, the vertical accuracy reached 0.043 m.With increasing or decreasing symmetry, the vertical accuracy was reduced.By contrast, the planimetric accuracy was hardly affected by the stereo forward/backward angle.

Symmetry and Vertical Accuracy
First, we verified the relationship between symmetry and vertical accuracy.During the initial simulation conditions, the attitude systematic error was set to omega =5″, while the attitude stability and random error were set to 0. Under different symmetry conditions, the vertical accuracies of  stereo image production according to the forward-and backward-view cameras are shown in Table 3 and Figure 9.At symmetry condition value 1, at which the pitch angles of the forward and backward cameras were symmetrical, the vertical accuracy reached 0.043 m.With increasing or decreasing symmetry, the vertical accuracy was reduced.By contrast, the planimetric accuracy was hardly affected by the stereo forward/backward angle.Assuming the base:height ratio was consistent, with forward/backward cameras at −22°/+22° and a base:height ratio of 0.8, we simulated the forward camera with a 5° interval and calculated the corresponding angle of the backward camera.The attitude systematic error was omega =5″, while attitude stability and random error were 0. The vertical accuracies under different symmetries were determined according to the production of stereo images by the forward/backward camera (Table 4 and Figure 10).Again, we found that the vertical accuracy of stereo images was highest under symmetrical mode.With the reduction of symmetry, the vertical accuracy decreased, and the symmetric stereo was unable to improve the planimetric accuracy.Assuming the base:height ratio was consistent, with forward/backward cameras at −22 • /+22 • and a base:height ratio of 0.8, we simulated the forward camera with a 5 • interval and calculated the corresponding angle of the backward camera.The attitude systematic error was omega =5", while attitude stability and random error were 0. The vertical accuracies under different symmetries were determined according to the production of stereo images by the forward/backward camera (Table 4 and Figure 10).Again, we found that the vertical accuracy of stereo images was highest under symmetrical mode.With the reduction of symmetry, the vertical accuracy decreased, and the symmetric stereo was unable to improve the planimetric accuracy.The simulation results show that, in line with the foregoing theoretical derivation, the symmetric mode is the optimal stereo observation mode when attitude systematic errors exist.In actual forward/backward camera stereo images, even in symmetry stereo mode, the attitude errors at different times are different (i.e., the attitude error is not systematic).This will lead to a decrease in vertical accuracy, which requires matrix CCD image processing.

Promotion of Matrix CCD Images for Vertical Accuracy
As shown in Figure 8, adjacent images from the small matrix CCD were matched to acquire the homonymous points.When all of these points were put into the regional network adjustment mode, we were able to obtain the attitude compensation matrix for every image.Then, we compared the vertical accuracies of two stereo products: (1) products produced by forward intersection of forward/backward images through the original simulation attitude; and (2) products produced by forward intersection of forward/backward images through a modified attitude, which are compensated for by matrix .Our experiment was performed to evaluate the vertical accuracy of the forward intersection products under the following simulation conditions: attitude stability (the three axis errors were 0.0001°/s or 0.0005°/s), attitude system error (the three axis errors were 5″ or 10″), and attitude random error (the three axis errors were 1″or 2″).Vertical accuracy results are given after small matrix CCD image modification (Figure 11).
By calculating the root-mean-square error (RMSE) of the intersection point accuracy, we aquire the following Table 5:  The simulation results show that, in line with the foregoing theoretical derivation, the symmetric mode is the optimal stereo observation mode when attitude systematic errors exist.In actual forward/backward camera stereo images, even in symmetry stereo mode, the attitude errors at different times are different (i.e., the attitude error is not systematic).This will lead to a decrease in vertical accuracy, which requires matrix CCD image processing.

Promotion of Matrix CCD Images for Vertical Accuracy
As shown in Figure 8, adjacent images from the small matrix CCD were matched to acquire the homonymous points.When all of these points were put into the regional network adjustment mode, we were able to obtain the attitude compensation matrix Ru for every image.Then, we compared the vertical accuracies of two stereo products: (1) products produced by forward intersection of forward/backward images through the original simulation attitude; and (2) products produced by forward intersection of forward/backward images through a modified attitude, which are compensated for by matrix Ru.
Our experiment was performed to evaluate the vertical accuracy of the forward intersection products under the following simulation conditions: attitude stability (the three axis errors were 0.0001 • /s or 0.0005 • /s), attitude system error (the three axis errors were 5" or 10"), and attitude random error (the three axis errors were 1"or 2").Vertical accuracy results are given after small matrix CCD image modification (Figure 11).
By calculating the root-mean-square error (RMSE) of the intersection point accuracy, we aquire the following Table 5:  Based on Figure 11 and Table 5, we have the following results.
1.For the original RMSE accuracy of (Figure 11a,b), the attitude stability was 0.0005°/s and attitude random error was 2″.When the systematic error increased from 5″ to 10″, the vertical accuracy only changed from 5.566 to 4.757 m, but the planimetric error increased from 17.978 to 37.393 m.When the systematic error doubled, the vertical accuracy was almost the same (Figure 11d,e).This confirms that under the symmetric stereo mode, attitude systematic error has little influence on vertical accuracy.However, the symmetric stereo mode cannot offset attitude Based on Figure 11 and Table 5, we have the following results.

1.
For the original RMSE accuracy of (Figure 11a,b), the attitude stability was 0.0005 • /s and attitude random error was 2".When the systematic error increased from 5" to 10", the vertical accuracy only changed from 5.566 to 4.757 m, but the planimetric error increased from 17.978 to 37.393 m.When the systematic error doubled, the vertical accuracy was almost the same (Figure 11d,e).This confirms that under the symmetric stereo mode, attitude systematic error has little influence on vertical accuracy.However, the symmetric stereo mode cannot offset attitude systematic error for planimetric accuracy.After matrix CCD processing, the vertical accuracy increased from 5.566 to 1.662 m (Figure 11a), and from 2.826 to 0.395 m (Figure 11d), while the planimetric accuracy did not change significantly.These results confirm that, under the symmetric stereo condition, images from the small matrix CCD can effectively restore the relative accuracy of attitude, but are unable to improve the absolute accuracy of attitude.In summary, vertical accuracy can be increased greatly by using matrix CCD processing, but planimetric accuracy changes little.

2.
Compared with the original accuracy (Figure 1a,c), under attitude stability of 0.0005 • /s and systematic error of 5", and as attitude random error changed from 2" to 1", vertical accuracy rose from 5.566 to 3.952 m.After small matrix CCD image processing, vertical accuracy increased to 1.662 m.The same situation can be seen in Figure 11d,f, which verifies that under symmetric stereo conditions, vertical accuracy is affected by the random error of the measurement (including the attitude low-frequency drift that removed the systematic part).After matrix CCD processing, the relative accuracy of the attitude was increased, and the vertical accuracy was also improved.

3.
By comparing Figure 11a-c with Figure 11d-f, the relative accuracy of the attitude can be seen to have increased indirectly with increasing attitude stability, such that the vertical accuracy without GCPs improved under symmetric stereo conditions.When attitude stability is high, restraining the influence of random errors and low-frequency drift of attitude through matrix CCD can result in vertical accuracy of better than 0.4 m without GCPs, which is suitable for large-scale mapping.

Conclusions
Based on the stereo vertical accuracy deduction of a linear push-broom satellite, we have developed a method for surveying and mapping assisted by a small matrix CCD camera.The results show that this method greatly improves vertical accuracy without the use of GCPs.Simulation parameters (i.e., orbit and attitude) were set according to those of the ZY3 satellite.At the same time, forward/backward linear array images and matrix CCD images were simulated.The symmetric stereo and matrix CCD processing effects were then verified.The experimental results show the following.

1.
Under the symmetric stereo mode, attitude systematic errors can be canceled out, which guarantees a high vertical accuracy.

2.
Attitude stability and attitude random errors are the main error sources for vertical accuracy under the symmetric stereo mode.When these are decreased, the vertical accuracy increases.

3.
The relative accuracy of the attitude can be improved using matrix CCD image processing.By using the improved attitude in stereo image production, the vertical accuracy without GCPs can be significantly improved.

Figure 2 .
Figure 2. Elevation errors at the stereo intersection.The solid line shows the true intersection; the dashed line shows the intersection with error.

Figure 2 .
Figure 2. Elevation errors at the stereo intersection.The solid line shows the true intersection; the dashed line shows the intersection with error.

Figure 2 .
Figure 2. Elevation errors at the stereo intersection.The solid line shows the true intersection; the dashed line shows the intersection with error.

Figure 4 .
Figure 4. Attitude low-frequency drift errors of the on-orbit satellite.

Figure 5 .
Figure 5. Imaging schematic of a small matrix charge-coupled device (CCD) and linear array CCD.

Figure 4 .
Figure 4. Attitude low-frequency drift errors of the on-orbit satellite.

Figure 4 .
Figure 4. Attitude low-frequency drift errors of the on-orbit satellite.

Figure 5 .Figure 5 .
Figure 5. Imaging schematic of a small matrix charge-coupled device (CCD) and linear array CCD.Figure 5. Imaging schematic of a small matrix charge-coupled device (CCD) and linear array CCD.

Figure 6 .
Figure 6.Schematic of the intersection of homonymous points.

Figure 6 .
Figure 6.Schematic of the intersection of homonymous points.

Figure 8
Figure8shows the simulation results for the forward and backward view of the linear array CCD images and small matrix CCD images.

Figure 8 .
Figure 8. Simulation image schematic of the linear array charge-coupled device (CCD) and the small matrix CCD.

Figure 8 .
Figure 8. Simulation image schematic of the linear array charge-coupled device (CCD) and the small matrix CCD.

Figure 10 .
Figure 10.Vertical accuracy under the same base:height ratio.

Figure 10 .
Figure 10.Vertical accuracy under the same base:height ratio.

Table 1 .
Orbit and attitude parameters.

Table 3 .
Simulation results for symmetry and vertical accuracy.

Table 3 .
Simulation results for symmetry and vertical accuracy.
Forward/Backward Camera Angle Symmetry Planimetric Accuracy/m Vertical Accuracy/m

Table 4 .
Vertical accuracy simulation of stereo images under the same base:height ratio.

Table 4 .
Vertical accuracy simulation of stereo images under the same base:height ratio.