Detection of Degraded Star Observation Using Singular Values for Improved Attitude Determination

This study introduces an innovative approach aimed at enhancing the accuracy of attitude determination through the computation of star observation quality. The proposed algorithm stems from the inherent invariance of singular values under attitude transformations, leveraging the concept of assessing error magnitude through the deviation of singular values. Quantization becomes imperative to employ this error magnitude as a weighting factor within the attitude determination process. To fulfill this purpose, this study applies p-value hypothesis testing to calculate quantized error levels. Simulation results validate that the calculated weights derived from the proposed algorithm lead to a discernible enhancement in attitude determination performance.


Introduction
Accurate attitude determination is an essential requirement for the success of a spacecraft in achieving the specified mission objectives [1].Since even a minor deviation in the spacecraft's attitude can lead to significant consequences, more precise attitude determination is being emphasized to mitigate the risk of potential bottlenecks that could limit the success of space missions.Among the various attitude determination sensors available to satisfy this requirement, star trackers stand out as key devices that provide the highest level of accuracy in on-orbit attitude determination [2,3].
After the star is projected onto the image plane of the star tracker, it goes through a series of workflows: centroid extraction, star identification from a star catalog, and attitude determination [4].For precise attitude determination, the extracted center position of the star from the image must be accurate [5].However, with the recent advancements in maneuvering performance, many satellites have gained the potential for agile maneuvers [6].These maneuvers cause image smearing during the exposure time of the star tracker [7], leading to a decrease in the accuracy of the centroid error [8].This effect is heightened by increased maneuvering speeds [9].Such a phenomenon is particularly noticeable in dim stars with lower signal-to-noise ratios when compared to bright stars [10].Several studies have been proposed to attain comparable performance in such situations [11][12][13][14].However, there inevitably exists a degradation in the accuracy of the centroid error, with the magnitude of this degradation varying from star to star based on their respective signal-to-noise ratios [15].
The signal-to-noise ratio is related to the brightness of the star.Nevertheless, owing to the ambiguity in magnitude information [16,17], it should not be depended upon to determine the extent of the error.Alongside signal-to-noise ratio, noise characteristics such as distortion, image degradation due to temperature, and starlight deflection can also introduce varying degrees of noise among different stars.Until now, several studies have presented strategies for utilizing star sensors in situations involving degradation.However, the predominant focus has been on strengthening the robust detection [18] and identification [8,19] of stars.As a result, the need arises for a method to quantify the varying levels of accuracy among stars and to incorporate this information into attitude determination for enhanced precision.
Accurate attitude determination also necessitates a robust star identification process.The occurrence of false matches during identification is another factor that diminishes the accuracy and reliability of attitude determination.While research in star identification has made significant strides, recent advancements have introduced algorithms that employ pattern-based feature extraction methods [20].It is widely acknowledged that patternbased recognition algorithms perform optimally when a substantial number of stars are present within an image frame, ensuring the distinctiveness of each pattern [21,22].While a larger quantity stars can indeed enhance the stability of matching performance, the inclusion of dim stars can intensify the centroid error due to their inherently low signalto-noise ratios.Hence, it becomes crucial to achieve not only precise centroid error but also robust star identification to ensure accurate attitude determination.Nonetheless, these two prerequisites present conflicting challenges, which unavoidably result in variations in accuracy among individual stars.
Therefore, this study proposes an innovative approach to assess the quality of star observations preceding the attitude determination, with the aim of enhancing accuracy.The primary contribution lies in the computation of individual star observation quality by utilizing the singular values among stars.Singular values have also been adopted in star identification techniques [23,24] and calibration [25] due to their insensitivity to attitude changes, making them suitable as reference values for determining the degree of error in star observations.To quantify this, we introduce p-value hypothesis testing to derive a quantized error level.In the end, this value can serve as a weighting factor in the attitude determination process, thereby contributing to the improvement of accuracy.
The organization of this paper can be summarized as follows.Section 2 furnishes the requisite background information relevant to this study.Moving ahead, Section 3 outlines the procedure for computing the quality level based on singular values, which constitutes the fundamental essence of this research.In Section 4, we validate the accuracy of the proposed method in quantifying the degree of error through simulation results and further analyze its practicality by applying it to attitude determinations.Lastly, Section 5 concludes by summarizing the main findings of the research presented in this paper.

Invariance of Singular Values
In this section, our objective is to offer concise insights to readers who might not be acquainted with the attitude-invariant properties of singular values and attitude determination.Singular value decomposition techniques [26] have found extensive application in domains like signal and image processing [27], as well as attitude estimation [28], enabling the extraction of valuable insights from correlation matrices.In fact, this methodology is renowned as one of the most resilient algorithms in the realm of numerical algebra [29].
A vector observation in a star tracker represents the relationship between the unit vector r i of the i th star, presented in the inertial frame, and the unit vector b i presented in the body frame, as illustrated in Figure 1.This relationship can be mathematically expressed through a matrix operation that rotates an arbitrary vector from the inertial frame to the body frame, which is defined as where C B I represents the directional cosine matrix (DCM) containing the attitude information in the body frame relative to the inertial frame.The subscript 'I' and the superscript 'B' signify the inertial and body frames, respectively.Notably, this matrix is characterized by its property of being a unitary matrix that conforms to the following expression: where  × represents an identity matrix of size 3 × 3.If N stars are observed from the star tracker, these operations can be extended using  .The vector matrices for  stars observed in the inertial frame and the body frame are denoted as  × and  × , respectively, with subscripts indicating that the matrices have sizes of 3 × N.
The observation vector matrices  × and  × are defined as follows: These  × and  × matrices can be factorized through singular value decomposition.
where  and  are 3 × 3 orthogonal matrices, representing the normalized eigenvectors of  ×  × and  ×  × , respectively. and  are 3 × N diagonal matrices, with their components being the singular values  and  of  × and  × . and  are N×N orthogonal matrices, denoting the normalized eigenvectors of  ×  × and  ×  × , respectively.Since a star's unit vector signifies its orientation in threedimensional space, if more than three stars are observed ( ≥ 3), both matrices  × and  × will have a rank of 3. Therefore, when observing more than three stars, the computation of three non-zero singular values will always be feasible.Notably, this matrix is characterized by its property of being a unitary matrix that conforms to the following expression: where I 3×3 represents an identity matrix of size 3 × 3.If N stars are observed from the star tracker, these operations can be extended using C B I .The vector matrices for N stars observed in the inertial frame and the body frame are denoted as R 3×N and B 3×N , respectively, with subscripts indicating that the matrices have sizes of 3 × N.
The observation vector matrices R 3×N and B 3×N are defined as follows: These R 3×N and B 3×N matrices can be factorized through singular value decomposition.
where P r and P b are 3 × 3 orthogonal matrices, representing the normalized eigenvectors of R 3×N R 3×N T and B 3×N B 3×N T , respectively.D r and D b are 3 × N diagonal matrices, with their components being the singular values σ r i and σ b i of R 3×N and B 3×N .Q r and Q b are N×N orthogonal matrices, denoting the normalized eigenvectors of R 3×N T R 3×N and B 3×N T B 3×N , respectively.Since a star's unit vector signifies its orientation in threedimensional space, if more than three stars are observed (N ≥ 3), both matrices R 3×N and B 3×N will have a rank of 3. Therefore, when observing more than three stars, the computation of three non-zero singular values will always be feasible.
By multiplying Equation ( 4) with the transpose of B 3×N , we can reformulate it as Equation (9).This transformation exhibits the characteristics of a similarity transformation, a consequence of the unitary property of the directional cosine matrix C B I .
The relationship between the eigenvalues σ r i and σ b i is conserved throughout the similarity transformation in Equation (9), allowing us to deduce or This suggests that the singular values of R 3×N and B 3×N remain unaltered, regardless of any shifts in coordinates resulting from attitude transformations.

Attitude Determination
Attitude can be determined using N known nonlinear vectors, which are obtained from the vectors in the inertial frame along with their corresponding observations in the body frame [30].This indicates the N star observation vector matrices R 3×N and B 3×N mentioned earlier.Attitude determination plays a crucial role in satellite control and adopts an algebraic approach for reconstructing attitude.The primary objective of attitude determination algorithms is to minimize a loss function, known as Wahba's problem, as illustrated in the equation below [31].
Here, k i represents the confidence level of the i th observation with a non-negative weight.To date, one of the most widely adopted methods for addressing Wahba's problem is the QUEST (QUaternion ESTimator) algorithm [32].Originating from the Davenport q-method [33], which offered an initial solution to Wahba's problem, the QUEST algorithm has the capacity to identify the optimal solution for Wahba's problem when presented with two or more vector observations.
The loss function enables us to formulate the subsequent gain function.As such, it becomes a problem of finding the optimal solution that maximizes this gain function.
The gain function is derived through the following trace operation: Here, 'tr' represents a trace operation.The matrix M and the normalized weight w i are defined as The attitude C B I can be represented using the corresponding quaternion Q = q 0 , q T T as follows.Here, [q] × refers to skew-symmetric matrices with the vector part of the quater- nion Q.
By substituting the above expression into the previous one, the gain function is summarized in the subsequent bilinear form: where K is the following matrix with a size of 4 × 4.
The optimal quaternion Q opt , obtained via the QUEST algorithm, corresponds to the eigenvector coupled with the largest eigenvalue of the matrix K.
In terms of w i , it is important to note that existing works solely mention the constraint that the set of weights sums to one, without any loss of generality.Therefore, we embark on our approach by introducing the concept of enhancing attitude accuracy by reflecting the individual confidence of the observations into the weights w i .
In summary, the invariance of singular values to attitude variations can serve as a reference value.The proposed work suggests an innovative approach by computing the divergence between the singular value derived from observed stars and employing it as a weighting factor in the attitude determination process, thereby enhancing accuracy.This novel aspect stands as the primary contribution of the proposed work, and the following section provides a detailed explanation of this approach.

Variation of Singular Values
Before proceeding with the calculation of the quality value, we examined whether singular values indeed remain unaltered with respect to attitude variations and assessed the magnitude of changes induced in singular values when there is image noise.For the simulations, synthesized star images are generated using the well-known Hipparcos catalog as reference star catalog [34].For each simulation, random attitude is provided to emulate the star tracker's boresight direction to view various directions within the celestial sphere.
The field of view, pixel size, and resolution of the star tracker are assumed to be 12 degrees, 10 −5 m, and 512 by 512, respectively.We introduced an error in the center position of the star projected onto the image plane and varied its magnitude to examine the effects on the singular values.The introduced error follows a random Gaussian noise distribution, and Figure 2 displays the average change after performing 10,000 simulations for each error size.distribution, and Figure 2 displays the average change after performing 10,000 simulations for each error size.
The three calculated singular values are ordered, with the first singular value defined as the one with the largest magnitude, followed by the remaining singular values in descending order of their magnitudes.It is noticeable that the change in singular value tends to increase linearly with the magnitude of the noise.This result suggests that quality can be evaluated by analyzing the error size introduced into the observations based on the differences in singular values.
Moreover, it is noticeable that the alteration in the first singular value is significantly smaller compared to the other two singular values.It is crucial to understand the covariance of singular value changes under normal sensor conditions, as it is a key factor in calculating the quality value, which will be discussed in a later section.In this study, the error under normal star tracker conditions is assumed to be 10 arc-seconds, and the changes and standard deviations of each singular value at that time are summarized in Table 1.The three calculated singular values are ordered, with the first singular value defined as the one with the largest magnitude, followed by the remaining singular values in descending order of their magnitudes.It is noticeable that the change in singular value tends to increase linearly with the magnitude of the noise.This result suggests that quality can be evaluated by analyzing the error size introduced into the observations based on the differences in singular values.
Moreover, it is noticeable that the alteration in the first singular value is significantly smaller compared to the other two singular values.It is crucial to understand the covariance of singular value changes under normal sensor conditions, as it is a key factor in calculating the quality value, which will be discussed in a later section.In this study, the error under normal star tracker conditions is assumed to be 10 arc-seconds, and the changes and standard deviations of each singular value at that time are summarized in Table 1.

Quality Value Calculation
The remaining question is how to quantify the quality of each observation.In star trackers, we can attribute noise to centroid errors resulting from satellite maneuvers, temperature fluctuations, distortion, and more.Bias can be linked to star deflection and false matches in the identification algorithm.This study suggests utilizing the singular value to quantify both error and bias, leveraging its inherent invariance to attitude changes in an ideal scenario.
As outlined in the background, a minimum of three stars is necessary to calculate the singular value.The subsequent question is to identify which stars, employed in this calculation, might exhibit error and bias.To address this, we specifically choose three stars, the minimum requisite for calculating the singular value, and create a subset as follows: where U denotes the set comprising all star observations: The singular value computed for a specific subset can be represented in the following manner: with v SV i and b SV i signifying the noise and bias associated with the i th singular value calculated from the subset sub j .Calculating the true value of the singular value is achievable using database information for a recognized star.When we rewrite the preceding expression as depicted below, the difference between the observed value and the true singular value becomes expressible solely as a function of noise and bias, irrespective of their specific values.
Therefore, the procedure for evaluating the error level can be characterized as a problem of quantifying the extent of bias within the singular value difference SV di f f sub j of the subgroup sub j , as defined below.
To address this issue, we establish the null hypothesis H 0 and the alternative hypothesis H 1 and introduce statistical hypothesis testing.Hypothesis testing is a significance test that determines whether to accept or reject the null hypothesis.H 0 represents the ideal scenario where the singular value derived from H 0 is devoid of errors and biases.
H 0 = (VB) sub j = 0 (32) where By utilizing the error covariance value in the normal situation of the star tracker, SV di f f sub j is normalized as follows: Since Γ SV di f f sub j represents the dot product of identical vectors with consideration for variance, it follows a chi-square distribution, and R sv is the error covariance of the singular value in the normal situation.
The p-value represents the probability of observing an event of a particular magnitude under the assumption that the null hypothesis is true.It is defined by P sub j as follows: F 0 stands for the cumulative probability function of the null hypothesis.When the Γ SV di f f sub j value is equal to 0, it completely matches the null hypothesis, resulting in a p-value of 1.Therefore, the p-value, denoted as P sub j , indicates the likelihood of how closely the subgroup sub j belongs to H 0 and ranges between 0 and 1.A higher P sub j implies stronger support for H 0 , indicating higher quality for the observations in sub j .Consequently, the individual quality value q s i for the i th star can be computed by aggregating the p-values from all subgroups containing that star.
where n sub j represents the number of all subsets containing the i th star s i and N is the total number of observed stars.The calculated q s i has the same meaning as k i in Equation (18).Therefore, it can be substituted into Equation (22) to calculate the normalized weight of the attitude determination process.This adjustment allows us to assign more weight to observations with higher quality, ultimately leading to improved attitude accuracy.

Quality Value Verification
To verify the performance of the proposed approach, the initial step involves checking that the quality value calculation is carried out correctly as intended.Subsequently, the effectiveness of the calculated quality value when applied to the attitude determination process is assessed.As previously mentioned, it is assumed that the star tracker would exhibit a standard error of 10 arc-seconds (1σ) under normal operating conditions, as discussed earlier.The R sv value used in the quality value calculation is based on this assumption and made use of the values provided in Table 1.
Several scenarios are then considered in which larger than normal errors are intentionally added to specific stars, and Table 2 summarizes them.In all scenarios, six stars are initially recognized by star tracker.Scenario 1 involves the introduction of errors of 25 arc-seconds and 60 arc-seconds to the 2nd and 6th stars, respectively.In Scenario 2, errors of 35 arc-seconds and 70 arc-seconds are added to the 2nd and 6th stars, respectively.Quality values are calculated for 100 s in both Scenarios 1 and 2, confirming that varying degrees of error yield the correct quality values.In Scenario 3, the quality value is computed over a 200-s duration, with the error size being altered midway, and the outcomes are observed.Implemented errors in the scenario is assumed to follow a random Gaussian distribution.To ensure a reliable acquisition of stars and facilitate the analysis of simulation results, it is assumed that the true attitude of the star tracker is maintained in a stationary state.The star tracker's update rate is 1Hz, and the simulation time is in the time-scaled values of this rate.For instance, if 100 s of simulation time elapse, it means the acquisition of 100 image frames.Each image frame consistently captures six stars, and the identities of these stars in the catalog are summarized in Table 3. Figures 3 and 4 illustrate the simulation outcomes for Scenarios 1 and 2, respectively.These results affirm the effective performance of the proposed algorithm in calculating the quality value, which accurately reflects the magnitude of errors present in the observed stars.Notably, the quality values of the 6th star in Scenario 1 and Scenario 2, both of which contain errors exceeding 60 arc-seconds, significantly larger than the normal situation, are observed to be nearly zero.These results suggest that when incorrect star observations arise due to sensor malfunction or high rate slewing motion, the quality value computation can effectively filter out these erroneous observations before entering the attitude determination process.Furthermore, the errors introduced to the 2nd star in Scenario 1 and Scenario 2 exhibit slight variations.Upon comparing these two cases, it becomes evident that the 2nd star in Scenario 2, characterized by a larger error, possesses a relatively smaller quality value.The results of Scenario 3 are presented in Figure 5.It can be observed that the error level in the observation changes after 100 s, and the corresponding shift in quality levels is accurately reflected.From the simulation results, it can be inferred that the proposed approach is capable of effectively accommodating the existing errors in the observations and can serve as a valuable weighting factor.The outcomes across various scenarios indicate that incorrect stars exhibit quality values comparable to those of normal stars.However, this resemblance can be attributed to the assumption that errors assigned to incorrect stars follow a random Gaussian distribution.Within the simulation process, instances should arise where the error is approximately 10 arc-seconds.In such cases, the algorithm naturally calculates a quality value similar to that of stars in standard situations.To effectively illustrate whether the proposed work accurately reflects the magnitude of errors, the results of quality value calculations are verified based on the error size.The outcomes across various scenarios indicate that incorrect stars exhibit quality values comparable to those of normal stars.However, this resemblance can be attributed to the assumption that errors assigned to incorrect stars follow a random Gaussian distribution.Within the simulation process, instances should arise where the error is approximately 10 arc-seconds.In such cases, the algorithm naturally calculates a quality value similar to that of stars in standard situations.To effectively illustrate whether the proposed work accurately reflects the magnitude of errors, the results of quality value calculations are verified based on the error size.Figure 6 presents the trend of quality values in relation to the error size in Scenario 1.
As depicted in Figure 6, the quality value exhibits a clear linear decrease with increasing error.An additional noteworthy point from the results is that stars with the same error size do not consistently share the same quality value.The proposed methodology utilizes changes in singular values associated with errors to assess the magnitude of the error.However, not all stars show an identical magnitude of singular value change for a given error size.While the variation in singular values due to errors generally displays a similar level of fluctuation on average, individual stars exhibit some variability.As observed in the figure, there is a subtle variation in the quality value even for the same error size.
to the assumption that errors assigned to incorrect stars follow a random Gaussian distribution.Within the simulation process, instances should arise where the error is approximately 10 arc-seconds.In such cases, the algorithm naturally calculates a quality value similar to that of stars in standard situations.To effectively illustrate whether the proposed work accurately reflects the magnitude of errors, the results of quality value calculations are verified based on the error size.Figure 6 presents the trend of quality values in relation to the error size in Scenario 1.The current significance of this research lies in its ability to assess and quantify errors that deviate from the normal range.For instance, issues such as the failure of star identification and the deflection of starlight can induce errors in the form of biases.Simulation results in Scenario 1, where the error form is changed from a random Gaussian to a bias, are illustrated in Figure 7.In contrast to prior findings, it is noteworthy that the quality value displays a distinct shifted pattern based on the magnitude of the error.Therefore, from the perspective of fault detection and isolation, the proposed work can significantly contribute to improving attitude determination.As depicted in Figure 6, the quality value exhibits a clear linear decrease with increasing error.An additional noteworthy point from the results is that stars with the same error size do not consistently share the same quality value.The proposed methodology utilizes changes in singular values associated with errors to assess the magnitude of the error.However, not all stars show an identical magnitude of singular value change for a given error size.While the variation in singular values due to errors generally displays a similar level of fluctuation on average, individual stars exhibit some variability.As observed in the figure, there is a subtle variation in the quality value even for the same error size.
The current significance of this research lies in its ability to assess and quantify errors that deviate from the normal range.For instance, issues such as the failure of star identification and the deflection of starlight can induce errors in the form of biases.Simulation results in Scenario 1, where the error form is changed from a random Gaussian to a bias, are illustrated in Figure 7.In contrast to prior findings, it is noteworthy that the quality value displays a distinct shifted pattern based on the magnitude of the error.Therefore, from the perspective of fault detection and isolation, the proposed work can significantly contribute to improving attitude determination.

Effects Depending on the Number of Abnormal Stars
The proposed method has proven effective in accurately calculating the quality value through a comparison of singular value variations among recognized stars.However, the key factor closely tied to the performance of the proposed method is the ratio between the number of stars within the normal error range and the number of abnormal stars.This ratio may vary from frame to frame.To validate this, additional simulations are conducted, introducing an increased number of stars with significant errors.
In Scenarios 4 and 5, similar to earlier scenarios, it is assumed that a total of six stars

Effects Depending on the Number of Abnormal Stars
The proposed method has proven effective in accurately calculating the quality value through a comparison of singular value variations among recognized stars.However, the key factor closely tied to the performance of the proposed method is the ratio between the number of stars within the normal error range and the number of abnormal stars.This ratio may vary from frame to frame.To validate this, additional simulations are conducted, introducing an increased number of stars with significant errors.
In Scenarios 4 and 5, similar to earlier scenarios, it is assumed that a total of six stars are recognized.However, in each scenario, three and four stars are assumed to have abnormal errors, respectively.Table 4 summarizes the errors associated with abnormal stars, and, for ease of comparison with Figure 7, these errors are presented in terms of bias types.In Scenario 4, with the increase in the number of abnormal stars to three, it indicates a situation where only half of the stars in the image frame exhibit normal errors.The results of Scenario 4 are represented in Figure 8. Comparing it to Figure 7, where there are two abnormal stars, it becomes evident that the distinction in quality values is less pronounced.This underscores that the algorithm performs more effectively when the proportion of normal stars within the image frame is higher.
Sensors 2024, 24, x FOR PEER REVIEW 13 of 18 abnormal errors, respectively.Table 4 summarizes the errors associated with abnormal stars, and, for ease of comparison with Figure 7, these errors are presented in terms of bias types.In Scenario 4, with the increase in the number of abnormal stars to three, it indicates a situation where only half of the stars in the image frame exhibit normal errors.The results of Scenario 4 are represented in Figure 8. Comparing it to Figure 7, where there are two abnormal stars, it becomes evident that the distinction in quality values is less pronounced.This underscores that the algorithm performs more effectively when the proportion of normal stars within the image frame is higher.Additionally, the outcomes of Scenario 5, where the number of abnormal stars is four, are depicted in Figure 9.In this case, the image frame becomes more dominated by the presence of abnormal stars than normal stars.The 4th and 6th stars, characterized by larger errors, exhibit notably lower quality values than the others.However, Additionally, the outcomes of Scenario 5, where the number of abnormal stars is four, are depicted in Figure 9.In this case, the image frame becomes more dominated by the presence of abnormal stars than normal stars.The 4th and 6th stars, characterized by larger errors, exhibit notably lower quality values than the others.However, distinguishing a meaningful difference in quality values for the 1st and 3rd stars, which have normal errors when compared to stars #2 and #5, , proves challenging.
Sensors 2024, 24, x FOR PEER REVIEW 14 of 18 distinguishing a meaningful difference in quality values for the 1st and 3rd stars, which have normal errors when compared to stars #2 and #5, , proves challenging.These outcomes arise from the methodology of computing singular values, which involves grouping three stars rather than assessing individual stars.The proposed approach determines the quality value for a specific star by averaging the assessed quality values of all subgroups containing that star of interest.Consequently, in scenarios where the number of abnormal stars dominates within an image frame, identifying the specific star contributing to the singular value error within that group becomes difficult.This underscores a limitation inherent in the approach utilizing singular values as a reference.

Attitude Determination
In this section, the difference in attitude determination accuracy will be compared when the calculated quality value is applied as a weighting factor for each observation versus when it is not applied.For comparison, Scenario 3 previously used to examine the quality value is employed.Additionally, in contrast to the scenarios, an ideal case is introduced where all stars maintain a normal condition, making comparison more straightforward.Attitude determination utilizes the QUEST algorithm as described in the preceding section, with the outcomes presented in Figure 10.
The results indicate that reflecting the error level is effective in enhancing positional accuracy.A crucial aspect of the proposed work is that even in scenarios with incorrect stars, it helps approach the positional accuracy of the ideal case.The process of attitude determination unfolds in three key stages: Sensor Measurement (Vector Observation), Observation Error Level Assessment, Attitude Determination.The primary goal of the presented work is to enhance attitude determination accuracy.Specifically, the focus is on evaluating the quality value of vector observations used in attitude determination.This These outcomes arise from the methodology of computing singular values, which involves grouping three stars rather than assessing individual stars.The proposed approach determines the quality value for a specific star by averaging the assessed quality values of all subgroups containing that star of interest.Consequently, in scenarios where the number of abnormal stars dominates within an image frame, identifying the specific star contributing to the singular value error within that group becomes difficult.This underscores a limitation inherent in the approach utilizing singular values as a reference.

Attitude Determination
In this section, the difference in attitude determination accuracy will be compared when the calculated quality value is applied as a weighting factor for each observation versus when it is not applied.For comparison, Scenario 3 previously used to examine the quality value is employed.Additionally, in contrast to the scenarios, an ideal case is introduced where all stars maintain a normal condition, making comparison more straightforward.Attitude determination utilizes the QUEST algorithm as described in the preceding section, with the outcomes presented in Figure 10.
proaches like the least square method and QUEST, as well as dynamic methods like the the Extended Kalman Filter (EKF) that integrates gyro data.The pre-evaluated error level can find application in the cost function of the QUEST algorithm or serve as a weighting factor in the EKF update process.Consequently, the presented work holds a significant advantage, being universally applicable and capable of enhancing performance, irrespective of the attitude determination method employed.Additional simulation is carried out to observe the trends in algorithm behavior concerning error magnitude.Assuming errors exist in two of the six stars, mirroring Scenario 1, a total of 100 Monte Carlo simulations are executed for each error size.The error magnitude was progressively increased, and the average results were analyzed.Figure 11 exclusively illustrates results up to 50 arc-seconds, offering a detailed analysis of the trend in the algorithm's accuracy improvement with respect to error size.The noticeable divergence in attitude determination accuracy becomes evident when the error size exceeds 10 arc-seconds.
The observed trend can be directly attributed to the value of  used in the quality value calculation.This particular parameter carries significant weight within the algorithm, representing the error of the star sensor under normal conditions.In this study, it is assumed to be 10 arc-seconds.This assumption guides the proposed algorithm to assign a lower weighting factor to star observations with errors larger than those typically encountered.
Ensuring the accuracy of the algorithm requires prior knowledge of errors in the sensor's typical operating conditions.However, changes in sensor performance can manifest during mission operations.Short-term fluctuations may arise from temperature variations in the space environment, while long-term degradation can occur as hardware approaches the end of its life.Therefore, instead of keeping the value applied to  constant, adapting it based on the situation is deemed to significantly enhance algorithm performance.The results indicate that reflecting the error level is effective in enhancing positional accuracy.A crucial aspect of the proposed work is that even in scenarios with incorrect stars, it helps approach the positional accuracy of the ideal case.The process of attitude determination unfolds in three key stages: Sensor Measurement (Vector Observation), Observation Error Level Assessment, Attitude Determination.The primary goal of the presented work is to enhance attitude determination accuracy.Specifically, the focus is on evaluating the quality value of vector observations used in attitude determination.This aligns more closely with the preliminary task of assessing observation error levels rather than committing to a specific attitude determination method.
Various methods are employed for attitude determination, encompassing static approaches like the least square method and QUEST, as well as dynamic methods like the the Extended Kalman Filter (EKF) that integrates gyro data.The pre-evaluated error level can find application in the cost function of the QUEST algorithm or serve as a weighting factor in the EKF update process.Consequently, the presented work holds a significant advantage, being universally applicable and capable of enhancing performance, irrespective of the attitude determination method employed.
Additional simulation is carried out to observe the trends in algorithm behavior concerning error magnitude.Assuming errors exist in two of the six stars, mirroring Scenario 1, a total of 100 Monte Carlo simulations are executed for each error size.The error magnitude was progressively increased, and the average results were analyzed.Figure 11 exclusively illustrates results up to 50 arc-seconds, offering a detailed analysis of the trend in the algorithm's accuracy improvement with respect to error size.The noticeable divergence in attitude determination accuracy becomes evident when the error size exceeds 10 arc-seconds.
The observed trend can be directly attributed to the value of R sv used in the quality value calculation.This particular parameter carries significant weight within the algorithm, representing the error of the star sensor under normal conditions.In this study, it is assumed to be 10 arc-seconds.This assumption guides the proposed algorithm to assign a lower weighting factor to star observations with errors larger than those typically encountered.
Ensuring the accuracy of the algorithm requires prior knowledge of errors in the sensor's typical operating conditions.However, changes in sensor performance can manifest during mission operations.Short-term fluctuations may arise from temperature variations in the space environment, while long-term degradation can occur as hardware approaches the end of its life.Therefore, instead of keeping the value applied to R sv constant, adapting it based on the situation is deemed to significantly enhance algorithm performance.In the final analysis, errors were increased up to 150 arc-seconds, and the results are presented in Figure 12.The results demonstrate that when the proposed algorithm is applied as the weighting factor for observations, more precise attitude determination can be achieved.A noteworthy point is that as the error present in an observation increases, the confidence level of that observation decreases, making it a negligible factor in the attitude determination.As a result, utilizing only observations that have higher quality leads to significantly improved attitude accuracy.While the performance improvement is not substantial under normal circumstances, the results reveal significant enhancements as the error value increases.However, the normal sensor behavior cannot always be guaranteed during operations; continuous errors or sudden anomalies like image spikes may occur.In such situations, it is believed that the proposed algorithm can play a pivotal role in enhancing the stability required for the spacecraft's mission.In the final analysis, errors were increased up to 150 arc-seconds, and the results are presented in Figure 12.The results demonstrate that when the proposed algorithm is applied as the weighting factor for observations, more precise attitude determination can be achieved.A noteworthy point is that as the error present in an observation increases, the confidence level of that observation decreases, making it a negligible factor in the attitude determination.As a result, utilizing only observations that have higher quality leads to significantly improved attitude accuracy.While the performance improvement is not substantial under normal circumstances, the results reveal significant enhancements as the error value increases.However, the normal sensor behavior cannot always be guaranteed during operations; continuous errors or sudden anomalies like image spikes may occur.In such situations, it is believed that the proposed algorithm can play a pivotal role in enhancing the stability required for the spacecraft's mission.In the final analysis, errors were increased up to 150 arc-seconds, and the results are presented in Figure 12.The results demonstrate that when the proposed algorithm is applied as the weighting factor for observations, more precise attitude determination can be achieved.A noteworthy point is that as the error present in an observation increases, the confidence level of that observation decreases, making it a negligible factor in the attitude determination.As a result, utilizing only observations that have higher quality leads to significantly improved attitude accuracy.While the performance improvement is not substantial under normal circumstances, the results reveal significant enhancements as the error value increases.However, the normal sensor behavior cannot always be guaranteed during operations; continuous errors or sudden anomalies like image spikes may occur.In such situations, it is believed that the proposed algorithm can play a pivotal role in enhancing the stability required for the spacecraft's mission.

Conclusions
In this study, a novel method is proposed to enhance attitude determination accuracy using star trackers.The key innovation involves the utilization of singular values, which remain invariant despite attitude changes, as reference values for assessing the magnitude of errors introduced in observations.The calculated singular value difference is quantified through hypothesis testing and is subsequently incorporated as an observational weight in the attitude determination process to improve attitude accuracy.Simulations confirm that the algorithm correctly evaluates the quality level.A noteworthy result is that errors attributable to sensor failures are rarely used as observations, which can enhance the reliability of attitude determination.Therefore, it is anticipated that the proposed algorithm can be applied to missions requiring high-precision attitude, as well as contribute to the stability of spacecraft systems.

Figure 1 .
Figure 1.A set of two stars seen from two different reference frames: (a) an inertial reference frame; (b) the body reference frame.

Figure 1 .
Figure 1.A set of two stars seen from two different reference frames: (a) an inertial reference frame; (b) the body reference frame.

∼
S r and ∼ S b represent 3 × 3 diagonal matrices with the squares of the three eigenvalues, σ 2 r i and σ 2 b i

Figure 2 .
Figure 2. Variation of singular values at each noise level.

Figure 2 .
Figure 2. Variation of singular values at each noise level.

Figure 6
presents the trend of quality values in relation to the error size in Scenario 1.

Figure 6 .
Figure 6.Quality value trend on error size.

Figure 6 .
Figure 6.Quality value trend on error size.

Figure 7 .
Figure 7. Quality value calculation results with bias type noise.

Figure 7 .
Figure 7. Quality value calculation results with bias type noise.

Figure 10 .
Figure 10.Comparing results of the attitude determination accuracy under Scenario 3.

Figure 10 .
Figure 10.Comparing results of the attitude determination accuracy under Scenario 3.

Figure 12 .Figure 11 .
Figure 12.Comparing results of the attitude determination accuracy up to 150 arc-seconds.

Figure 12 .Figure 12 .
Figure 12.Comparing results of the attitude determination accuracy up to 150 arc-seconds.

Table 1 .
Variation and standard deviation of singular values.

Table 1 .
Variation and standard deviation of singular values.

Table 2 .
Scenarios for quality value verification.

Table 3 .
Identity of image stars utilized in simulation.

Table 4 .
Scenarios for verifying the effects based on the number of abnormal stars.

Table 4 .
Scenarios for verifying the effects based on the number of abnormal stars.