[发明专利]基于点扩散函数波前修正的星点像质心计算方法

权利要求书

1.一种基于点扩散函数波前修正的星点像质心计算方法，包括：

步骤101、采集多帧星图；

步骤102、在所述星图中，找到每颗星点像峰值所在像素，并以星点像峰值所在像素为中心选择一个矩形区域，这个矩形区域包括星点像绝大部分能量，这个矩形区域由测量精度需求决定，矩形区域行列像素数为奇数；

步骤103、在步骤102所得到的每颗星点像所在的矩形区域内，用光子权重平均的方法计算每颗星点像的k阶矩；

步骤104、根据步骤103所得到的k阶矩以及一个与所述k阶矩和峰值像素坐标的差成正比的偏差，得到星点像的真实质心坐标。

2.根据权利要求1所述的基于点扩散函数波前修正的星点像质心计算方法，其特征在于，在步骤101中，还包括对所采集的多帧星图进行预处理，所述预处理包括去除暗噪声、本底噪声，并进行平场修正。

3.根据权利要求1或2所述的基于点扩散函数波前修正的星点像质心计算方法，其特征在于，在步骤101中，所采集的星图为真实星图，或在实验室利用光学成像模拟的星图。

4.根据权利要求1或2所述的基于点扩散函数波前修正的星点像质心计算方法，其特征在于，在步骤103中，用光子权重平均的方法计算每颗星点像的k阶矩的计算公式如下：

rg=Σi=1mΣj=1nIijki/Σi=1mΣj=1nIijk]]>

cg=Σi=1mΣj=1nIijkj/Σi=1mΣj=1nIijk]]>

其中，k为自然数，(r_g，c_g)为光子权重平均计算出的k阶矩，r表示行，c表示列；探测器阵列为m×n，(i，j)为像素阵列的行列指标，I_ij为第i行、第j列的灰度值。