[发明专利]一种自适应光学成像中的点扩散函数估计方法

权利要求书

1.一种自适应光学成像中的点扩散函数估计方法，其特征在于，包括以下步骤：

1）指定样条函数的次数d及刻画光滑程度的连续性阶数r，确定波前区域内三角剖分的形式；

2）建立斜率测量模型的偏微分方程并确定方程中各参数矩阵的值；

3）根据具体的参数值确定变换矩阵D；

4）结合相邻样条函数间的光滑拼接条件确定重构矩阵R；

5）利用斜率测量数据计算各二元样条函数的系数向量；

6）根据系数重构出畸变波前φ(x,y)；

7）由重构的降质波前去估计大气瞬时点扩散函数h(x,y)。

2.根据权利要求1所述的自适应光学成像中的点扩散函数估计方法，其特征在于：所述步骤1）中，根据微透镜阵列排布方式确定三角剖分的形式，确定三角剖分的形式时采用基于几何布局的网格单元直接剖分方法。

3.根据权利要求1所述的自适应光学成像中的点扩散函数估计方法，其特征在于：所述步骤2）中，斜率测量模型的偏微分方程的建立方法具体为：

首先将波前表示为：

φ(x)≈sdr(b(x))=Σ|λ|=dcλBλd(b(x))=Bd(b(x))·c]]>

其m阶方向导数可表示为Dumsdr(b(x))=d!(d-m)!Bd-mPud,d-m(a)·c,]]>其中，为二元单纯形样条，d为样条函数次数，r为连续性阶数，b(x)为一个三角剖分上笛卡尔坐标到面积坐标的变换，为d次Bernstein多项式基函数，λ=(λ₀,λ₁,λ₂)为多重符号且其长度为|λ|=λ₀+λ₁+λ₂=d，c_λ为多项式系数，B^d(b(x))和c分别为和c_λ的矩阵形式，u为笛卡尔坐标下方向导数的方向向量，a为与u相对应的面积坐标下的方向坐标，为分块对角矩阵，每对角块为对应于a和u的de Casteljau矩阵；其次根据相邻三角剖分间的r阶光滑拼接条件，每条边都有个约束条件，从而构成约束方程Hc=0，H为由E条边构成的E*Z行约束矩阵，且H的每行都包含一个单独的连续性约束，这样，波前传感器的斜率测量偏微分方程即可表示为：

σu=dBd-1Pud,d-1c+n]]>

Hc=0

其中，σ_u表示波前传感器斜率测量矩阵，n表示测量噪声及模型误差。