Research into Image Point Positioning Compensation of a High-Orbit Thermal Infrared Staring Camera

: The sight of the high-orbit thermal infrared staring camera is concentrated and sensitive to temperature changes, and therefore its image point positioning is key to ensuring the geometric quality of the thermal infrared staring image and eliminating the errors of the imaging system in orbit. According to the geometric imaging characteristics of the high-orbit thermal infrared staring camera, internal and external positioning compensation models are proposed in this paper based on the two-dimensional pointing angle of temperature change, and four kinds of experimental schemes of image point positioning compensation are designed. The experimental results show that the method proposed in this paper has a good compensation effect on the image point positioning error.


Introduction
Infrared thermal imaging plays an increasingly important role in military and civilian fields because of its good concealment, strong anti-interference, strong target recognition ability and all-weather work.An infrared thermal imager uses photoelectric technology to detect the infrared radiation emitted by the object in a passive way [1], calculates the temperature of each point on the surface of the object, displays different temperatures in different colors, and converts them into images and graphics that can be distinguished by human vision.An infrared thermal imager can break through human visual barriers and detect objects in a completely dark environment.It can detect objects even in the presence of smoke and dust, and does not require light source illumination, so it can be used 24 h per day.Compared with low-orbit satellites, high-orbit satellites have a larger observation range and incomparable temporal resolution advantages due to their high orbital altitude, and can complete more complex observation tasks, such as mobile target search and the tracking and periodic monitoring of targets.As an important link in the space-based thermal infrared observation system, high-orbit thermal infrared detectors are used for tactical intelligence collection, field situational awareness tasks, forest fire positioning and so on, so their imaging quality is closely related to the observation quality.In space-based early warning, the geosynchronous orbit satellite is mainly used for the detection and monitoring of the active section of ballistic missiles.It is equipped with an infrared scanning detector and a high-resolution infrared staring detector.The working mechanism is that the infrared scanning detector detects the tail flame of the ballistic missile.After discovering the ballistic missile target, it is transferred to the infrared staring detector to stare and track the target, determine the active section trajectory of the target, and achieving the accurate detection of the target.Therefore, it is necessary to study the imaging of the geosynchronous orbit early warning satellite.The on-orbit imaging process of the high-orbit camera, similar to other types of optical cameras, is actually a process in which multiple on-board loads such as attitude measurement, orbit measurement, and time synchronization measurement work together.In this process, there are various external measurement results errors.In addition, due to the complex operating environment of the geosynchronous satellite, the camera itself also has a certain degree of internal distortion.Therefore, the geometric processing of the high-orbit image is very important.
In recent years, high-precision geometry processing technologies for middle and low orbit satellites and aerial images have become more developed [2], including geometric imaging models, in-orbit geometric calibration models and sensor rectifying models [3][4][5].Xiu, J. [6] developed the line-of-sight (LOS) algorithm and control flow and image motion compensation based on the imaging optical path.Kornus [7] and others conducted an in-depth exploration on the sensor laboratory calibration and the in-orbit rectifying of the MOMS-2P satellite, and for each CCD array they built a corresponding in-orbit geometric rectifying model, mainly correcting the internal geometric distortion of the camera, including the deviation of the image principal point, the deviation of the focal length, and the bending and rotating of the CCD, and finally a uniform solution was obtained through the method of bundle adjustment.Gachet [8] and others carried out the in-orbit geometric calibration of a SPOT5 HRG camera and HRS camera, respectively, through high-precision 50 cm high-resolution reference data provided by the calibration site on the ground.They, respectively, used the calibration model of quintic polynomial fitting to compensate for the errors in the directions alongside the orbit and perpendicular to the orbit, and this could ultimately reduce the distortion error inside the camera to less than 0.1 pixel.Cai, Y. [9] established different positioning models and analyzed the calculation principle of the aerial remote sensing image positioning algorithm.Jacobsen [10] and others set 15 additional parameters according to the characteristics of India's IRS-1C satellite camera and verified the impact of each additional parameter on the positioning accuracy by designing experiments with different combinations of additional parameters.Japan's ALOS [11] satellite used the self-rectifying regional network adjustment with additional parameters to carry out the overall calibration.Su Wenbo [12] studied the in-orbit geometric calibration of the space linear array CCD sensor, established the detection model of the external orientation elements with the help of SPOT5 HRS image auxiliary data, and calibrated the space linear array CCD sensor using the self-rectifying bundle adjustment to achieve good calibration accuracy.Wang Chongyang [13] proposed a simplified line-of-sight vector correction model in the form of a non-pointing angle for in-orbit geometric calibration and carried out relevant experiments with the help of simulation data, proving that the model did not need to solve the satellite earth distance, and with a simplified calculation process, it had stronger applicability.
Most of the above research work is aimed at the in-orbit geometric calibration of the low-orbit remote sensing linear array satellites to the ground.At present, there has been little research into the error compensation of high-orbit infrared cameras.For the geosynchronous orbit satellite, the solar illumination is uneven, and the satellite temperature changes dramatically [14], with the variation in laws of its temperature of annual long periodicity and short diurnal periodicity [15].The infrared camera collects the temperature information of sensitive targets in the geosynchronous orbit, and due to its high-orbit height, the concentrated imaging line of sight, and severe temperature changes, the geometric imaging model needs to be studied [16].The geosynchronous orbit satellite is affected by the drastic change in temperature, which will cause the relative installation relationship between the star sensor and the early warning camera to change.The change in the installation angle will be directly reflected in the error of attitude measurement.Therefore, we have carried out a detailed study on the influence of temperature on the camera distortion error.
In addition, most of the objects observed by this type of satellite show point features, and there is almost no ground control point information in the process of real-time geo-metric processing, which brings great difficulties regarding the accuracy of the motion parameter estimation of high-value targets [17].Therefore, according to the geometric imaging characteristics and the error sources of the high-orbit thermal infrared staring camera, this paper proposes a geometric imaging model based on the two-dimensional pointing angle of temperature change.It can provide effective help to realize the twodimensional trajectory reconstruction and provide more accurate observation trajectories for target parameter estimation.

Compensation Scheme
In this paper, we construct an error compensation model for the image point positioning of early warning area array cameras according to the two-dimensional pointing angle model, including an external parameter error compensation model and an internal parameter error compensation model.
The external parameter compensation model is as follows: is the probe coordinate corresponding to the image point on the image, ψ x (x, y), ψ y (x, y) are the sight directions of the x direction and the y direction, respectively, in the camera coordinate system, (X S , Y S , Z S ) is the position coordinate of the satellite in the earth fixed coordinate system, (X, Y, Z) is the position coordinate of the object point corresponding to the image point in the earth fixed coordinate system, R J2000 WGS84 is the conversion matrix from the earth fixed coordinate system to the J2000 coordinate system, R CAM I NS (α, β, θ) (t) refers to the time-varying conversion matrix between the star sensor and the camera under different temperature conditions, and (α, β, θ) refers to the relative installation angle between the star sensor and the camera, which is also the compensation object, and λ is the scale coefficient.
The calibration results of the installation angle at different temperatures are calculated using the simulation data samples.We establish the variation curve of the installation angle and time, and then apply it to the compensation of the external parameters.Finally, the image point positioning error before and after compensation is directly calculated to verify the accuracy of the external parameter error compensation.
The internal parameters compensation model is as follows: in which xa 0 ∼ xa 9 , ya 0 ∼ ya 9 are internal compensation parameters.
For the internal parameter compensation, the two-dimensional pointing angle corresponding to the probe on the area array camera can be calculated by using the calculated internal calibration parameters, and the deviation between the two-dimensional pointing angle of the sample probe and the two-dimensional pointing angle after adding the distortion parameters can be counted, and the pointing angle after calibration can be displayed intuitively.Or the two-dimensional pointing angle of the probe element can be calculated according to the calculated internal calibration parameters, and compared with the image point coordinates calculated according to the original internal parameters to compare the image point positioning errors before and after calibration.

Calibration Principle
To solve the problem of the motion parameter estimation error of high-value point targets caused by the lack of ground control points in the process of real-time geometric processing [18], this paper proposes a calibration method of virtual control points based on image space uniformity, with the image points of the target projected onto the Earth's surface through the photography center to obtain the corresponding projection point sequence of the Earth's surface, so as to obtain the virtual control points, as shown in Figure 1.For the internal parameter compensation, the two-dimensional pointing angle corresponding to the probe on the area array camera can be calculated by using the calculated internal calibration parameters, and the deviation between the two-dimensional pointing angle of the sample probe and the two-dimensional pointing angle after adding the distortion parameters can be counted, and the pointing angle after calibration can be displayed intuitively.Or the two-dimensional pointing angle of the probe element can be calculated according to the calculated internal calibration parameters, and compared with the image point coordinates calculated according to the original internal parameters to compare the image point positioning errors before and after calibration.

Calibration Principle
To solve the problem of the motion parameter estimation error of high-value point targets caused by the lack of ground control points in the process of real-time geometric processing [18], this paper proposes a calibration method of virtual control points based on image space uniformity, with the image points of the target projected onto the Earth's surface through the photography center to obtain the corresponding projection point sequence of the Earth's surface, so as to obtain the virtual control points, as shown in Figure 1.At the observation time, the image points of the virtual control point sequence are completely coincident with the image points of the target point at that time, so a series of virtual control points are generated by using the intersection points obtained.Given the At the observation time, the image points of the virtual control point sequence are completely coincident with the image points of the target point at that time, so a series of virtual control points are generated by using the intersection points obtained.Given the coordinates (x, y) of the image plane of the target and the internal and external orientation elements at the observation time, the spatial coordinates (X, Y, Z) of the corresponding virtual control points in the earth fixed coordinate system can be obtained from the geometric imaging model: For the generated virtual control image, N generated control points are selected as the calibration points, with the object space coordinates corresponding to the calibration points set as (X i , Y i , Z i ) and the corresponding image point coordinates set as (x i , y i ), i = 1, 2, 3, • • •, N, and then the calibration model is decomposed into the pointing angle error along the direction of x and the pointing angle error along the direction of y according to the two-dimensional pointing angle model: The initial values of the extrinsic calibration parameters are determined according to the initial installation matrix of the camera, and the initial values of the intrinsic calibration parameters are determined according to the parameters of the theoretical two-dimensional pointing angle model when the camera has no internal distortion.
The current value of the intrinsic calibration parameter W N is regarded as the true value, and the extrinsic calibration parameter M W is regarded as the unknown parameter to be determined.For each virtual calibration control point, the above equation is linearized to establish the corresponding error equation: When solving the intrinsic calibration parameters, the extrinsic calibration parameter M W obtained above is regarded as the true value, and the intrinsic calibration parameter W N as the unknown parameter to be determined.For each virtual calibration control point, the above equation is linearized to establish the corresponding error equation:

Step-by-Step Compensation Scheme
Due to the strong correlation between the internal and external parameters, if the internal and external parameters are compensated at the same time, an ill-conditioned equation will be generated.Accordingly, this paper adopts the strategy of step-by-step error compensation, alternately solving the internal and external parameter values of the error equation, so as to carry out more accurate error compensation for image point positioning.According to the different order of solving parameters, the compensation schemes can be divided into the following two kinds: internal compensation before external compensation, and external compensation before internal compensation.The scheme of "the external first and then the internal" compensates the external parameters of the camera, updates the internal parameter compensation model in real time using the external parameters after compensation, and then compensates the internal parameters using the updated internal parameter compensation model.On the other hand, the scheme of "the internal first and then the external" compensates the camera's internal parameters first, updates the external parameter compensation model in real time using the compensated internal parameters, and then compensates the external parameters using the updated external parameter compensation model.
In order to better analyze the accuracy of parameter compensation, in the abovementioned external-internal compensation scheme, after the first step of external parameter compensation, the statistical accuracy is only the accuracy of the external parameter compensation.Similarly, in the internal and external compensation scheme, the accuracy calculation is carried out after the internal parameter compensation, which is the accuracy of the internal parameter compensation.At this time, there are four kinds of compensation schemes (as shown in Figure 2) actually established: first outside and then inside, first inside and then outside, only inside compensation, and only outside compensation.
Appl.Sci.2023, 13, 9688 7 of 16 In order to better analyze the accuracy of parameter compensation, in the abovementioned external-internal compensation scheme, after the first step of external parameter compensation, the statistical accuracy is only the accuracy of the external parameter compensation.Similarly, in the internal and external compensation scheme, the accuracy calculation is carried out after the internal parameter compensation, which is the accuracy of the internal parameter compensation.At this time, there are four kinds of compensation schemes (as shown in Figure 2) actually established: first outside and then inside, first inside and then outside, only inside compensation, and only outside compensation.The specific compensation process is shown in the following steps: 1.The spatial motion model of the target and the strict imaging geometric model of the area array camera are established.The three-dimensional motion trajectory of the generated target space is observed and the imaging simulation results of the motion trajectory on the two-dimensional image are obtained.2. In the two-dimensional imaging process, three-dimensional errors are added to the attitude parameters and orbital position parameters, and distortion errors are added to the camera parameters.The calculated position is used as the positioning result before error compensation.3. The grid points with uniform and regular distribution are selected on the camera array.The pixel coordinates (0, 0) in the upper left corner are taken as the starting grid points, and the pixel coordinates (10,000, 10,000) in the lower right corner are taken as the ending grid points.The interval between the grid points in the X direction and the Y direction is 100 pixels.4. According to the pixel size design parameters, the pixel coordinates of these grid points are transformed into the coordinates in the camera coordinate system, and the two-dimensional pointing angle in the camera coordinate system is calculated as the initial coordinate value.5.For the image point positioning error caused by the internal orientation elements, the actual two-dimensional pointing angle coordinates of each grid point in the camera The specific compensation process is shown in the following steps: 1.
The spatial motion model of the target and the strict imaging geometric model of the area array camera are established.The three-dimensional motion trajectory of the generated target space is observed and the imaging simulation results of the motion trajectory on the two-dimensional image are obtained.2.
In the two-dimensional imaging process, three-dimensional errors are added to the attitude parameters and orbital position parameters, and distortion errors are added to the camera parameters.The calculated position is used as the positioning result before error compensation.

3.
The grid points with uniform and regular distribution are selected on the camera array.The pixel coordinates (0, 0) in the upper left corner are taken as the starting grid points, and the pixel coordinates (10,000, 10,000) in the lower right corner are taken as the ending grid points.The interval between the grid points in the X direction and the Y direction is 100 pixels.

4.
According to the pixel size design parameters, the pixel coordinates of these grid points are transformed into the coordinates in the camera coordinate system, and the two-dimensional pointing angle in the camera coordinate system is calculated as the initial coordinate value.5.
For the image point positioning error caused by the internal orientation elements, the actual two-dimensional pointing angle coordinates of each grid point in the camera coordinate system are calculated, and the difference between the actual two-dimensional pointing angle and the initial pointing angle is used to reflect the distortion inside the camera, and the internal parameter compensation model is constructed to calculate the internal calibration parameters.6.
For the image point positioning error caused by the external orientation elements, the relative installation angle change caused by the temperature change is used to compensate the external parameters, and the external parameter compensation model is constructed to calculate the external calibration parameters.7.
After calculating the internal and external calibration parameters, the error compensation is carried out according to the two schemes proposed in this paper.The positioning results before and after error compensation are compared with the theoretical values to obtain the optimal compensation scheme.
For external parameter compensation, the image point positioning error before and after compensation is directly calculated.For internal parameter compensation, the twodimensional pointing angle corresponding to the probe on the area array camera is calculated with the solved intrinsic calibration parameters, and the deviation between the two-dimensional pointing angle of the probe and the two-dimensional pointing angle after adding distortion parameters is counted.For the temperature error, the calibration results of the installation angle at different temperatures are calculated, the change curve of the installation angle and time is established, polynomial fitting is performed, and then it is applied to the compensation of external parameters.

Two-Dimensional Pointing Angle Experiment
The orbital height of the geosynchronous orbit satellite studied in this paper is 36,000 km.Although the three-axis stability for the research object of this paper is the geosynchronous orbit satellite, the ambient temperature changes drastically [14].The drastic temperature difference will cause the installation bracket of the camera to deform violently, and the installation angle of the camera will also change accordingly, resulting in the same effect as the attitude error.For the image point positioning error, we collected the observation data of different seasons and different times of the day in a year and established the error empirical model of the camera installation angle related to temperature and time (as shown in Figure 3).
According to the camera's internal distortion parameters (as shown in Table 1), assuming that there is no external parameter error, the actual two-dimensional pointing angle coordinates of each grid point in the camera coordinate system are calculated, and the difference between the two-dimensional pointing angle and the initial coordinate value is calculated, as shown in Figure 4.It can be seen that after the distortion parameter is added, the difference of the two-dimensional probe pointing angle of the camera can represent the geometric distortion of the camera.chronous orbit satellite, the ambient temperature changes drastically [14].The drastic temperature difference will cause the installation bracket of the camera to deform violently, and the installation angle of the camera will also change accordingly, resulting in the same effect as the attitude error.For the image point positioning error, we collected the observation data of different seasons and different times of the day in a year and established the error empirical model of the camera installation angle related to temperature and time (as shown in Figure 3).According to the camera's internal distortion parameters (as shown in Table 1), assuming that there is no external parameter error, the actual two-dimensional pointing angle coordinates of each grid point in the camera coordinate system are calculated, and the difference between the two-dimensional pointing angle and the initial coordinate value is calculated, as shown in Figure 4.It can be seen that after the distortion parameter is added, the difference of the two-dimensional probe pointing angle of the camera can represent the geometric distortion of the camera.According to the internal distortion parameters of the designed camera, it is assumed that there is no external parameter error; that is, there is only internal geometric distortion error.The actual two-dimensional pointing angle coordinates of each grid point in the camera coordinate system are calculated, and the difference between the actual two-dimensional pointing angle and the initial pointing angle is used to reflect the distortion inside the camera.According to the internal parameter solution method introduced in Section 2.1, the internal error compensation is performed on the distorted two-dimensional pointing angle, and the difference between the compensated two-dimensional pointing angle and the initial pointing angle is obtained (as shown in Figure 5).It can be clearly seen that the interpolation of the two-dimensional pointing angle after compensation is reduced, and the fitting accuracy can reach the millimeter level, which verifies the effectiveness of the internal error compensation scheme proposed in this paper.According to the internal distortion parameters of the designed camera, it is assumed that there is no external parameter error; that is, there is only internal geometric distortion error.The actual two-dimensional pointing angle coordinates of each grid point in the camera coordinate system are calculated, and the difference between the actual twodimensional pointing angle and the initial pointing angle is used to reflect the distortion inside the camera.According to the internal parameter solution method introduced in Section 2.1, the internal error compensation is performed on the distorted two-dimensional pointing angle, and the difference between the compensated two-dimensional pointing angle and the initial pointing angle is obtained (as shown in Figure 5).It can be clearly seen that the interpolation of the two-dimensional pointing angle after compensation is reduced, and the fitting accuracy can reach the millimeter level, which verifies the effectiveness of the internal error compensation scheme proposed in this paper.
Figures 4 and 5 show the effectiveness of the two-dimensional pointing angle model proposed in this paper in solving the internal compensation parameters.Figures 4 and 5 show the effectiveness of the two-dimensional pointing angle model proposed in this paper in solving the internal compensation parameters.

Internal and External Parameter Error Compensation Experiment
The relative installation angle error (Table 2) is increased.The probe pointing angle of the virtual grid point in the camera coordinate system is calculated and then converted to the body coordinate system through the setting of relative installation angle to obtain the pointing angle deviation in the body coordinate system, as shown in Figure 6.

Internal and External Parameter Error Compensation Experiment
The relative installation angle error (Table 2) is increased.The probe pointing angle of the virtual grid point in the camera coordinate system is calculated and then converted to the body coordinate system through the setting of relative installation angle to obtain the pointing angle deviation in the body coordinate system, as shown in Figure 6.Figures 4 and 5 show the effectiveness of the two-dimensional pointing angle model proposed in this paper in solving the internal compensation parameters.

Internal and External Parameter Error Compensation Experiment
The relative installation angle error (Table 2) is increased.The probe pointing angle of the virtual grid point in the camera coordinate system is calculated and then converted to the body coordinate system through the setting of relative installation angle to obtain the pointing angle deviation in the body coordinate system, as shown in Figure 6.The internal and external parameters, respectively, calculated by the four schemes, are used to calculate the pointing angle of the compensated grid point in the body coordinates, and then the difference is found with the pointing angle calculated by the real internal and external parameters, as shown in Figures 7-10.The internal and external compensation parameters calculated from the above four schemes are used for the compensation calculation of the image point positioning of the target track.The positioning error after compensation is shown in Table 3.The internal and external parameters, respectively, calculated by the four schemes, are used to calculate the pointing angle of the compensated grid point in the body coordinates, and then the difference is found with the pointing angle calculated by the real internal and external parameters, as shown in Figures 7-10.The internal and external compensation parameters calculated from the above four schemes are used for the compensation calculation of the image point positioning of the target track.The positioning error after compensation is shown in Table 3.The internal and external parameters, respectively, calculated by the four schemes, are used to calculate the pointing angle of the compensated grid point in the body coordinates, and then the difference is found with the pointing angle calculated by the real internal and external parameters, as shown in Figures 7-10.The internal and external compensation parameters calculated from the above four schemes are used for the compensation calculation of the image point positioning of the target track.The positioning error after compensation is shown in Table 3.In the process of two-dimensional imaging, three-direction errors are added to the attitude parameters and orbital position parameters, and distortion errors are added to In the process of two-dimensional imaging, three-direction errors are added to the attitude parameters and orbital position parameters, and distortion errors are added to In the process of two-dimensional imaging, three-direction errors are added to the attitude parameters and orbital position parameters, and distortion errors are added to In the process of two-dimensional imaging, three-direction errors are added to the attitude parameters and orbital position parameters, and distortion errors are added to the camera parameters as the positioning results before error compensation.The positioning results before and after error compensation are compared with the theoretical values to obtain the optimal compensation scheme (as shown in Table 3).
It can be seen from the experimental results that the error compensation accuracy of the four error compensation schemes, namely, the internal first and then the external, the external first and then the internal, the external only, and the internal only, decreases sequentially.The difference between the two error compensation schemes, namely, the internal first and then the external, and the external first and then the internal, is not significant.This shows that after multiple iterations, the internal and external parameters tend to be stable, but the scheme in which the internal compensation is carried out first, followed by the external, is still relatively better than the one in which the external takes place before the internal.This is due to the different methods of solving the compensation error, and it can be concluded that it is better to finally conduct the external parameter compensation scheme.
In addition, the error compensation effect of the external scheme alone is much better than that of the internal scheme alone, and this indicates that there is still a large external parameter error after the internal parameter compensation, and the positioning error does not decrease but increases.This is because when the internal and external errors exist, there is a correlation, which will offset some of the errors.After the internal parameter compensation, most of the positioning errors caused by distortion are eliminated.On the other hand, the error caused by the external parameters is more obvious, so the external parameters need to be compensated.However, after the external parameter compensation, it can be seen that the external parameter compensation will absorb some of the internal parameter errors; that is, there will be partial substitution between the internal and external parameters.

Temperature Installation Angle Experiment
According to the image point positioning error compensation process designed above, two schemes are used to compensate the image point positioning under different temperature conditions; that is, only external parameter compensation and internal and external compensation are performed.At different times of the day, under the temperature condition, the same target trajectory is staring imaging, and the camera distortion is the same each time.The specific distortion parameters are shown in Table 1.The external relative installation angle parameters of each experiment are calculated, and a difference is made with the original relative installation angle to obtain the corresponding curve of the installation angle error related with time and compare it with the original added relative installation angle error.The compensation results of the two schemes under different temperature conditions are compared with the installation angle of the original simulation, and the results are shown in Figure 11.
Although temperature has a great effect on the relative installation angle, the change rule of the relative installation angle and time can be recovered after compensation, and the compensation scheme of "the internal first and then the external" is better than that of the external only.It can be found that even when the relative installation angle error is large, the two schemes can better calculate the true relative installation angle, and this can also verify the correctness and effectiveness of the two compensation process schemes designed above.Although temperature has a great effect on the relative installation angle, the change rule of the relative installation angle and time can be recovered after compensation, and the compensation scheme of "the internal first and then the external" is better than that of the external only.It can be found that even when the relative installation angle error is large, the two schemes can better calculate the true relative installation angle, and this can also verify the correctness and effectiveness of the two compensation process schemes designed above.

Conclusions
According to the periodicity of the temperature of the high-orbit thermal infrared staring camera, a two-dimensional internal and external positioning compensation model based on temperature change is proposed in this paper.In this paper, four kinds of experimental schemes of image point positioning compensation were designed: only internal parameter error compensation, only external parameter error compensation, first inside and outside, and first outside and inside.The internal and external parameters calculated by the four schemes were used to calculate the pointing angle of the compensated grid points in the body coordinate, and the pointing angle calculated by the real internal and external parameters was different (Figures 7-10), and the internal and external compensation parameters were calculated, respectively.The compensation calculation of the target trajectory image point positioning was carried out.The positioning error after compensation was obtained, and it was compared with the positioning result before error correction and the ideal value (as shown in Table 3).The experiment proves that the parameter error compensation scheme is the optimal solution.The experimental results show that the compensation accuracy was the highest when internal compensation was per- Although temperature has a great effect on the relative installation angle, the change rule of the relative installation angle and time can be recovered after compensation, and the compensation scheme of "the internal first and then the external" is better than that of the external only.It can be found that even when the relative installation angle error is large, the two schemes can better calculate the true relative installation angle, and this can also verify the correctness and effectiveness of the two compensation process schemes designed above.

Conclusions
According to the periodicity of the temperature of the high-orbit thermal infrared staring camera, a two-dimensional internal and external positioning compensation model based on temperature change is proposed in this paper.In this paper, four kinds of experimental schemes of image point positioning compensation were designed: only internal parameter error compensation, only external parameter error compensation, first inside and outside, and first outside and inside.The internal and external parameters calculated by the four schemes were used to calculate the pointing angle of the compensated grid points in the body coordinate, and the pointing angle calculated by the real internal and external parameters was different (Figures 7-10), and the internal and external compensation parameters were calculated, respectively.The compensation calculation of the target trajectory image point positioning was carried out.The positioning error after compensation was obtained, and it was compared with the positioning result before error correction and the ideal value (as shown in Table 3).The experiment proves that the parameter error compensation scheme is the optimal solution.The experimental results show that the compensation accuracy was the highest when internal compensation was per-

Conclusions
According to the periodicity of the temperature of the high-orbit thermal infrared staring camera, a two-dimensional internal and external positioning compensation model based on temperature change is proposed in this paper.In this paper, four kinds of experimental schemes of image point positioning compensation were designed: only internal parameter error compensation, only external parameter error compensation, first inside and outside, and first outside and inside.The internal and external parameters calculated by the four schemes were used to calculate the pointing angle of the compensated grid points in the body coordinate, and the pointing angle calculated by the real internal and external parameters was different (Figures 7-10), and the internal and external compensation parameters were calculated, respectively.The compensation calculation of the target trajectory image point positioning was carried out.The positioning error after compensation was obtained, and it was compared with the positioning result before error correction and the ideal value (as shown in Table 3).The experiment proves that the parameter error compensation scheme is the optimal solution.The experimental results show that the compensation accuracy was the highest when internal compensation was performed first, followed by external compensation.The phase angle of the external parameter compensation under different temperature conditions was consistent with the added installation error.In this paper, two schemes were used to compensate the image point positioning under different temperature conditions, namely, only external parameter compensation and internal and external compensation (as shown in Figure 11).The compensation results were compared with the installation angle of the original simulation.The two schemes could better recover the change rule of the relative installation angle and time, and the compensation scheme of internal and external compensation was better than that of external parameter compensation only.Even when the relative installation angle error was large, the two schemes could better calculate the real relative installation angle.
Despite the success of our approach, we acknowledge that there are still some limitations and challenges that need to be addressed in future research.For instance, in the process of solving the internal and external compensation parameters, due to the redundant observation virtual control points, the least squares evaluation method was adopted to solve the compensation parameters, but the method needs to give the initial value of the estimated parameters in the process, and in the process of establishing the error equation, it needs to be Taylor-expanded and linearized, and there is a certain error in the theory.Therefore, in follow-up research work, it is necessary to establish a compensation parameter solution method that does not require a linearization process in order to improve the accuracy and practicability of parameter compensation.In the experimental process of internal and external calibration compensation, although this paper proposes a two-dimensional pointing angle model to replace the internal distortion model of the early warning area array camera, it fails to fully solve the correlation between internal and external parameters, resulting in unsatisfactory solution results.Therefore, in the process of on-orbit application, a more preferred parameter compensation model than the two-dimensional pointing angle model should be found to minimize the influence of strong correlation between internal and external parameters.
In general, the internal and external compensation model proposed in this paper had a good compensation effect on the target image point positioning error, which provides a certain reference for its on-orbit application.External parameter compensation under different temperature conditions was carried out, and the installation angle and time were established.The change model was consistent with the phase angle of the added installation error, which could also prove the validity and correctness of the internal and external parameter compensation model proposed in this paper.The compensation of the target image point positioning error caused by the existence of various errors in the imaging process was successfully carried out.The results of this paper can provide accurate trajectory reconstruction for the subsequent motion parameter estimation of high-value point targets.This paper focuses on the error compensation method of internal and external parameters of an area array camera, and the main error components were internal geometric distortion and installation angle error.Therefore, the error compensation method of internal and external parameters in the imaging model proposed in this paper can be used for all orbit altitudes and types of satellites.

Figure 1 .
Figure 1.Schematic diagram of ground virtual control points.

Figure 1 .
Figure 1.Schematic diagram of ground virtual control points.

Figure 2 .
Figure 2. Design of internal and external parameter compensation process.(a) Internal compensation before external compensation; (b) external compensation before internal compensation.

Figure 2 .
Figure 2. Design of internal and external parameter compensation process.(a) Internal compensation before external compensation; (b) external compensation before internal compensation.

Figure 3 .
Figure 3. Curve of relative installation angle and time change.Figure 3. Curve of relative installation angle and time change.

Figure 3 .
Figure 3. Curve of relative installation angle and time change.Figure 3. Curve of relative installation angle and time change.

Figure 4 .
Figure 4. Pointing angle difference diagram after camera distortion.(a) Pointing angle difference in X direction; (b) pointing angle difference in Y direction.

Figure 4 .
Figure 4. Pointing angle difference diagram after camera distortion.(a) Pointing angle difference in X direction; (b) pointing angle difference in Y direction.

Figure 5 .
Figure 5. Pointing angle difference diagram after internal parameter compensation.(a) Pointing angle difference in X direction; (b) pointing angle difference in Y direction.

Table 2 .
Relative installation angle error parameters.

Figure 5 .
Figure 5. Pointing angle difference diagram after internal parameter compensation.(a) Pointing angle difference in X direction; (b) pointing angle difference in Y direction.

Figure 5 .
Figure 5. Pointing angle difference diagram after internal parameter compensation.(a) Pointing angle difference in X direction; (b) pointing angle difference in Y direction.

Figure 7 .
Figure 7. Pointing angle difference diagram only after internal compensation scheme.(a) Pointing angle difference in X direction; (b) pointing angle difference in Y direction.

Figure 7 .
Figure 7. Pointing angle difference diagram only after internal compensation scheme.(a) Pointing angle difference in X direction; (b) pointing angle difference in Y direction.Figure 7. Pointing angle difference diagram only after internal compensation scheme.(a) Pointing angle difference in X direction; (b) pointing angle difference in Y direction.

Figure 7 .
Figure 7. Pointing angle difference diagram only after internal compensation scheme.(a) Pointing angle difference in X direction; (b) pointing angle difference in Y direction.Figure 7. Pointing angle difference diagram only after internal compensation scheme.(a) Pointing angle difference in X direction; (b) pointing angle difference in Y direction.

Figure 8 .Figure 9 .Figure 10 .
Figure 8. Pointing angle difference diagram only after external compensation scheme.(a) Pointing angle difference in X direction; (b) pointing angle difference in Y direction.

Figure 8 .Figure 8 .Figure 9 .Figure 10 .
Figure 8. Pointing angle difference diagram only after external compensation scheme.(a) Pointing angle difference in X direction; (b) pointing angle difference in Y direction.

Figure 9 .Figure 8 .Figure 9 .Figure 10 .
Figure 9. Pointing angle difference diagram of the compensation scheme of the internal first and then the external.(a) Pointing angle difference in X direction; (b) pointing angle difference in Y direction.

Figure 10 .
Figure 10.Pointing angle difference diagram of the compensation scheme of the external first and then the internal.(a) Pointing angle difference in X direction; (b) pointing angle difference in Y direction.

Figure 11 .
Figure 11.The relationship between the relative installation angle and time change of the two schemes.(a) Compensation result of relative installation angle α ; (b) compensation result of rela- tive installation angle β ; (c) compensation result of relative installation angle θ .

Figure 11 .
Figure 11.The relationship between the relative installation angle and time change of the two schemes.(a) Compensation result of relative installation angle α ; (b) compensation result of rela- tive installation angle β ; (c) compensation result of relative installation angle θ .

Figure 11 .
Figure 11.The relationship between the relative installation angle and time change of the two schemes.(a) Compensation result of relative installation angle α; (b) compensation result of relative installation angle β; (c) compensation result of relative installation angle θ.

Table 1 .
Target image simulation parameters.

Table 1 .
Target image simulation parameters.

Table 2 .
Relative installation angle error parameters.

Table 2 .
Relative installation angle error parameters.

Table 3 .
Positioning error after four compensation schemes.