CN105629696B - 一种基于迭代去噪收缩阈值算法的数字全息重构方法 - Google Patents
一种基于迭代去噪收缩阈值算法的数字全息重构方法 Download PDFInfo
- Publication number
- CN105629696B CN105629696B CN201610021891.8A CN201610021891A CN105629696B CN 105629696 B CN105629696 B CN 105629696B CN 201610021891 A CN201610021891 A CN 201610021891A CN 105629696 B CN105629696 B CN 105629696B
- Authority
- CN
- China
- Prior art keywords
- hologram
- algorithm
- delta
- matrix
- light
- 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.)
- Expired - Fee Related
Links
- 238000004422 calculation algorithm Methods 0.000 title claims abstract description 64
- 238000000034 method Methods 0.000 title claims abstract description 34
- 239000011159 matrix material Substances 0.000 claims abstract description 38
- 230000010363 phase shift Effects 0.000 claims abstract description 19
- 238000012545 processing Methods 0.000 claims abstract description 6
- 241000350158 Prioria balsamifera Species 0.000 claims description 6
- 230000000694 effects Effects 0.000 claims description 5
- 238000005259 measurement Methods 0.000 claims description 4
- 238000005070 sampling Methods 0.000 claims description 4
- 238000005457 optimization Methods 0.000 claims description 3
- ONCZDRURRATYFI-QTCHDTBASA-N methyl (2z)-2-methoxyimino-2-[2-[[(e)-1-[3-(trifluoromethyl)phenyl]ethylideneamino]oxymethyl]phenyl]acetate Chemical compound CO\N=C(/C(=O)OC)C1=CC=CC=C1CO\N=C(/C)C1=CC=CC(C(F)(F)F)=C1 ONCZDRURRATYFI-QTCHDTBASA-N 0.000 claims 3
- 238000004364 calculation method Methods 0.000 abstract description 5
- 230000008602 contraction Effects 0.000 abstract 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 abstract 1
- 230000006870 function Effects 0.000 description 19
- 238000001093 holography Methods 0.000 description 10
- 230000006835 compression Effects 0.000 description 6
- 238000007906 compression Methods 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 5
- 230000003287 optical effect Effects 0.000 description 5
- 238000003384 imaging method Methods 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 238000010587 phase diagram Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 241000208140 Acer Species 0.000 description 1
- 241000272168 Laridae Species 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 244000005700 microbiome Species 0.000 description 1
- 230000008447 perception Effects 0.000 description 1
- GGCZERPQGJTIQP-UHFFFAOYSA-N sodium;9,10-dioxoanthracene-2-sulfonic acid Chemical compound [Na+].C1=CC=C2C(=O)C3=CC(S(=O)(=O)O)=CC=C3C(=O)C2=C1 GGCZERPQGJTIQP-UHFFFAOYSA-N 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 238000002834 transmittance Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
- G03H—HOLOGRAPHIC PROCESSES OR APPARATUS
- G03H1/00—Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
- G03H1/04—Processes or apparatus for producing holograms
- G03H1/08—Synthesising holograms, i.e. holograms synthesized from objects or objects from holograms
- G03H1/0808—Methods of numerical synthesis, e.g. coherent ray tracing [CRT], diffraction specific
-
- G—PHYSICS
- G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
- G03H—HOLOGRAPHIC PROCESSES OR APPARATUS
- G03H1/00—Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
- G03H1/04—Processes or apparatus for producing holograms
- G03H1/08—Synthesising holograms, i.e. holograms synthesized from objects or objects from holograms
- G03H1/0866—Digital holographic imaging, i.e. synthesizing holobjects from holograms
-
- G—PHYSICS
- G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
- G03H—HOLOGRAPHIC PROCESSES OR APPARATUS
- G03H1/00—Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
- G03H1/04—Processes or apparatus for producing holograms
- G03H1/08—Synthesising holograms, i.e. holograms synthesized from objects or objects from holograms
- G03H1/0808—Methods of numerical synthesis, e.g. coherent ray tracing [CRT], diffraction specific
- G03H2001/0816—Iterative algorithms
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Holo Graphy (AREA)
Abstract
本发明公开了一种基于迭代去噪收缩阈值算法的数字全息重构方法:1、读入原始图像,其振幅为O0(ξ,η),利用此原始图像的物光信息获得衍射光复振幅O(x,y);2、由菲涅尔衍射得到的衍射光复振幅O(x,y)与加入了0和π/2相移量后的参照光R(x,y)相干涉到达全息干板上的光波复振幅分别为U1(x,y)和U2(x,y),经过数值处理形成两幅全息图I(1)和I(2);3、对两幅全息图I(1)和I(2)进行叠加生成一幅相移全息图构造传感矩阵A,使得其中Φ为已知的测量矩阵,R表示稀疏基矩阵,y表示得到的观测数据;4、利用迭代去噪收缩阈值算法(IDNST)求解获得重构原始物体的图像本发明的IDNST算法在(TwIST)算法的基础上引入阈值和正则化参数的双收缩,信噪比得以提高,并且加快了收敛速度,提高了重构精度,使得再现质量达到更优水平。
Description
技术领域
本发明涉及数字全息技术领域,特别涉及一种基于迭代去噪收缩阈值算法的数字全息重构方法。
背景技术
压缩感知技术应用于数字全息技术中,和传统的光学全息术相比,具有以下优点:1)是由于用电子成像器件代替全息干板来记录全息图,所需的曝光时间可以很短,而没有化学银盐干版的湿处理过程,可连续记录运动物体的各个瞬间过程,有利于实现实时在线全息记录;2)是数字计算得到的再现图样可通过计算机进行定量分析。数字再现全息图得到的是物场的复振幅分布,并可通过计算得到强度和相位分布,易于实现两个或多个全息图的相加减、增减背景图像、叠加图像[1]等操作,而这些在光学全息中很难做到,由此也显现了数字全息技术更具实用性;3)将压缩感知理论与数字全息图再现相结合,消除了传统算法对数字全息重构过程中出现的0极像、共轭像干扰[2,3]的问题,对实现计算全息实时性具有重要的意义。
在后续工作中针对稀疏表示中还可以考虑冗余字典上的稀疏分解理论,进一步减少测量值提高重构效果。目前国内外已经有一系列的工作针对基于压缩感知的数字全息成像问题展开。文献[4]研究了2种不同变换域对全息图稀疏化的影响,根据压缩感知理论直接获取图像的压缩表示,并且利用傅里叶稀疏域下使用正交匹配追踪(OMP)重构算法能获得1.5%的压缩率。文献[5]将压缩感知理论与数字全息图再现相结合提出基于全变差的重构算法,并应用于数字全息图压缩感知全息图重建,消除了传统算法对数字全息重构过程中会出现0级像、共轭像干扰的问题。文献[6]将压缩感知应用在相移全息图中,通过利用可变密度采样结构将更多信号积聚在中心位置,很少的信号分布在全息图的边缘,利用全变分求解最优值,从而获得再现像。文献[7]验证了压缩感知技术从低维数据中产生高维图像的重要性,并且阐明了利用压缩感知2D全息数据中重构出一幅3D层析图。文献[8]详细阐明了基于全变分的两步迭代收缩阈值的重构算法,尤其说明了对于病态的问题,相比IST算法,此法对于目标函数求解最优化值有很快的收敛速度。文献[9]介绍了基于压缩感知的次曝光同轴全息结构,能够捕捉微生物的行为动态,相比在线数字全息术,此方法只需要一次曝光特性。文献[10]主要将压缩感知技术应用在获得多个视图全息投影,并且通过对生成的全息图二次抽样重构3D场景。
现有技术中,数字全息图像重构过程中,收敛速度慢,重构精度低。
参考文献
[1]T.-C.Poon and P.P.Banerjee,Contemporary Optical Image Processingwith MATLAB(Elsevier,2001).
[2]Yair Rivenson and Adrian Stern.”Conditions for practicingCompressive Fresnel holography,”Opt.Lett.36,3365-3367(2010).
[3]V.Mico,J.Garc1a,Z.Zalevsky,and B.Javidi,“Phase-shifting Gaborholography,”Opt.Express 17,1492–1494(2009).
[4]李科,李军.基于压缩传感的全息图压缩研究[J].华南师范大学学报,2012,44(4):61-65.
[5]简献忠,周海,乔静远,李莹,王佳.基于全变差重构算法的数字全息研究[J].激光技术,
2014,38(2):236-239.
[6]Yair Rivenson,Adrian Sterna and Bahram,Javidi.”Compressive FresnelHolography,”J.Display Technol,vol.6,NO.10,pp.506-509,Oct.2010.
[7]David J.Brady,Kerkil Choi,Daniel L.Marks,Ryoichi Horisaki andSehoon Lim.”Compressive Holography,”Opt.Express 17(15),13040-13049.
[8]J.M.Bioucas-Dias and M.A.T.Figueiredo,“A new TwIST:two-stepiterative shrinkage/thresholding algorithms for image restoration,”IEEETrans.Image Process.16,2992–3004(2007).
[9]Y.Rivenson,A.Stern and B.Javidi,“Improved depth resolution bysingle-exposure in-line compressive holography,”Appl.Opt.52,A223–A229(2013).
[10]Yair Rivenson,Adrian Stern and Joseph Rosen.”Compressive multipleview projection incoherent holography,”Opt.Express 19(7),6109-6118(2011).
发明内容
本发明的目的在于提供一种基于迭代去噪收缩阈值算法的数字全息重构方法;在两步迭代收缩阈值算法(TwIST)的基础上,引入阈值和正则化参数的双收缩,提出了迭代去噪收缩阈值(IDNST)算法;利用前两次迭代的值和不断更新的参数t以及不断收缩的正则化参数来获得新的迭代值,提高了算法的重构精度,缩短了迭代时间,并且加强了抗噪能力。本发明提出的重构算法IDNS,加快了收敛速度,提高了图像的重构精度,使得再现质量达到更优水平。
为了实现上述目的,本发明采用如下技术方案:
一种基于迭代去噪收缩阈值算法的数字全息重构方法,包括以下步骤:
步骤一:读入原始图像,其振幅为O0(ξ,η),利用此原始图像的物光信息进行菲涅尔衍射得到距离孔径平面为z的观察平面上任意一点P(x,y)的衍射光复振幅O(x,y);
步骤二:由菲涅尔衍射得到的衍射光复振幅O(x,y)与加入了0和π/2相移量后的参照光R(x,y)相干涉到达全息干板上的光波复振幅分别为U1(x,y)和U2(x,y),经过数值处理形成两幅全息图I(1)和I(2);
步骤三:对两幅全息图I(1)和I(2)进行叠加得到相移全息图构造传感矩阵A,使得y=ΦI=ΦRO'=AO',其中Φ为已知的测量矩阵,R表示稀疏基矩阵,y表示得到的观测数据;
步骤四:利用迭代去噪收缩阈值算法重构出原始图像
进一步的,步骤一中:
其中,λ是光的波长,k是波数,z为物体所在平面与全息面的距离。
进一步的,步骤二具体包括以下步骤:
步骤2.1:衍射光O(x,y)与参考光R(x,y)相干涉在全息记录平面xy上所叠加的总光场分布U(x,y)为:
U(x,y)=O(x,y)+R(x,y)
步骤2.2:利用全息记录平面xy上的总光场U(x,y)计算形成的全息图光强分布为:
I=U(x,y)·U*(x,y)=(O+R)(O*+R*)
=|O|2+|R|2+O·R*+O*·R
式中I代表全息图I(1)和I(2);上式中U(x,y)为U1(x,y)时,计算获得全息图I(1);U(x,y)为U2(x,y)时,计算获得全息图I(2);
步骤2.1中,在参考光R(x,y)的相移量为0时,获得U1(x,y);在参考光R(x,y)的相移量为π/2时,获得U2(x,y)。
进一步的,步骤三中:传感矩阵A:
步骤3.1:将步骤二中得到的两幅全息图I(1)和I(2)叠加得
步骤3.2:根据步骤二,CCD上的光强分布I(x,y)具体为
式中,α表示CCD量子系数,I0(x,y)表示CCD背景噪音。式中第一项(α+I0)表示全息图的直流分量,可以忽略;第二项Re表示实部,2αRe(o*r)表示物体干涉形成的全息图,其中o和r分别为衍射光和参考光的振幅大小;第三项α(o*r)(o*r)*表示物体之间对干涉的影响;
步骤3.3:α(o*r)(o*r)*<<2αRe(o*r),因而可以忽略。因而全息图可以简化为:
式中F、F-1分别表示二维傅里叶变换和傅里叶逆变换。
步骤3.4:假设CCD的分辨力为Nx×Ny,像素大小为Δx×Δy,被记录的目标空间o(x0,y0)被划分为步长为Δx×Δy的Nx×Ny个采样区间,对上式进行离散化,得
步骤3.5:为了应用压缩感知理论,将目标空间二维矩阵与全息图化为一维向量,定义则上式简化为
I′=T2DWQBO′=RO′
其中,和O'分别为式中和的简写;式中B=F2D是大小为(Nx×Ny)×(Nx×Ny)二维分块对角阵,表示二维离散傅里叶变换;BO'对应式中F[o(mΔx,nΔy,qΔz)];Q为(Nx×Ny)×(Nx×Ny)的对角矩阵,其[(n-1)×Nx+m]为F[r(kΔx-mΔx,lΔy-nΔy)]矩阵第n行m列元素;F[r(kΔx-mΔx,lΔy-nΔy)]表示菲涅尔变换核r(mΔx,nΔy)的傅里叶变换;W=I1,为(Nx×Ny)×(Nx×Ny)单位矩阵;T2D大小为(Nx×Ny)×(Nx×Ny),表示二维傅里叶逆变换;R=T2DWQB为压缩感知理论的稀疏矩阵,大小为(Nx×Ny)×(Nx×Ny),则A=ΦR=ΦT2DWQB。
进一步的,步骤四具体包括以下步骤:
步骤4.1:初始化;初始图像O'0=0,初始两步迭代次数TwIST_iters=0,迭代次数t1=1,最大迭代次数Mt=1000,算法终止条件值tolA=10-3,正则化收缩因子u=0.9,计数器s=1,以及各种临时参数O'1=O'0,O'2=O'0;
步骤4.2:定义目标函数为
根据步骤三得
其中A是传感矩阵,y表示得到的观测数据,I'代表是步骤三获得的叠加生成的相移全息图;τ代表正则化参数,ΘTV(I')为总变分函数,其公式如下:
和代表对于图像像素值I'的一阶水平和垂直差分操作
步骤4.3:在最优化问题中将式转化为最小化式
ψλ(O')=arg min{X}
并计算O'0对应的目标函数X(O'0),并且赋值pref=X(O'0);
步骤4.4:对O'0进行去噪处理,O'1=Γλ(O'0),其中,
Γλ(O')=Ψλ(O'+ΦT(y-ΦO'))
Ψλ为去噪软收缩函数,在此引入迭代因子t,则
从而
并计算O'k+1对应的目标函数X(O'k+1);
步骤4.5:判断算法是否已经进行过两步迭代;
如果TwIST_iters≠0不成立,则计算O'的目标函数值X(O');之后再判断条件
X(O'k)>X(O'k+1)是否成立;
若X(O'k)>X(O'k+1)成立,则s=s×2;并判断条件s>1000是否成立,若成立,则终止程序,最终输出估计值O'k,若不成立,则赋值TwIST_iters=0并转至步骤4.4;
若X(O'k)>X(O'k+1)不成立,则赋值TwIST_iters=TwIST_iters+1,并执行步骤4.6;
如果TwIST_iters≠0成立,则两步估计O'k+1=(1-α)O'k-1+(α-β)O'k+βΓλ(O'k),同时计算O'k+1的目标函数值X(O'),之后再判断X(O')>pref是否成立;若成立则赋值TwIST_iters=0,并转至步骤4.4;若不成立则赋值TwIST_iters=TwIST_iters+1,O'=O'k+1,并执行步骤4.6;
步骤4.6:判断算法终止条件;判断条件迭代次数s>1000是否成立,若成立则终止算法,输出最终估计值O';若不成立则赋值O'k=O'k-1,O'k+1=O'k,并将迭代次数加1,并计算算法终止条件并赋值X(O'k)=X(O'k-1);
步骤4.7:判断C(O'k,O'k-1)≤tolA和t>Mt是否成立:
若都不成立,则收缩正则化参数,τ=τ×u,转至步骤4.4;
只要有一条件成立,则终止算法;最终重构出原始图像
相对于现有技术,本发明具有以下有益效果:
为了解决全息图像数据在传输过程中占用大量内存并在一定程度上增加设计成本的问题,本发明方法基于压缩感知对全息图压缩采样再通过通信信道传输后重构恢复出原始图像。
因此首先由物体光场分布的复振幅O0(ξ,η)利用菲涅尔衍射,得到距离孔径平面为z的观察平面衍射光场中任意一点P(x,y)的复振幅O(x,y);由得到的物光波衍射光O(x,y)与加入了0和π/2相移量后的参照光R(x,y)相干涉到达全息干板上的光波复振幅分别为U1(x,y)和U2(x,y),经过数值处理形成两幅全息图I(1)和I(2);对两幅全息图I(1)和I(2)进行叠加生成一幅相移全息图构造传感矩阵A,使得其中Φ为已知的测量矩阵,R表示稀疏基矩阵,y表示得到的观测数据;引入总变分函数ΘTV(x)去掉噪声影响,同时利用迭代去噪收缩阈值(IDNST)重构算法求解获得最后的复原图像本发明在实际应用中大量减少占用内存,实用性强。仿真结果表明,本发明所提出的基于迭代去噪收缩阈值压缩全息方法,能够高概率的恢复出原始图像。
本发明用电子成像器件代替全息干板来记录全息图,所需的曝光时间可以很短,同时将压缩感知理论与数字全息图再现相结合,消除了传统算法对数字全息重构过程中出现的0极像、共轭像干扰的问题,对实现计算全息实时性具有重要的意义;另一方面运用我们所提出的基于迭代去噪收缩阈值(IDNST)重构算法成功地由二维的编码像复原出了原始图像的数据,通过引入阈值和正则化参数的双收缩减小了迭代运行时长,提高了算法的重构精度和抗噪能力,根据图8所示,与TwIST算法相比重构误差明显减小,重构效果达到更优。
附图说明
图1(a)为原始图像,图1(b)和图1(c)为利用原有的TwIST算法生成的全息图1和全息图2;
图2(a)为原始图像,图2(b)和图2(c)为利用本发明IDNST算法生成的全息图1和全息图2;
图3为本发明方法的流程图;
图4为IDNST重构算法流程图;
图5(a)、(b)和(c)为利用IDNST算法重建光的幅值与相位图和复原的图像;
图6(a)、(b)和(c)为利用原有的TwIST算法重建光的幅值与相位图和复原的图像;
图7为采用本发明IDNST算法与原有TwIST算法的信噪比对比图;
图8为采用本发明IDNST算法与原有TwIST算法的均方误差对比图。
图9为数字全息记录与再现的直观立体图。
具体实施方式
一种基于迭代去噪收缩阈值算法的数字全息重构方法,包括以下步骤:
步骤一:由已知物体所在平面,物体自身的物光波由该平面衍射出,设该平面为x0y0平面,物体光场分布的振幅O0(ξ,η)利用菲涅尔衍射,根据惠更斯-菲涅尔子波干涉原理得到距离孔径平面为z的观察平面上任意一点P(x,y)的衍射光复振幅O(x,y),并将其二次相位因子展开,得到菲涅尔衍射的傅里叶变换表达式;读入原始图像,其振幅为O0(ξ,η),利用此原始图像的物光信息进行菲涅尔衍射得到在观测平面上的复振幅O(x,y);
步骤二:由得到的物光波衍射光复振幅O(x,y)与加入了0和π/2相移量后的参照光R(x,y)相干涉到达全息干板上的光波复振幅分别为U1(x,y)和U2(x,y),经过数值处理形成两幅全息图I(1)和I(2);
步骤三:对两幅全息图I(1)和I(2)进行叠加生成一幅相移全息图构造传感矩阵A,使得其中Φ为已知的测量矩阵,R表示稀疏基矩阵,y表示得到的观测数据;
步骤四:引入总变分函数ΘTV(x)去掉随机噪声对重构图像的影响,同时利用基于迭代去噪收缩阈值(IDNST)重构算法求解获得重构出原始图像
步骤一中物体物光波的振幅O0(ξ,η)利用菲涅尔衍射,得到距离孔径平面为z的观察平面上任意一点P(x,y)的衍射光复振幅O(x,y),具体包括以下步骤:
步骤1.1:衍射光场在观察平面中任意一点的复振幅表示成在孔径平面S上的面积积分形式:
其中,O0(ξ,η)是自身物体物光波的振幅,O(x,y)是观察平面衍射光的复振幅,K(θ)是倾斜因子,λ是光的波长,k是波数,z是物体所在平面与全息面的距离,r是物体所在平面上的点到全息面上的点的距离。
步骤1.2:在基尔霍夫边界条件下,孔径平面上孔径S以外的区域内光场为0,因此可以将上式积分扩展到无穷。在傍轴条件下K(θ)近似为1,再利用二项式近似
从而可以得到菲涅尔衍射公式
步骤1.3:根据傅立叶变换的定义,可以将上式表达为
其中称是参与记录过程的系统传递函数。
步骤二中分别在参考光R(x,y)中加入了0和π/2相移量后与衍射光复振幅O(x,y)干涉形成两幅同轴全息图I(1)和I(2)的具体步骤如下:
步骤2.1:在参考光R(x,y)的相移量分别为0和π/2时,则衍射光O(x,y)与参考光R(x,y)相干涉在全息记录平面xy上所叠加的总光场分布U(x,y)(相移量为0时通过下式计算获得U1(x,y),相移量为π/2时通过下式计算获得U2(x,y))为:
U(x,y)=O(x,y)+R(x,y)
步骤2.2:利用全息记录平面xy上的总光场U(x,y)计算形成的全息图光强分布为:
I(x,y)=U(x,y)·U*(x,y)=(O+R)(O*+R*)
=|O|2+|R|2+O·R*+O*·R
式中I(x,y)代表全息图I(1)和I(2)(上式中U(x,y)为U1(x,y)时,计算获得全息图I(1);U(x,y)为U2(x,y)时,计算获得全息图I(2))。上式中,O和R为O(x,y)和R(x,y)的简写。第三项和第四项包含了参考光波与物光波相互干涉后所形成干涉条纹的信息,干涉条纹不仅包含了物光波的振幅信息同时也将相位信息记录在了其中。
步骤2.3:在整个全息图物光波的波前记录过程中,按照一定的数量关系使得入射光光强能够透过全息记录平面,则称该数量关系为记录平面的透过率函数:
TH(x,y)=T0+β·I(x,y)
其中T0和β都是常数,所以在分析中,我们将无关紧要的T0和β省略。
步骤三中构造传感矩阵A,具体步骤如下:
步骤3.1:将步骤二中得到的两幅全息图I(1)和I(2)叠加得
步骤3.2:根据步骤二,CCD上的光强分布具体为
式中,α表示CCD量子系数,I0(x,y)表示CCD背景噪音。式中第一项(α+I0)表示全息图的直流分量,可以忽略;第二项Re表示实部,2αRe(o*r)表示物体干涉形成的全息图,其中o和r分别为衍射光和参考光的振幅大小;第三项α(o*r)(o*r)*表示物体之间对干涉的影响;
步骤3.3:α(o*r)(o*r)*<<2αRe(o*r),因而可以忽略。因而全息图可以简化为:
式中F、F-1分别表示二维傅里叶变换和傅里叶逆变换。
步骤3.4:假设CCD的分辨力为Nx×Ny,像素大小为Δx×Δy,被记录的目标空间o(x0,y0)被划分为步长为Δx×Δy的Nx×Ny个采样区间,对上式进行离散化,得
步骤3.5:为了应用压缩感知理论,将目标空间二维矩阵与全息图化为一维向量,定义则上式简化为
I′=T2DWQBO′=RO′
其中,和O'分别为式中和的简写。式中B=F2D是大小为(Nx×Ny)×(Nx×Ny)二维分块对角阵,表示二维离散傅里叶变换;BO'对应式中F[o(mΔx,nΔy,qΔz)];Q为(Nx×Ny)×(Nx×Ny)的对角矩阵,其[(n-1)×Nx+m]为F[r(kΔx-mΔx,lΔy-nΔy)]矩阵第n行m列元素;F[r(kΔx-mΔx,lΔy-nΔy)]表示菲涅尔变换核r(mΔx,nΔy)的傅里叶变换;W=I1,为(Nx×Ny)×(Nx×Ny)单位矩阵;T2D大小为(Nx×Ny)×(Nx×Ny),表示二维傅里叶逆变换;R=T2DWQB为压缩感知理论的稀疏矩阵,大小为(Nx×Ny)×(Nx×Ny),则A=ΦR=ΦT2DWQB。
步骤四中利用IDNST算法重构出两幅全息图的具体步骤如下:
步骤4.1:初始化;初始图像O'0=0,初始两步迭代次数TwIST_iters=0,迭代次数t1=1,最大迭代次数Mt=1000,算法终止条件值tolA=10-3,正则化收缩因子u=0.9,计数器s=1,以及各种临时参数O'1=O'0,O'2=O'0;
步骤4.2:定义目标函数为
根据步骤三得
其中A是传感矩阵,y表示得到的观测数据,I'代表是步骤二获得的待重构的数据全息图;τ代表正则化参数,ΘTV(I')为总变分函数,其公式如下:
和代表对于图像像素值I'的一阶水平和垂直差分操作
步骤4.3:在最优化问题中将式转化为最小化式
ψλ(O')=arg min{X}
并计算O'0对应的目标函数X(O'0),并且赋值pref=X(O'0);
步骤4.4:对O'0进行去噪处理,O'1=Γλ(O'0),其中,
Γλ(O')=Ψλ(O'+ΦT(y-ΦO'))
Ψλ为去噪软收缩函数,在此引入迭代因子t,则
从而
并计算O'k+1对应的目标函数X(O'k+1);
步骤4.5:判断算法是否已经进行过两步迭代;
如果TwIST_iters≠0不成立,则计算O'的目标函数值X(O');之后再判断条件X(O'k)>X(O'k+1)是否成立;
若X(O'k)>X(O'k+1)成立,则s=s×2;并判断条件s>1000是否成立,若成立,则终止程序,最终输出估计值O'k,若不成立,则赋值TwIST_iters=0并转至步骤4.4;
若X(O'k)>X(O'k+1)不成立,则赋值TwIST_iters=TwIST_iters+1,并执行步骤4.6;
如果TwIST_iters≠0成立,则两步估计O'k+1=(1-α)O'k-1+(α-β)O'k+βΓλ(O'k),同时计算O'k+1的目标函数值X(O'),之后再判断X(O')>pref是否成立;若成立则赋值TwIST_iters=0,并转至步骤4.4;若不成立则赋值TwIST_iters=TwIST_iters+1,O'=O'k+1,并执行步骤4.6;
步骤4.6:判断算法终止条件;判断条件迭代次数s>1000是否成立,若成立则终止算法,输出最终估计值O';若不成立则赋值O'k=O'k-1,O'k+1=O'k,并将迭代次数加1,并计算算法终止条件并赋值X(O'k)=X(O'k-1);
步骤4.7:判断C(O'k,O'k-1)≤tolA和t>Mt是否成立:
若都不成立,则收缩正则化参数,τ=τ×u,转至步骤4.4;
只要有一条件成立,则终止算法;最终重构出原始图像
请参阅图1至图4所示,本发明利用PC平台仿真代替干涉仪实验装置获取全息图。硬件为acer AL1916W 3GHz、内存3G,软件为Windows 7旗舰版,MATLAB R2014a版本。仿真中光源采用He-Ne激光器,其波长为632.8nm;原始图像像素为256pixel×256pixel。在IDNST重构算法中,参数设定:在编程过程中,根据经验得从而得α=ρ2+1和β=2α/λ1+λm。将图5和图6的(c)图进行对照可以看出在相同的初始化条件下本文提出的IDNST算法恢复效果更佳。根据图7和图8,可以看出IDNST算法迭代时间简短,信噪比更大并且均方误差更小。
Claims (4)
1.一种基于迭代去噪收缩阈值算法的数字全息重构方法,其特征在于,包括以下步骤:
步骤一:读入原始图像,其振幅为O0(ξ,η),利用此原始图像的物光信息进行菲涅尔衍射得到距离孔径平面为z的观察平面上任意一点P(x,y)的衍射光复振幅O(x,y);
步骤二:由菲涅尔衍射得到的衍射光复振幅O(x,y)与加入了0和π/2相移量后的参照光R(x,y)相干涉到达全息干板上的光波复振幅分别为U1(x,y)和U2(x,y),经过数值处理形成两幅全息图I(1)和I(2);
步骤三:对两幅全息图I(1)和I(2)进行叠加生成一幅相移全息图构造传感矩阵A,使得其中Φ为已知的测量矩阵,R表示稀疏基矩阵,y表示得到的观测数据;
步骤四:利用迭代去噪收缩阈值算法重构出原始图像
步骤四具体包括以下步骤:
步骤4.1:初始化;初始图像O'0=0,初始两步迭代次数TwIST_iters=0,迭代次数t1=1,最大迭代次数Mt=1000,算法终止条件值tolA=10-3,正则化收缩因子u=0.9,计数器s=1,以及各种临时参数O'1=O'0,O'2=O'0;
步骤4.2:定义目标函数为
<mrow>
<mi>X</mi>
<mrow>
<mo>(</mo>
<mi>I</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mo>|</mo>
<mo>|</mo>
<mi>y</mi>
<mo>-</mo>
<mi>&Phi;</mi>
<mi>I</mi>
<mo>|</mo>
<msubsup>
<mo>|</mo>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msub>
<mi>&tau;&Theta;</mi>
<mrow>
<mi>T</mi>
<mi>V</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msup>
<mi>I</mi>
<mo>&prime;</mo>
</msup>
<mo>)</mo>
</mrow>
</mrow>
根据步骤三得
<mrow>
<mi>X</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>I</mi>
<mo>&prime;</mo>
</msup>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mo>|</mo>
<mo>|</mo>
<mi>y</mi>
<mo>-</mo>
<msup>
<mi>AO</mi>
<mo>&prime;</mo>
</msup>
<mo>|</mo>
<msubsup>
<mo>|</mo>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msub>
<mi>&tau;&Theta;</mi>
<mrow>
<mi>T</mi>
<mi>V</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msup>
<mi>I</mi>
<mo>&prime;</mo>
</msup>
<mo>)</mo>
</mrow>
</mrow>
其中A是传感矩阵,y表示得到的观测数据,I'代表是步骤三获得的叠加生成的相移全息图;τ代表正则化参数,ΘTV(I')为总变分函数,其公式如下:
<mrow>
<msub>
<mi>&Theta;</mi>
<mrow>
<mi>T</mi>
<mi>V</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munder>
<mo>&Sigma;</mo>
<mi>i</mi>
</munder>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>&Delta;</mi>
<mi>i</mi>
<mi>h</mi>
</msubsup>
<msup>
<mi>I</mi>
<mo>&prime;</mo>
</msup>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>&Delta;</mi>
<mi>i</mi>
<mi>v</mi>
</msubsup>
<msup>
<mi>I</mi>
<mo>&prime;</mo>
</msup>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mrow>
和代表对于图像像素值I'的一阶水平和垂直差分操作
<mrow>
<msubsup>
<mi>&Delta;</mi>
<mi>i</mi>
<mi>h</mi>
</msubsup>
<msup>
<mi>I</mi>
<mo>&prime;</mo>
</msup>
<mo>=</mo>
<msub>
<msup>
<mi>I</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<msup>
<mi>I</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msubsup>
<mi>&Delta;</mi>
<mi>i</mi>
<mi>v</mi>
</msubsup>
<msup>
<mi>I</mi>
<mo>&prime;</mo>
</msup>
<mo>=</mo>
<msub>
<msup>
<mi>I</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<msup>
<mi>I</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
</msub>
</mrow>
步骤4.3:在最优化问题中将式转化为最小化式
ψλ(O')=argmin{X}
并计算O'0对应的目标函数X(O'0),并且赋值pref=X(O'0);
步骤4.4:对O'0进行去噪处理,O'1=Γλ(O'0),其中,
Γλ(O')=Ψλ(O'+ΦT(y-ΦO'))
Ψλ为去噪软收缩函数,在此引入迭代因子t,则
<mrow>
<msub>
<mi>t</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mn>1</mn>
<mo>+</mo>
<msqrt>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mn>4</mn>
<msubsup>
<mi>t</mi>
<mi>k</mi>
<mn>2</mn>
</msubsup>
</mrow>
</msqrt>
</mrow>
<mn>2</mn>
</mfrac>
</mrow>
从而
<mrow>
<msub>
<msup>
<mi>O</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>=</mo>
<msub>
<msup>
<mi>O</mi>
<mo>&prime;</mo>
</msup>
<mi>k</mi>
</msub>
<mo>+</mo>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>t</mi>
<mi>k</mi>
</msub>
<mo>-</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>t</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<msup>
<mi>O</mi>
<mo>&prime;</mo>
</msup>
<mi>k</mi>
</msub>
<mo>-</mo>
<msub>
<msup>
<mi>O</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
并计算O'k+1对应的目标函数X(O'k+1);
步骤4.5:判断算法是否已经进行过两步迭代;
如果TwIST_iters≠0不成立,则计算O'的目标函数值X(O');之后再判断条件X(O'k)>X(O'k+1)是否成立;
若X(O'k)>X(O'k+1)成立,则s=s×2;并判断条件s>1000是否成立,若成立,则终止程序,最终输出估计值O'k,若不成立,则赋值TwIST_iters=0并转至步骤4.4;
若X(O'k)>X(O'k+1)不成立,则赋值TwIST_iters=TwIST_iters+1,并执行步骤4.6;
如果TwIST_iters≠0成立,则两步估计O'k+1=(1-α)O'k-1+(α-β)O'k+βΓλ(O'k),同时计算O'k+1的目标函数值X(O'),之后再判断X(O')>pref是否成立;若成立则赋值TwIST_iters=0,并转至步骤4.4;若不成立则赋值TwIST_iters=TwIST_iters+1,O'=O'k+1,并执行步骤4.6;
步骤4.6:判断算法终止条件;判断条件迭代次数s>1000是否成立,若成立则终止算法,输出最终估计值O';若不成立则赋值O'k=O'k-1,O'k+1=O'k,并将迭代次数加1,并计算算法终止条件并赋值X(O'k)=X(O'k-1);
步骤4.7:判断C(O'k,O'k-1)≤tolA和t>Mt是否成立:
若都不成立,则收缩正则化参数,τ=τ×u,转至步骤4.4;
只要有一条件成立,则终止算法;最终重构出原始图像
2.根据权利要求1所述的一种基于迭代去噪收缩阈值算法的数字全息重构方法,其特征在于,步骤一中:
<mrow>
<mi>O</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mi>i</mi>
<mi>&lambda;</mi>
<mi>z</mi>
</mrow>
</mfrac>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mi>k</mi>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mi>exp</mi>
<mo>&lsqb;</mo>
<mfrac>
<mrow>
<mi>i</mi>
<mi>&pi;</mi>
</mrow>
<mrow>
<mn>2</mn>
<mi>z</mi>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<msup>
<mi>x</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mi>y</mi>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mi>F</mi>
<mo>{</mo>
<msub>
<mi>O</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>&xi;</mi>
<mo>,</mo>
<mi>&eta;</mi>
<mo>)</mo>
</mrow>
<mi>exp</mi>
<mo>&lsqb;</mo>
<mfrac>
<mrow>
<mi>i</mi>
<mi>k</mi>
</mrow>
<mrow>
<mn>2</mn>
<mi>z</mi>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<msup>
<mi>&xi;</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mi>&eta;</mi>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>}</mo>
</mrow>
其中,λ是光的波长,k是波数。
3.根据权利要求1所述的一种基于迭代去噪收缩阈值算法的数字全息重构方法,其特征在于,步骤二具体包括以下步骤:
步骤2.1:衍射光O(x,y)与参考光R(x,y)相干涉在全息记录平面xy上所叠加的总光场分布U(x,y)为:
U(x,y)=O(x,y)+R(x,y)
步骤2.2:利用全息记录平面xy上的总光场U(x,y)计算形成的全息图光强分布为:
I=U(x,y)·U*(x,y)=(O+R)(O*+R*)
=|O|2+|R|2+O·R*+O*·R
式中I代表全息图I(1)和I(2);上式中U(x,y)为U1(x,y)时,计算获得全息图I(1);U(x,y)为U2(x,y)时,计算获得全息图I(2);
步骤2.1中,在参考光R(x,y)的相移量为0时,获得U1(x,y);在参考光R(x,y)的相移量为π/2时,获得U2(x,y)。
4.根据权利要求1所述的一种基于迭代去噪收缩阈值算法的数字全息重构方法,其特征在于,步骤三中:构造传感矩阵A的步骤为:
步骤3.1:将步骤二中得到的两幅全息图I(1)和I(2)叠加得
<mrow>
<mover>
<mi>I</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>2</mn>
<mi>r</mi>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>I</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>iI</mi>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
步骤3.2:根据步骤二,CCD上的光强分布具体为
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mover>
<mi>I</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>&alpha;</mi>
<mo>|</mo>
<mi>U</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>I</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mi>&alpha;</mi>
<mo>+</mo>
<msub>
<mi>I</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mn>2</mn>
<mi>&alpha;</mi>
<mi>Re</mi>
<mrow>
<mo>(</mo>
<mi>o</mi>
<mo>*</mo>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>&alpha;</mi>
<mrow>
<mo>(</mo>
<mi>o</mi>
<mo>*</mo>
<mi>r</mi>
<mo>)</mo>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<mi>o</mi>
<mo>*</mo>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mo>*</mo>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
式中,α表示CCD量子系数,I0(x,y)表示CCD背景噪音;式中第一项(α+I0)表示全息图的直流分量;第二项Re表示实部,2αRe(o*r)表示物体干涉形成的全息图,其中o和r分别为衍射光和参照光的振幅大小;第三项α(o*r)(o*r)*表示物体之间对干涉的影响;
步骤3.3:α(o*r)(o*r)*2αRe(o*r),将全息图简化为:
<mrow>
<mover>
<mi>I</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>Re</mi>
<mrow>
<mo>(</mo>
<mi>o</mi>
<mo>*</mo>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>Re</mi>
<mo>{</mo>
<msup>
<mi>F</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>&lsqb;</mo>
<mi>F</mi>
<mrow>
<mo>(</mo>
<mi>o</mi>
<mo>)</mo>
</mrow>
<mi>F</mi>
<mrow>
<mo>(</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>}</mo>
</mrow>
式中F、F-1分别表示二维傅里叶变换和傅里叶逆变换;
步骤3.4:CCD的分辨力为Nx×Ny,像素大小为Δx×Δy,被记录的目标空间o(x0,y0)被划分为步长为Δx×Δy的Nx×Ny个采样区间,对上式进行离散化,得
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mover>
<mi>I</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mi>&Delta;</mi>
<mi>x</mi>
<mo>,</mo>
<mi>l</mi>
<mi>&Delta;</mi>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mi>n</mi>
<msub>
<mi>N</mi>
<mi>y</mi>
</msub>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mi>m</mi>
<msub>
<mi>N</mi>
<mi>x</mi>
</msub>
</munderover>
<mi>o</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mi>&Delta;</mi>
<mi>x</mi>
<mo>,</mo>
<mi>n</mi>
<mi>&Delta;</mi>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>*</mo>
<mi>r</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mi>&Delta;</mi>
<mi>x</mi>
<mo>-</mo>
<mi>m</mi>
<mi>&Delta;</mi>
<mi>x</mi>
<mo>,</mo>
<mi>l</mi>
<mi>&Delta;</mi>
<mi>y</mi>
<mo>-</mo>
<mi>n</mi>
<mi>&Delta;</mi>
<mi>y</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<msup>
<mi>F</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>{</mo>
<mi>F</mi>
<mo>&lsqb;</mo>
<mi>o</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mi>&Delta;</mi>
<mi>x</mi>
<mo>,</mo>
<mi>n</mi>
<mi>&Delta;</mi>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mi>F</mi>
<mo>&lsqb;</mo>
<mi>r</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mi>&Delta;</mi>
<mi>x</mi>
<mo>-</mo>
<mi>m</mi>
<mi>&Delta;</mi>
<mi>x</mi>
<mo>,</mo>
<mi>l</mi>
<mi>&Delta;</mi>
<mi>y</mi>
<mo>-</mo>
<mi>n</mi>
<mi>&Delta;</mi>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>}</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
步骤3.5:将目标空间二维矩阵与全息图化为一维向量,定义 则上式简化为
I'=T2DWQBO'=RO'
其中,和O'分别为式中和的简写;B=F2D是大小为(Nx×Ny)×(Nx×Ny)二维分块对角阵,表示二维离散傅里叶变换;BO'对应式中F[o(mΔx,nΔy,qΔz)];Q为(Nx×Ny)×(Nx×Ny)的对角矩阵,其[(n-1)×Nx+m]为F[r(kΔx-mΔx,lΔy-nΔy)]矩阵第n行m列元素;F[r(kΔx-mΔx,lΔy-nΔy)]表示菲涅尔变换核r(mΔx,nΔy)的傅里叶变换;W=I1,为(Nx×Ny)×(Nx×Ny)单位矩阵;T2D大小为(Nx×Ny)×(Nx×Ny),表示二维傅里叶逆变换;R=T2DWQB为压缩感知理论的稀疏矩阵,大小为(Nx×Ny)×(Nx×Ny),则传感矩阵A=ΦR=ΦT2DWQB。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610021891.8A CN105629696B (zh) | 2016-01-13 | 2016-01-13 | 一种基于迭代去噪收缩阈值算法的数字全息重构方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610021891.8A CN105629696B (zh) | 2016-01-13 | 2016-01-13 | 一种基于迭代去噪收缩阈值算法的数字全息重构方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105629696A CN105629696A (zh) | 2016-06-01 |
CN105629696B true CN105629696B (zh) | 2018-04-17 |
Family
ID=56044769
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610021891.8A Expired - Fee Related CN105629696B (zh) | 2016-01-13 | 2016-01-13 | 一种基于迭代去噪收缩阈值算法的数字全息重构方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105629696B (zh) |
Families Citing this family (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR102669967B1 (ko) * | 2016-11-24 | 2024-05-29 | 삼성전자주식회사 | 홀로그램 생성 방법 및 홀로그램 생성 장치 |
CN108765575A (zh) * | 2018-02-24 | 2018-11-06 | 石化盈科信息技术有限责任公司 | 一种基于ar的工业设备图鉴展示方法和系统 |
CN108646538B (zh) * | 2018-04-24 | 2020-11-20 | 安徽大学 | 一种单次曝光复振幅物体全息重建方法、设备及系统 |
CN112236725B (zh) * | 2018-08-23 | 2022-06-14 | 杜尔利塔斯有限公司 | 全息图计算的方法 |
CN109270492B (zh) * | 2018-09-27 | 2022-11-25 | 重庆大学 | 一种用于大全息距离的正则化参数选取方法 |
CN109581849B (zh) * | 2019-01-04 | 2020-10-16 | 中国工程物理研究院激光聚变研究中心 | 一种同轴全息重建方法及系统 |
CN111273533B (zh) * | 2019-09-27 | 2021-08-06 | 中国工程物理研究院激光聚变研究中心 | 一种同轴数字全息自动聚焦方法及系统 |
CN111897197B (zh) * | 2020-08-18 | 2021-11-16 | 四川大学 | 基于双相位编码的傅里叶相位全息图生成方法 |
CN112180707B (zh) * | 2020-09-28 | 2021-11-02 | 四川大学 | 基于球面自衍射模型的球面纯相位全息图生成方法 |
CN112666129B (zh) * | 2020-12-14 | 2023-03-31 | 西安邮电大学 | 一种考虑折射率差异的三波长相干衍射成像方法 |
CN113917819B (zh) * | 2021-10-13 | 2023-03-21 | 中国工程物理研究院激光聚变研究中心 | 基于菲涅尔掩膜的非相干三维全息分层重构方法 |
CN114237000B (zh) * | 2021-12-15 | 2023-05-23 | 中国工程物理研究院激光聚变研究中心 | 一种离轴数字全息优化重建方法及系统 |
CN114442313B (zh) * | 2021-12-16 | 2022-11-04 | 南京大学 | 一种基于多次迭代的光学超晶格优化设计方法 |
CN114413794B (zh) * | 2022-01-29 | 2023-09-22 | 中国工程物理研究院激光聚变研究中心 | 大口径kdp晶体最佳相位匹配角测量系统及其测量方法 |
CN115100300B (zh) * | 2022-05-26 | 2024-07-16 | 北京大学 | 一种基于注意力量子启发式神经网络的计算机全息图压缩方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102597891A (zh) * | 2009-10-05 | 2012-07-18 | 原子能与替代能源委员会 | 叠加可见图像和合成全息图 |
EP2024792B1 (en) * | 2006-05-11 | 2015-04-22 | Cambridge Enterprise Limited | Phase retrieval and phase hologram synthesis |
CN104952034A (zh) * | 2014-03-27 | 2015-09-30 | 香港城市大学 | 复杂全息图向相位全息图的转换 |
US9208533B2 (en) * | 2010-09-08 | 2015-12-08 | Commissariat A L'energie Atomique Et Aux Energies Alternatives | Method for concealing a synthetic hologram in a binary image |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2006090320A1 (en) * | 2005-02-23 | 2006-08-31 | Lyncee Tec S.A. | Wave front sensing method and apparatus |
-
2016
- 2016-01-13 CN CN201610021891.8A patent/CN105629696B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2024792B1 (en) * | 2006-05-11 | 2015-04-22 | Cambridge Enterprise Limited | Phase retrieval and phase hologram synthesis |
CN102597891A (zh) * | 2009-10-05 | 2012-07-18 | 原子能与替代能源委员会 | 叠加可见图像和合成全息图 |
US9208533B2 (en) * | 2010-09-08 | 2015-12-08 | Commissariat A L'energie Atomique Et Aux Energies Alternatives | Method for concealing a synthetic hologram in a binary image |
CN104952034A (zh) * | 2014-03-27 | 2015-09-30 | 香港城市大学 | 复杂全息图向相位全息图的转换 |
Also Published As
Publication number | Publication date |
---|---|
CN105629696A (zh) | 2016-06-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105629696B (zh) | 一种基于迭代去噪收缩阈值算法的数字全息重构方法 | |
Huang et al. | Bicubic interpolation and extrapolation iteration method for high resolution digital holographic reconstruction | |
CN108416723B (zh) | 基于全变分正则化和变量分裂的无透镜成像快速重构方法 | |
CN110675326B (zh) | 基于U-Net网络的计算鬼成像重建恢复的方法 | |
Shi et al. | Deep prior-based sparse representation model for diffraction imaging: A plug-and-play method | |
Yu et al. | PDNet: A lightweight deep convolutional neural network for InSAR phase denoising | |
CN105184295B (zh) | 一种基于小波变换与连通域的全息扫描空间距离提取方法 | |
CN118071866B (zh) | 一种稀疏数字全息图像重建方法 | |
CN108646538B (zh) | 一种单次曝光复振幅物体全息重建方法、设备及系统 | |
Xia et al. | Comparative analysis for combination of unwrapping and de-noising of phase data with high speckle decorrelation noise | |
Yang et al. | Dynamic proximal unrolling network for compressive imaging | |
Madsen et al. | On-axis digital holographic microscopy: Current trends and algorithms | |
Katkovnik et al. | Sparse approximations of phase and amplitude for wave field reconstruction from noisy data | |
Lu et al. | Complex-valued speckle effect and its suppression for high quality of phase unwrapping reconstruction in coherent digital holographic microscopy | |
Katkovnik et al. | Multiwavelength surface contouring from phase-coded noisy diffraction patterns: wavelength-division optical setup | |
Liang et al. | Multi-fidelity and learning-regularization for single image super resolution | |
Luo et al. | Fast transmission of computer-generated hologram with compressed sensing and quantum-inspired neural network | |
JP2007502445A (ja) | デジタル・ホログラフィにおける画像の再構成での空間解像度を修正する方法 | |
CN111397733A (zh) | 一种单/多帧快照式光谱成像方法、系统及介质 | |
Gemechu | Sparse Regularization Based on Orthogonal Tensor Dictionary Learning for Inverse Problems | |
Xu et al. | Interferogram blind denoising using deep residual learning for phase-shifting interferometry | |
Bo | Deep learning approach for computer-generated holography | |
Reddy et al. | Robust off-axis digital holographic imaging system in compressive sensing framework | |
CN113448232B (zh) | 三维分层目标压缩全息术的测量矩阵降维方法 | |
Migukin et al. | Advanced multi-plane phase retrieval using graphic processing unit: augmented Lagrangian technique with sparse regularization |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180417 Termination date: 20210113 |