CN108873323B - 一种实现边缘增强成像的方法及系统 - Google Patents

一种实现边缘增强成像的方法及系统 Download PDF

Info

Publication number
CN108873323B
CN108873323B CN201810744412.4A CN201810744412A CN108873323B CN 108873323 B CN108873323 B CN 108873323B CN 201810744412 A CN201810744412 A CN 201810744412A CN 108873323 B CN108873323 B CN 108873323B
Authority
CN
China
Prior art keywords
imaging
filter
fourier transform
fourier
gaussian
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
Application number
CN201810744412.4A
Other languages
English (en)
Other versions
CN108873323A (zh
Inventor
姚海南
王飞
蔡阳健
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Suzhou University
Original Assignee
Suzhou University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Suzhou University filed Critical Suzhou University
Priority to CN201810744412.4A priority Critical patent/CN108873323B/zh
Publication of CN108873323A publication Critical patent/CN108873323A/zh
Application granted granted Critical
Publication of CN108873323B publication Critical patent/CN108873323B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G02OPTICS
    • G02BOPTICAL ELEMENTS, SYSTEMS OR APPARATUS
    • G02B27/00Optical systems or apparatus not provided for by any of the groups G02B1/00 - G02B26/00, G02B30/00
    • GPHYSICS
    • G02OPTICS
    • G02BOPTICAL ELEMENTS, SYSTEMS OR APPARATUS
    • G02B27/00Optical systems or apparatus not provided for by any of the groups G02B1/00 - G02B26/00, G02B30/00
    • G02B27/09Beam shaping, e.g. changing the cross-sectional area, not otherwise provided for
    • G02B27/0927Systems for changing the beam intensity distribution, e.g. Gaussian to top-hat
    • GPHYSICS
    • G02OPTICS
    • G02BOPTICAL ELEMENTS, SYSTEMS OR APPARATUS
    • G02B27/00Optical systems or apparatus not provided for by any of the groups G02B1/00 - G02B26/00, G02B30/00
    • G02B27/30Collimators

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Optics & Photonics (AREA)
  • Image Processing (AREA)

Abstract

本发明公开了一种实现边缘增强成像的方法及系统,通过将扩束和准直后的高斯光束投射到物体上,接着,对透过物体的光束进行第一傅里叶变换,对经过第一傅里叶变换后的光束进行调制,并对调制后的光束进行第二傅里叶变换后成像,其中,对光束进行调制的滤波函数的得到过程如下:先结合振幅和相位对边缘增强成像的调制作用,通过将四个高斯函数叠加得到点扩散函数,满足其主瓣周围没有多余旁瓣,并通过利用傅里叶变换,对点扩散函数进行反向计算,从而得到的,可以在第二傅里叶透镜的傅里叶平面实现分辨率高而且各向同性的边缘增强成像。

Description

一种实现边缘增强成像的方法及系统
技术领域
本发明涉及边缘增强成像技术领域,特别涉及一种实现边缘增强成像的方法及系统。
背景技术
边缘增强成像技术在成像领域有着很大的应用,1942年,Zernike首次实现了相位对比成像,从那之后,大量的工作都用来研究这种技术,Marr和Torre等人在理论上对边缘增强成像有很大的贡献。目前常用来实现边缘增强成像的技术主要有①数值空间微分成像②微分干涉对比度成像③希尔伯特转换滤波成像。
2018年,朱腾峰利用表面等离激元结构在实验中实现了空间微分的边缘增强成像,但这种数字成像技术在物体没有明显的特征的情况下不能很好的实现相位对。关于微分干涉对比度成像在实现起来就相对简单一点,可以利用空间光调制器来简化光路进而完成成像操作,不过这种技术的成像结果是各向异性的。而希尔伯特转换滤波成像则是利用里空间滤波的思想,在4f系统的频谱面上加上滤波片来实现边缘增强成像,想要实现希尔伯特转换,最实用的方法是利用螺旋位相板。
在光学微波领域,螺旋位相板被大量的使用来重构生物标本边缘的振幅和相位信息。随着技术的发展,方向可选的边缘增强成像可以通过矢量光学滤波片、分数阶涡旋滤波片和相移涡旋滤波片来实现。同时涡旋透镜在边缘增强成像中也有很大的应用价值。但是我们发现由于传统的涡旋位相板的点扩散函数的主瓣两侧存在大量的多余旁瓣,这种现象会让成像结果在边缘处出现衍射噪声,使结果不均匀,而且这种效应会随着拓扑和的增加而变得严重。为了解决这个问题从而提高成像质量,拉盖尔高斯滤波片、贝塞尔型滤波片、艾利螺旋位相滤波片被设计出来用于抑制多余旁瓣,从而得到均匀的、分辨率高的边缘成像结果。目前边缘增强成像技术在红外照明、生物成像、天文观测、指纹鉴定和远距离遥感等领域都有重大的应用前景。
发明内容
针对现有技术的不足,本发明目的之一在于提供一种实现边缘增强成像的滤波片。采用如下技术方案:
一种实现边缘增强成像的方法,其包括:
产生高斯光束;
对所述高斯光束进行扩束和准直,并将扩束和准直后的高斯光束投射到物体上;
对透过物体的光束进行第一傅里叶变换;
对经过第一傅里叶变换后的光束进行调制;
对调制后的光束进行第二傅里叶变换后成像;
其中,对光束进行调制的滤波片的滤波函数
Figure GDA0002222638530000021
为:
Figure GDA0002222638530000022
Figure GDA0002222638530000023
(其中,
Figure GDA0002222638530000024
代表频谱面的径向坐标,FT代表傅里叶变换,代表频谱面上半径为R的圆形孔,该圆形孔为定义域,定义域内值为1,定义域外值为0,h(x,y)为滤波函数
Figure GDA0002222638530000026
的点扩散函数,其由四个高斯函数叠加得到,(x,y)代表成像面的坐标,ω0是高斯束腰,d0是任意常数)。
作为本发明的进一步改进,对所述高斯光束进行扩束和准直,具体包括:通过扩束镜对高斯光束进行扩束和准直。
作为本发明的进一步改进,所述第一傅里叶透镜和第二傅里叶透镜的焦距相等。
本发明目的之二在于提供一种实现边缘增强成像的系统。采用如下技术方案:
一种实现边缘增强成像的系统,其包括:
光源,用于产生高斯光束;
扩束准直装置,用于对所述高斯光束进行扩束和准直,并将扩束和准直后的高斯光束投射到物体上;
第一傅里叶透镜,用于对透过物体的光束进行第一傅里叶变换;
滤波片,用于对经过第一傅里叶变换后的光束进行调制;
第二傅里叶透镜,用于对调制后的光束进行第二傅里叶变换,在傅里叶平面处成像;
其中,所述滤波片的滤波函数
Figure GDA0002222638530000031
为:
Figure GDA0002222638530000032
Figure GDA0002222638530000033
(其中,代表频谱面的径向坐标,FT代表傅里叶变换,
Figure GDA0002222638530000035
代表频谱面上半径为R的圆形孔,该圆形孔为定义域,定义域内值为1,定义域外值为0,h(x,y)为滤波函数
Figure GDA0002222638530000036
的点扩散函数,其由四个高斯函数叠加得到,(x,y)代表成像面的坐标,ω0是高斯束腰,d0是任意常数)。
作为本发明的进一步改进,所述第一傅里叶透镜和第二傅里叶透镜的焦距相等。
作为本发明的进一步改进,所述扩束准直装置为扩束镜。
作为本发明的进一步改进,还包括光束分析仪,所述光束分析仪设于所述第二傅里叶透镜的傅里叶平面。
本发明的有益效果:
本发明公开了一种实现边缘增强成像的方法及系统,通过将扩束和准直后的高斯光束投射到物体上,接着,对透过物体的光束进行第一傅里叶变换,对经过第一傅里叶变换后的光束进行调制,并对调制后的光束进行第二傅里叶变换后成像,其中,对光束进行调制的滤波函数的得到过程如下:先结合振幅和相位对边缘增强成像的调制作用,通过将四个高斯函数叠加得到点扩散函数,满足其主瓣周围没有多余旁瓣,并通过利用傅里叶变换,对点扩散函数进行反向计算,从而得到滤波片的滤波函数,可以在第二傅里叶透镜的傅里叶平面实现分辨率高而且各向同性的边缘增强成像。
上述说明仅是本发明技术方案的概述,为了能够更清楚了解本发明的技术手段,而可依照说明书的内容予以实施,并且为了让本发明的上述和其他目的、特征和优点能够更明显易懂,以下特举较佳实施例,并配合附图,详细说明如下。
附图说明
图1是本发明实施例中实现边缘增强成像的系统的示意图;
图2是本发明实施例中4f系统的示意图;
图3(a)和3(b)分别是本发明实施例中滤波函数
Figure GDA0002222638530000041
(涡旋函数)的点扩散函数实部的二维平面分布和径向值在径向方向上的截面分布图;图3(c)和3(d)分别是本发明实施例中滤波函数
Figure GDA0002222638530000042
(涡旋函数)的点扩散函数虚部的二维平面分布和径向值在径向方向上的截面分布图;
图4(a)和4(b)分别是本发明实施例中滤波函数
Figure GDA0002222638530000043
(Bessel函数)的点扩散函数虚部的二维平面分布和径向值在径向方向上的截面分布图;
图5(a)和5(b)分别是本发明实施例中滤波函数
Figure GDA0002222638530000044
的点扩散函数实部的二维平面分布和径向值在径向方向上的截面分布图;图5(c)和5(d)分别是本发明实施例中滤波函数
Figure GDA0002222638530000045
的点扩散函数虚部的二维平面分布和径向值在径向方向上的截面分布图;
图6是本发明实施例中振幅物体的理论和实验对比图;
图7是本发明实施例中相位物体的理论和实验对比图。
附图说明:1、光源;2、扩束镜;3、物体;4、第一傅里叶透镜;5、滤波片;6、第二傅里叶透镜;7、光束分析仪。
具体实施方式
下面结合附图和具体实施例对本发明作进一步说明,以使本领域的技术人员可以更好地理解本发明并能予以实施,但所举实施例不作为对本发明的限定。
本发明实施例中的实现边缘增强成像的方法,包括以下步骤:
S11、产生高斯光束;
S12、对所述高斯光束进行扩束和准直,并将扩束和准直后的高斯光束投射到物体上;
具体的,通过扩束镜对高斯光束进行扩束和准直。
S13、对透过物体的光束进行第一傅里叶变换;
具体的,通过第一傅里叶透镜对透过物体的光束进行第一傅里叶变换。
S14、对经过第一傅里叶变换后的光束进行调制;
具体的,利用空间光调制器加载的滤波片对光束进行调制,其中,该滤波片的滤波函数
Figure GDA0002222638530000051
为:
Figure GDA0002222638530000052
Figure GDA0002222638530000053
(其中,
Figure GDA0002222638530000054
代表频谱面的径向坐标,FT代表傅里叶变换,
Figure GDA0002222638530000055
代表频谱面上半径为R的圆形孔,该圆形孔为定义域,定义域内值为1,定义域外值为0,h(x,y)为滤波函数
Figure GDA0002222638530000056
的点扩散函数,其由四个高斯函数叠加得到,(x,y)代表成像面的坐标,ω0是高斯束腰,d0是任意常数)。
S15、对调制后的光束进行第二傅里叶变换后成像;
具体的,通过第二傅里叶透镜对对调制后的光束进行第二傅里叶变换,在第二傅里叶透镜的后焦平面成像,即可实现边缘增强成像。
其中,第一傅里叶透镜和第二傅里叶透镜的焦距相等,物体设于第一傅里叶透镜的前焦平面,滤波片设于第一傅里叶透镜的后焦平面上,同时也是第二傅里叶透镜的前焦平面。
如图1所示,为本发明实施例中实现边缘增强成像的系统,该系统包括:光源1、扩束镜2、物体3、第一傅里叶透镜4、滤波片5、第二傅里叶透镜6、光束分析仪7,第一傅里叶透镜4和第二傅里叶透镜6的焦距相等,物体3设于第一傅里叶透镜4的前焦平面,滤波片5设于第一傅里叶透镜4的后焦平面上,同时也是第二傅里叶透镜6的前焦平面,光束分析仪7设于第二傅里叶透镜6的后焦平面。
光源1用于产生高斯光束,扩束镜2用于对所述高斯光束进行扩束和准直,并将扩束和准直后的高斯光束投射到物体3上,第一傅里叶透镜4用于对透过物体3的光束进行第一傅里叶变换,滤波片5用于对经过第一傅里叶变换后的光束进行调制,第二傅里叶透镜6用于对调制后的光束进行第二傅里叶变换,在傅里叶平面(第二傅里叶透镜6的后焦平面)处成像,光束分析仪7用于拍摄光强信息。
其中,滤波片5的滤波函数
Figure GDA0002222638530000061
为:
Figure GDA0002222638530000062
Figure GDA0002222638530000063
(其中,
Figure GDA0002222638530000064
代表频谱面的径向坐标,FT代表傅里叶变换,
Figure GDA0002222638530000065
代表频谱面上半径为R的圆形孔,该圆形孔为定义域,定义域内值为1,定义域外值为0,h(x,y)为滤波函数
Figure GDA0002222638530000066
的点扩散函数,其由四个高斯函数叠加得到,(x,y)代表成像面的坐标,ω0是高斯束腰,d0是任意常数)。
如图2所示,为本实施例中4f系统的示意图,4f系统为本发明实现边缘增强成像的方法及系统理论基础。L1和L2是傅里叶薄透镜,它们的焦距都是f。在图2中,(x0,y0),(u,v)和(x,y)分别代表了入射面、频谱面和成像面的笛卡尔坐标。这里假设入射面上物体的复振幅是g(x0,y0),在频谱面上的滤波片的函数是H(u,v),成像面上的出射光场的复振幅为根据4f系统的光路计算,我们可以将出射场描述成如下形式:
Figure GDA0002222638530000072
其中,h(x,y)是滤波函数H(u,v)的傅里叶变换,符号
Figure GDA0002222638530000073
代表卷积算符)。根据卷积的理论,可以得到:如果想要实现边缘加强成像,滤波函数需要两个主要的性质,一个是能够使物体内部变暗,另一个是可以让物体的边缘变亮。一般情况下,为了满足第一个性质,滤波函数需要满足下列条件:
Figure GDA0002222638530000074
该公式中省去了一些常数项,从中可以得到只要滤波函数的中心等于零就可以让物体的中间变暗。而另一个性质就比较难以满足,因此我们首先单独分析滤波函数的相位和振幅对成像结果的影响,为此,我们找到了一个纯相位的滤波函数(涡旋函数)和一个纯振幅的滤波函数(Bessel函数),他们的表达式分别为:
Figure GDA0002222638530000075
Figure GDA0002222638530000076
其中,
Figure GDA0002222638530000077
代表频谱面的径向坐标,
Figure GDA0002222638530000078
代表一个半径为R的圆孔孔径,Jl表示lth阶第一类Bessel函数,kr是模式常数。通过傅里叶变换,这两个滤波函数的点扩散函数可以分别近似的表示成如下形式:
Figure GDA0002222638530000081
Figure GDA0002222638530000082
其中,(r,θ)表示成像面的径向坐标,λ表示波长,J0代表0th阶第一类Bessel函数。
为了更直观的看出滤波函数对成像结果的影响,通过数值模拟画出滤波函数的点扩散函数分布图。
如图3所示,图3(a)和3(b)分别是滤波函数(涡旋函数)的点扩散函数实部的二维平面分布和径向值在径向方向上的截面分布图;图3(c)和3(d)分别是滤波函数
Figure GDA0002222638530000084
(涡旋函数)的点扩散函数虚部的二维平面分布和径向值在径向方向上的截面分布图;其中R=700mm,f=400mm,PSF表示点扩散函数。
如图4所示,图4(a)和4(b)分别是滤波函数
Figure GDA0002222638530000085
(Bessel函数)的点扩散函数虚部的二维平面分布和径向值在径向方向上的截面分布图;其中,l=1,kr=1.1mm-1
通过图3和图4,再结合卷积的成像理论,我们发现在滤波函数
Figure GDA0002222638530000086
(涡旋函数)的点扩散函数的实部和虚部中都存在一个主极大瓣和一个主极小瓣,在主极大瓣和主极小瓣的周围会出现很多多余的旁瓣,所以我们可以预测利用滤波函数
Figure GDA0002222638530000087
(涡旋函数)进行的边缘加强成像可以在所有方向上实现,但是在成像的边缘会出现一些污点,使得边缘成像是不均匀的。而对于滤波函数
Figure GDA0002222638530000088
(Bessel函数)来说,它的点扩散函数只有虚部,其存在两个主极大瓣和一个主极小瓣,而且两个主极大瓣和一个主极小瓣差别明显,但在主极大瓣和主极小瓣的旁边只存在一些很微小的旁瓣,衍射噪声被很好地抑制住了,所以对于滤波函数
Figure GDA0002222638530000089
(Bessel函数)来说,它的边缘成像也可以在全部方向上实现,成像结果的质量比较好,而且是各向同性的,但是在边缘处会出现两个边缘图像。
通过前面的分析,我们发现纯相位的滤波函数和纯振幅的滤波函数都可以实现边缘增强成像,但是在边缘增强成像中都存在一些缺陷,这说明滤波函数的振幅和相位对边缘增强成像都有一定的影响。
因此为了得到更好的成像质量,消除纯相位的滤波函数和纯振幅的滤波函数的成像缺陷,结合相位和振幅对边缘加强成像的调制效应,通过计算滤波函数的点扩散函数,确保点扩散函数的实部和虚部满足只存在一个主极大瓣和一个主极小瓣,而且在主极大瓣和主极小瓣周围的多余旁瓣可以被全部抑制,那么该点扩散函数对应的滤波函数即能够实现分辨率高而且各向同性的边缘增强成像,实现同时对振幅和相位进行调制。本发明实施例中的实现边缘增强成像提出滤波片的设计方法包括以下步骤:
步骤21、将四个高斯函数叠加得到滤波函数的点扩散函数;具体的,点扩散函数h(x,y)为:
Figure GDA0002222638530000091
(其中,(x,y)代表成像面的坐标,ω0是高斯束腰,d0是任意常数);
由于高斯函数是最简单的滤波函数,一个高斯函数只有一个极值,两个高斯光束叠加会出现一个极大值和一个极小值,因此,通过将两个纯实数的高斯函数与两个纯虚部的高斯函数叠加可以保证点扩散函数的实部和虚部都只存在一个主极大瓣和主极小瓣。
步骤22、利用傅里叶变换,对上述点扩散函数进行反向计算,得到滤波函数。
该滤波函数为:
(其中,代表频谱面的径向坐标,FT代表傅里叶变换,
Figure GDA0002222638530000102
代表频谱面上半径为R的圆形孔,该圆形孔为定义域,定义域内值为1,定义域外值为0)。
由于频谱面上圆形孔的影响,滤波函数
Figure GDA0002222638530000104
实际的点扩散函数在径向坐标下可以表示成:
如图5所示,图5(a)和5(b)分别是滤波函数
Figure GDA0002222638530000106
的点扩散函数实部的二维平面分布和径向值在径向方向上的截面分布图;图5(c)和5(d)分别是滤波函数
Figure GDA0002222638530000107
的点扩散函数虚部的二维平面分布和径向值在径向方向上的截面分布图;其中,ω0=d0=26mm。
通过图5,可以发现滤波函数
Figure GDA0002222638530000108
点扩散函数的实部和虚部中多余的旁瓣(衍射噪声)已经完全被抑制了,所以对于滤波函数
Figure GDA0002222638530000109
来说,它可以在4f成像系统中实现各个方向的边缘加强成像,而且是各向同性的,成像均匀分布,它的成像质量要比滤波函数
Figure GDA00022226385300001010
(涡旋函数)和滤波函数
Figure GDA00022226385300001011
(Bessel函数)要好很多,非常好的消除了成像中存在的缺陷,可以得到一个效果更好、分辨率更高的边缘像。
如图6所示,为本发明实施例中振幅物体的理论和实验对比图。基于上述实现边缘增强成像的方法及系统。其中,振幅物体是一个简单的圆形孔(半径是7mm),第一行是理论结果,第二行是实验结果,(a)和(e)是物体的照片,(b)和(f)是经过滤波片(滤波函数为
Figure GDA00022226385300001012
)后的像,(c)和(g)是经过滤波片(滤波函数为
Figure GDA00022226385300001013
)后的像,(d)和(h)是经过滤波片(滤波函数为)后的像。
如图7所示,为本发明实施例中相位物体的理论和实验对比图。基于上述实现边缘增强成像的方法及系统。其中,相位物体是一个熊猫(相位变化为0-π),第一行是理论结果,第二行是实验结果,(a)和(e)是物体的照片,(b)和(f)是经过滤波片(滤波函数为
Figure GDA0002222638530000111
)后的像,(c)和(g)是经过滤波片(滤波函数为)后的像,(d)和(h)是经过滤波片(滤波函数为
Figure GDA0002222638530000113
)后的像。
从图6和7可以发现,在实施例中使用滤波函数为
Figure GDA0002222638530000114
的滤波片可以同时对物体的振幅和相位进行调制,使用滤波函数为
Figure GDA0002222638530000115
的滤波片在本实施例中可以实现分辨率高而且各向同性的边缘增强成像,其相对传统的成像方法,在成像质量上有了很大的提高,消除了成像缺陷。
本发明的有益效果:
本发明公开了一种实现边缘增强成像的方法及系统,通过将扩束和准直后的高斯光束投射到物体上,接着,对透过物体的光束进行第一傅里叶变换,对经过第一傅里叶变换后的光束进行调制,并对调制后的光束进行第二傅里叶变换后成像,其中,对光束进行调制的滤波函数的得到过程如下:先结合振幅和相位对边缘增强成像的调制作用,通过将四个高斯函数叠加得到点扩散函数,满足其主瓣周围没有多余旁瓣,并通过利用傅里叶变换,对点扩散函数进行反向计算,从而得到滤波片的滤波函数,可以在第二傅里叶透镜的傅里叶平面实现分辨率高而且各向同性的边缘增强成像。
以上实施例仅是为充分说明本发明而所举的较佳的实施例,本发明的保护范围不限于此。本技术领域的技术人员在本发明基础上所作的等同替代或变换,均在本发明的保护范围之内。本发明的保护范围以权利要求书为准。

Claims (6)

1.一种实现边缘增强成像的方法,其特征在于,包括:
产生高斯光束;
对所述高斯光束进行扩束和准直,并将扩束和准直后的高斯光束投射到物体上;
对透过物体的光束进行第一傅里叶变换;
对经过第一傅里叶变换后的光束进行调制;
对调制后的光束进行第二傅里叶变换后成像;
其中,对光束进行调制的滤波片的滤波函数
Figure FDA0002222638520000011
为:
Figure FDA0002222638520000012
Figure FDA0002222638520000013
其中,代表频谱面的径向坐标,FT代表傅里叶变换,
Figure FDA0002222638520000015
代表频谱面上半径为R的圆形孔,该圆形孔为定义域,定义域内值为1,定义域外值为0,h(x,y)为滤波函数的点扩散函数,其由四个高斯函数叠加得到,(x,y)代表成像面的坐标,ω0是高斯束腰,d0是任意常数。
2.如权利要求1所述的实现边缘增强成像的方法,其特征在于,对所述高斯光束进行扩束和准直,具体包括:通过扩束镜对高斯光束进行扩束和准直。
3.一种实现边缘增强成像的系统,其特征在于,包括:
光源,用于产生高斯光束;
扩束准直装置,用于对所述高斯光束进行扩束和准直,并将扩束和准直后的高斯光束投射到物体上;
第一傅里叶透镜,用于对透过物体的光束进行第一傅里叶变换;
滤波片,用于对经过第一傅里叶变换后的光束进行调制;
第二傅里叶透镜,用于对调制后的光束进行第二傅里叶变换,在傅里叶平面处成像;
其中,所述滤波片的滤波函数为:
Figure FDA0002222638520000022
Figure FDA0002222638520000023
其中,
Figure FDA0002222638520000024
代表频谱面的径向坐标,FT代表傅里叶变换,
Figure FDA0002222638520000025
代表频谱面上半径为R的圆形孔,该圆形孔为定义域,定义域内值为1,定义域外值为0,h(x,y)为滤波函数
Figure FDA0002222638520000026
的点扩散函数,其由四个高斯函数叠加得到,(x,y)代表成像面的坐标,ω0是高斯束腰,d0是任意常数。
4.如权利要求3所述的实现边缘增强成像的系统,其特征在于,所述第一傅里叶透镜和第二傅里叶透镜的焦距相等。
5.如权利要求3所述的实现边缘增强成像的系统,其特征在于,所述扩束准直装置为扩束镜。
6.如权利要求3所述的实现边缘增强成像的系统,其特征在于,还包括光束分析仪,所述光束分析仪设于所述第二傅里叶透镜的傅里叶平面。
CN201810744412.4A 2018-07-09 2018-07-09 一种实现边缘增强成像的方法及系统 Active CN108873323B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810744412.4A CN108873323B (zh) 2018-07-09 2018-07-09 一种实现边缘增强成像的方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810744412.4A CN108873323B (zh) 2018-07-09 2018-07-09 一种实现边缘增强成像的方法及系统

Publications (2)

Publication Number Publication Date
CN108873323A CN108873323A (zh) 2018-11-23
CN108873323B true CN108873323B (zh) 2020-01-10

Family

ID=64300077

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810744412.4A Active CN108873323B (zh) 2018-07-09 2018-07-09 一种实现边缘增强成像的方法及系统

Country Status (1)

Country Link
CN (1) CN108873323B (zh)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109683461B (zh) * 2019-01-24 2020-11-10 杭州光粒科技有限公司 基于光场渲染的全息图生成方法、系统、存储介质及近眼ar全息三维显示系统
CN109828381B (zh) * 2019-04-01 2021-01-22 北京航空航天大学 一种基于光擦写频域滤波板的边缘增强成像方法
CN111179181B (zh) * 2019-11-15 2023-05-09 南京大学 一种基于拉盖尔高斯模式的图像边缘增强和复制装置
CN114429549A (zh) * 2020-10-28 2022-05-03 北京小米移动软件有限公司 图像边缘提取方法、图像边缘提取装置及存储介质
CN112304246B (zh) * 2020-12-21 2021-03-23 苏州大学 基于空间相干结构调控的光学成像系统及成像方法
CN113362256B (zh) * 2021-07-02 2023-02-17 中国船舶重工集团公司第七二四研究所 一种基于分数傅里叶变换的高分辨图像边缘增强方法
CN113487637B (zh) * 2021-07-05 2024-06-25 南京邮电大学 基于叠加螺旋相位滤波器的多方向边缘检测方法
CN113885196A (zh) * 2021-09-17 2022-01-04 常州北邮新一代信息技术研究院有限公司 基于光学成像系统的图像边缘增强方法、装置
CN113804625A (zh) * 2021-09-30 2021-12-17 常州北邮新一代信息技术研究院有限公司 自动跟踪细胞成像方法及系统
CN114217454B (zh) * 2022-02-22 2022-06-10 华中科技大学 基于衍射光学元件的空间频谱调制器件的设计与实现方法
CN117274287B (zh) * 2023-08-31 2024-06-18 哈尔滨理工大学 一种基于无干涉编码孔径相关全息术的边缘检测方法

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0567558A (ja) * 1991-09-06 1993-03-19 Nikon Corp 露光方法
CN101794024A (zh) * 2010-04-02 2010-08-04 上海理工大学 一种产生柱矢量光束的装置及方法
CN203232201U (zh) * 2013-01-18 2013-10-09 苏州大学 一种部分相干艾里光束的产生装置
GB2536276B (en) * 2015-03-12 2019-11-13 Powerphotonic Ltd Field mapper
CN105225243B (zh) * 2015-10-15 2017-10-10 徐德明 一种可抗噪声的图像边缘检测方法
CN105182547B (zh) * 2015-10-19 2017-07-11 山东师范大学 一种基于双折射偏振分束器产生矢量光束的方法及装置
CN105607267B (zh) * 2016-03-07 2017-11-17 东南大学 一种生成无衍射针型光场的装置
CN106444048B (zh) * 2016-09-06 2018-09-11 苏州大学 一种高斯谢尔模光束的整形方法及整形系统
CN106644105B (zh) * 2016-10-26 2019-04-30 深圳大学 基于双螺旋点扩散函数的波前传感器、探测方法及系统
CN106441084B (zh) * 2016-11-21 2019-02-01 深圳大学 基于微全息阵列的波前传感器、波前探测方法及系统
CN107621701B (zh) * 2017-09-07 2023-08-25 苏州大学 产生双指数贝塞尔高斯光束的方法及系统
CN108254917A (zh) * 2018-01-29 2018-07-06 哈尔滨工业大学 一种基于分数阶光子轨道角动量的边缘增强成像装置及方法

Also Published As

Publication number Publication date
CN108873323A (zh) 2018-11-23

Similar Documents

Publication Publication Date Title
CN108873323B (zh) 一种实现边缘增强成像的方法及系统
US11644791B2 (en) Holographic imaging device and data processing method therefor
CN104102017B (zh) 基于结构光照的随机散射光学超衍射极限成像系统及方法
Kashter et al. Resolving images by blurring: superresolution method with a scattering mask between the observed objects and the hologram recorder
CN113487637B (zh) 基于叠加螺旋相位滤波器的多方向边缘检测方法
CN109212749B (zh) 一种实现边缘增强成像的滤波片及其设计方法
CN105759440A (zh) 基于结构光照的随机散射光学超衍射极限成像方法
Sharma et al. Directional edge enhancement using superposed vortex filter
CN105043543A (zh) 一种可调控超瑞利散斑场的制作装置及其制作方法
Augustin et al. Mask responses for single-pixel terahertz imaging
CN104345626A (zh) 一种离轴数字全息波前记录和再现方法及实施装置
Li et al. Fast quantitative phase imaging based on Kramers-Kronig relations in space domain
Weigel et al. Widefield microscopy with infinite depth of field and enhanced lateral resolution based on an image inverting interferometer
Bu et al. Tunable edge enhancement by higher-order spiral Fresnel incoherent correlation holography system
Li et al. Scalar and vectorial vortex filtering based on geometric phase modulation with a q-plate
Li et al. Far‐Field Synthetic Aperture Imaging via Fourier Ptychography with Quasi‐Plane Wave Illumination
Kim et al. Phase Differences Averaging (PDA) Method for Reducing the Phase Error in Digital Holographic Microscopy (DHM)
Grebenyuk et al. An off-axis digital holographic microscope with quasimonochromatic partially spatially coherent illumination in transmission
Gorevoy et al. Multi-wavelength off-axis digital holographic microscopy with broadly tunable low-coherent sources: theory, performance and limitations
US20230236404A1 (en) System for quantitative differential phase contrast microscopy with isotropic transfer function
He et al. Edge enhancement in three-dimensional vortex imaging based on FINCH by Bessel-like spiral phase modulation
JP2022527827A (ja) 構造化照明顕微鏡法を使用した高度な試料画像化
CN110133849A (zh) 一种点扩散函数重建方法
Wei et al. Prior-free mixed speckle simplex separation strategy for multi-object imaging through thin scattering media beyond the optical memory effect
Xu et al. Implementing the edge enhancement with vortex filter in both linear and nonlinear optics

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant