CN103714516B - 一种自适应光学成像中的点扩散函数估计方法 - Google Patents
一种自适应光学成像中的点扩散函数估计方法 Download PDFInfo
- Publication number
- CN103714516B CN103714516B CN201310651641.9A CN201310651641A CN103714516B CN 103714516 B CN103714516 B CN 103714516B CN 201310651641 A CN201310651641 A CN 201310651641A CN 103714516 B CN103714516 B CN 103714516B
- Authority
- CN
- China
- Prior art keywords
- matrix
- wavefront
- function
- spline
- point spread
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 54
- 230000003044 adaptive effect Effects 0.000 title claims abstract description 11
- 238000012634 optical imaging Methods 0.000 title claims abstract description 10
- 239000011159 matrix material Substances 0.000 claims abstract description 67
- 238000005259 measurement Methods 0.000 claims abstract description 44
- 230000009466 transformation Effects 0.000 claims abstract description 17
- 210000001747 pupil Anatomy 0.000 claims abstract description 9
- 239000013598 vector Substances 0.000 claims description 31
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000009792 diffusion process Methods 0.000 claims description 3
- 230000017105 transposition Effects 0.000 claims description 3
- 238000003384 imaging method Methods 0.000 abstract description 4
- 230000007547 defect Effects 0.000 abstract description 2
- 230000000593 degrading effect Effects 0.000 abstract 1
- 238000005192 partition Methods 0.000 abstract 1
- 230000003287 optical effect Effects 0.000 description 6
- 238000004422 calculation algorithm Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000012937 correction Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
Landscapes
- Length Measuring Devices By Optical Means (AREA)
- Image Processing (AREA)
Abstract
本发明公开了一种自适应光学成像中的点扩散函数估计方法。该发明基于研究多元样条的B网方法,首先将经过大气降质的畸变波前近似表示为单纯形剖分上的二元样条,进而得到波前传感器斜率测量数据与样条函数系数之间的数学关系;其次,结合相邻三角剖分上样条函数的光滑拼接条件,利用斜率测量数据通过约束最小二乘估计样条函数系数进而重构畸变波前;最后将得到的畸变波前结合透镜光瞳函数,即可估计出经大气湍流成像的瞬时点扩散函数。本发明与传统的Zernike模式法估计点扩散函数相比克服了变换矩阵秩不完备的缺陷,提高了估计精度;并且适用于中心被遮挡的波前传感器及微透镜阵列的多种几何排布的情况。
Description
技术领域
本发明属于空间目标探测识别技术领域,涉及到空间目标自适应光学图像复原理论与方法,具体涉及到一种自适应光学成像中的点扩散函数估计方法。
背景技术
根据空间目标光学成像探测的特性,自适应光学图像复原方法主要可分为两类:盲解卷积和波前解卷积。对于盲解卷积,一个相对准确的点扩散函数估计值作为迭代算法的初始估计将会加快收敛、减少计算量;对于波前解卷积,点扩散函数估计值与真实值之间的偏差更是会直接影响复原的效果。因此快速准确的点扩散函数估计方法对于自适应光学图像的复原来讲意义显著。
点扩散函数与波前存在着理论上的对应关系,在自适应光学***中用波前传感器探测波前信息。目前使用H-S型波前传感器的自适应光学***居多,所获得数据是波前相位的斜率采样值。通过波前重构可以在有噪声及干扰的情况下,利用有限离散观测信息恢复出完整的波前,进而重构短曝光图像的瞬时点扩散函数。波前重构的基本方法大致上可以分为两大类:区域法和模式法。由于Zernike多项式各模式的像差意义明显,故Zernike模式法在面向图像复原的波前重构方法中一直得到广泛的应用。Zernike模式法最终归结为求解线性方程组,但由于Zernike函数偏导数的不完全正交性以及在有限采样点上函数的非正交性,在求最小二乘解时会出现变换矩阵D的秩不完备、DTD奇异或接近奇异的情况,这将使梯度矢量的测量误差被放大,导致波前重构结果精度降低,进而影响点扩散函数的估计精度;其次,许多大型望远镜的中心常常被遮挡,若仍用Zernike多项式展开波前,因为展开域是环形的,所以Zernike展开不再具有正交性,仍应用此法则需要十分繁琐的数学推导与计算,这些因素都限制了Zernike模式法的应用效力。另外,虽然近年来在大型自适应光学设备快速发展的推动下不断涌现出新的面向校正控制的波前重构算法,但把这些方法进一步的应用到点扩散函数估计直至图像复原中的效力还并未得到有效的研究和验证。
发明内容
本发明的目的在于克服现有成像点扩散函数估计方法的不足,提出一种自适应光学成像中的点扩散函数估计方法,该发明基于研究多元样条的B网方法,将经过大气降质的畸变波前近似表示为单纯形剖分上的二元样条,避免了求最小二乘解时不完备或奇异矩阵的出现,提高了估计精度,并可以应对多种波前传感器微透镜阵列的几何排布方式。
为了达到上述目的,本发明所采用的技术方案包括以下步骤:
1)指定样条函数的次数d及刻画光滑程度的连续性阶数r,确定波前区域内三角剖分的形式;
2)建立斜率测量模型的偏微分方程并确定方程中各参数矩阵的值;
3)根据具体的参数值确定变换矩阵D;
4)结合相邻样条函数间的光滑拼接条件确定重构矩阵R;
5)利用斜率测量数据计算各二元样条函数的系数向量
6)根据系数重构出畸变波前φ(x,y);
7)由重构的降质波前去估计大气瞬时点扩散函数h(x,y)。
所述步骤1)中,根据微透镜阵列排布方式确定三角剖分的形式,确定三角剖分的形式时采用基于几何布局的网格单元直接剖分方法。
所述步骤2)中,斜率测量模型的偏微分方程的建立方法具体为:
首先将波前表示为:
其m阶方向导数可表示为其中,为二元单纯形样条,d为样条函数次数,r为连续性阶数,b(x)为一个三角剖分上笛卡尔坐标到面积坐标的变换,为d次Bernstein多项式基函数,λ=(λ0,λ1,λ2)为多重符号且其长度为|λ|=λ0+λ1+λ2=d,cλ为多项式系数,Bd(b(x))和c分别为和cλ的矩阵形式,u为笛卡尔坐标下方向导数的方向向量,a为与u相对应的面积坐标下的方向坐标,为分块对角矩阵,每对角块为对应于a和u的deCasteljau矩阵;其次根据相邻三角剖分间的r阶光滑拼接条件,每条边都有个约束条件,从而构成约束方程Hc=0,H为由E条边构成的E*Z行约束矩阵,且H的每行都包含一个单独的连续性约束,这样,波前传感器的斜率测量偏微分方程即可表示为:
Hc=0
其中,σu表示波前传感器斜率测量矩阵,n表示测量噪声及模型误差。
相邻三角剖分间的r阶光滑拼接条件描述如下:
设S(b(x))和分别是定义在相邻三角形t=<v1,v2,v3>和上的d次多项式,vi(i=1,2,3)为三角形的顶点,cλ和分别是S(b(x))和关于t和的贝塞尔坐标,则S(b(x))和之间的r阶光滑拼接条件为:
其中,是关于t的面积坐标,λs=(s,λ1,λ2),λ0=(0,λ1,λ2),λ1+λ2=d-s。
所述步骤3)中,确定变换矩阵D的方法具体为:
步骤2)建立了一个约束最小二乘问题σu=Dc+n,Hc=0;与Zernike模式法相对应,通过对比此约束最小二乘问题与斜率测量偏微分方程,即可将变换矩阵D表示为:
所述步骤4)中,确定重构矩阵R的方法具体为:
将由斜率测量数据去重构样条函数系数向量的矩阵定义为重构矩阵,结合步骤1)和步骤2),通过求解约束最小二乘问题,多项式系数向量c的显式解表达式为:
重构矩阵R,则表示为:
R=(DTD)-1DT-(DTD)-1HT[H(DTD)-1HT]-1(DTD)-1DT
其中,T表示矩阵转置,σ为波前斜率测量向量,此时由重构矩阵所确定的样条函数系数是在约束Hc=0下的最优线性无偏估计。
所述步骤5)中,根据步骤4),重构矩阵是由斜率测量数据去重构样条函数系数向量的矩阵,二元样条函数的系数向量直接由重构矩阵与波前斜率测量向量相乘得到,则其计算公式为:
所述步骤6)中,畸变波前φ(x,y)的重构方法为:
根据步骤2),已将畸变波前表示为单纯形剖分上的二元样条,步骤5)又得到了样条函数的系数向量而d次Bernstein多项式基函数Bd(b(x))是确定的,将二者相乘即可得到畸变波前的估计量即:
所述步骤7)中,估计大气瞬时点扩散函数h(x,y)的方法具体为:
已知透镜光瞳函数为A,根据步骤6),由傅里叶光学理论,结合透镜光瞳函数A,得到其波前对应的点扩散函数为:
其中,FT-1表示傅里叶逆变换,
与现有技术相比,本发明具有以下有益效果:
本发明在重构波前时满足变换矩阵D满秩、DTD非奇异的条件,不会放大斜率测量误差,并且基于二元样条的方法可以选择使用非线性基函数,这都促使波前重构精度的提高,进而可以提高点扩散函数的估计精度;另外,本发明不受波前传感器微透镜阵列几何排布方式的限制,对圆形或方形布局都有效,并且还适用于中心被遮挡的情况,方法的通用性较强。
附图说明
图1为本发明的原理示意图;
图2为本发明的具体流程图;
图3为本发明一种典型54子孔径波前传感器微透镜阵列排布图;
图4为本发明针对图3波前传感器做的一种三角剖分示例图。
具体实施方式
下面结合附图和实施例对本发明作进一步详细的描述:
参见图1和图2,本发明包括以下步骤:
1)指定样条函数的次数d及刻画光滑程度的连续性阶数r,确定波前区域内三角剖分的形式;
根据微透镜阵列排布方式确定三角剖分的形式,确定三角剖分的形式时采用基于几何布局的网格单元直接剖分方法。
2)建立斜率测量模型的偏微分方程并确定方程中各参数矩阵的值;
斜率测量模型的偏微分方程的建立方法具体为:
首先将波前表示为:
其m阶方向导数可表示为其中,为二元单纯形样条,d为样条函数次数,r为连续性阶数,b(x)为一个三角剖分上笛卡尔坐标到面积坐标的变换,为d次Bernstein多项式基函数,λ=(λ0,λ1,λ2)为多重符号且其长度为|λ|=λ0+λ1+λ2=d,cλ为多项式系数,Bd(b(x))和c分别为和cλ的矩阵形式,u为笛卡尔坐标下方向导数的方向向量,a为与u相对应的面积坐标下的方向坐标,为分块对角矩阵,每对角块为对应于a和u的deCasteljau矩阵;其次根据相邻三角剖分间的r阶光滑拼接条件,每条边都有个约束条件,从而构成约束方程Hc=0,H为由E条边构成的E*Z行约束矩阵,且H的每行都包含一个单独的连续性约束,这样,波前传感器的斜率测量偏微分方程即可表示为:
Hc=0
其中,σu表示波前传感器斜率测量矩阵,n表示测量噪声及模型误差。
相邻三角剖分间的r阶光滑拼接条件描述如下:
设S(b(x))和分别是定义在相邻三角形t=<v1,v2,v3>和上的d次多项式,vi(i=1,2,3)为三角形的顶点,cλ和分别是S(b(x))和关于t和的贝塞尔坐标,则S(b(x))和之间的r阶光滑拼接条件为:
其中,是关于t的面积坐标,λs=(s,λ1,λ2),λ0=(0,λ1,λ2),λ1+λ2=d-s。
3)根据具体的参数值确定变换矩阵D;
确定变换矩阵D的方法具体为:
步骤2)建立了一个约束最小二乘问题σu=Dc+n,Hc=0;与Zernike模式法相对应,通过对比此约束最小二乘问题与斜率测量偏微分方程,即可将变换矩阵D表示为:
4)结合相邻样条函数间的光滑拼接条件确定重构矩阵R;
确定重构矩阵R的方法具体为:
将由斜率测量数据去重构样条函数系数向量的矩阵定义为重构矩阵,结合步骤1)和步骤2),通过求解约束最小二乘问题,多项式系数向量c的显式解表达式为:
这样重构矩阵R即可表示为:
R=(DTD)-1DT-(DTD)-1HT[H(DTD)-1HT]-1(DTD)-1DT
其中,T表示矩阵转置,σ为波前斜率测量向量,此时由重构矩阵所确定的样条函数系数是在约束Hc=0下的最优线性无偏估计。
5)利用斜率测量数据计算各二元样条函数的系数向量
根据步骤4),重构矩阵是由斜率测量数据去重构样条函数系数向量的矩阵,二元样条函数的系数向量可直接由重构矩阵与波前斜率测量向量相乘得到,则其计算公式为:
6)根据系数重构出畸变波前φ(x,y);
畸变波前φ(x,y)的重构方法为:
根据步骤2),已将畸变波前表示为单纯形剖分上的二元样条,步骤5)又得到了样条函数的系数向量而d次Bernstein多项式基函数Bd(b(x))是确定的,这样将二者相乘即可得到畸变波前的估计量即:
7)由重构的降质波前去估计大气瞬时点扩散函数h(x,y)。
估计大气瞬时点扩散函数h(x,y)的方法具体为:
已知透镜光瞳函数为A,根据步骤6),由傅里叶光学理论,结合透镜光瞳函数A,则其波前对应的点扩散函数为:
其中,FT-1表示傅里叶逆变换,
本发明的原理在于:
本发明基于研究多元样条的B网方法,首先将经过大气降质的畸变波前近似表示为单纯形剖分上的二元样条,进而得到波前传感器斜率测量数据与样条函数系数之间的数学关系;其次,结合相邻三角剖分上样条函数的光滑拼接条件,利用斜率测量数据由约束最小二乘估计样条函数系数进而重构畸变波前;最后,畸变波前结合透镜光瞳函数,即可估计出大气湍流的成像点扩散函数。
实施例:
为使本发明的目的、技术方案、优点等更加清楚明白,以下结合具体实例,并参照附图,对本发明做进一步详细说明。本发明是一种基于多元样条的自适应光学成像点扩散函数估计方法,图2为本发明估计点扩散函数的具体流程,以下分别做详细介绍:
1)初始化阶段。依次为,首先依据波前传感器微透镜阵列的布局方式确定波前单纯形区域上三角剖分的形式,这里以54子孔径波前传感器为例,如图3所示,此正六边形微透镜阵列的布局方式本质上是三角的,直观上即可以依据不同的样条函数次数确定出几种三角剖分,不需要复杂的剖分算法,但对于此种波前传感器来讲,传统的基于有限差分的波前重构方法将失效,并且由于其中心被遮挡,Zernike模式法更是失去了正交性,使用时需要十分繁琐的数学修正,此时本发明的优点就显现出来;
其次确定二元样条函数的次数d及相邻三角剖分间光滑拼接的连续性阶数r,如图4所示,针对图3的布局,54个小圆圈表示微透镜阵列的位置,粗实线区域被剖分为3个三角形,整个环形波前区域共分为6个粗实线区域18个三角形,确定样条函数的次数d=2,连续性阶数r=1;
再次,依据以上的描述建立波前斜率测量偏微分方程并确定方程参数矩阵,波前近似表示为:
其m阶方向导数为其中为d次Bernstein多项式基函数,cλ为多项式系数,相应的粗体为其矩阵形式,u为笛卡尔坐标下方向导数的方向向量,a为与u相对应的面积坐标下的方向坐标,为分块对角矩阵,每对角块为对应于a和u的deCasteljau矩阵;其次根据相邻三角剖分间的r阶光滑拼接条件每条边都有个约束条件,从而构成约束方程Hc=0,H为由E条边构成的E*Z行的约束矩阵。这样,波前传感器的斜率测量偏微分方程即可表示为:
其中σu表示波前传感器斜率测量值,n表示测量噪声及模型误差。模型建立好其相应的Bd-1、与H阵也已随即确定;
2)确定变换矩阵D。与Zernike模式法类似,在此也将样条函数系数与斜率测量数据之间的变换关系用变换矩阵D来描述。约束最小二乘问题的模型为σu=Dc+n;Hc=0,结合步骤1)则有:
与Zernike模式法不同的是这里的DTD是非奇异的,这是由于的构成方式以及加入了剖分之间光滑拼接条件约束所直接决定的,从而不会放大斜率数据的测量噪声。
3)计算重构矩阵R。将由斜率测量数据去重构样条函数系数向量的矩阵定义为重构矩阵。结合步骤1)和步骤2),由约束最小二乘问题的显式解表达式:
R=(DTD)-1DT-(DTD)-1HT[H(DTD)-1HT]-1(DTD)-1DT(4)
此时由重构矩阵所确定的样条函数系数是在约束Hc=0下的最优线性无偏估计。
4)计算样条函数系数向量由步骤3)可得:
5)重构畸变波前φ(x,y)。由步骤1):
6)估计点扩散函数h(x,y)。在步骤5)的基础上,由傅里叶光学理论,结合透镜光瞳函数A,其波前对应的点扩散函数为:
其中FT-1表示傅里叶逆变换。至此,则完成了成像点扩散函数的估计。
以上实施例仅用以说明本发明而并非限制本发明所描述的技术方案;因此尽管本说明书参照上述的实施例对本发明已进行了详细的说明,但是本领域的普通技术人员应当理解,仍然可以对本发明进行修改或等同替换;而一切不脱离发明的精神和范围的技术方案及其改进,均应涵盖在本发明的权利要求范围之内。
Claims (3)
1.一种自适应光学成像中的点扩散函数估计方法,其特征在于,包括以下步骤:
1)指定样条函数的次数d及刻画光滑程度的连续性阶数r,确定波前区域内三角剖分的形式;
2)建立斜率测量模型的偏微分方程并确定方程中各参数矩阵的值;
首先将波前表示为:
其m阶方向导数可表示为 其中,为二元单纯形样条,d为样条函数次数,r为连续性阶数,b(x)为一个三角剖分上笛卡尔坐标到面积坐标的变换,为d次Bernstein多项式基函数,λ=(λ0,λ1,λ2)为多重符号且其长度为|λ|=λ0+λ1+λ2=d,cλ为多项式系数,Bd(b(x))和c分别为和cλ的矩阵形式,u为笛卡尔坐标下方向导数的方向向量,a为与u相对应的面积坐标下的方向坐标,为分块对角矩阵,每对角块为对应于a和u的deCasteljau矩阵;其次根据相邻三角剖分间的r阶光滑拼接条件,每条边都有个约束条件,从而构成约束方程Hc=0,H为由E条边构成的E*Z行约束矩阵,且H的每行都包含一个单独的连续性约束,这样,波前传感器的斜率测量偏微分方程即可表示为:
Hc=0
其中,σu表示波前传感器斜率测量矩阵,n表示测量噪声及模型误差;
3)根据具体的参数值确定变换矩阵D,具体方法为:
步骤2)建立了一个约束最小二乘问题σu=Dc+n,Hc=0;与Zernike模式法相对应,通过对比此约束最小二乘问题与斜率测量偏微分方程,即可将变换矩阵D表示为:
4)结合相邻样条函数间的光滑拼接条件确定重构矩阵R;
确定重构矩阵R的方法具体为:
将由斜率测量数据去重构样条函数系数向量的矩阵定义为重构矩阵,结合步骤1)和步骤2),通过求解约束最小二乘问题,多项式系数向量c的显式解表达式为:
重构矩阵R,则表示为:
R=(DΤD)-1DΤ-(DΤD)-1HT[H(DΤD)-1HT]-1(DΤD)-1DΤ
其中,T表示矩阵转置,σ为波前斜率测量向量,此时由重构矩阵所确定的样条函数系数是在约束Hc=0下的最优线性无偏估计;
5)利用斜率测量数据计算各二元样条函数的系数向量
根据步骤4),重构矩阵是由斜率测量数据去重构样条函数系数向量的矩阵,二元样条函数的系数向量直接由重构矩阵与波前斜率测量向量相乘得到,则其计算公式为:
6)根据系数重构出畸变波前φ(x,y);
畸变波前φ(x,y)的重构方法为:
根据步骤2),已将畸变波前表示为单纯形剖分上的二元样条,步骤5)又得到了样条函数的系数向量而d次Bernstein多项式基函数Bd(b(x))是确定的,将二者相乘即可得到畸变波前的估计量即:
7)由重构的降质波前去估计大气瞬时点扩散函数h(x,y);
估计大气瞬时点扩散函数h(x,y)的方法具体为:
已知透镜光瞳函数为A,根据步骤6),由傅里叶光学理论,结合透镜光瞳函数A,得到其波前对应的点扩散函数为:
其中,FT-1表示傅里叶逆变换,
2.根据权利要求1所述的自适应光学成像中的点扩散函数估计方法,其特征在于:所述步骤1)中,根据微透镜阵列排布方式确定三角剖分的形式,确定三角剖分的形式时采用基于几何布局的网格单元直接剖分方法。
3.根据权利要求1所述的自适应光学成像中的点扩散函数估计方法,其特征在于:相邻三角剖分间的r阶光滑拼接条件描述如下:
设S(b(x))和分别是定义在相邻三角形t=<ν1,ν2,ν3>和上的d次多项式,νi(i=1,2,3)为三角形的顶点,cλ和分别是S(b(x))和关于t和的贝塞尔坐标,则S(b(x))和之间的r阶光滑拼接条件为:
其中,是关于t的面积坐标,λs=(s,λ1,λ2),λ0=(0,λ1,λ2),λ1+λ2=d-s。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310651641.9A CN103714516B (zh) | 2013-12-03 | 2013-12-03 | 一种自适应光学成像中的点扩散函数估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310651641.9A CN103714516B (zh) | 2013-12-03 | 2013-12-03 | 一种自适应光学成像中的点扩散函数估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103714516A CN103714516A (zh) | 2014-04-09 |
CN103714516B true CN103714516B (zh) | 2016-05-25 |
Family
ID=50407460
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310651641.9A Active CN103714516B (zh) | 2013-12-03 | 2013-12-03 | 一种自适应光学成像中的点扩散函数估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103714516B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106845024B (zh) * | 2017-03-02 | 2020-04-21 | 哈尔滨工业大学 | 一种基于波前反演的光学卫星在轨成像仿真方法 |
DE102019204575B3 (de) * | 2019-04-01 | 2020-08-06 | Carl Zeiss Smt Gmbh | Verfahren, Vorrichtung und Computerprogramm zum Bestimmen einer Wellenfront eines massebehafteten Teilchenstrahls |
CN111127473B (zh) * | 2019-11-07 | 2023-11-07 | 北京理工大学 | 一种基于优化策略的高速高精度区域分块波前重构方法 |
CN112905952A (zh) * | 2021-02-09 | 2021-06-04 | 南京信息工程大学 | 一种用于任意孔径光学元件波前梯度数据重构方法 |
CN113296175B (zh) * | 2021-05-25 | 2022-02-25 | 北京理工大学 | 一种加工具有多个数值孔径的微透镜阵列的方法 |
CN114970339A (zh) * | 2022-05-20 | 2022-08-30 | 大连理工大学 | 数据驱动识别偏微分方程的序列奇异值过滤方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102717305A (zh) * | 2012-06-15 | 2012-10-10 | 天津大学 | 一种光学自由曲面原位测量方法 |
-
2013
- 2013-12-03 CN CN201310651641.9A patent/CN103714516B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102717305A (zh) * | 2012-06-15 | 2012-10-10 | 天津大学 | 一种光学自由曲面原位测量方法 |
Non-Patent Citations (5)
Title |
---|
《Deconvolution from wave-front sensing: a new technique for compensating turbulence-degraded images》;J.Primot等;《J. Opt. Soc. Am. A》;19900930;第7卷(第9期);全文 * |
《自适应光学***的实时模式复原算法》;李新阳等;《强激光与粒子束》;20020131;第14卷(第1期);全文 * |
曾明等.《基于部分岭估计的Zernike模式法人眼波前像差重建》.《光电工程》.2013,第40卷(第7期),全文. * |
赵金宇.《采用波前测量图像恢复技术的误差分析》.《红外与激光工程》.2011,第40卷(第9期),全文. * |
陈波等.《面向图像复原的广义岭估计模式法波前重构》.《红外与激光工程》.2012,第41卷(第2期),全文. * |
Also Published As
Publication number | Publication date |
---|---|
CN103714516A (zh) | 2014-04-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103714516B (zh) | 一种自适应光学成像中的点扩散函数估计方法 | |
CN103310427B (zh) | 图像超分辨率与像质增强方法 | |
JP7406853B2 (ja) | コンピュータにより実行される法線ベクトルに基づく点群平滑化フィルタリング方法 | |
CN103258328B (zh) | 一种宽视场镜头的畸变中心定位方法 | |
CN103164846B (zh) | 一种基于转向核估计的红外条纹非均匀性校正方法 | |
CN103761712B (zh) | 基于自适应光学***点扩散函数重建的图像盲卷积方法 | |
JP2013502151A5 (zh) | ||
CN109752715B (zh) | 一种sar数据全散射体探测方法及装置 | |
CN104881866A (zh) | 一种扩展小孔成像模型的鱼眼相机矫正与标定方法 | |
Huang et al. | Shape reconstruction from gradient data in an arbitrarily-shaped aperture by iterative discrete cosine transforms in Southwell configuration | |
CN112629678B (zh) | 一种通用形状无衍射迭代计算的快速相位恢复方法 | |
CN112665529B (zh) | 基于条纹密度区域分割和校正的物体三维形貌测量方法 | |
CN108507564B (zh) | 一种基于点扩散函数拟合的星敏感器质心定位方法 | |
CN106707283B (zh) | 基于无味信息滤波的相位展开算法 | |
CN108230251A (zh) | 组合式图像恢复方法及装置 | |
CN104933713A (zh) | 一种利用边缘分析的图像mtf估计方法 | |
KR101726771B1 (ko) | 스트렐 비율을 이용한 변형거울의 평활화 방법 및 장치 | |
Ceruso et al. | Reconstructing wavefront phase from measurements of its slope, an adaptive neural network based approach | |
Singh et al. | Wavefront fitting and comparison of camera aberrations using Zernike circle polynomials | |
Rodríguez | Online self-camera orientation based on laser metrology and computer algorithms | |
US9582887B2 (en) | Methods and apparatus for determining field of view dependent depth map correction values | |
Rais et al. | Improving the accuracy of a Shack-Hartmann wavefront sensor on extended scenes | |
CN114329872A (zh) | 一种光学遥感相机椭率稳定度的计算仿真方法及*** | |
Idier et al. | Statistical behavior of joint least-square estimation in the phase diversity context | |
JP4984668B2 (ja) | ソーナーシステムおよびその位相誤差補正方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |