#### Abstract

A novel method was presented for increasing the accuracy of subpixel centroid estimation for smearing star image. Model of the smearing trajectory of smearing star was built. It helped to study the analytical form of the errors, caused by image smearing, for centroid estimation. In the algorithm, the errors were estimated with accuracy and used to revise the centroid processed by CoM (centre of mass). Simulations have been run to study the effect of angular rates, integration time, and actual position of star on the accuracy of centroid estimation. Results were presented which suggested that the proposed algorithm had a precision better than 1/10 of a pixel when the angular rate was up to 3.0 deg/s.

#### 1. Introduction

A star sensor is a device which outputs pose for its carrier by measuring the directions of stars in the FOV (field of view). Those directions are compared with star directions stored in star catalog. Then pose of the carrier can be determined. Being one of the most accurate means for pose determination with accuracies down to arc seconds [1, 2], star sensor has been used in a variety of spacecrafts and experiments with its increasing accuracies.

In static condition angular rates of the star sensor are low, usually less than 0.3 deg/s. Intensity of a star imaged is often considered following 2D Gaussian distribution in engineering. Adopting interpolation or Gaussian centroid, the processed centroid can be better than 1/10 of a pixel [3]. Contrarily, angular rates are high in dynamic condition. Imaged stars are nonideal that their gray intensity does not follow 2D Gaussian distribution. Accuracy of the processed centroid might be worse than several pixels.

Many works have been done on processing centroid for smeared star. Xiaojuan and Xinlong [4] have described how the three-axis angular rates affect the centre determination of star spots. Weina et al. [5] have described a method to restore the smeared star image before processing centroid, while only the smear caused by velocities of cross-sight direction was analyzed. Image deconvolution could increase the ability to detect stars for a sensor while stars in the deconvolved image present a nonusual aspect (with several peaks) [6]. A gyrocontrol platform might be designed to eliminate the effect caused by angular movement of the spacecraft. But this scheme resulted in increased size, weight, and complexity of the star sensor. Adjusting the integration time to the variation of the angular rates can weaken the effect on processing centroid, while smear is inevitable and errors caused by the smearing still cannot be neglected. In this paper, it is hoped that the question will be resolved with a better way.

To overcome the difficulties arisen from smeared star image, a novel and effective method was presented in this paper, aiming at obtaining better centroid. In our case, the smearing trajectory of star was deduced firstly. Then, within the allowed precision it was approximated by Taylor’s series. After that, analytical expressions of errors, resulting from image smearing when CoM was applied, were deduced. The errors are related to the unknown actual position. Combining the relation between centroid by CoM and the actual position, the errors could be expressed by centroid processed by CoM. Lastly, a novel and efficient algorithm was proposed to obtain better centroid for smearing star image.

This paper consists of eight sections. Section 1 makes a brief introduction of the background and significance of the study. After this introduction, the mechanism of operation of a typical star sensor and its imaging geometrical model are outlined. Section 3 introduces the trajectory of smearing star. Section 4 further studies the mathematical model of the smearing star and its trajectory is approximated by polynomial from engineering perspective. In Section 5, errors of the processed centroid by CoM for smeared star are studied. Main steps of our method are drawn in Section 6. Simulations are constructed to verify the feasibility and efficiency of our method in Section 7. Lastly, conclusions are drawn in Section 8.

#### 2. Performance of Star Sensor

A star sensor usually consists of sensor system and data-processing system. It captures star images with its electronic camera. Afterward, a series of sequential implementations, including star image pretreatment, processing centroid, identification, and star pattern recognition are constructed. After that, it outputs the pose of the carrier.

The pose of carrier is often parameterized by a direct cosine matrix which satisfies the following equation: where and denote the unit vectors in the body-fixed frame and the inertial based coordinates system, respectively. For pose determination, the star sensor observes directional vectors of stars, while are derived from processed centroids. And then it obtains the reference vectors from star catalog by identifying star constellations. However, centroids without high precision would result in a worse and decrease the success rate of identifying star constellations. It is suggested that the much more accurate centroids might lay a solid foundation for star sensor to output more accurate pose of the spacecraft [1].

The imaging geometrical model of star sensor is drawn in Figure 1. All the analysis after this section is based on the coordinate system shown in it. , , and are the angular rates around -axis, *-*axis, and -axis, respectively. is the focal length and denotes the position of a star on the focal plane.

#### 3. Smearing Trajectory

In static condition, a star sensor obtains dispersive star spot of 3 by 3 to 7 by 7 pixels by defocusing method. Maximum angular rate that a star sensor can endure in order to output pose with required precision might be estimated briefly [7]. However, errors caused by image smearing cannot be ignored when the angular rate is two or three times the maximum one.

Dynamic variation of coordinate is related to and (Figure 1). and are the two components of the -axis line velocity. The synthesized line velocity about time can be described as follows:

In the same way, another differential expression about time could also be obtained:

Continuing to evaluate the second-order derivative of (2), a differential equation set can be derived:

Then, a differential equation about can be obtained:

Solve (5) and function can be expressed as

Similarly, function can be expressed as , , , and are constant coefficients of real numbers relying on initial conditions of the smeared star, when , , , , and . Those constant coefficients can be deduced when the initial conditions are combined with (6) and (7):

#### 4. Characteristics of Smearing Star

Trajectory of a smearing star has already been deduced above. Gray value of the smeared star in coordinate can be calculated in theory [8]: where and is the integration time. and are components of the trajectory on axis and axis, respectively.

When parameter is removed the smearing trajectory can be expressed only by parameters and : The trajectory is along a certain circle, its centre locates at and its radius is of .

Centre of a star almost keeps motionless on the focal plane during the integration time that is the case in static condition. However, its centre moves along the trajectory in dynamic condition which makes the star a trail imaged on the focal plane.

In reality, a star sensor does not output the pose during the two phases of propelled flight and entry because of the body vibration conducted from the engine and high speed spin of the spacecraft; its only functioning period is the coast phase. In this period, angular rates of the spacecraft can be regarded as constants during its integration time, which is as short as 200 ms, 100 ms, or less. This hypothesis is strictly valid for the case of midcourse of a missile in space, when the engine shuts down and no control moment is exerted.

Assuming that the initial condition of a given smeared star is known, its smearing trajectory can be expanded as the first-order Taylor series considering that is only a small amount: When the high-order terms are omitted, the trajectory can be approximated from engineering perspective:

This polynomial approximation plays a vital role in the following discussion and will not bring about notable errors when one of the main silicon array detectors is applied. In Section 7.1 the accuracy of this approximation will be analyzed through numerical simulation method.

#### 5. Errors of Centroid by CoM

Prior to the computation of centroid, pixels of smeared stars are extracted crudely from the noisy background. Then, size of centroid window for each detected star is determined and the centroid can be calculated by denotes the size of the centroid window, , . and represent coordinate of gray on axis and axis, respectively. denotes the gray value in and is the mass centre of the star calculated by CoM [3].

is the actual position of the given smeared star; it is the position where the star located at the beginning of the integration time and where the star should be in static condition. Relationship between and can be built when the errors resulting from image smearing on axis and axis are noted as and , respectively:

Movement of the centre of smearing star can be served as uniform motion in straight line approximately during the short integration time. can be taken as the location where centre of the star is when time is exactly at half of the integration time [7]. Figure 2 shows the position relation instructed above. Namely, is substituted by in (15) and relationship between and can be built:

**(a)**

**(b)**

is set up constantly and equation set (16) has the unique solution

Analytical expressions of the errors caused by image smearing for centroid estimation are deduced withhigh-order term (17). Four things become apparent from the results: the variation of and can affect the degree of the smear; a detector with suitable pixel size should be chosen to increase the performance of the star sensor in dynamic condition; position of a star in the FOV is also a weighty element, usually ignored, and might also affect its smearing degree; angular rate of the spacecraft might vary all the time during the flight and also contribute to the smearing of the star image as would be discussed later.

In the following section, a novel centroid method developed from the basic CoM will be proposed. It might process better centroid when star sensor works in high dynamic condition.

#### 6. Processing Centroid

is a small amount and is set up all the time. The high-order terms (in (17)) might be omitted from engineering perspective. Then the errors resulting from image smearing might be estimated: and denote the estimated values of and , respectively.

Although and cannot be ignored, they both can be estimated with high accuracy. In (18), and are sure for star sensor in real time. Usually, a star sensor has been calibrated on ground before it is launched with a spacecraft, so is sure in advance. Even though may drift away from the calibrated value, it can be estimated with high precision [9]. In addition, the angular rates are available by the embedded gyroscope or can be estimated by* Kalman filtering* [10–12].

A flowchart of the proposed centroid-estimation algorithm is shown in Figure 3. The subpixel centroid is calculated in the following way. Grays of the stars are extracted from star image firstly. At the same time, the number of stars detected is counted. Star tracker will not output the pose of the spacecraft when the number of detected stars is less than 3.* Num_star* denotes the number of stars detected in the FOV. Afterwards, centre of each star is estimated by CoM. Then is used to estimate the error caused by image smearing for each detected star. Lastly, is revised by .

#### 7. Results and Analysis

Digital values of light intensity for a smeared star can be calculated (9). Firstly, a mimic image of gray values resulting from the rotation of the spacecraft is obtained for a star. Afterwards, all the other stars in the FOV are simulated. Then, a smeared star image can be formed and all the smeared stars in it have an exact position. Furthermore, angular rates and are changeable and star images with different smearing degree can be simulated which might help to validate the efficiency of the proposed method.

All the simulations and experiments discussed herein are based on MATLAB 2010 in PC. Several parameters are set constantly. is set to 60 mm; spread of the PSF is 0.6 of a pixel, ensuring that more than 95% energy of the star is restricted within a 3 by 3 pixels area when it is imaged in static condition and pixel size is 15 by 15 m^{2}.

##### 7.1. Precision Analysis of Polynomial Approximation

Trajectory of a smearing star can be deduced when is zero in the same way as it was in Section 3. Namely, It shows that makes the trajectory a small section of an arc while and can only drive the star moving along a straight line strictly. As proves that is the factor which contributes to the errors resulting from the polynomial approximation (12).

Maximum value of the difference between analytical form of the trajectory and its approximated form (20) might depict the errors caused by the overlook of the high-order term. Namely, where and represent the error caused by the overlook of the high-order terms.

A simulation was run. was set to 100 ms. and both kept 3.0 deg/s. was the only factor changeable, varying from 1.0 deg/s to 20.0 deg/s. The smeared star located at (150.0000, 150.0000), which might help to eliminate the effect caused by varied position of the star.

Both and increased with the increase of and reached their peaks when increased to 20.0 deg/s (Figure 4). Actually, is far lower than 20.0 deg/s during functioning period of the star sensor. However, the peak was still lower than 1/20 of a pixel for such an incredible range, demonstrating that the polynomial approximations of the smearing trajectory are suitable from engineering perspective.

##### 7.2. Effects of Vital Elements on Processing Centroid

Smearing degree of an imaged star in dynamic condition is determined by many parameters. Actual position and integration time are important elements. In this section we will further study how those elements have effect on the accuracy of processed centroid.

Stars in different positions of the FOV will smear with difference. A simulation was constructed to study the relationship. All the angular rates kept 3.0 deg/s. was set to (50.0000, 50.0000), (100.0000, 100.0000), (150.0000, 150.0000), (200.0000, 200.0000), (250.0000, 250.0000), and (300.0000, 300.0000) successively, meaning that the distance between the simulated star and axis (Figure 1) became longer. CoM was used to process centroid and the calculated and were viewed in Figure 5.

The results indicate that the farther a star is away from axis the more serious the smearing of star will be. and alike reach about 4.5 pixels when the actual position is set at (300.0000, 300.0000). Hence, stars near the centre of the focal plane will have a better accuracy when CoM is applied.

Another simulation was run to study the accuracy when varied. It varied from 20 ms to 180 ms which met the range of the integration time that a star sensor spent to output the star image. Actual position was set at (150.0000, 150.0000) and all the angular rates were 1.0 deg/s.

Results of this simulation analysis are shown in Figure 6. The variation of and indicates that the bias errors, the absolute values of and , increase when becomes longer. For a certain star, the smearing length increases and gray intensity grows thinner. That is to say, energy of stars has the probability of being swallowed by the noisy background. Crude extraction of stars from the background will be a head-scratching problem. Clearly, should be suitable, ensuring the efficient accumulation of star energy by detector and decimating the effect of the noisy background on estimating centroid simultaneity.

##### 7.3. Estimate Centroid

The bias errors may reach several pixels when CoM is applied to a smeared star. A star sensor can hardly output the pose of the carrier in that situation. In this section, our method would be applied to smeared stars. Comparison between our method and CoM will be drawn to verify the feasibility of our method.

*Case 1. *Actual positions of smeared stars all were ; increased from 0.5 deg/s to 3.0 deg/s while and kept 1.0 deg/s. In the presence of white Gaussian noise, the SNR (signal to noise ratio) was 60 dB indicating a low noise level.

, , , , and the estimated centroid are listed in Table 1. With the increase of , keeps about −3.4 pixels while increases and reaches about 10.5 pixels when increases to 3.0 deg/s. However, our method has the accuracy better than 1/10 of a pixel even when increases to 3.0 deg/s which is the peak set in the simulation. Meanwhile, we can see that increases while keeps almost constant with the increase of . The variation of might cause the variation of while it might not be attributed to the variation of , just as it was described in (3).

*Case 2. * of smeared stars were set to (150.0000, 150.0000); increased from 0.5 deg/s to 3.0 deg/s while and were 1.0 deg/s constantly. The noise level was 60 dB.

When increases from 1.0 deg/s to 3.0 deg/s keeps about 3.4 pixels. But decreases from −3.4 pixels to −10.4 pixels, meaning that bias error in axis increases (Table 2). What also draws our attention is that the variation of might cause the variation of while it might not be attributed to the variation of (2). Our method shows that the results are better than 1/10 of a pixel even when reaches about 3.0 deg/s.

*Case 3. *Actual positions are all set to (150.0000, 150.0000); increases from 0.5 deg/s to 3.0 deg/s while the other two angular rates are 1.0 deg/s constantly. SNR was 60 dB.

The increase of has the least influence on the precision of estimation centroid (Table 3), because the focal length and the CCD array size are not of the same order of the magnitude, with the former one much higher than the latter one. When the three angular rates are at the same level contributes the least to the alteration of the line velocities in axis and axis directions ((2) and (3)).

In summary, smearing of star makes the obtained centroid worse than it should be in static condition. Precision of CoM for smeared star is even worse than 1.0 pixel, which will certainly result in the failure of pattern recognition of star image. However, our method provides a concise and efficient way to process the centroid. Revising the centroid, processed by CoM, with the corresponding errors estimated by (18), accuracy of the processed centroid can be better than 1/10 of a pixel for smearing stars.

#### 8. Conclusion

The smearing trajectory of star on the focal plane has been modeled. Within the allowed precision, the smearing trajectory has been approximated by polynomial from engineering prospective. Effects of actual position of star, integration time, pixel size, focal length, and angular rates on the accuracy of subpixel centroid estimation have been modeled in analytical form.

Also, simulations have been run to further study the effects of integration time, actual position, and angular rates on the accuracy of processed centroid. When the integration time becomes longer the length of smearing also gets longer. It is suggested that the estimated centroids of stars which are near the centre of the focal plane are proved to be less affected by the rotation of the star sensor around axis. Also, it is proved that among the three angular rates has the least effect on the bias errors; can only contribute to and can only contribute to .

By revising the centroid (processed by CoM) with the estimated errors, it can be better than 1/10 of a pixel in high dynamic condition, validated by simulation method. Moreover, our method will be useful in the field where high accuracy of pose measurement is needed, based on point targets, in high dynamic condition, such as adaptive optics, spatial acquisition systems for optical links, and laser triangulation sensor.

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgments

The authors thank Xiangdong Zhou and Wenming Zhang for their technical support. The support for this work was provided by the National Key Basic Research Program of China Grant no. 2014CB744200. Special thanks are to the anonymous referees for their thorough review and useful comments.