CN106406065A - 一种全息再现图像的零级像消除方法 - Google Patents
一种全息再现图像的零级像消除方法 Download PDFInfo
- Publication number
- CN106406065A CN106406065A CN201610849069.0A CN201610849069A CN106406065A CN 106406065 A CN106406065 A CN 106406065A CN 201610849069 A CN201610849069 A CN 201610849069A CN 106406065 A CN106406065 A CN 106406065A
- Authority
- CN
- China
- Prior art keywords
- image
- zero
- order
- contourlet
- order image
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 34
- 230000008030 elimination Effects 0.000 title abstract description 7
- 238000003379 elimination reaction Methods 0.000 title abstract description 7
- 238000012545 processing Methods 0.000 claims abstract description 7
- 230000009466 transformation Effects 0.000 claims description 27
- 230000005540 biological transmission Effects 0.000 claims description 4
- 239000000203 mixture Substances 0.000 claims description 4
- 230000008569 process Effects 0.000 claims description 3
- 230000000694 effects Effects 0.000 abstract description 5
- 230000000717 retained effect Effects 0.000 abstract description 2
- 230000000007 visual effect Effects 0.000 abstract 1
- 238000001914 filtration Methods 0.000 description 7
- 238000000354 decomposition reaction Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 239000000686 essence Substances 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000010363 phase shift Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000002203 pretreatment Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000011426 transformation method 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/22—Processes or apparatus for obtaining an optical image from holograms
- G03H1/2286—Particular reconstruction light ; Beam properties
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Holo Graphy (AREA)
Abstract
本发明公开了一种全息再现图像的零级像消除方法,涉及全息技术领域,利用Contourlet变换去除零级像,包括以下步骤:1:利用M激光器得到所需要处理的数字全息图像;2:读入数字全息图像,进行Contourlet变换,得到低频系数和高频系数;3:在数字全息再现像中,通过对Contourlet分解,分离出全息图样的直流成分,从而达到去除噪声的目的;4:对低频系数和高频系数进行Contourlet变换逆变换,得到去除零级像后的数字全息图像;5:对处理过后的全息图进行再现像处理,得到结果。本发明从主观视觉效果上能够清楚地发现其良好的零级像去除能力,更好地保留了物像的信息,从能量角度进行对比也发现该方法效果更优。
Description
技术领域
本发明涉及全息技术领域,尤其是一种全息再现图像的零级像消除方法。
背景技术
消除零级像的传统方法需要移相器材以及其他辅助设备,且对周围的环境要求极高,十分繁琐。随着数字图像处理技术的发展,越来越多的学者关注于利用图像处理技术去除零级像,提出了许多消除零级像的方法。如零级像的空域滤波预处理消除法,该方法只需记录一幅全息图像,不需要增加其他的实验装置,直接对全息图像在空域进行预处理,便可消除零级像。例如利用频域滤波的零级像消除方法,该方法用一个有限脉冲响应滤波器,滤除零级像的部分,从而消除零级衍射斑对再现像的影响。例如利用小波变换消除零级像的方法,通过小波变换将数字全息图变换到小波域,去除低频部分,再利用小波逆变换重现全息图。尽管这些方法极大地改善了全息再现像的像质,但是他们在图像增强质量、再现速度等方面还有所欠缺,如空域滤波法再现像较为模糊,零级像消除不彻底。而频域滤波法需要多次采用傅里叶变换和频谱滤波,再现速度变慢,另外窗函数的选择也相对苛刻,全息图的细节成分也容易被滤除。小波变换则在图像的细节表达上有所保留,方向性不足,去除全息图的低频成分,会有损部分真实像的信息。
发明内容
本发明提出的一种全息再现图像的零级像消除方法,有更好的图像增强效果,不仅零级像去除得更为彻底,真实像也得到了一定程度的增强。
本发明的技术方案是这样实现的:
一种全息再现图像的零级像消除方法,利用Contourlet变换去除零级像,该方法包括以下步骤:
步骤1:利用M激光器发出的波长为λ=6300nm,波束k=2π/λ,菲涅尔传输距离z=0.3m的光波得到所需要处理的数字全息图像;
步骤2:读入数字全息图像,确定Contourlet分解的层数,进行Contourlet变换,得到低频系数和高频系数d0,d1,1,……d1,j;d2,1,……d2,j,……di,1,……di,j;
步骤3:在数字全息再现像中,零级像是低通成分,与直流成分相应,即d0,可通过对Contourlet分解,分离出全息图样的直流成分,从而达到去除噪声的目的;
步骤4:对低频系数和高频系数进行Contourlet变换逆变换,得到去除零级像后的数字全息图像;
步骤5:对处理过后的全息图进行再现像处理,得到结果。
本发明通过提供的全息再现图像的零级像消除方法,其有益效果在于:从主观视觉效果上能够清楚地发现其良好的零级像去除能力,更好地保留了物像的信息,从能量角度进行对比也发现该方法效果更优。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为Contourlet变换原理图;
图2为Contourlet变换滤波器组结构图;
图3为Contourlet频率分解原理图;
图4为三维文字实物零级像消除效果图;
图5为三维骰子零级像消除效果图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
Contourlet变换的基本原理是首先利用拉普拉斯塔形滤波器(LaplacianPyramid,LP)获取图像的奇异点,再利用方向滤波器组(Directional Filter Bank,DFB)对方向相同的奇异点滤波,得到高频方向上的子带系数。即Contourlet变换的基本原理可以概括为子带滤波和方向滤波,如图1所示。
Contourlet变换的方向滤波器组结构如图2所示,图像经金字塔形方向滤波器组(Pyramidal Direction Filter Bank,PDFB)滤波,可得到一个低频子图像和高频子图像。尺度的多少和分解的方向数成正比,即尺度越大得到的方向数也越多。Contourlet变换频域分解原理如图3所示。
图像经Contourlet变换后得到的系数非常零散,绝大多数得到的Contourlet变换系数幅值等于0,图像的轮廓和边缘周围的系数幅值一般较大,这些幅值较大的系数存在一定的关联性和继承性。
零级像主要集中在数字全息频谱的低频部分,并且占据了其绝大部分能量。消除零级像的最直接的方法就是采用频域滤波,直接滤除低频部分,但是这种方法在消除零级像的同时也会损失部分真实像的信息。小波很好地弥补了这一缺陷,但是小波方向选择性表达不足,分解后对数字全息的细节纹理信息表达不足,可能出现模糊真实像的情况。所以,本申请选用有更好的方向选择性的Contourlet变换去除零级像。该方法步骤如下:
步骤1:利用M激光器发出的波长为λ=6300nm,波束k=2π/λ,菲涅尔传输距离z=0.3m的光波得到所需要处理的数字全息图像;
步骤2:读入数字全息图像,确定Contourlet分解的层数,进行Contourlet变换,得到低频系数和高频系数d0,d1,1,……d1,j;d2,1,……d2,j,……di,1,……di,j;
步骤3:在数字全息再现像中,零级像是低通成分,与直流成分相应,即d0,可通过对Contourlet分解,分离出全息图样的直流成分,从而达到去除噪声的目的;
步骤4:对低频系数和高频系数进行Contourlet变换逆变换,得到去除零级像后的数字全息图像;
步骤5:对处理过后的全息图进行再现像处理,得到结果。
本实施例使用的全息图是利用Matlab计算仿真得到的,其中激光器λ=6300nm,波束k=2π/λ,菲涅尔传输距离z=0.3m。记录的物体是三维文字实物和三维骰子图像。
图4、5中,(a)是数字全息再现的直接结果;(b)为空域滤波法再现像结果;(c)为对全息图进行二阶巴特沃斯滤波器消除零级像后的再现像结果;(d)是对图像进行小波变换滤去低频成分后重建的再现像结果;(e)是对图像进行两层Contourlet变换,滤去低频成分后的重建结果。从结果中我们不难看出如果不对全息图像进行预处理,零级像会占据很大的能量,并且由于CCD板的面上的凹凸不平会产生隶属于乘性噪声的散斑,使成像质量非常不好。从实验结果上看空域滤波法没有很好地消除零级像;巴特沃斯滤波法似乎比较好地消除了零级像,但从结果中不难发现,由于高频部分没有相对增强,造成了再现像的对比度不强,像质改善得不够;小波变换法使得像质有了明显提高,能够在消除零级像的同时较好地提高真实物像像质,Contourlet变换法相比小波变换法显得再现像的像质更清晰,零级像消除效果更好。
本申请实验除主观视觉观察外,还利用处理前后零级像所在区域的能量得到的数字全息图像的衰减程度来评价零级像去除效果,如表1所示。
表1区域的能量衰减程度
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (1)
1.一种全息再现图像的零级像消除方法,利用Contourlet变换去除零级像,其特征在于,该方法包括以下步骤:
步骤1:利用M激光器发出的波长为λ=6300nm,波束k=2π/λ,菲涅尔传输距离z=0.3m的光波得到所需要处理的数字全息图像;
步骤2:读入数字全息图像,确定Contourlet分解的层数,进行Contourlet变换,得到低频系数和高频系数d0,d1,1,……d1,j;d2,1,……d2,j,……di,1,……di,j;
步骤3:在数字全息再现像中,零级像是低通成分,与直流成分相应,即d0,可通过对Contourlet分解,分离出全息图样的直流成分,从而达到去除噪声的目的;
步骤4:对低频系数和高频系数进行Contourlet变换逆变换,得到去除零级像后的数字全息图像;
步骤5:对处理过后的全息图进行再现像处理,得到结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610849069.0A CN106406065A (zh) | 2016-09-24 | 2016-09-24 | 一种全息再现图像的零级像消除方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610849069.0A CN106406065A (zh) | 2016-09-24 | 2016-09-24 | 一种全息再现图像的零级像消除方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106406065A true CN106406065A (zh) | 2017-02-15 |
Family
ID=57997515
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610849069.0A Pending CN106406065A (zh) | 2016-09-24 | 2016-09-24 | 一种全息再现图像的零级像消除方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106406065A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109445262A (zh) * | 2018-12-28 | 2019-03-08 | 天津利信科技有限公司 | 全息数据的存储方法及装置 |
CN115065817A (zh) * | 2022-05-24 | 2022-09-16 | 昆明理工大学 | 一种全息图压缩方法、编码器和全息再现模块 |
-
2016
- 2016-09-24 CN CN201610849069.0A patent/CN106406065A/zh active Pending
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109445262A (zh) * | 2018-12-28 | 2019-03-08 | 天津利信科技有限公司 | 全息数据的存储方法及装置 |
CN115065817A (zh) * | 2022-05-24 | 2022-09-16 | 昆明理工大学 | 一种全息图压缩方法、编码器和全息再现模块 |
CN115065817B (zh) * | 2022-05-24 | 2024-01-12 | 昆明理工大学 | 一种全息图压缩方法、编码器和全息再现模块 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Knaus et al. | Dual-domain image denoising | |
Mohideen et al. | Image de-noising using discrete wavelet transform | |
CN103854264B (zh) | 一种基于改进型阈值函数的小波变换图像去噪方法 | |
CN101482617B (zh) | 基于非下采样轮廓波的合成孔径雷达图像去噪方法 | |
CN104159094B (zh) | 一种改善光学扫描全息层析成像效果的方法 | |
JP5348145B2 (ja) | 画像処理装置および画像処理プログラム | |
CN103020916B (zh) | 一种二维希尔伯特变换和bemd结合的图像去噪方法 | |
CN102930512A (zh) | 基于HSV色彩空间结合Retinex的水下图像增强方法 | |
WO2007114363A1 (ja) | 画像処理方法 | |
CN103399348A (zh) | 基于Shearlet变换的地震信号去噪方法 | |
CN106406065A (zh) | 一种全息再现图像的零级像消除方法 | |
CN107895354A (zh) | 一种基于非下采样Shearlet变换的高分三SAR图像降斑方法 | |
CN103489157A (zh) | 一种增强合成孔径声纳干涉图质量的滤波方法及系统 | |
CN101540039A (zh) | 单帧图象超分辨方法 | |
CN105976340A (zh) | 改进的基于小波分解的旋滤波算法 | |
CN110175959B (zh) | 一种台风云图增强方法 | |
CN104156925B (zh) | 对超声图像去除散斑和边界增强的处理方法及系统 | |
CN109242806A (zh) | 一种基于高斯核函数的小波域值去噪方法 | |
CN102289800A (zh) | 基于Treelet的Contourlet域图像去噪方法 | |
Tayade et al. | Medical image denoising and enhancement using DTCWT and Wiener filter | |
CN106408537A (zh) | 一种裸眼3d再现图像散斑抑制方法 | |
Cho et al. | Laplace operation-based DC noise reduction in digital holography | |
KR101734419B1 (ko) | 이미지 복원 장치 및 방법 | |
Jeng et al. | Application of multiresolution analysis in removing ground-penetrating radar noise | |
Xuhui et al. | Ultrasound image denoising based on the contourlet transform and anisotropic diffusion |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
CB02 | Change of applicant information | ||
CB02 | Change of applicant information |
Address after: 100041 Beijing, Shijingshan District Xing Xing street, building 30, No. 3, building 2, C-0058 Applicant after: Beijing micro cloud Software Co., Ltd. Address before: 100041 Beijing, Shijingshan District Xing Xing street, building 30, No. 3, building 2, C-0058 Applicant before: Micro light capital investment management (Beijing) Co., Ltd. |
|
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20170215 |