[发明专利]一种图像处理中构建光滑点扩散函数的方法

权利要求书

1.一种图像处理中构建光滑点扩散函数的方法，其特征在于，包括：

步骤1、利用期望最大化主成分分析法来拟合PSF，并获取第j阶主成分P_j，k以及主成分系数C_i，j；

i是正整数，表示点源序号；

k是正整数，表示第i点源的第k像素，即像素序号；

j是非负整数，表示主成分的阶；

P_j，k是第j阶主成分里的第k像素上对应的值；

步骤2、利用Moffatlets基函数F_l，k拟合各阶主成分P_j，k并确定基函数系数D_j，l；

步骤3、将第j阶更新为第j+1阶，重复步骤1、2，直到残差符合设定阈值。

2.根据权利要求1所述的一种图像处理中构建光滑点扩散函数的方法，其特征在于，所述步骤1中，假设一系列点光源的亮度为O_i，k，即观测点光源的PSF，是焦平面上的一个二维密度轮廓。

3.根据权利要求2所述的一种图像处理中构建光滑点扩散函数的方法，其特征在于，所述步骤1中，EMPCA中任何一个光源的二维空间分布用各阶主成分的和式表示，对于第i个点源的亮度I_i，k，表示为：

N是用来构建PSF的主成分的最大阶数；

C_i，j是主成分系数矩阵，表示第i点源的第j阶主成分上的系数。

4.根据权利要求3所述的一种图像处理中构建光滑点扩散函数的方法，其特征在于，步骤2中，用Moffatlets基函数拟合各阶主成分P_j，k，得到光滑的PSF；

对于一个主成分P_j，k，用Moffatlets基函数来表示：

M是用来拟合第j阶主成分的基函数最大阶数；

D_j，l是基函数系数，表示为拟合第j阶主成分，所使用的第l阶基函数对应的系数；

F_l，k是第l阶基函数在第k像素上对应的值；

l是非负整数，表示基函数的阶。

5.根据权利要求4所述的一种图像处理中构建光滑点扩散函数的方法，其特征在于，步骤2中，通过最小二乘法构建关于D_j，l的线性方程组χ²函数，求解可得D_j，l；

所述χ²函数为：

W_i，k是权重函数。

6.根据权利要求5所述的一种图像处理中构建光滑点扩散函数的方法，其特征在于，对于W_i，k，选取本地噪音。