CN117274287A - 一种基于无干涉编码孔径相关全息术的边缘检测方法 - Google Patents
一种基于无干涉编码孔径相关全息术的边缘检测方法 Download PDFInfo
- Publication number
- CN117274287A CN117274287A CN202311120341.8A CN202311120341A CN117274287A CN 117274287 A CN117274287 A CN 117274287A CN 202311120341 A CN202311120341 A CN 202311120341A CN 117274287 A CN117274287 A CN 117274287A
- Authority
- CN
- China
- Prior art keywords
- image
- amplitude
- spread function
- reconstruction
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 39
- 238000003708 edge detection Methods 0.000 title claims abstract description 36
- 238000001093 holography Methods 0.000 title claims abstract description 19
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 28
- 238000001228 spectrum Methods 0.000 claims abstract description 26
- 230000008569 process Effects 0.000 claims description 10
- 238000001914 filtration Methods 0.000 claims description 8
- 230000001419 dependent effect Effects 0.000 claims description 7
- 238000005286 illumination Methods 0.000 claims description 3
- 230000010287 polarization Effects 0.000 claims description 3
- 230000008901 benefit Effects 0.000 abstract description 6
- 230000004069 differentiation Effects 0.000 abstract description 6
- 238000004364 calculation method Methods 0.000 abstract description 3
- 238000012634 optical imaging Methods 0.000 abstract description 3
- 238000003384 imaging method Methods 0.000 description 8
- 230000003287 optical effect Effects 0.000 description 7
- 230000003595 spectral effect Effects 0.000 description 4
- 230000006872 improvement Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000003909 pattern recognition Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 238000012984 biological imaging Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000000739 chaotic effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/13—Edge detection
-
- 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/0005—Adaptation of holography to specific applications
-
- 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/0443—Digital holography, i.e. recording holograms with digital recording means
-
- 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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/20—Image enhancement or restoration using local operators
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/168—Segmentation; Edge detection involving transform domain methods
-
- 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/0005—Adaptation of holography to specific applications
- G03H2001/0033—Adaptation of holography to specific applications in hologrammetry for measuring or analysing
- G03H2001/0038—Adaptation of holography to specific applications in hologrammetry for measuring or analysing analogue or digital holobjects
-
- 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/0443—Digital holography, i.e. recording holograms with digital recording means
- G03H2001/0454—Arrangement for recovering hologram complex amplitude
- G03H2001/0456—Spatial heterodyne, i.e. filtering a Fourier transform of the off-axis record
-
- 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
- G03H2001/0883—Reconstruction aspect, e.g. numerical focusing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20172—Image enhancement details
- G06T2207/20192—Edge enhancement; Edge preservation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Computing Systems (AREA)
- Holo Graphy (AREA)
Abstract
本发明公开了一种基于无干涉编码孔径相关全息术的边缘检测方法,涉及计算光学成像领域。通过非线性重建NLR(Non‑Linear Reconstruction)算法,将点扩展函数与物体全息图图像重建的互相关计算转换到频域,并将转换后的点扩展函数频谱振幅与物体全息图频谱振幅分别取不同的幂次;引入Laguerre‑Gaussian振幅型叠加涡旋相位滤波器,与改变后的点扩展函数频谱振幅和物体全息图频谱相乘后,再傅里叶逆变换实现边缘图像重建。本发明在非线性重建NLR算法中引入了Laguerre‑Gaussian振幅型叠加涡旋相位滤波器,实现了无干涉编码孔径相关全息系统的边缘检测,通过调整频谱振幅幂次、滤波器拓扑因子和初始相位能够获得各项同性、各向异性和二阶微分的边缘增强图像。本方法与传统非相干全息系统的边缘检测方法相比能够实现方向可控的各向异性和图像的二阶微分,并且具有参数动态可调的优点。
Description
技术领域
本发明涉及一种基于非相干全息术的边缘检测方法,属于计算光学成像领域,尤其涉及一种基于无干涉编码孔径相关全息术的边缘检测应用领域。
背景技术
边缘检测作为一项重要的目标检测技术,在特征识别、天文观测、生物成像和X射线成像等领域有着广泛的应用。与计算机数字边缘检测方法相比,光学边缘检测具有速度快、能耗低、并行处理和信息容量大等优点。光学边缘检测主要是借助4f空间滤波成像系统,通过在频谱平面上加载涡旋相位滤波器进行图像径向对称的希尔伯特变换,基于希尔伯特变换的涡旋滤波由于对物体复折射率梯度的敏感性,可以实现图像的边缘对比度增强。现阶段光学边缘检测技术发展主要取决于光模拟滤波算法的改进和滤波器性能的提升。
非相干全息术具有高分辨率、成本低、无散斑噪声和成像速度快等优点。非相干全息术成像步骤主要分为光学记录和图像重建。菲涅尔非相干相关数字全息技术FINCH(Fresnel Incoherent Correlation Holography)利用非相干光源的自干涉完成全息记录,通过菲涅耳反向传播完成图像重建。在后续研究中,一种无需双波干涉的无干涉编码孔径相关全息I-COACH(Interferenceless Coded Aperture Correlation Holography)术被提出。I-COACH系统只有单一的光束被混沌相位掩膜散射,通过图像传感器分别记录系统的点扩散函数PSF(Point Spread Function)和物体全息图,再通过匹配滤波器、逆滤波器和非线性重建NLR(Non-Linear Reconstruction)算法等各种方法重建图像。
2021年提出了一种基于非相干全息术的计算光学成像技术,并且分别基于无干涉成像系统和FINCH成像系统实现了图像的边缘检测,但是这种边缘检测方法仅能实现一阶微分的各项同性边缘检测。在本发明中提出了一种基于无干涉编码孔径相关全息术的边缘检测方法,通过在NLR算法引入Laguerre-Gaussian振幅型叠加涡旋相位滤波器,在I-COACH系统中实现了各项同性、各向异性和二阶微分的图像边缘检测同时具有动态可调的优点。
发明内容
一种基于无干涉编码孔径相关全息术的边缘检测方法,本发明通过在NLR算法引入高阶拉盖尔-高斯振幅涡旋相位滤波器,实现了各项同性边缘检测和图像二阶微分,通过在NLR算法中叠加两个不同相位的涡旋相位滤波器实现了可调的各向异性边缘增强。
本方法包括以下三个步骤:S1,点扩散函数和物体全息图记录;S2,在NLR算法中引入Laguerre-Gaussian振幅型叠加涡旋相位滤波器对记录的点扩散函数和物体全息图进行滤波处理;S3,调整参数,实现各项同性、各向异性和二阶微分边缘图像重建。
S1:编码相位掩膜CPM(Coded Phase Mask)由Gerchberg-Saxton算法合成,CPM起到了散射调制光的作用,在图像传感器上产生散斑状的强度图。
定义垂直光束传播方向的物平面上任意一物点的像素点坐标为(xs,ys),垂直光束传播方向的CPM平面上像素点坐标为(x1,y1),垂直光束传播方向的图像传感器平面上像素点坐标为(x2,y2)。
假设物平面有一物点振幅为/>其经过CPM掩膜调制后在图像传感器平面强度分布为
其中,表示图像传感器的坐标,/>zs为物平面到CPM距离,f0为CPM前面透镜的焦距,透镜到CPM的距离近似看作为0,λ为光波长,/>为随机CPM的相位,/>为CPM平面的坐标(x1,y1),/>zh为CPM平面到图像传感器平面的距离。
考虑到采用非相干光照明,物体形成的强度分布可以视为PSF强度的叠加。任意一个二维物体产生的强度分布可以表示为
其中,an为与物点相关的常数。
S2:振幅涡旋相位滤波器是一种幅值带通滤波器,是由一个振幅项乘以涡流相位滤波器,可以滤除低频和部分高频噪声,因此本发明选择了Laguerre-Gaussian振幅型涡旋相位滤波器表示为
其中,(r,θ)为频域极坐标,ω1是控制最大振幅位置的参数,l为涡旋滤波器的拓扑因子。
通过叠加两个不同相位的涡旋滤波器实现各向异性边缘增强,叠加涡旋滤波器由下式表示
S2=exp(il1(θ+β))+cexp(-il2(θ+β))
其中,l1和l2为两个拓扑因子,β为初始相位,c为权重因子。因子c决定了正、负涡旋的权重比,控制了边缘增强功率。
Laguerre-Gaussian振幅型叠加涡旋相位滤波器如下式表示
无干涉编码孔径相关全息系统的重建过程实际上是一个光学模式识别问题,NLR重建算法本质其实是空域互相关重建算法转换到了频域进行计算,在NLR算法中首先需要将IOBJ以及IPSF变换到频域如下表示
其中,为傅里叶变换。
将和/>分别改变m和n的幂次,在NLR重建算法中/>将作为重构函数,
作为目标函数。重建图像可由下式表示
其中,为傅里叶逆变换
选择幂次m和n可分别用于调整目标的功率谱和重构函数,当m≠1时对物体光谱大小的影响使得多点目标的重建过程是非线性的,重建过程中的这种非线性提高了重建图像的信噪比,而不会失去线性相关性的移位不变性。
为了实现边缘检测的功能,通过Laguerre-Gaussian振幅型叠加涡旋相位滤波器在NLR算法中对和/>进行滤波处理,如下式表示
于是得到具有边缘检测功能的NLR图像重建算法,如下式表示
S3:当权重因子c=0,拓扑因子l1=3,频谱振幅幂次m=-1和n=0.8时,实现图像的各向同性的边缘图像重建。当权重因子c=0,拓扑因子l1=4,频谱振幅幂次m=-1和n=0.5时,实现图像的各向同性二阶微分的边缘图像重建,二阶微分边缘重建的图像对精细的边缘更为敏感,产生的双线边缘结构更为精细,图像具有更高的对比度。
各向异性边缘增强的方向可以通过调整初始相位β和加权因子c来控制,权重因子c=1,拓扑因子l1=l2=3,频谱振幅幂次m=-1和n=0.8,当β=0、和/>时,分别得到了包括水平方向(H)、垂直方向(V)、对角线方向(D)和反对角线方向(AD)方向的边缘图像。
本方法与传统非相干全息系统的边缘检测方法相比能够实现方向可控的各向异性和图像的二阶微分,并且具有参数动态可调的优点。
附图说明
图1无干涉编码孔径相关全息系统示意图。
图2一种基于无干涉编码孔径相关全息术的边缘检测方法流程图。
图3各项同性边缘图像重建结果:(a)原图,(b)各项同性边缘检测结果。
图4二阶微分边缘图像重建结果:(a)原图,(b)二阶微分边缘检测结果。
图5各向异性边缘图像重建结果:(a)原图,(b)水平方向,(c)垂直方向,(d)对角线方向,(e)反对角线方向。
附图标记说明:
1、单色LED,2、第一透镜,3、目标物体,4、第二透镜,5、偏振片,6、相位型空间光调制器PSLM(Phase spatial light modulator),7、图像传感器。
具体实施方式
为了更好地解释本发明的实施过程,下面将以一个实施例对本发明作进一步详细说明,但本发明不限于该实施例。
实施例
首先给出Laguerre-Gaussian振幅型叠加涡旋相位滤波器在频域的表达式。
Laguerre-Gaussian振幅型涡旋相位滤波器如下式表示
其中,(r,θ)为频域极坐标,ω1是控制最大振幅位置的参数,l为涡旋滤波器的拓扑因子。
涡旋相位延伸了将一维希尔伯特变换到径向空间,保证任意径向方向的相位差为π,得到各向同性边缘增强。为了得到各项异性边缘检测,通过叠加两个不同相位的涡旋滤波器实现,叠加涡旋滤波器由下式表示
S2=exp(il1(θ+β))+cexp(-il2(θ+β))
其中,l1和l2为拓扑因子,β为初始相位,c为权重因子。因子c决定了正、负涡旋的权重比,控制了边缘增强功率。
于是,Laguerre-Gaussian振幅型叠加涡旋相位滤波器如下式表示
定义垂直光束传播方向的物平面上任意物点的像素点坐标为(xs,ys),垂直光束传播方向的CPM平面上像素点坐标为(x1,y1),垂直光束传播方向的图像传感器平面上像素点坐标为(x2,y2)。
单色LED(1)中心波长为532nm,单色LED(1)发出的单色非相干光经过第一透镜(2)聚焦到目标物体(3),单色LED(1)到第一透镜(2)距离为z1,第一透镜(2)到目标物体(3)的距离为z2,第一透镜(2)焦距为f1。z1、z2和f1满足透镜成像公式z1 -1+z2 -1=f1 -1。目标物体(3)首先使用了25μm的针孔,近似看作在目标平面有一物点振幅为/>其在第二透镜(4)所在平面的复振幅可以表示为
其中,zs为物点到第二透镜(4)的距离,λ为光源中心波长,/>zs的大小等于第二透镜(4)的焦距f0。经过第二透镜(4)后携带目标物体(3)信息的光束被准直,因此第二透镜(4)到PSLM(6)的距离可以近似看作0,zs也为物点到PSLM(6)的距离。
第二透镜(4)出射的光束经过偏振片(5)后入射到PSLM(6),偏振片(5)与PSLM(6)的调制偏振方向一致,PSLM(6)上加载了由Gerchberg-Saxton算法合成的CPM,经过CPM掩膜调制后在图像传感器上的平面强度分布为
其中,表示图像传感器的坐标(x2,y2),/> 为随机CPM的相位,/>为CPM平面的坐标(x1,y1),/>zh为CPM平面到图像传感器(7)的距离。
此时图像传感器(7)上记录的为PSF强度图,反映了一个物点经过系统的强度分布。
在系统目标物体(3)处同一位置将针孔换成目标图像,记录物体全息图,考虑到采用非相干光照明,物体形成的强度分布可以视为PSF强度的叠加。任意一个二维物体产生的强度分布可以表示为
其中,an为与物点相关的常数。
无干涉编码孔径相关全息术系统的重建过程实际上是一个光学模式识别问题,通过IOBJ与IPSF的相互关联重建图像,本质仍是一种互相关算法,而NLR算法本质其实是互相关重建算法转换到了频域进行计算,在NLR算法中首先需要将IOBJ以及IPSF变换到频域如下表示
将和/>·分别改变m和n的幂次,重建图像可由下式表示
在NLR重建算法中将·作为重构函数,/>作·为目标函数,选择参数m和n可分别用于调整目标的功率谱和重构函数,当m≠1时对物体光谱大小的影响使得多点目标的重建过程是非线性的。重建过程中的这种非线性提高了重建图像的信噪比,而不会失去线性相关性的移位不变性。
为了实现边缘检测的功能,通过Laguerre-Gaussian振幅型叠加涡旋相位滤波器在NLR算法中对·
和·进·行滤波处理,如下式表示
得到具有边缘检测功能的NLR图像重建算法,如下式表示
由于每个对象点的位置都在的相位中编码,因此保持了移位不变性的性质,而参数m仅对/>的幅度而不是相位进行操作,理论上,产生最清晰图像点的m和n的值是满足方程m+n=0的值。不同的实验环境产生不同的最优参数,不遵循这个方程。在逆滤波器(m或n=-1)和匹配滤波器(m或n=1)之间的范围内寻找最佳成像的参数m和n对。
当权重因子c=0,拓扑因子l1=3,频谱振幅幂次m=-1和n=0.8时,实现图像的各向同性的边缘图像重建。如图3所示,图3(a)为原图,图3(b)为各向同性边缘图像重建结果。
当权重因子c=0,拓扑因子l1=4,频谱振幅幂次m=-1和n=0.5时,实现图像的各向同性二阶微分的边缘图像重建。如图4所示,图4(a)为原图,图4(b)为各向同性二阶微分边缘图像重建结果,需要注意的是二阶微分计算的是一阶微分的变化率,二阶微分边缘重建的图像对精细的边缘更为敏感,因此通常会产生更窄或更精确的边缘,如图4(b)所示可以看到圆形图案存在两条更为精细的边缘,图像具有更高的对比度,这在图像处理中具有重要意义。
在任意方向上的各向异性边缘检测也可以通过在频域中调整滤波器的初始相位角β来实现,各向异性边缘增强的方向可以通过调整初始相位β和加权因子c来控制。
当权重因子c=1、拓扑因子l1=l2=3,频谱振幅幂次m=-1和n=0.8,实现图像的各向异性的边缘图像重建。为了更好的观察,选择了圆形图案被用作示例来说明方向增强效应,如图5所示,图5(a)为原图,当β=0、和/>时,分别得到了包括图5(b)水平方向(H)、图5(c)垂直方向(V)、图5(d)对角线方向(D)和图5(e)反对角线方向(AD)的边缘图像。
本方法与传统非相干全息系统的边缘检测方法相比能够实现方向可控的各向异性和图像的二阶微分,并且具有参数动态可调的优点。
以上所述仅为本发明的较佳实施案例,不代表本发明的保护范围。凡在本发明的原则之内所做的任何修改、等同替换、改进等,均包含在本发明的保护范围。
Claims (5)
1.一种基于无干涉编码孔径相关全息术的边缘检测方法,其特征在于,所述方法采用具体无干涉编码孔径相关全息术的装置包括单色LED、第一透镜、物体、第二透镜、偏振片、相位型空间光调制器PSLM(Phase spatial light modulator)、图像传感器,其中:
所述单色LED发出的单色非相干光经过所述第一透镜聚焦到所述物体上,所述物体首先使用了25μm的针孔,近似看作在所述物体平面有一物点;所述物体到所述第二透镜的距离大小等于所述第二透镜的焦距,经过所述第二透镜后光束被准直;准直后的光束经过所述偏振片后入射到所述PSLM,所述偏振片偏振方向与所述PSLM的调制偏振轴对齐,所述PSLM上加载了编码相位掩膜CPM(Coded Phase Mask);经所述PSLM反射的光束入到所述图像传感器,所述图像传感器此时记录的图像作为点扩散函数;仅将所述物体处的针孔换成目标图像,光束传播过程与所述物体使用针孔时完全相同,此时所述图像传感器上记录的图像作为物体全息图;
所述方法包括以下三个步骤:S1,点扩散函数和物体全息图记录;S2,在非线性重建NLR(Non-Linear Reconstruction)算法中引入Laguerre-Gaussian振幅型叠加涡旋相位滤波器对记录的点扩散函数和物体全息图进行滤波处理;S3,调整参数,实现各项同性、各向异性和二阶微分边缘图像重建;
定义垂直光束传播方向的物体平面上任意一物点的像素点坐标为(xs,ys),垂直光束传播方向的CPM平面上像素点坐标为(x1,y1),垂直光束传播方向的图像传感器平面上像素点坐标为(x2,y2);
S1:记录所述点扩散函数和物体全息图,在无干涉编码孔径相关全息系统的所述物体平面依次放置针孔和目标图像,针孔位置和目标图像位置处在同一平面,在所述图像传感器上分别记录了针孔的强度分布IPSF和目标图像的强度分布IOBJ;
S2:通过在所述NLR算法中引入所述Laguerre-Gaussian振幅型叠加涡旋相位滤波器对记录的所述点扩散函数频谱和物体全息图频谱进行滤波处理。在所述NLR算法中,首先需要改变所述物体全息图频谱和点扩散函数频谱,所述物体全息图频谱振幅改变m次幂,相位不变,点扩散函数频谱振幅改变n次幂,相位取负;将所述Laguerre-Gaussian振幅型叠加涡旋相位滤波器、改变后的全息图频谱和改变后的点扩散函数频谱三者相乘,再傅里叶逆变换重建图像边缘;
S3:选择最佳图像边缘重建的点扩散函数频谱振幅幂次m和物体全息图频谱振幅幂次n,同时调整滤波器拓扑因子l、权重因子c和初始相位β,获得无干涉编码孔径相关全息系统的各项同性,各向异性和二阶微分边缘图像。
2.如权利要求1所述的一种基于无干涉编码孔径相关全息术的边缘检测方法,其特征在于:在所述S1中,在所述物体平面放置针孔记录系统的点扩展函数,近似看做所述物平面有一物点其振幅为/>其表达式为
式中,表示所述图像传感器的坐标(x2,y2),/>zs为所述物体平面到所述CPM平面的距离,λ为所述单色LED的中心波长,/> f0为所述第二透镜的焦距,所述第二透镜到所述CPM平面的距离近似看作为0,/> 为所述CPM的相位,/>为所述CPM平面的坐标(x1,y1),zh为所述CPM平面到所述图像传感器平面的距离;
在所述物体平面的位置将针孔替换为目标图像,所述图像传感器记录物体全息图,其表达式为
式中,an为与物点相关的常数。因为使用非相干光照明,任意一个二维物体产生的强度分布可以视为点扩展函数强度的叠加。
3.如权利要求1所述的一种基于无干涉编码孔径相关全息术的边缘检测方法,其特征在于:在所述S2中,所述Laguerre-Gaussian振幅型叠加涡旋相位滤波器在频域的表达式为
式中,(r,θ)为频域极坐标,ω1是控制最大振幅位置的参数,l1和l2为两个拓扑因子,β为初始相位,c为权重因子。初始相位β控制了各项异性边缘增强的方向,权重因子c决定了正、负涡旋的权重比,控制了边缘增强功率。
4.如权利要求1所述的一种基于无干涉编码孔径相关全息术的边缘检测方法,其特征在于:在所述S2中,通过所述Laguerre-Gaussian振幅型叠加涡旋相位滤波器对改变后的点扩散函数频谱和物体全息图频谱进行滤波处理,其表达式为
式中, 为傅里叶变换。
对S取傅里叶逆变换,得到具有边缘检测功能的NLR图像重建算法,表达式如下
式中,为傅里叶逆变换。
5.如权利要求1所述的一种基于无干涉编码孔径相关全息术的边缘检测方法,其特征在于:在所述S3中,选择最佳图像边缘重建的点扩散函数频谱振幅幂次m和物体全息图频谱振幅幂次n,调整拓扑因子l、权重因子c和初始相位β,得到各项同性、各向异性和二阶微分边缘图像重建。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311120341.8A CN117274287B (zh) | 2023-08-31 | 2023-08-31 | 一种基于无干涉编码孔径相关全息术的边缘检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311120341.8A CN117274287B (zh) | 2023-08-31 | 2023-08-31 | 一种基于无干涉编码孔径相关全息术的边缘检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN117274287A true CN117274287A (zh) | 2023-12-22 |
CN117274287B CN117274287B (zh) | 2024-06-18 |
Family
ID=89215087
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311120341.8A Active CN117274287B (zh) | 2023-08-31 | 2023-08-31 | 一种基于无干涉编码孔径相关全息术的边缘检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117274287B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN118096802A (zh) * | 2024-01-02 | 2024-05-28 | 哈尔滨理工大学 | 一种基于并行任意阶拓扑光学微分的图像处理方法及装置 |
Citations (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2007116365A2 (en) * | 2006-04-07 | 2007-10-18 | Ecole Polytechnique Federale De Lausanne (Epfl) | Method and apparatus to measure and compute the amplitude point spread function and associated parameters of a coherent optical imaging system |
US20120116703A1 (en) * | 2009-04-24 | 2012-05-10 | Nicolas Pavillon | Method and apparatus for enhanced spatial bandwidth wavefronts reconstructed from digital interferograms or holograms |
CN103984103A (zh) * | 2014-02-18 | 2014-08-13 | 上海大学 | 一种利用计算全息图产生涡旋光束的方法 |
CN105204311A (zh) * | 2015-07-06 | 2015-12-30 | 电子科技大学 | 一种基于高斯切趾的光学扫描全息边缘检测方法 |
CN106485764A (zh) * | 2016-11-02 | 2017-03-08 | 中国科学技术大学 | Mri图像的快速精确重建方法 |
CN106896498A (zh) * | 2017-03-08 | 2017-06-27 | 华南师范大学 | 拉盖尔高斯涡旋光束的产生/三维重建装置及方法 |
US20180284025A1 (en) * | 2017-03-31 | 2018-10-04 | Richard Gozali | Oam microscope for edge enhancement of biomedical and condensed matter samples and objects |
CN108873323A (zh) * | 2018-07-09 | 2018-11-23 | 苏州大学 | 一种实现边缘增强成像的方法及系统 |
CN110488397A (zh) * | 2019-07-25 | 2019-11-22 | 西安电子科技大学 | 复合透镜、可编码单帧成像装置及成像方法 |
CN111307279A (zh) * | 2019-11-19 | 2020-06-19 | 中国人民解放军战略支援部队航天工程大学 | 一种基于三步相移法的涡旋光模式检测方法 |
US20200378745A1 (en) * | 2019-05-28 | 2020-12-03 | University Of Central Florida Research Foundation, Inc. | Iterative Optical Diffraction Tomography (iODT) Method and Applications |
WO2021093259A1 (zh) * | 2019-11-14 | 2021-05-20 | 深圳大学 | 一种任意奇点光束阶数检测装置及方法 |
CN113487637A (zh) * | 2021-07-05 | 2021-10-08 | 南京邮电大学 | 基于叠加螺旋相位滤波器的多方向边缘检测方法 |
CN113885196A (zh) * | 2021-09-17 | 2022-01-04 | 常州北邮新一代信息技术研究院有限公司 | 基于光学成像系统的图像边缘增强方法、装置 |
CN114019690A (zh) * | 2021-11-04 | 2022-02-08 | 浙江大学 | 产生任意阶光学涡旋阵列和带缺陷有限光晶格的光学系统 |
CN114429549A (zh) * | 2020-10-28 | 2022-05-03 | 北京小米移动软件有限公司 | 图像边缘提取方法、图像边缘提取装置及存储介质 |
CN114488548A (zh) * | 2022-04-18 | 2022-05-13 | 华南师范大学 | 一种高操作性光镊的产生方法及系统 |
CN114578576A (zh) * | 2022-05-05 | 2022-06-03 | 华南师范大学 | 一种旋角、旋速和聚焦可控的光学扳手的制备方法及系统 |
CN114754674A (zh) * | 2022-03-29 | 2022-07-15 | 中国人民解放军战略支援部队航天工程大学 | 一种基于共轭叠加态涡旋光的物体旋转中心定位方法 |
WO2023279424A1 (zh) * | 2021-07-08 | 2023-01-12 | 苏州大学 | 一种多边形部分相干涡旋光束产生系统及方法 |
CN115984314A (zh) * | 2022-11-25 | 2023-04-18 | 哈尔滨理工大学 | 一种基于计算全息二阶微分的图像边缘检测方法及系统 |
CN116051584A (zh) * | 2023-01-13 | 2023-05-02 | 哈尔滨理工大学 | 一种目标边缘红外场景生成与高效率边缘检测方法及装置 |
CN116485829A (zh) * | 2023-04-24 | 2023-07-25 | 哈尔滨理工大学 | 基于复振幅滤波拓扑空间微分运算的图像边缘检测方法、装置及复振幅调制模块 |
-
2023
- 2023-08-31 CN CN202311120341.8A patent/CN117274287B/zh active Active
Patent Citations (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2007116365A2 (en) * | 2006-04-07 | 2007-10-18 | Ecole Polytechnique Federale De Lausanne (Epfl) | Method and apparatus to measure and compute the amplitude point spread function and associated parameters of a coherent optical imaging system |
US20120116703A1 (en) * | 2009-04-24 | 2012-05-10 | Nicolas Pavillon | Method and apparatus for enhanced spatial bandwidth wavefronts reconstructed from digital interferograms or holograms |
CN103984103A (zh) * | 2014-02-18 | 2014-08-13 | 上海大学 | 一种利用计算全息图产生涡旋光束的方法 |
CN105204311A (zh) * | 2015-07-06 | 2015-12-30 | 电子科技大学 | 一种基于高斯切趾的光学扫描全息边缘检测方法 |
CN106485764A (zh) * | 2016-11-02 | 2017-03-08 | 中国科学技术大学 | Mri图像的快速精确重建方法 |
CN106896498A (zh) * | 2017-03-08 | 2017-06-27 | 华南师范大学 | 拉盖尔高斯涡旋光束的产生/三维重建装置及方法 |
US20180284025A1 (en) * | 2017-03-31 | 2018-10-04 | Richard Gozali | Oam microscope for edge enhancement of biomedical and condensed matter samples and objects |
CN108873323A (zh) * | 2018-07-09 | 2018-11-23 | 苏州大学 | 一种实现边缘增强成像的方法及系统 |
US20200378745A1 (en) * | 2019-05-28 | 2020-12-03 | University Of Central Florida Research Foundation, Inc. | Iterative Optical Diffraction Tomography (iODT) Method and Applications |
CN110488397A (zh) * | 2019-07-25 | 2019-11-22 | 西安电子科技大学 | 复合透镜、可编码单帧成像装置及成像方法 |
WO2021093259A1 (zh) * | 2019-11-14 | 2021-05-20 | 深圳大学 | 一种任意奇点光束阶数检测装置及方法 |
CN111307279A (zh) * | 2019-11-19 | 2020-06-19 | 中国人民解放军战略支援部队航天工程大学 | 一种基于三步相移法的涡旋光模式检测方法 |
CN114429549A (zh) * | 2020-10-28 | 2022-05-03 | 北京小米移动软件有限公司 | 图像边缘提取方法、图像边缘提取装置及存储介质 |
CN113487637A (zh) * | 2021-07-05 | 2021-10-08 | 南京邮电大学 | 基于叠加螺旋相位滤波器的多方向边缘检测方法 |
WO2023279424A1 (zh) * | 2021-07-08 | 2023-01-12 | 苏州大学 | 一种多边形部分相干涡旋光束产生系统及方法 |
CN113885196A (zh) * | 2021-09-17 | 2022-01-04 | 常州北邮新一代信息技术研究院有限公司 | 基于光学成像系统的图像边缘增强方法、装置 |
CN114019690A (zh) * | 2021-11-04 | 2022-02-08 | 浙江大学 | 产生任意阶光学涡旋阵列和带缺陷有限光晶格的光学系统 |
CN114754674A (zh) * | 2022-03-29 | 2022-07-15 | 中国人民解放军战略支援部队航天工程大学 | 一种基于共轭叠加态涡旋光的物体旋转中心定位方法 |
CN114488548A (zh) * | 2022-04-18 | 2022-05-13 | 华南师范大学 | 一种高操作性光镊的产生方法及系统 |
CN114578576A (zh) * | 2022-05-05 | 2022-06-03 | 华南师范大学 | 一种旋角、旋速和聚焦可控的光学扳手的制备方法及系统 |
CN115984314A (zh) * | 2022-11-25 | 2023-04-18 | 哈尔滨理工大学 | 一种基于计算全息二阶微分的图像边缘检测方法及系统 |
CN116051584A (zh) * | 2023-01-13 | 2023-05-02 | 哈尔滨理工大学 | 一种目标边缘红外场景生成与高效率边缘检测方法及装置 |
CN116485829A (zh) * | 2023-04-24 | 2023-07-25 | 哈尔滨理工大学 | 基于复振幅滤波拓扑空间微分运算的图像边缘检测方法、装置及复振幅调制模块 |
Non-Patent Citations (12)
Title |
---|
A. FASSBENDER: "Broadband Laguerre-Gaussian metasurfaces and direct phase mapping (Conference Presentation)", PHYSICS, ENGINEERING, 23 May 2018 (2018-05-23) * |
ABDERRAHMEN TRICHILI: "Encoding information using Laguerre Gaussian modes", SPIE OPTICAL ENGINEERING, 25 August 2015 (2015-08-25) * |
S. SINGH: "Tailoring Large Asymmetric Laguerre–Gaussian Beam Array Using Computer-Generated Holography", COMPUTER SCIENCE, 24 February 2023 (2023-02-24) * |
刘义东;王建东;刘普生;: "一种提高信息传输密度的新方法研究", 实验科学与技术, no. 05, 28 October 2013 (2013-10-28) * |
卜远壮: "基于多模Laguerre-Gaussian光束相位调制的FINCH系统边缘检测应用研究", CNKI硕士电子期刊库, 15 May 2022 (2022-05-15) * |
李志焕;吴福全;陈君;崔祥霞;: "基于新型螺旋相位滤波器提高图像边缘增强对比度的方法", 泰山学院学报, no. 06, 25 November 2012 (2012-11-25) * |
汪慧超;胡阿健;陈培锋;: "空间光调制器产生拉盖尔-高斯光束方法研究", 激光技术, no. 03, 25 May 2017 (2017-05-25) * |
王林;袁操今;聂守平;李重光;张慧力;赵应春;张秀英;冯少彤;: "数字全息术测定涡旋光束拓扑电荷数", 物理学报, no. 24, 23 December 2014 (2014-12-23) * |
范全平;王少义;魏来;杨祖华;张强强;陈勇;巫殷忠;曹磊峰;: "一种产生边缘增强图像的螺旋波带片编码成像技术", 激光与光电子学进展, no. 03, 30 September 2017 (2017-09-30) * |
邹丽;王乐;张士兵;赵生妹;: "基于波前校正的轨道角动量复用通信系统抗干扰研究", 通信学报, no. 10, 25 October 2015 (2015-10-25) * |
陈君;崔祥霞;肖静;杨兆华;刘婷婷;: "基于径向希尔伯特变换的图像边缘增强特性研究", 光子学报, no. 03, 15 March 2011 (2011-03-15) * |
高衍;付时尧;张世坤;王庆;胡新奇;: "高阶涡旋光束生成技术研究", 光学技术, no. 02, 15 March 2017 (2017-03-15) * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN118096802A (zh) * | 2024-01-02 | 2024-05-28 | 哈尔滨理工大学 | 一种基于并行任意阶拓扑光学微分的图像处理方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN117274287B (zh) | 2024-06-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104102017B (zh) | 基于结构光照的随机散射光学超衍射极限成像系统及方法 | |
Jesacher et al. | Shadow effects in spiral phase contrast microscopy | |
CN108873323B (zh) | 一种实现边缘增强成像的方法及系统 | |
CN105204311B (zh) | 一种基于高斯切趾的光学扫描全息边缘检测方法 | |
JP6424313B2 (ja) | ホログラフィック顕微鏡および高分解能ホログラム画像用のデータ処理方法 | |
CN105759440A (zh) | 基于结构光照的随机散射光学超衍射极限成像方法 | |
Kostencka et al. | Autofocusing method for tilted image plane detection in digital holographic microscopy | |
CN113487637B (zh) | 基于叠加螺旋相位滤波器的多方向边缘检测方法 | |
WO2006072581A1 (en) | Spiral phase contrast imaging in microscopy | |
CN115984314B (zh) | 一种基于计算全息二阶微分的图像边缘检测方法及系统 | |
CN104345626A (zh) | 一种离轴数字全息波前记录和再现方法及实施装置 | |
CN109164691A (zh) | 一种实现透过散射介质聚焦的离轴数字全息相位共轭方法 | |
CN115308185A (zh) | 高速无伪影超分辨结构光照明显微成像方法 | |
Bu et al. | Tunable edge enhancement by higher-order spiral Fresnel incoherent correlation holography system | |
CN114324245A (zh) | 基于部分相干结构光照明的定量相位显微装置和方法 | |
CN113376993B (zh) | 编码孔径相关全息术的虚拟点扩散函数记录方法 | |
Wang et al. | Zero-order term suppression in off-axis holography based on deep learning method | |
CN117274287B (zh) | 一种基于无干涉编码孔径相关全息术的边缘检测方法 | |
CN109613712A (zh) | 一种采用轴棱锥相位分布和方位角随机调制产生无衍射散斑的方法及其实现装置 | |
CN113504717A (zh) | 基于时空聚焦的均匀全息双光子显微系统 | |
Johnson et al. | Analysis of a joint transform correlator using a phase-only spatial light modulator | |
Zhang et al. | Computational incoherent edge enhancement imaging based on self-interference digital holography | |
Palacios et al. | Methods of Fourier optics in digital holographic microscopy | |
Kumar et al. | Recognition of three-dimensional objects using joint fractional correlator and nonlinear joint fractional correlator with the help of digital Fresnel holography: a comparative study | |
Feng et al. | Volume holographic wavelet correlation processor |
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 |