CN103714516A - 一种自适应光学成像中的点扩散函数估计方法 - Google Patents
一种自适应光学成像中的点扩散函数估计方法 Download PDFInfo
- Publication number
- CN103714516A CN103714516A CN201310651641.9A CN201310651641A CN103714516A CN 103714516 A CN103714516 A CN 103714516A CN 201310651641 A CN201310651641 A CN 201310651641A CN 103714516 A CN103714516 A CN 103714516A
- Authority
- CN
- China
- Prior art keywords
- matrix
- wavefront
- point spread
- estimation method
- spread function
- 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.)
- Granted
Links
Images
Landscapes
- Length Measuring Devices By Optical Means (AREA)
- Image Processing (AREA)
Abstract
本发明公开了一种自适应光学成像中的点扩散函数估计方法。该发明基于研究多元样条的B网方法,首先将经过大气降质的畸变波前近似表示为单纯形剖分上的二元样条,进而得到波前传感器斜率测量数据与样条函数系数之间的数学关系;其次,结合相邻三角剖分上样条函数的光滑拼接条件,利用斜率测量数据通过约束最小二乘估计估计样条函数系数进而重构畸变波前;最后将得到的畸变波前结合透镜光瞳函数,即可估计出经大气湍流成像的瞬时点扩散函数。本发明与传统的Zernike模式法估计点扩散函数相比克服了变换矩阵秩不完备的缺陷,提高了估计精度;并且适用于中心被遮挡的波前传感器及微透镜阵列的多种几何排布的情况。
Description
技术领域
本发明属于空间目标探测识别技术领域,涉及到空间目标自适应光学图像复原理论与方法,具体涉及到一种自适应光学成像中的点扩散函数估计方法。
背景技术
根据空间目标光学成像探测的特性,自适应光学图像复原方法主要可分为两类:盲解卷积和波前解卷积。对于盲解卷积,一个相对准确的点扩散函数估计值作为迭代算法的初始估计将会加快收敛、减少计算量;对于波前解卷积,点扩散函数估计值与真实值之间的偏差更是会直接影响复原的效果。因此快速准确的点扩散函数估计方法对于自适应光学图像的复原来讲意义显著。
点扩散函数与波前存在着理论上的对应关系,在自适应光学***中用波前传感器探测波前信息。目前使用H-S型波前传感器的自适应光学***居多,所获得数据是波前相位的斜率采样值。通过波前重构可以在有噪声及干扰的情况下,利用有限离散观测信息恢复出完整的波前,进而重构短曝光图像的瞬时点扩散函数。波前重构的基本方法大致上可以分为两大类:区域法和模式法。由于Zernike多项式各模式的像差意义明显,故Zernike模式法在面向图像复原的波前重构方法中一直得到广泛的应用。Zernike模式法最终归结为求解线性方程组,但由于Zernike函数偏导数的不完全正交性以及在有限采样点上函数的非正交性,在求最小二乘解时会出现变换矩阵D的秩不完备、DΤD奇异或接近奇异的情况,这将使梯度矢量的测量误差被放大,导致波前重构结果精度降低,进而影响点扩散函数的估计精度;其次,许多大型望远镜的中心常常被遮挡,若仍用Zernike多项式展开波前,因为展开域是环形的,所以Zernike展开不再具有正交性,仍应用此法则需要十分繁琐的数学推导与计算,这些因素都限制了Zernike模式法的应用效力。另外,虽然近年来在大型自适应光学设备快速发展的推动下不断涌现出新的面向校正控制的波前重构算法,但把这些方法进一步的应用到点扩散函数估计直至图像复原中的效力还并未得到有效的研究和验证。
发明内容
本发明的目的在于克服现有成像点扩散函数估计方法的不足,提出一种自适应光学成像中的点扩散函数估计方法,该发明基于研究多元样条的B网方法,将经过大气降质的畸变波前近似表示为单纯形剖分上的二元样条,避免了求最小二乘解时不完备或奇异矩阵的出现,提高了估计精度,并可以应对多种波前传感器微透镜阵列的几何排布方式。
为了达到上述目的,本发明所采用的技术方案包括以下步骤:
1)指定样条函数的次数d及刻画光滑程度的连续性阶数r,确定波前区域内三角剖分的形式;
2)建立斜率测量模型的偏微分方程并确定方程中各参数矩阵的值;
3)根据具体的参数值确定变换矩阵D;
4)结合相邻样条函数间的光滑拼接条件确定重构矩阵R;
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的de Casteljau矩阵;其次根据相邻三角剖分间的r阶光滑拼接条件,每条边都有个约束条件,从而构成约束方程Hc=0,H为由E条边构成的E*Z行约束矩阵,且H的每行都包含一个单独的连续性约束,这样,波前传感器的斜率测量偏微分方程即可表示为:
Hc=0
其中,σu表示波前传感器斜率测量矩阵,n表示测量噪声及模型误差。
相邻三角剖分间的r阶光滑拼接条件描述如下:
所述步骤3)中,确定变换矩阵D的方法具体为:
步骤2)建立了一个约束最小二乘问题σu=Dc+n,Hc=0;与Zernike模式法相对应,通过对比此约束最小二乘问题与斜率测量偏微分方程,即可将变换矩阵D表示为:
所述步骤4)中,确定重构矩阵R的方法具体为:
将由斜率测量数据去重构样条函数系数向量的矩阵定义为重构矩阵,结合步骤1)和步骤2),通过求解约束最小二乘问题,多项式系数向量c的显式解表达式为:
重构矩阵R,则表示为:
R=(DΤD)-1DΤ-(DΤD)-1HT[H(DΤD)-1HT]-1(DΤD)-1DΤ
其中,Τ表示矩阵转置,σ为波前斜率测量向量,此时由重构矩阵所确定的样条函数系数是在约束Hc=0下的最优线性无偏估计。
所述步骤5)中,根据步骤4),重构矩阵是由斜率测量数据去重构样条函数系数向量的矩阵,二元样条函数的系数向量直接由重构矩阵与波前斜率测量向量相乘得到,则其计算公式为:
所述步骤6)中,畸变波前φ(x,y)的重构方法为:
所述步骤7)中,估计大气瞬时点扩散函数h(x,y)的方法具体为:
已知透镜光瞳函数为A,根据步骤6),由傅里叶光学理论,结合透镜光瞳函数A,得到其波前对应的点扩散函数为:
与现有技术相比,本发明具有以下有益效果:
本发明在重构波前时满足变换矩阵D满秩、DΤD非奇异的条件,不会放大斜率测量误差,并且基于二元样条的方法可以选择使用非线性基函数,这都促使波前重构精度的提高,进而可以提高点扩散函数的估计精度;另外,本发明不受波前传感器微透镜阵列几何排布方式的限制,对圆形或方形布局都有效,并且还适用于中心被遮挡的情况,方法的通用性较强。
附图说明
图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的de Casteljau矩阵;其次根据相邻三角剖分间的r阶光滑拼接条件,每条边都有个约束条件,从而构成约束方程Hc=0,H为由E条边构成的E?Z行约束矩阵,且H的每行都包含一个单独的连续性约束,这样,波前传感器的斜率测量偏微分方程即可表示为:
Hc=0
其中,σu表示波前传感器斜率测量矩阵,n表示测量噪声及模型误差。
相邻三角剖分间的r阶光滑拼接条件描述如下:
设S(b(x))和分别是定义在相邻三角形t=<ν1,ν2,ν3>和上的d次多项式,νi(i=1,2,3)为三角形的顶点,cλ和分别是S(b(x))和关于t和的贝塞尔坐标,则S(b(x))和之间的r阶光滑拼接条件为:
3)根据具体的参数值确定变换矩阵D;
确定变换矩阵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Τ
其中,Τ表示矩阵转置,σ为波前斜率测量向量,此时由重构矩阵所确定的样条函数系数是在约束Hc=0下的最优线性无偏估计。
5)利用斜率测量数据计算各二元样条函数的系数向量?c;
畸变波前φ(x,y)的重构方法为:
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的de Casteljau矩阵;其次根据相邻三角剖分间的r阶光滑拼接条件每条边都有个约束条件,从而构成约束方程Hc=0,H为由E条边构成的E*Z行的约束矩阵。这样,波前传感器的斜率测量偏微分方程即可表示为:
Hc=0
其中σu表示波前传感器斜率测量值,n表示测量噪声及模型误差。模型建立好其相应的Bd-1、与H阵也已随即确定;
2)确定变换矩阵D。与Zernike模式法类似,在此也将样条函数系数与斜率测量数据之间的变换关系用变换矩阵D来描述。约束最小二乘问题的模型为σu=Dc+n;Hc=0,结合步骤1)则有:
3)计算重构矩阵R。将由斜率测量数据去重构样条函数系数向量的矩阵定义为重构矩阵。结合步骤1)和步骤2),由约束最小二乘问题的显式解表达式:
R=(DΤD)-1DΤ-(DΤD)-1HT[H(DΤD)-1HT]-1(DΤD)-1DΤ (4)
此时由重构矩阵所确定的样条函数系数是在约束Hc=0下的最优线性无偏估计。
5)重构畸变波前φ(x,y)。由步骤1):
6)估计点扩散函数h(x,y)。在步骤5)的基础上,由傅里叶光学理论,结合透镜光瞳函数A,其波前对应的点扩散函数为:
其中FT-1表示傅里叶逆变换。至此,则完成了成像点扩散函数的估计。
以上实施例仅用以说明本发明而并非限制本发明所描述的技术方案;因此尽管本说明书参照上述的实施例对本发明已进行了详细的说明,但是本领域的普通技术人员应当理解,仍然可以对本发明进行修改或等同替换;而一切不脱离发明的精神和范围的技术方案及其改进,均应涵盖在本发明的权利要求范围之内。
Claims (9)
2.根据权利要求1所述的自适应光学成像中的点扩散函数估计方法,其特征在于:所述步骤1)中,根据微透镜阵列排布方式确定三角剖分的形式,确定三角剖分的形式时采用基于几何布局的网格单元直接剖分方法。
3.根据权利要求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的de Casteljau矩阵;其次根据相邻三角剖分间的r阶光滑拼接条件,每条边都有个约束条件,从而构成约束方程Hc=0,H为由E条边构成的E*Z行约束矩阵,且H的每行都包含一个单独的连续性约束,这样,波前传感器的斜率测量偏微分方程即可表示为:
Hc=0
其中,σu表示波前传感器斜率测量矩阵,n表示测量噪声及模型误差。
5.根据权利要求3所述的自适应光学成像中的点扩散函数估计方法,其特征在于:所述步骤3)中,确定变换矩阵D的方法具体为:
步骤2)建立了一个约束最小二乘问题σu=Dc+n,Hc=0;与Zernike模式法相对应,通过对比此约束最小二乘问题与斜率测量偏微分方程,即可将变换矩阵D表示为:
6.根据权利要求1所述的自适应光学成像中的点扩散函数估计方法,其特征在于:所述步骤4)中,确定重构矩阵R的方法具体为:
将由斜率测量数据去重构样条函数系数向量的矩阵定义为重构矩阵,结合步骤1)和步骤2),通过求解约束最小二乘问题,多项式系数向量c的显式解表达式为:
重构矩阵R,则表示为:
R=(DΤD)-1DΤ-(DΤD)-1HT[H(DΤD)-1HT]-1(DΤD)-1DΤ
其中,Τ表示矩阵转置,σ为波前斜率测量向量,此时由重构矩阵所确定的样条函数系数是在约束Hc=0下的最优线性无偏估计。
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 true CN103714516A (zh) | 2014-04-09 |
CN103714516B 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) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106845024A (zh) * | 2017-03-02 | 2017-06-13 | 哈尔滨工业大学 | 一种基于波前反演的光学卫星在轨成像仿真方法 |
CN111127473A (zh) * | 2019-11-07 | 2020-05-08 | 北京理工大学 | 一种基于优化策略的高速高精度区域分块波前重构方法 |
CN112905952A (zh) * | 2021-02-09 | 2021-06-04 | 南京信息工程大学 | 一种用于任意孔径光学元件波前梯度数据重构方法 |
CN113296175A (zh) * | 2021-05-25 | 2021-08-24 | 北京理工大学 | 一种加工具有多个数值孔径的微透镜阵列的方法 |
TWI749498B (zh) * | 2019-04-01 | 2021-12-11 | 德商卡爾蔡司Smt有限公司 | 確定質量粒子束波前的方法和裝置 |
CN116504341A (zh) * | 2022-05-20 | 2023-07-28 | 大连理工大学 | 数据驱动识别偏微分方程的序列奇异值过滤方法 |
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 |
---|
J.PRIMOT等: "《Deconvolution from wave-front sensing: a new technique for compensating turbulence-degraded images》", 《J. OPT. SOC. AM. A》 * |
曾明等: "《基于部分岭估计的Zernike模式法人眼波前像差重建》", 《光电工程》 * |
李新阳等: "《自适应光学***的实时模式复原算法》", 《强激光与粒子束》 * |
赵金宇: "《采用波前测量图像恢复技术的误差分析》", 《红外与激光工程》 * |
陈波等: "《面向图像复原的广义岭估计模式法波前重构》", 《红外与激光工程》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106845024A (zh) * | 2017-03-02 | 2017-06-13 | 哈尔滨工业大学 | 一种基于波前反演的光学卫星在轨成像仿真方法 |
CN106845024B (zh) * | 2017-03-02 | 2020-04-21 | 哈尔滨工业大学 | 一种基于波前反演的光学卫星在轨成像仿真方法 |
TWI749498B (zh) * | 2019-04-01 | 2021-12-11 | 德商卡爾蔡司Smt有限公司 | 確定質量粒子束波前的方法和裝置 |
US11355311B2 (en) | 2019-04-01 | 2022-06-07 | Carl Zeiss Smt Gmbh | Method and apparatus for determining a wavefront of a massive particle beam |
CN111127473A (zh) * | 2019-11-07 | 2020-05-08 | 北京理工大学 | 一种基于优化策略的高速高精度区域分块波前重构方法 |
CN111127473B (zh) * | 2019-11-07 | 2023-11-07 | 北京理工大学 | 一种基于优化策略的高速高精度区域分块波前重构方法 |
CN112905952A (zh) * | 2021-02-09 | 2021-06-04 | 南京信息工程大学 | 一种用于任意孔径光学元件波前梯度数据重构方法 |
CN113296175A (zh) * | 2021-05-25 | 2021-08-24 | 北京理工大学 | 一种加工具有多个数值孔径的微透镜阵列的方法 |
CN113296175B (zh) * | 2021-05-25 | 2022-02-25 | 北京理工大学 | 一种加工具有多个数值孔径的微透镜阵列的方法 |
CN116504341A (zh) * | 2022-05-20 | 2023-07-28 | 大连理工大学 | 数据驱动识别偏微分方程的序列奇异值过滤方法 |
CN116504341B (zh) * | 2022-05-20 | 2023-11-07 | 大连理工大学 | 数据驱动识别偏微分方程的序列奇异值过滤方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103714516B (zh) | 2016-05-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103714516A (zh) | 一种自适应光学成像中的点扩散函数估计方法 | |
JP7406853B2 (ja) | コンピュータにより実行される法線ベクトルに基づく点群平滑化フィルタリング方法 | |
CN104459696A (zh) | 一种基于平地相位的sar干涉基线精确估计方法 | |
CN102208108B (zh) | 摄像机大视场高精度快速现场全局标定方法 | |
CN103258328B (zh) | 一种宽视场镜头的畸变中心定位方法 | |
CN103310427B (zh) | 图像超分辨率与像质增强方法 | |
CN104182982A (zh) | 双目立体视觉摄像机标定参数的整体优化方法 | |
CN109752715B (zh) | 一种sar数据全散射体探测方法及装置 | |
CN107014313B (zh) | 基于s变换脊值的加权最小二乘相位展开的方法及*** | |
WO2020240720A1 (ja) | 合成開口レーダの信号処理装置および信号処理方法 | |
CN109376422B (zh) | 一种均匀圆阵优化设计评估方法和装置 | |
CN112629678B (zh) | 一种通用形状无衍射迭代计算的快速相位恢复方法 | |
CN107544069B (zh) | 基于平面近似模型的多基线相位解缠绕方法 | |
CN112665529B (zh) | 基于条纹密度区域分割和校正的物体三维形貌测量方法 | |
CN104915549A (zh) | 一种基于半参数模型的InSAR干涉相位解缠方法 | |
CN108230251A (zh) | 组合式图像恢复方法及装置 | |
CN107402074B (zh) | 一种基于光学稀疏子孔径区域内标准正交多项式计算进行重构波前的方法 | |
CN113532308B (zh) | 数字图像相关中带初值的岭回归应变测量方法 | |
Ceruso et al. | Reconstructing wavefront phase from measurements of its slope, an adaptive neural network based approach | |
CN115507737A (zh) | 时序InSAR干涉像对的选取方法及装置 | |
Guo et al. | Sequential deconvolution from wave-front sensing using bivariate simplex splines | |
Vyas et al. | Advanced methods for improving the efficiency of a Shack Hartmann wavefront sensor | |
CN103674815A (zh) | 一种基于传感器特征的气溶胶mtf监测方法 | |
CN106228526A (zh) | 一种基于中频的衍射极限模糊图像盲复原方法 | |
CN102322826A (zh) | 一种物体三维表面数据的改进的测量方法 |
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 |