Self-Calibration Method Based on Surface Micromaching of Light Transceiver Focal Plane for Optical Camera

In remote sensing photogrammetric applications, inner orientation parameter (IOP) calibration of remote sensing camera is a prerequisite for determining image position. However, achieving such a calibration without temporal and spatial limitations remains a crucial but unresolved issue to date. The accuracy of IOP calibration methods of a remote sensing camera determines the performance of image positioning. In this paper, we propose a high-accuracy self-calibration method without temporal and spatial limitations for remote sensing cameras. Our method is based on an auto-collimating dichroic filter combined with a surface micromachining (SM) point-source focal plane. The proposed method can autonomously complete IOP calibration without the need of outside reference targets. The SM procedure is used to manufacture a light transceiver focal plane, which integrates with point sources, a splitter, and a complementary metal oxide semiconductor sensor. A dichroic filter is used to fabricate an auto-collimation light reflection element. The dichroic filter, splitter, and SM point-source focal plane are integrated into a camera to perform an integrated self-calibration. Experimental measurements confirm the effectiveness and convenience of the proposed method. Moreover, the method can achieve micrometer-level precision and can satisfactorily complete real-time calibration without temporal or spatial limitations.


Introduction
High-resolution earth observation applications are essential in many fields, such as mapping, environment monitoring, and resource exploration.Such applications require high-resolution images and high-accuracy image positioning determinations [1][2][3].However, high-accuracy and high-resolution optical imaging payloads need high attitude control precision, orientation, and transfer matrix of attitude [4][5][6][7].The imaging payloads' positioning accuracy determines the spacecraft performance, which is crucial for completing space missions.The GeoEye-1 satellite can collect panchromatic imagery with a resolution of 0.46 m and multispectral imagery with a resolution of 1.84 m; the system can also precisely locate an object within 5 m of its true location on the surface of the Earth with CE90 (Circular Error of 90%) [8], owing to the high-performance electro-optical camera developed by the International Telephone & Telegraph Space Systems Division [9].The Pleiades-1B can acquire panchromatic (0.5 m) and multispectral (2 m) images with high geo-location accuracy (3 m @ CE90) and high resolution [10].The system also employs a push-broom imaging concept [11].NAOMI [12] is a product line of high-resolution optical cameras applied to SPOT-6 satellite and is based on the same telescope concept implementing one or several focal plane units in the same camera.The high-accuracy space smart payload system integrated with attitude and position (SSPIAP) can provide high resolution (0.5 m) and accurate location (3.7 m @ CE90) of images [13].The system also integrates a high-resolution remote camera and miniature attitude-and position-sensitive devices, such as star tracker, micro-electro-mechanical system gyroscope, and global position system, into a new smart payload system.Inner orientation parameter (IOP) calibration is a prerequisite for image position determination of a camera.Given that remote sensing optical camera is much needed for high-accuracy position determination, its accuracy calibration must be high.
To complete the high-accuracy positioning of remote sensing camera, high imaging quality and high calibration accuracy of IOPs are needed.Numerous technologies have been proposed to improve the imaging quality of remote sensing cameras [14,15].Specifically, a special optical hardware, such as wave front sensor [16], can determine the quality of wave fronts in front of detector array to correct dynamic optical aberrations that blur images.Advanced technologies can ensure high imaging quality for remote sensing cameras.Camera calibration methods have been proposed for different applications, such as IOP calibration methods for space optical camera [17][18][19][20][21][22][23][24][25].Furthermore, laboratory measuring angle methods are widely used for IOP calibration of remote sensing cameras [26][27][28][29][30][31].In this type of simple and micrometer-level calibration method, an uncalibrated optical camera is placed on a precision turntable in capturing the parallel lights of star points emitted by a collimator.Computer vision techniques [14,32,33] are used for camera calibration based on three dimensional (3D) coordinates of several given and image points; this type of method uses multiple views to calibrate intrinsic parameters, distortions, and image deformations.In 2001, multiple two dimensional (2D) diffraction grids are generated using crossed-phase diffractive optical elements [34].The equally spaced dots can be used to calibrate wide-angle geometric cameras for photogrammetry.Moreover, these methods can be precisely calibrated in a laboratory; however, they may drift to a large degree because of vibration during complicated launches and on-orbit working conditions.Furthermore, laboratory calibration approaches generally need a calibrated reference object.To mediate this problem, several self-calibration methods that do not require a calibrated reference object have been developed for calculating IOPs [35][36][37][38][39][40][41].Self-calibration methods can use unfamiliar scenes and motions to calibrate a camera; moreover, these methods use constraints among the system parameters for calibration.Nevertheless, these self-calibration methods exhibit computational complexity, heavy computation, and nontrivial solutions of equations.Moreover, calibration accuracy cannot be guaranteed and its robustness is low.
Traditional methods for on-orbit calibration of remote sensing cameras use ground control point (GCP) or stars control point (SCP) methods [42,43].However, the availability and access to GCPs is sometimes difficult.Moreover, operational SCP implementation is highly restrictive because acquiring star images wastes a long orbit portion before completing the attitude change of a space camera pointing to between Earth and space.Therefore, self-calibrating bundle adjustment is ideal for mediating control point limitations for on-orbit camera calibrations [44].Lichti et al. [45] compared three geometric self-calibration methods for ranged cameras, and the self-calibration bundle adjustment is found to be slightly superior.However, the self-calibrating bundle adjustment method presents long computation times [46].Moreover, 180 • satellite maneuvers have been used to calibrate on-orbit cameras.Self-calibrating bundle adjustment needs a standard ground calibration field [47].Delvit et al. [48] proposed an auto-reverse method during the commissioning phase.This method is efficient and does not require external reference data; however, its operational implementation is highly restrictive because acquiring a same-site image pair wastes along orbit portion before completing the alignment to ground projection of scan-line on ground velocity.Therefore, existing on-ground and on-orbit methods present common characteristics with the aid of external targets and are limited by space and time.
We develop an efficient self-calibration method for optical cameras.The proposed method considers the temporal and spatial limitations when conducting the calibration function based on surface micromachining (SM) of light transceiver focal plane.This method is expected to be highly accurate, easy to operate, require few auxiliary equipment, and replicate real orbit conditions.The proposed auto-collimating self-calibration method is based on a bi-plane dichroic filter element and an SM point-source focal plane (SM-PSFP).The on-orbit camera can perform integrated self-calibration without temporal and spatial limitations and is therefore a potential calibration standard.Moreover, this method can verify the functions of the on-orbit status of optical camera and can be used on-ground and on-orbit.The rest of the paper is organized as follows.The principle of the proposed method is introduced in Sections 2 and 3, and the experimental results are presented in Section 4.

Proposed Self-Calibration Method
Positioning without GCPs is among the key technical problems for remote sensing photogrammetry.The IKonos-2 satellite can reach a positioning accuracy of 15 m without GCPs [49].The WorldView-2 satellite can reach a positioning accuracy of 6.5 m without GCPs [50].The GeoEye-1 satellite can provide a positioning accuracy of 4 m without GCPs [51].The SSPIAP also adopts a positioning method without GCPs and requires a positioning accuracy of 5 m.Given that a camera generally uses forward and back-looking images to complete the digital mapping, the geometric performance of the camera must be accurate.To extract highly accurate topographical information from two overlapping strips of images, a camera must provide highly accurate IOPs.After adjusting the optical camera, the IOPs deviate from the ideal values specified by the design, manufacture, and assembly [52][53][54].Therefore, accurate IOP calibration is necessary.The positioning accuracy (without GCPs) of a camera mainly depends on the accuracies of satellite station positioning, attitude measurement, image point measurement, and IOP measurement.Moreover, the error distribution of a remote sensing camera is obtained by the positioning algorithm.Without GCPs, the error distribution equation can be built based on the position of homologous pixel points of forward-and back-looking images, IOP, and a series of coordinate transforms from camera to ground.The error distribution can generally be divided into two aspects: (1) camera imaging with accuracy depending on pixel point measurement error and IOP measurement accuracy; and (2) position and attitude accuracies at forwardand back-looking imaging times.Combined with device performance of satellite, error distribution equations can calculate the optimal error distribution.For example, a positioning accuracy of 5 m can obtained by a camera if the distribution of its primary errors exhibits the following characteristics: (1) an attitude determination accuracy within 10 (angular seconds); (2) a precise orbit determination within 0.2 m; (3) an angle calibration accuracy between the star tracker and the optical camera within 5 ; (4) a camera lens distortion calibration accuracy within 5 µm; (5) a principal distance calibration accuracy within 50 µm; and (6) a principal point calibration accuracy within a third of a pixel.After adjusting the camera, all of the parameters need to be calibrated.In this paper, we mainly discuss the on-orbit method of the IOPs, including principal distance and principal point, when a camera is working on a satellite.We propose an integrated self-calibration method that can be used on-ground and on-orbit.
A bi-plane auto-collimation dichroic filter (BADF) and an SM-PSFP are integrated into a camera to calibrate IOPs.The principle of the proposed auto-collimating calibration of IOPs for camera is shown in Figure 1.The optical system of the camera is a co-axial, three-mirror structure, including a primary mirror, a secondary mirror, a tertiary mirror, a folding mirror, a focusing mirror, and a focal plane.When the remote sensing camera works on satellite, the lights of ground target reflection and radiation pass through the dichroic filter and optical system.These lights are then sensed by the charge-coupled device (CCD) sensor to complete on-orbit imaging.In Figure 1a denotes the auto-collimation dichroic filter, which is plated on a bi-plane substrate located near the secondary mirror to fabricate an auto-collimation light reflection element.In Figure 1b denotes the focal plane of the remote sensing camera.In the focal plane B, multiple CCD sensors use a stagger arrangement strategy (Figure 1b).Two conjugated SM-PSFPs are installed on the interlace area on focal plane assembly B.Moreover, C denotes one of the SM-PSFPs.In the interlace area, the point-source photomask, beam splitter prism, and image sensor are integrated into SM-PSFP C (Figure 1d).When the LED of the SM-PSFP is lighted, the lights pass through the photomask to produce multiple point sources.The splitter prism turns the lights of the point sources by 90°.In other words, the point sources outgo lights from the sensor surface of the SM-PSFP on the focal plane.The gray and color scales of an equivalent point source are shown in Figure 1d.These scales show a Gaussian distribution.Figure 1c shows the equivalent path.The auto-collimating light path includes two point sources, a camera lens, an auto-collimation dichroic filter, and an image detector.According to the optical path of camera, the SM point sources are auto-collimating lights when they pass through the focusing, folding, tertiary, primary, and secondary mirrors of camera lenses.Moreover, the auto-collimating lights are reflected by the auto-collimation dichroic filter and return to the camera lens.Outgoing lights of the camera are incidents on the focal plane detector.
In the proposed method, lighting elements, a mask, beam splitters, and a detection module are integrated into an SM focal plane (SMFP), which is introduced into the optical system of the camera.These lights are then sensed by the charge-coupled device (CCD) sensor to complete on-orbit imaging.In Figure 1a denotes the auto-collimation dichroic filter, which is plated on a bi-plane substrate located near the secondary mirror to fabricate an auto-collimation light reflection element.In Figure 1b denotes the focal plane of the remote sensing camera.In the focal plane B, multiple CCD sensors use a stagger arrangement strategy (Figure 1b).Two conjugated SM-PSFPs are installed on the interlace area on focal plane assembly B.Moreover, C denotes one of the SM-PSFPs.In the interlace area, the point-source photomask, beam splitter prism, and image sensor are integrated into SM-PSFP C (Figure 1d).When the LED of the SM-PSFP is lighted, the lights pass through the photomask to produce multiple point sources.The splitter prism turns the lights of the point sources by 90 • .In other words, the point sources outgo lights from the sensor surface of the SM-PSFP on the focal plane.The gray and color scales of an equivalent point source are shown in Figure 1d.These scales show a Gaussian distribution.Figure 1c shows the equivalent path.The auto-collimating light path includes two point sources, a camera lens, an auto-collimation dichroic filter, and an image detector.According to the optical path of camera, the SM point sources are auto-collimating lights when they pass through the focusing, folding, tertiary, primary, and secondary mirrors of camera lenses.Moreover, the auto-collimating lights are reflected by the auto-collimation dichroic filter and return to the camera lens.Outgoing lights of the camera are incidents on the focal plane detector.
In the proposed method, lighting elements, a mask, beam splitters, and a detection module are integrated into an SM focal plane (SMFP), which is introduced into the optical system of the camera.The SMFP is installed on the interleaving area between two CCD sensors of the focal plane assemblies.A dichroic filter element is also integrated into the optical system of the camera.The calibration optical path can be designed based on the optical system of the camera.The dichroic filter element can be installed on the supporting truss of secondary mirror of the optical system.Based on the size of the truss of the secondary mirror, the size of a dichroic filter element can be determined without extra light obstruction.The reflection angle of the dichroic filter element is designed based on the field of view of the camera.The proposed method can be summarized in four steps.Figure 2 shows the logical chart of the proposed method.The details of these steps are discussed below.The SMFP is installed on the interleaving area between two CCD sensors of the focal plane assemblies.A dichroic filter element is also integrated into the optical system of the camera.The calibration optical path can be designed based on the optical system of the camera.The dichroic filter element can be installed on the supporting truss of secondary mirror of the optical system.Based on the size of the truss of the secondary mirror, the size of a dichroic filter element can be determined without extra light obstruction.The reflection angle of the dichroic filter element is designed based on the field of view of the camera.The proposed method can be summarized in four steps.Figure 2 shows the logical chart of the proposed method.The details of these steps are discussed below.In the first step, an auto-collimation dichroic filter is plated on a bi-plane glass substrate.The dichroic filter can selectively admit light in a small range of bands while reflecting light in other bands.Figure 3 shows the reflection and transmission ratios of the dichroic filter for the visible LED.The dichroic filter allows one of the wavelength bands to pass through while reflecting other wavelength bands.On the basis of the details of the dichroic filter, the bands of SM point sources can be determined.The fabricated dichroic filter element is installed on the truss of secondary mirror of the optical camera; the filter does not obstruct the light reflected from and radiated off a target when the camera is working on a satellite.The lights of the SM point sources are reflected by the dichroic filter to complete an on-orbit calibration.Moreover, the angle between the two planes of the bi-plane dichroic filter element is α, which is set by the user.The user-defined angle determines the calibration resolution of principal distance of camera.The calibration resolution can be expressed as where ds is the centroid extraction accuracy of SM point images.Following the requirements of calibration accuracy for the remote sensing camera, Equation (1) can determine the angle of the two planes of the dichroic filter element.For example, the extraction accuracy of our centroid extraction algorithm is less than 0.1 pixels and the pixel size of the image sensor is 6 µm.When α ≥ 0.35°, the designed calibration accuracy is less than 0.05 mm.Therefore, the angle between the two planes of the dichroic filter element must be larger than 0.35°.In the first step, an auto-collimation dichroic filter is plated on a bi-plane glass substrate.The dichroic filter can selectively admit light in a small range of bands while reflecting light in other bands.Figure 3 shows the reflection and transmission ratios of the dichroic filter for the visible LED.The dichroic filter allows one of the wavelength bands to pass through while reflecting other wavelength bands.On the basis of the details of the dichroic filter, the bands of SM point sources can be determined.The fabricated dichroic filter element is installed on the truss of secondary mirror of the optical camera; the filter does not obstruct the light reflected from and radiated off a target when the camera is working on a satellite.The lights of the SM point sources are reflected by the dichroic filter to complete an on-orbit calibration.Moreover, the angle between the two planes of the bi-plane dichroic filter element is α, which is set by the user.The user-defined angle determines the calibration resolution of principal distance of camera.The calibration resolution can be expressed as where ds is the centroid extraction accuracy of SM point images.Following the requirements of calibration accuracy for the remote sensing camera, Equation (1) can determine the angle of the two planes of the dichroic filter element.For example, the extraction accuracy of our centroid extraction algorithm is less than 0.1 pixels and the pixel size of the image sensor is 6 µm.When α ≥ 0.35 • , the designed calibration accuracy is less than 0.05 mm.Therefore, the angle between the two planes of the dichroic filter element must be larger than 0.35 • .In the second step, we fabricate an SM-PSFP and then integrate it into the optical focal plane of camera.The SM point sources are installed on the focal plane.This arrangement indicates that the monitoring optical path is auto-collimating.Thus, principal distance and principal point can be monitored as needed.To ensure a sufficiently small size and low power consumption, we use SM procedures to fabricate a point-source focal plane.The point-source focal plane is mainly composed of the fabricated mask by the SM process, a housing, and an electrical system.The assembly of the point-source focal plane is shown in Figure 4.The fabrication process of the mask can be summarized into two steps: (1) chromium, gold, and tantalum materials are plated on a specified glass substrate; and (2) photoetching is performed on the metal layers to obtain several small apertures.In the second step, we fabricate an SM-PSFP and then integrate it into the optical focal plane of camera.The SM point sources are installed on the focal plane.This arrangement indicates that the monitoring optical path is auto-collimating.Thus, principal distance and principal point can be monitored as needed.To ensure a sufficiently small size and low power consumption, we use SM procedures to fabricate a point-source focal plane.The point-source focal plane is mainly composed of the fabricated mask by the SM process, a housing, and an electrical system.The assembly of the point-source focal plane is shown in Figure 4.The fabrication process of the mask can be summarized into two steps: (1) chromium, gold, and tantalum materials are plated on a specified glass substrate; and (2) photoetching is performed on the metal layers to obtain several small apertures.In the second step, we fabricate an SM-PSFP and then integrate it into the optical focal plane of camera.The SM point sources are installed on the focal plane.This arrangement indicates that the monitoring optical path is auto-collimating.Thus, principal distance and principal point can be monitored as needed.To ensure a sufficiently small size and low power consumption, we use SM procedures to fabricate a point-source focal plane.The point-source focal plane is mainly composed of the fabricated mask by the SM process, a housing, and an electrical system.The assembly of the point-source focal plane is shown in Figure 4.The fabrication process of the mask can be summarized into two steps: (1) chromium, gold, and tantalum materials are plated on a specified glass substrate; and (2) photoetching is performed on the metal layers to obtain several small apertures.In the third step, light emitting diodes (LEDs), an image sensor, a beam splitter prism, and a mask are packaged into a point-source focal plane.The integrated packaging with a complementary metal oxide semiconductor (CMOS) sensor is shown in Figure 5. LEDs are installed under the mask.Moreover, light can pass through the etched apertures.Given that other parts are covered by a chromium layer, they block the incidental light.The beam splitter can emit lights from point sources to the surface of the CMOS sensor and receive lights into the CMOS sensor.In this way, the light can transceive on the same focal plane.The SM point sources can be controlled by the camera controller in real time.In the imaging mode of the camera, the LEDs can be switched off without effecting the image.The SM light sources are installed on the side of the image sensors.For the optical remote sensor, several image sensors are butted together into one greater image sensor [55][56][57], and the SM focal planes are generally placed on the butting area.We use two SMFPs to build a differential calibration structure.In the first focal plane, one point source emits red light and the other emits blue light.
In the third step, light emitting diodes (LEDs), an image sensor, a beam splitter prism, and a mask are packaged into a point-source focal plane.The integrated packaging with a complementary metal oxide semiconductor (CMOS) sensor is shown in Figure 5. LEDs are installed under the mask.Moreover, light can pass through the etched apertures.Given that other parts are covered by a chromium layer, they block the incidental light.The beam splitter can emit lights from point sources to the surface of the CMOS sensor and receive lights into the CMOS sensor.In this way, the light can transceive on the same focal plane.The SM point sources can be controlled by the camera controller in real time.In the imaging mode of the camera, the LEDs can be switched off without effecting the image.The SM light sources are installed on the side of the image sensors.For the optical remote sensor, several image sensors are butted together into one greater image sensor [55][56][57], and the SM focal planes are generally placed on the butting area.We use two SMFPs to build a differential calibration structure.In the first focal plane, one point source emits red light and the other emits blue light.Two SMFPs are used to build a differential calibration structure.Based on the differential SMFPs, we build a mathematical calibration model.Following the geometrical relationships between point source positions and their images, the mathematical calibration equation is repeatedly solved to calculate the IOPs of the camera.
In the fourth and final step, a camera controller manipulates the SM light sources, and the CMOS senses the SM source images.A centroid extraction algorithm processes input images to extract the star point positions.In our method, we use the centroid extraction algorithm based on mathematical morphology to capture the position information.Finally, the principal distance and principal point can be calculated using the mathematical calibration model.

Calibration Calculation Based on Differential SMFPs
We adopt two SMFPs to build a differential structure for calculating calibration.Two conjugate SMFPs are symmetrically installed on the focal planes of optical system of the camera.In the red SMFP, lights emitted by point sources pass through beam splitters and into the optical system of the Two SMFPs are used to build a differential calibration structure.Based on the differential SMFPs, we build a mathematical calibration model.Following the geometrical relationships between point source positions and their images, the mathematical calibration equation is repeatedly solved to calculate the IOPs of the camera.
In the fourth and final step, a camera controller manipulates the SM light sources, and the CMOS senses the SM source images.A centroid extraction algorithm processes input images to extract the star point positions.In our method, we use the centroid extraction algorithm based on mathematical morphology to capture the position information.Finally, the principal distance and principal point can be calculated using the mathematical calibration model.

Calibration Calculation Based on Differential SMFPs
We adopt two SMFPs to build a differential structure for calculating calibration.Two conjugate SMFPs are symmetrically installed on the focal planes of optical system of the camera.In the red SMFP, lights emitted by point sources pass through beam splitters and into the optical system of the camera.
Lights from the primary mirror become collimated light because the SMFP is located on the focal plane of the camera.The collimated lights are reflected by one of the dichroic filter planes into the primary mirror and through the other mirrors of the camera system.Thereafter, lights from the focusing mirror pass through the source prism and approach the CMOS of the red SMFP.The principle of the blue SMFP is the same with that of the red focal plane.
We consider an ideal optical system.The imaging relationship of the equivalent optical paths of self-calibration in an ideal optical system is shown in Figure 6.The auto-collimating light path of S 1 is composed of an SM point source S 1 , a camera optical system, an auto-collimation dichroic filter, and an image detector.We use OXYZ to represent the camera coordinate system.In the OXYZ, the origin is the center of the optical system of the remote sensing camera.The origin of the coordinate system could be taken as a point O(0 0 0) T .We define z-axis as the optical axis, and its positive direction is from left to right; y-axis is perpendicular to the z-axis, and its positive direction is from the bottom going up; x-axis is perpendicular to the paper surface, and its positive direction is vertically shooting out from the paper.The image plane is denoted by I and locates f units along optic axis, where f is the focal length of optical system.The SM point sources S 1 and S 2 and receiver are located on the focal plane of optical camera I.The position vectors of two point sources S 1 and S 2 are defined as T , where h 1 and h 2 are the distances between the SM point sources S 1 and S 2 and optical axis, respectively.The distance between the SM point source S 2 and optical axis is h 2 .
The normal vectors of the direction vector of dichroic filter are defined as n 1 ( 0 −sinα cosα ) T and n 2 ( 0 sinα cosα ) T , where α is the rotated angle of dichroic filter around the x-axis; that is, the angle between dichroic filter and y-axis.The direction vector of two point sources S 1 and S 2 are defined as l 1 ( 0 sinα −cosα ) T and l 2 ( 0 −sinα −cosα ) T .
Remote Sens. 2016, 8, 893 8 of 21 camera.Lights from the primary mirror become collimated light because the SMFP is located on the focal plane of the camera.The collimated lights are reflected by one of the dichroic filter planes into the primary mirror and through the other mirrors of the camera system.Thereafter, lights from the focusing mirror pass through the source prism and approach the CMOS of the red SMFP.The principle of the blue SMFP is the same with that of the red focal plane.
We consider an ideal optical system.The imaging relationship of the equivalent optical paths of self-calibration in an ideal optical system is shown in Figure 6.The auto-collimating light path of S1 is composed of an SM point source S1, a camera optical system, an auto-collimation dichroic filter, and an image detector.We use OXYZ to represent the camera coordinate system.In the OXYZ, the origin is the center of the optical system of the remote sensing camera.The origin of the coordinate system could be taken as a point O(0 0 0) T .We define z-axis as the optical axis, and its positive direction is from left to right; y-axis is perpendicular to the z-axis, and its positive direction is from the bottom going up; x-axis is perpendicular to the paper surface, and its positive direction is vertically shooting out from the paper.The image plane is denoted by I and locates f units along optic axis, where f is the focal length of optical system.The SM point sources S1 and S2 and receiver are located on the focal plane of optical camera I.The position vectors of two point sources S1 and S2 are defined as , where h1 and h2 are the distances between the SM point sources S1 and S2 and optical axis, respectively.The distance between the SM point source S2 and optical axis is h2.The normal vectors of the direction vector of dichroic filter are defined as According to the optical path of camera, the prime light ray emitted from point source S1 is autocollimating light when it passes through the camera lenses.Images of the two conjugate point sources are located on S1′ and S2′.The auto-collimating lights are reflected by the auto-collimation dichroic filter first and then return into the camera lens.Outgoing lights from the camera are incidents on the focal plane detector.Given that the angle between the prime light ray and optical axis is the same as the normal of dichroic filter, the optical path is auto-collimated; therefore, the position of SM point source and its image coincides.The intersection points of the incoming rays with the dichroic filter are: where ( , ) ⋅ ⋅ denotes the vector scalar products.The direction 1 l′  and 1 l  of the reflected ray is: According to the optical path of camera, the prime light ray emitted from point source S 1 is auto-collimating light when it passes through the camera lenses.Images of the two conjugate point sources are located on S 1 and S 2 .The auto-collimating lights are reflected by the auto-collimation dichroic filter first and then return into the camera lens.Outgoing lights from the camera are incidents on the focal plane detector.Given that the angle between the prime light ray and optical axis is the same as the normal of dichroic filter, the optical path is auto-collimated; therefore, the position of SM point source and its image coincides.The intersection points of the incoming rays with the dichroic filter are: where (•, •) denotes the vector scalar products.The direction l 1 and l 1 of the reflected ray is: Finally, we obtain the equation of the outgoing ray where k 1 and k 2 are the system parameters.Intersecting the ray with the image plane parallel to the XOY plane at the distance f from the origin and rewriting Equations ( 6) and ( 7) explicitly we get According to Equations ( 8) and ( 9), the equation of the optical path is as follows: where S i is the position vector of point source S i , I = 1 or 2. In an ideal optical system, the principal distance is its focal length.Let their respective images S 1 and S 2 in the focal plane coordinate system be located at x S 1 , y S 1 and x S 2 , y S 2 .The principal point can then be expressed as: When the optical system of a camera exhibits maladjustments, its prime planar position, principal point, and principal distance change accordingly.Considering the influences of different maladjustments of optical components on the first-order parameters of optical cameras, the change in the axial maladjustments of an optical component can vary the primary plane position and principal distance.Moreover, the deviation of an optical component from its original position in a vertical direction can vary the primary plane position, principal point, and principal distance.Furthermore, the slant of an optical component from its original position in a vertical direction can vary the primary plane position, principal point, and principal distance.During on-orbit imaging, a space camera estimates image quality and determines the adjustment operation of the focal length.Focusing can change the installed planes of the SM point sources and image sensor.The optical path when optical system of a camera presents maladjustments during focusing is shown in Figure 7.In Figure 7, yellow lines represent the original auto-collimating lights while red lines represent the auto-collimating lights of optical system with maladjustments.In the auto-collimating light path, the SMFPs are located on position I of optical system C without maladjustments.The lights emitted from the SM point sources S1 and S2 on the symmetrical SMFPs pass through the camera lens and dichroic filters F1 and F2, and these lights return to the positions of S1′ and S2′.The positions of S1′ and S2′ are coincidental with those of S1 and S2, respectively.The relative position vector of the image points of two SMFPs is defined as . When the optical camera is maladjusted, the SMFP is located on position I′ after focusing.In Figure 7, C′ represents the maladjusted optical system.The lights emitted from the symmetrical SM point sources S1 and S2 on I′ pass through the camera lens and dichroic filter, and these lights return to the positions of S1′ and S2′ on I′.The relative position vector of the image points of two SMFPs is defined as . Considering that the autocollimating light paths of S1 and S2 present the same spread rules, we first analyze the lights emitted from S1.The imaging relationship of S1 is shown in Figure 8.In Figure 8, yellow lines are original auto-collimating light paths.Red lines are optical paths of the maladjusted optical system.I′ is the focal plane of the maladjusted optical system.When the focal plane is located on position I for the optical system without maladjustments, the direction vector of outgoing light ray from S1 is 1 (0 . The origin of the system C could be taken as a point O(0 0 0) T .The position vector of the image point of S1 for the optical system without maladjustments In Figure 7, yellow lines represent the original auto-collimating lights while red lines represent the auto-collimating lights of optical system with maladjustments.In the auto-collimating light path, the SMFPs are located on position I of optical system C without maladjustments.The lights emitted from the SM point sources S 1 and S 2 on the symmetrical SMFPs pass through the camera lens and dichroic filters F 1 and F 2 , and these lights return to the positions of S 1 and S 2 .The positions of S 1 and S 2 are coincidental with those of S 1 and S 2 , respectively.The relative position vector of the image points of two SMFPs is defined as L ( 0 L 0 ) T .When the optical camera is maladjusted, the SMFP is located on position I after focusing.In Figure 7 In Figure 7, yellow lines represent the original auto-collimating lights while red lines represent the auto-collimating lights of optical system with maladjustments.In the auto-collimating light path, the SMFPs are located on position I of optical system C without maladjustments.The lights emitted from the SM point sources S1 and S2 on the symmetrical SMFPs pass through the camera lens and dichroic filters F1 and F2, and these lights return to the positions of S1′ and S2′.The positions of S1′ and S2′ are coincidental with those of S1 and S2, respectively.The relative position vector of the image points of two SMFPs is defined as . When the optical camera is maladjusted, the SMFP is located on position I′ after focusing.In Figure 7, C′ represents the maladjusted optical system.The lights emitted from the symmetrical SM point sources S1 and S2 on I′ pass through the camera lens and dichroic filter, and these lights return to the positions of S1′ and S2′ on I′.The relative position vector of the image points of two SMFPs is defined as . Considering that the autocollimating light paths of S1 and S2 present the same spread rules, we first analyze the lights emitted from S1.The imaging relationship of S1 is shown in Figure 8.In Figure 8, yellow lines are original auto-collimating light paths.Red lines are optical paths of the maladjusted optical system.I′ is the focal plane of the maladjusted optical system.When the focal plane is located on position I for the optical system without maladjustments, the direction vector of outgoing light ray from S1 is 1 (0 . The origin of the system C could be taken as a point O(0 0 0) T .The position vector of the image point of S1 for the optical system without maladjustments In Figure 8, yellow lines are original auto-collimating light paths.Red lines are optical paths of the maladjusted optical system.I is the focal plane the maladjusted optical system.When the focal plane is located on position I for the optical system without maladjustments, the direction vector of outgoing light ray from S 1 is d 1 ( 0 −sinα −cosα ) T .The origin of the system C could be taken as a point O(0 0 0) T .The position vector of the image point of S1 for the optical system without maladjustments → S d1 is: When the focal plane is located on position I for the optical system with maladjustments, the direction vector of outgoing light ray from S1 is d 1 ( 0 −sinα 1 −cosα 1 ) T .The position vector T of point source S1 can be calculated by Equation (10).For the optical system with maladjustments, the deviation vector between the principal point and its original optical axis is ∆ h 0 ∆h 0 T .In Figure 8, the origin of the system C could be taken as a point O'(0 ∆h 0)T.For optical system with long focal length, the distance ∆h is small and can be neglected.The relationship between point source S1 and direction vector d 1 can be expressed as The normal vector of the direction vector of dichroic filter is n F 1 ( 0 sinα cosα ) T .The direction vector of the incoming rays on dichroic filter is e ( 0 −sinα 1 −cosα 1 ) T .The angle of incidence between the incidence lights and dichroic filter plane is ε = α − α1.The angle between outgoing rays and z-axis is β = ε + α = 2α − α1.The direction vector of the outgoing rays from dichroic filter is Given that plane I is the focal plane of the current maladjusted focal plane, the lights emitted from point source S 1 are parallel when they pass through the camera system.Therefore, one of the rays reflected by the dichroic filter element passes through the center of the prime plane first and then return into S 1 .The direction vector of the outgoing rays from prime plane F 1 is c 1 ( 0 sinα 1 cosα 1 ) T .
The position vector of S 1 is S 1 ( 0 h 1 f ) T .The following relationship can be satisfied as Using the same method, the above imaging relationship of point source S 2 can be obtained.The variation vector of relative position of image points of two SMFPs can be expressed as where c 2 ( 0 sinα 2 cosα 2 ) T , which is the direction vector of the outgoing rays from primary plane T is the direction vector of outgoing light ray from S 2 when the focal plane is located on position I for the optical system without maladjustments.Substituting Equations ( 13)- (15) into Equation (16), a new equation is obtained as: Solving the above equation, the mathematical model of principal distance can be expressed as: where α and h i are the designed values.Using the same principle, the principal point can be expressed as: The variation ∆l is the distance in the centroid of point images of SMFP.The point images are captured by two CMOS sensors on SMFP.The centroid extraction algorithm is used to calculate the centroid of point images.On the basis of variation ∆l, principal distance and principal point can be calculated by Equations ( 18) and (19).

Experiment of Optical System
To verify the effectiveness of the proposed method, we use ZEMAX in simulating the calibration experiments.ZEMAX presents high accuracy of ray tracing.For a maladjusted optical system, the coordinate positions from rays with different fields of view (FOV) positions to the image plane can be traced accurately.In each experiment, the reference values of principal distance and principal point are first calculated.We use the on-ground calibration method based on measuring angles to calibrate the reference values of principal distance and principal point.The laboratory measuring angle method is widely used in calibrating IOPs of remote sensing cameras [20][21][22][23][24][25].In this method, an uncalibrated optical camera is placed on a precision turntable to capture the parallel lights of star points emitted by a collimator.At different rotation angles of the turntable, the camera can obtain multiple positions of star points and their images from different FOV positions.The recorded rotation angle and image point positions are used to build an imaginary geometrical equation between each star point and its image.According to the principle of camera distortion, an optical camera has its minimum distortion at its principal point location.A least squares method is used to solve the imaging geometrical equation for calculating principal distance and principal point.This method is simple and can achieve a micrometer-level calibration accuracy.Similar to the ground calibration method of principal distance and principal point, we use the least squares and multiple regression analyses to calculate the principal distance and principal point of each test and maladjusted system.Principal distance and principal point can be expressed as follows: where f is the principal distance, p x/y is the position of the principal point in the x direction or the y direction, i is the number of measurement points, W i is the measurement angle of the ith measurement point, and L i is the measurement height of the ith measurement point.
In the first step of experimentation, we use an ideal optical system to perform modeling because the ideal optical system does not exhibit aberrations and can adequately verify our method.In ZEMAX, the focal length and aperture diameter of our input optical system model are 4500 and 600 mm, respectively.Moreover, a mirror with a deflecting angle of 0.7 • is used to simulate the dichroic filter.The ideal optical path is shown in Figure 9.We set different maladjustments for S 2 and use two different methods to calculate the principal distance of different maladjusted optical systems.In our simulation, we use 10 reference points at different FOV positions to estimate principal distances and principal points for the on-ground method.We set several maladjustments of S2 along the z direction to test the proposed method.The deviations of the S2 mirror from its original position are from 0.010 mm to 0.05 mm.Then, we set maladjustments along the y direction and around the x-and y-axes.Along the y direction, we set two maladjustments of S2 to test the proposed method.The deviations of the S2 mirror from its original position are 0.020 mm along the x and y directions.Thereafter, we set several maladjustments of S2 around the x-axis to test the proposed method.The deviations of the S2 mirror from its original position are 30 s of arc around the x-and y-axes.The calculated results of the principal distance in the simulated conditions with different maladjustments of optical system are shown in Tables 1 and  2. In Table 2, Δl is the variation of relative positions of the two SMFPs (Equation ( 16)).Compared with the error in the on-ground method, the error in our method is less than 0.02 µm.The error is attributed to the round-off error, ignoring ∆h, and focusing.In the second step, we use an actual optical system in performing the modeling.The actual optical system presents aberration and is thus used to further verify our method.Our input optical system model parameters areas follows: the focal length is 8000 mm, the aperture diameter is 600 mm, the relative aperture is 1/13.3,FOV is 1.8°, and deflected FOV is 0.578°.Moreover, the co-axial, threemirror-anastigmatic optical system with deflected FOV is adopted.The optical system is composed We set several maladjustments of S 2 along the z direction to test the proposed method.The deviations of the S 2 mirror from its original position are from 0.010 mm to 0.05 mm.Then, we set maladjustments along the y direction and around the xand y-axes.Along the y direction, we set two maladjustments of S 2 to test the proposed method.The deviations of the S 2 mirror from its original position are 0.020 mm along the x and y directions.Thereafter, we set several maladjustments of S 2 around the x-axis to test the proposed method.The deviations of the S 2 mirror from its original position are 30 s of arc around the xand y-axes.The calculated results of the principal distance in the simulated conditions with different maladjustments of optical system are shown in Tables 1 and 2. In Table 2, ∆l is the variation of relative positions of the two SMFPs (Equation ( 16)).Compared with the error in the on-ground method, the error in our method is less than 0.02 µm.The error is attributed to the round-off error, ignoring ∆h, and focusing.In the second step, we use an actual optical system in performing the modeling.The actual optical system presents aberration and is thus used to further verify our method.Our input optical system model parameters areas follows: the focal length is 8000 mm, the aperture diameter is 600 mm, the relative aperture is 1/13.3,FOV is 1.8 • , and deflected FOV is 0.578 • .Moreover, the co-axial, three-mirror-anastigmatic optical system with deflected FOV is adopted.The optical system is composed of three aspherical and two spherical mirrors.These mirrors build an optical path of the secondary imagery.On the first image plane, a folding mirror is placed and rotated along the light direction by 90 • .Based on this optical system, we design a calibrated optical path.The dichroic filter that provides angle information is installed on the truss of the secondary mirror.Based on the size of the truss and the absence of extra obstruction, the size of the dichroic filter is designed to 40 mm × 40 mm under the condition of no additional blocking incident rays.The FOV and deflected angle of dichroic filter are designed based on the space camera.he parameters of the dichroic filter component are shown Table 3.The calibrated optical path designed based on the auto-collimating method and designed parameters is shown in Figure 10.The optical rays are reflected from the mirror, passes through the optical system twice, and then concentrates on the CMOS detector. of three aspherical and two spherical mirrors.These mirrors build an optical path of the secondary imagery.On the first image plane, a folding mirror is placed and rotated along the light direction by 90°.Based on this optical system, we design a calibrated optical path.The dichroic filter that provides angle information is installed on the truss of the secondary mirror.Based on the size of the truss and the absence of extra obstruction, the size of the dichroic filter is designed to 40 mm × 40 mm under the condition of no additional blocking incident rays.The FOV and deflected angle of dichroic filter are designed based on the space camera.he parameters of the dichroic filter component are shown Table 3.The calibrated optical path designed based on the auto-collimating method and designed parameters is shown in Figure 10.The optical rays are reflected from the mirror, passes through the optical system twice, and then concentrates on the CMOS detector.The optical system is composed of primary, secondary, tertiary, folding, and focusing mirrors.The maladjustments of the system can be attributed to axial maladjustment, vertical axis decentration, and slanting.By analyzing the maladjustments of the different elements of the optical system, we determine that the secondary and tertiary mirror maladjustments significantly affect the focal plane of the optical system.Given that the tertiary mirror is installed in the third mirror room, nearly no maladjustment occurs in the mirror.The highest effect is from the secondary mirror because it is installed on the truss of the camera (Figure 11).In actual on-orbit work, the truss of camera is easily changed by platform vibration.Therefore, the largest effect on the IOPs of the camera comes from the secondary mirror.Moreover, we set different maladjustments for the secondary mirror.Using the same approach, we set 10 reference points in a valid FOV and set maladjustments along the z direction.Furthermore, we set several maladjustments of S2 along the z direction to test the proposed method.The deviations of the S2 mirror from its original position are from 0.002 mm to −0.02 mm.We set maladjustments along the y direction and around the x-and y-axes.Along the y direction, we set two maladjustments of S2 to test the proposed method.The deviations of the S2 mirror from its original position are both 0.020 mm along the x and y directions.We set several maladjustments of S2 around the x-axis to test the proposed method.The deviations of the S2 mirror from its original position are both 5s of arc around the x-and y-axes.The calculated results of the principal distance in the simulation conditions with different maladjustments of optical system are shown in Tables 4  and 5. Compared with that in the on-ground method, the error in our method is less than 0.0148 mm.The optical system is composed of primary, secondary, tertiary, folding, and focusing mirrors.The maladjustments of the system can be attributed to axial maladjustment, vertical axis decentration, and slanting.By analyzing the maladjustments of the different elements of the optical system, we determine that the secondary and tertiary mirror maladjustments significantly affect the focal plane of the optical system.Given that the tertiary mirror is installed in the third mirror room, nearly no maladjustment occurs in the mirror.The highest effect is from the secondary mirror because it is installed on the truss of the camera (Figure 11).In actual on-orbit work, the truss of camera is easily changed by platform vibration.Therefore, the largest effect on the IOPs of the camera comes from the secondary mirror.Moreover, we set different maladjustments for the secondary mirror.Using the same approach, we set 10 reference points in a valid FOV and set maladjustments along the z direction.Furthermore, we set several maladjustments of S 2 along the z direction to test the proposed method.The deviations of the S 2 mirror from its original position are from 0.002 mm to −0.02 mm.We set maladjustments along the y direction and around the xand y-axes.Along the y direction, we set two maladjustments of S 2 to test the proposed method.The deviations of the S 2 mirror from its original position are both 0.020 mm along the x and y directions.We set several maladjustments of S 2 around the x-axis to test the proposed method.The deviations of the S 2 mirror from its original position are both 5s of arc around the xand y-axes.The calculated results of the principal distance in the simulation conditions with different maladjustments of optical system are shown in Tables 4 and 5. Compared with that in the on-ground method, the error in our method is less than 0.0148 mm.

Experiment of Setup in Lab
In the first step of experimentation, we discretely set up the experimental system components to verify the proposed method.Figure 12 shows the experimental calibration system.The system includes an optical camera, a reflective mirror, horizontal rotation adjustment, 2D displacement platform, an LED, and a CMOS.The F/ratio of the optical system is 11.8 and the aperture diameter is 127 mm.

Experiment of Setup in Lab
In the first step of experimentation, we discretely set up the experimental system components to verify the proposed method.Figure 12 shows the experimental calibration system.The system includes an optical camera, a reflective mirror, horizontal rotation adjustment, 2D displacement platform, an LED, and a CMOS.The F/ratio of the optical system is 11.8 and the aperture diameter is 127 mm.To verify the proposed method, we use the focusing knob to set three specific positions.We use a known line-pair method to calibrate the principal distance at three positions as reference values.We then perform our method along the three specific positions.Figure 13 shows the captured images and calibrated results at the three positions.The centroids of the three positions are (310.76,251.12), (242.43,320.70), and (124.76,378.99).Compared with the error in the on-ground method, the error in our method is less than 2 mm.The calibration error is large and mainly due to that the tooling fixtures do not achieve the requirements.Therefore, we use the SM method to further improve calibration accuracy.In the second step, we use SMFP to set up an experimental system for the integrated calibration with imaging to verify the proposed method.Figure 14 shows the experimental calibration system.The system includes an optical camera, an auto-collimating filter, a processing circuit, a collimator, an optical theodolite, and a high-accuracy turntable.The optical system is a co-axial Schmidt-Cassegrain optical system.The aperture diameter is 202.6 mm, the focal length is 1026 mm, and the F/ratio of the optical system is 10.The image sensor is a CMOS detector, and the image resolution is 1280 pixels × 1024 pixels.The SM point sources, beam splitter, and image sensor are installed on the focal plane.A collimator provided an infinite target for the test system.To calibrate the reference value, the three-axis turntable and collimator are adjusted evenly using a level.By adjusting the support tooling of the camera and using the benchmark prisms of the optical axis of camera, the visual axis of camera and collimator are moved to share a common shaft.The camera controller sets the imaging mode.Star points of the collimator are imaged on the target CMOS sensor.The turntable To verify the proposed method, we use the focusing knob to set three specific positions.We use a known line-pair method to calibrate the principal distance at three positions as reference values.We then perform our method along the three specific positions.Figure 13 shows the captured images and calibrated results at the three positions.The centroids of the three positions are (310.76,251.12), (242.43,320.70), and (124.76,378.99).Compared with the error in the on-ground method, the error in our method is less than 2 mm.The calibration error is large and mainly due to that the tooling fixtures do not achieve the requirements.Therefore, we use the SM method to further improve calibration accuracy.To verify the proposed method, we use the focusing knob to set three specific positions.We use a known line-pair method to calibrate the principal distance at three positions as reference values.We then perform our method along the three specific positions.Figure 13 shows the captured images and calibrated results at the three positions.The centroids of the three positions are (310.76,251.12), (242.43,320.70), and (124.76,378.99).Compared with the error in the on-ground method, the error in our method is less than 2 mm.The calibration error is large and mainly due to that the tooling fixtures do not achieve the requirements.Therefore, we use the SM method to further improve calibration accuracy.In the second step, we use SMFP to set up an experimental system for the integrated calibration with imaging to verify the proposed method.Figure 14 shows the experimental calibration system.The system includes an optical camera, an auto-collimating filter, a processing circuit, a collimator, an optical theodolite, and a high-accuracy turntable.The optical system is a co-axial Schmidt-Cassegrain optical system.The aperture diameter is 202.6 mm, the focal length is 1026 mm, and the F/ratio of the optical system is 10.The image sensor is a CMOS detector, and the image resolution is 1280 pixels × 1024 pixels.The SM point sources, beam splitter, and image sensor are installed on the focal plane.A collimator provided an infinite target for the test system.To calibrate the reference value, the three-axis turntable and collimator are adjusted evenly using a level.By adjusting the support tooling of the camera and using the benchmark prisms of the optical axis of camera, the visual axis of camera and collimator are moved to share a common shaft.The camera controller sets the imaging mode.Star points of the collimator are imaged on the target CMOS sensor.The turntable In the second step, we use SMFP to set up an experimental system for the integrated calibration with imaging to verify the proposed method.Figure 14 shows the experimental calibration system.The system includes an optical camera, an auto-collimating filter, a processing circuit, a collimator, an optical theodolite, and a high-accuracy turntable.The optical system is a co-axial Schmidt-Cassegrain optical system.The aperture diameter is 202.6 mm, the focal length is 1026 mm, and the F/ratio of the optical system is 10.The image sensor is a CMOS detector, and the image resolution is 1280 pixels × 1024 pixels.The SM point sources, beam splitter, and image sensor are installed on the focal plane.A collimator provided an infinite target for the test system.To calibrate the reference value, the three-axis turntable and collimator are adjusted evenly using a level.By adjusting the support tooling of the camera and using the benchmark prisms of the optical axis of camera, the visual axis of camera and collimator are moved to share a common shaft.The camera controller sets the imaging mode.Star points of the collimator are imaged on the target CMOS sensor.The turntable is revolved, and the rotation angle and capture image are recorded.The processing circuit output stars coordinates in real time.Based on the measured centroid position and rotation angle, the least squares and two-multiplication regression analyses are used to obtain the optimal estimation values of the IOPs.For our method, the camera controller sets the calibration mode and switches the SM point sources on or off.The SM point sources are lighted when the camera controller performs the turn-on command.Figure 15 shows the captured images when the two point sources are lighted.SM optical sources are an ideal point source and can be used in measuring internal parameters.Using the spot centroid algorithm, the position of images can be determined.Principal distance and principal point are calculated using Equations ( 18) and (19).We adjust the motion of the secondary mirror to simulate on-orbit maladjustments of the optical system.From analyzing the maladjustments of the different elements of the optical system, we observe that the secondary and primary mirror maladjustments affect the focal plane of the optical system.Given that the primary mirror is installed in the primary mirror room, it exhibits nearly no maladjustment.Figure 16 shows the principal distance variation under different maladjustments.We use the ground calibration method based on the least squares and multiple regression analyses as For our method, the camera controller sets the calibration mode and switches the SM point sources on or off.The SM point sources are lighted when the camera controller performs the turn-on command.Figure 15 shows the captured images when the two point sources are lighted.SM optical sources are an ideal point source and can be used in measuring internal parameters.Using the spot centroid algorithm, the position of images can be determined.Principal distance and principal point are calculated using Equations ( 18) and ( 19).For our method, the camera controller sets the calibration mode and switches the SM point sources on or off.The SM point sources are lighted when the camera controller performs the turn-on command.Figure 15 shows the captured images when the two point sources are lighted.SM optical sources are an ideal point source and can be used in measuring internal parameters.Using the spot centroid algorithm, the position of images can be determined.Principal distance and principal point are calculated using Equations ( 18) and (19).We adjust the motion of the secondary mirror to simulate on-orbit maladjustments of the optical system.From analyzing the maladjustments of the different elements of the optical system, we observe that the secondary and primary mirror maladjustments affect the focal plane of the optical system.Given that the primary mirror is installed in the primary mirror room, it exhibits nearly no maladjustment.Figure 16 shows the principal distance variation under different maladjustments.We use the ground calibration method based on the least squares and multiple regression analyses as We adjust the motion of the secondary mirror to simulate on-orbit maladjustments of the optical system.From analyzing the maladjustments of the different elements of the optical system, we observe that the secondary and primary mirror maladjustments affect the focal plane of the optical system.Given that the primary mirror is installed in the primary mirror room, it exhibits nearly no maladjustment.Figure 16 shows the principal distance variation under different maladjustments.We use the ground calibration method based on the least squares and multiple regression analyses as reference.The calibration deviation is less than 0.035 mm.The deviation is attributed to the self-method error between methods under the maladjusted condition.reference.The calibration deviation is less than 0.035 mm.The deviation is attributed to the selfmethod error between methods under the maladjusted condition.To test the monitoring accuracy of our method with an approximation to the on-orbit environment, our experiment is performed in a laboratory with constant temperature.The experimental turntable uses a gas-floating vibration isolation platform to avoid vibration disturbance.We use our method to monitor the variation in principal distance and principal point in the static case.We process tens of thousands of images to calculate the monitoring accuracy.The variations in principal distance and principal point position in 2000 real-time seconds are shown in Figure 17.Based on the statistical data in Figure 17, the mean square error formula is used to calculate the monitoring accuracy.The monitoring accuracy of the principal distance can reach 0.017 mm.Furthermore, the monitoring accuracy of the principal point can reach 0.005 and 0.0049 mm in the x and y directions, respectively.For example, the camera exhibits a 5 m image positioning accuracy; therefore, the camera requires less than 50 µm calibration accuracy of principal distance and principal point.Moreover, the monitoring accuracy can reach micrometer level and meet the mapping requirements of camera.To test the monitoring accuracy of our method with an approximation to the on-orbit environment, our experiment is performed in a laboratory with constant temperature.The experimental turntable uses a gas-floating vibration isolation platform to avoid vibration disturbance.We use our method to monitor the variation in principal distance and principal point in the static case.We process tens of thousands of images to calculate the monitoring accuracy.The variations in principal distance and principal point position in 2000 real-time seconds are shown in Figure 17.Based on the statistical data in Figure 17, the mean square error formula is used to calculate the monitoring accuracy.The monitoring accuracy of the principal distance can reach 0.017 mm.Furthermore, the monitoring accuracy of the principal point can reach 0.005 and 0.0049 mm in the x and y directions, respectively.For example, the camera exhibits a 5 m image positioning accuracy; therefore, the camera requires less than 50 µm calibration accuracy of principal distance and principal point.Moreover, the monitoring accuracy can reach micrometer level and meet the mapping requirements of camera.To test the monitoring accuracy of our method with an approximation to the on-orbit environment, our experiment is performed in a laboratory with constant temperature.The experimental turntable uses a gas-floating vibration isolation platform to avoid vibration disturbance.We use our method to monitor the variation in principal distance and principal point in the static case.We process tens of thousands of images to calculate the monitoring accuracy.The variations in principal distance and principal point position in 2000 real-time seconds are shown in Figure 17.Based on the statistical data in Figure 17, the mean square error formula is used to calculate the monitoring accuracy.The monitoring accuracy of the principal distance can reach 0.017 mm.Furthermore, the monitoring accuracy of the principal point can reach 0.005 and 0.0049 mm in the x and y directions, respectively.For example, the camera exhibits a 5 m image positioning accuracy; therefore, the camera requires less than 50 µm calibration accuracy of principal distance and principal point.Moreover, the monitoring accuracy can reach micrometer level and meet the mapping requirements of camera.

Conclusions
In this work, we propose the SMPF and BADF methods for constructing an on-orbit calibration method of camera with high accuracy.The BADF and SMPF methods are integrated into a camera.First, the point sources are installed on the focal plane, the SM method is used to fabricate point sources, and these point sources are packaged with the image sensor and beam splitter.Second, we integrate the auto-collimation dichroic filter into the optical system of the camera.Third, a mathematical model of IOPs is built based on a geometrical imaging model.Fourth, the centroid extraction algorithm is used to process images to extract the star point positions and calculate the IOPs.Finally, we use ZEMAX to simulate the proposed method and set up an experiment for verifying the feasibility of our method with micrometer-level monitoring accuracy.The proposed method can complete self-calibration without temporal and spatial limitations in real time.Our method can also be applied to improve the performance of other calibration methods.

Figure 1 .
Figure 1.Principle of self-calibration: (a) the internal structure of remote sensing camera; (b) the optical plane of remote sensing camera; (c) the equivalent light path of dichroic filter; and (d) the light transceiver focal plane for calibrating IOPs.

Figure 1 .
Figure 1.Principle of self-calibration: (a) the internal structure of remote sensing camera; (b) the optical plane of remote sensing camera; (c) the equivalent light path of dichroic filter; and (d) the light transceiver focal plane for calibrating inner orientation parameters (IOPs).

Figure 2 .
Figure 2. Logical chart of the proposed method.

Figure 2 .
Figure 2. Logical chart of the proposed method.

Figure 3 .
Figure 3. Principle of bi-plane dichroic filter and its filtering performance.

Figure 4 .
Figure 4. SM process for fabricating the point-source mask.

Figure 3 .
Figure 3. Principle of bi-plane dichroic filter and its filtering performance.

Figure 3 .
Figure 3. Principle of bi-plane dichroic filter and its filtering performance.

Figure 4 .
Figure 4. SM process for fabricating the point-source mask.

Figure 4 .
Figure 4. Surface micromachining (SM) process for fabricating the point-source mask.

Figure 5 .
Figure 5. Integrated packaging of the SM point-source focal plane.

Figure 5 .
Figure 5. Integrated packaging of the SM point-source focal plane.

Figure 6 .
Figure 6.Imaging relationship of an ideal optical system. α

Figure 7 .
Figure 7. Imaging relationship of actual optical system.

Figure 8 .
Figure 8. Optical path of one of the differential structures.

Figure 7 .
Figure 7. Imaging relationship of actual optical system.

Figure 8 .
Figure 8. Optical path of one of the differential structures.

Figure 8 .
Figure 8. Optical path of one of the differential structures.

Figure 9 .
Figure 9. Optical path of the ideal system.

Figure 9 .
Figure 9. Optical path of the ideal system.

Figure 10 .
Figure 10.Optical path simulation of the actual system.

Figure 10 .
Figure 10.Optical path simulation of the actual system.

Figure 13 .
Figure 13.Calibrated results of the three specific positions.

Figure 13 .
Figure 13.Calibrated results of the three specific positions.

Figure 13 .
Figure 13.Calibrated results of the three specific positions.
Remote Sens. 2016, 8, 893 17 of 21 is revolved, and the rotation angle and capture image are recorded.The processing circuit output stars coordinates in real time.Based on the measured centroid position and rotation angle, the least squares and two-multiplication regression analyses are used to obtain the optimal estimation values of the IOPs.

Figure 15 .
Figure 15.Captured image (a); and (b-d) gray scale, color scale, and 3D plots, respectively, of the same detailed view of one point.
Remote Sens. 2016, 8, 893 17 of 21 is revolved, and the rotation angle and capture image are recorded.The processing circuit output stars coordinates in real time.Based on the measured centroid position and rotation angle, the least squares and two-multiplication regression analyses are used to obtain the optimal estimation values of the IOPs.

Figure 15 .
Figure 15.Captured image (a); and (b-d) gray scale, color scale, and 3D plots, respectively, of the same detailed view of one point.

Figure 15 .
Figure 15.Captured image (a); and (b-d) gray scale, color scale, and 3D plots, respectively, of the same detailed view of one point.

Figure 16 .
Figure 16.(a) Calibrated results of reference method; (b) the calibrated values of this method; and (c) the deviation between two different methods.

Figure 16 .
Figure 16.(a) Calibrated results of reference method; (b) the calibrated values of this method; and (c) the deviation between two different methods.
Remote Sens. 2016, 8, 893 18 of 21reference.The calibration deviation is less than 0.035 mm.The deviation is attributed to the selfmethod error between methods under the maladjusted condition.

Figure 16 .
Figure 16.(a) Calibrated results of reference method; (b) the calibrated values of this method; and (c) the deviation between two different methods.

Table 1 .
Simulation calculation results using ground method.

Table 2 .
Simulation calculation results of the proposed method.

Table 1 .
Simulation calculation results using ground method.

Table 2 .
Simulation calculation results of the proposed method.

Table 3 .
Parameters of the dichroic filter component.

Table 3 .
Parameters of the dichroic filter component.

Table 4 .
Simulation calculation results using ground method.

Table 5 .
Simulation calculation results using the proposed method.

Table 4 .
Simulation calculation results using ground method.

Table 5 .
Simulation calculation results using the proposed method.