CN105761212A - 基于dmd编码的哈达码变换高灵敏成像的方法 - Google Patents
基于dmd编码的哈达码变换高灵敏成像的方法 Download PDFInfo
- Publication number
- CN105761212A CN105761212A CN201410779177.6A CN201410779177A CN105761212A CN 105761212 A CN105761212 A CN 105761212A CN 201410779177 A CN201410779177 A CN 201410779177A CN 105761212 A CN105761212 A CN 105761212A
- Authority
- CN
- China
- Prior art keywords
- coding
- image
- template
- dmd
- dimensional
- 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
Landscapes
- Image Processing (AREA)
Abstract
本发明公开了一种基于DMD编码的哈达码变换高灵敏成像的方法。基于探测器噪声独立于信号,且每次测量噪声也相互独立的假设,分析了在哈达码编码成像与经典成像中,噪声对图像信噪比的影响。推导出编码成像的信噪比提升与编码模板长度有关,为经典成像信噪比的倍。同时本发明采用分区编码的方式,减小了高分辨率图像的编码时间。实验结果表明,与经典成像方式相比,采用分区编码的哈达码变换成像方法明显的提高了图像的信噪比,同时可以在高分辨率图像条件下,缩短编码时间。
Description
技术领域
本发明属于弱信号探测技术领域,具体涉及一种基于DMD(数字微镜设备)编码的哈达码变换高灵敏成像的方法。
背景技术
微弱光电子成像极大的扩展了人眼的视觉范围,对微弱光信号的探测一直是天文观察、导弹预警(或告警)、卫星侦察、分子生物学、核物理等领域研究热点。
研究当前的微弱光成像方式主要是利用探测器单元对场景进行离散采样,在一定的积分时间内将光信号转换为电信号,单个探测器单元只对单一的场景单元探测。研究的重点与热点均在探测器光敏材料上,通过提高探测像元的灵敏度,降低暗电噪声来增强成像质量,但是这样的研制周期长,同时还有多方面的限制,在当前探测器探测水平的基础上研究一种新的成像方式有着积极现实的意义。
文献一(HadamardTransformOptics[M].MartinHarwit,NeilJ.A.Sloane)系统的阐述了HadamardTransform理论,并严格证明了该理论对弱信号探测的可行性。
文献二(A.Wuttig,R.Riesenberg.SensitiveHadamardtransformimagingspectrometerwithasimpleMEMS.ProceedingsofSPIE,Vol.48812003)是HadamardTransform理论在光谱仪领域应用的成功范例,也是该理论应用最多的地方。实践证明HadamardTransform理论的多通道探测优势,极大的提高了对弱信号探测的信噪比。但是该理论在直视成像领域的应用还未有具体的操作方法。
发明内容
本发明提出一种基于DMD编码的哈达码变换高灵敏成像的方法,目的在于弥补因探测器性能的限制,而无法实现对弱信号的有效探测。
实现本发明目的的技术方案为:一种基于DMD编码的哈达码变换高灵敏成像的方法,包括以下步骤:
步骤一:确定编码长度l值,哈达码一维编码模板长度l=4m+3,m=0,1,2……,其中m为自然数;
步骤二:确定单一二维编码模板分辨率l+1=i×j,其中i,j为编码模板的行和列;
步骤三:确定编码图像分辨率M×N=k2×(i×j),k=1,2,3……,其中M=k×i,N=k×j;
步骤四:采集编码图像光信号,左移步骤一中的一维编码模板,重复l次步骤二、三、四;
步骤五:对步骤四中采集的l个编码图像光信号进行哈达码逆变换处理,重建场景图像。
与现有技术相比,其显著优点:(1)可获取高信噪比图像;(2)解决了探测器非均匀性问题。
附图说明
图1是本发明方法系统框图。
图2是本发明方法中系统流程图。
图3是本发明方法中数字微镜原理图。
图4是本发明方法中DMD编码示意图。
图5是本发明方法中编码成像示意图。
图6是本发明方法中微光场景原图和常规直视成像图。
图7是本发明方法中微光场景图像不同编码复原对比图。
图8是本发明方法中红外场景原图和常规直视成像图。
图9是本发明方法中红外场景图像不同编码复原对比图。
图10是本发明方法中CT场景原图和常规直视成像图。
图11是本发明方法中CT场景图像不同编码复原对比图。
具体实施方式
本发明提出一种新的成像方式,利用哈达码变换具有多通道探测的优势,以DMD实现对空间像元的调制,采集相应数量的编码图像,最后重建出原始图像。同时采用编码的方式实现了采集时间与编码图像分辨率的最优化。利用这一成像方法有效的降低了对探测器的要求,相对的提高了探测灵敏度。
以下结合附图和具体实施例对本发明作进一步的详细说明。
哈达码变换(HadamardTransformHT)是统计学中的称量设计在光学中的应用。n个物体,分组称量所得各物体的重量,比一个一个单独称出的重量要准确。因此,如采用n个HT模板对试样信号进行调制,可得到个调制的信号,用检测器检测每一个调制信号的量值,n次测量后则可以通过HT把n次测得的调制信号还原成试样的信号。在常规测量中,检测器在每一时间间隔里只检测一个分辨单元的信号强度,而哈达码变换多通道检测技术在同一时间里却可以同时检测多个分辨单元里组合信号的总强度。在相同的实验条件下,经哈达码变换后,信噪比可提高倍。
结合图1,本发明提出了一种基于DMD编码的哈达码变换高灵敏成像的方法。该方法以HadamardTransform理论为基础。具体的实施方式如下:首先利用望远物镜将场景会聚到DMD(数字微镜),DMD在编码矩阵的驱动下对场景进行编码,准直物镜将编码后的反射光线准直,最后利用高倍率的成像物镜获取编码图像,采集相应数量的编码图像,利用Hadamard理论图像恢复算法复原场景图像。
DMD控制原理:DMD控制的目的是增大编码图像的分辨率的同时缩短成像时间。哈达码变换(HT)实现对图像的编码采集,首先将一维变换模板折叠成二维矩阵模板,折叠矩阵的大小决定了图像编码模板的大小,DMD像元尺寸是一定的,因此折叠矩阵的大小同时制约着编码图像的分辨率。以1023的编码长度为例,那么它折叠形成的二维矩阵为32×32(多出来的一个像元以0编码代替),实际上就是不编码。
如果这样的分辨率在某些场合不满足实际使用要求,本发明提供了一种控制方案,控制的过程如图3。首先根据使用要求和场合选择合适长度的编码,按照编码的长度设定的合适的长宽比,将一维编码折叠成二维矩阵编码。哈达码变化(HT)一维编码长度为素数,需要添加一位编码,使其可以实现二维矩阵。添加的编码在图像重建时,可利用相邻像素值求平均代替,其结果不影响整体图像质量。利用设定好的二维矩阵,首尾交替编码,就可以实现高分辨率的图像编码。
由图3可知折叠后的二维编码模板分辨率为Φ=[M×m,N×n],(m,n)为单个模板的分辨率。望远物镜将场景投影在DMD编码区域,编码模板中白色像元表示状态“1”,将光线反射进入成像光路。由于采用多个模板组合的方式,因此需要对DMD反射的光线进行准直,利用高倍率的成像物镜或者是低分辨率大像元的探测器(目的都是为了实现编码像元的会聚),充分实现多通道测量的优势。按照一维编码模板的长度,采集对应数量的编码图像,最后重建图像。
采用控制的多编码模板成像有别于单一模板的哈达码成像。主要的区别在于信噪比提高的计算上。采用单一模板的成像,编码图像由成像透镜汇聚成一点投影在探测器上,由哈达码理论知识可得,一维模板折叠前长度为n,则信噪比可提高倍。采用多模板成像时,探测器获取的图像不能为单点,应该是每一个模板对应一个点,否则将无法重建图像。图3中Φ=[M×m,N×n],理论上应该获取Θ=[M,N]个点,这时候也就是编码模板的低分辨率表示,这时候信噪比提高依然是理论值
图像重建步骤:
根据DMD控制成像过程如图4,建立图像复原模型。利用探测器获取L(一维模板的长度)幅编码图像,取出每幅图像中对应的像素y1(i,j),y2(i,j),…yl(i,j),令y(i,j)=[y1(i,j),y2(i,j)…yl(i,j)],哈达码变化矩阵为Ψl*l,求逆可得x(i,j)=Ψ-1×yT,为图3中单一模板所对应的编码区域,以此类推可以恢复原图像。由此可见本发明方法是利用增加采样时间和重构端的计算复杂度来实现高灵敏的成像。其效果可参见图5-11,实验结果表明,与经典成像方式相比,采用分区编码的哈达码变换成像方法明显的提高了图像的信噪比,同时可以在高分辨率图像条件下,缩短编码时间。
Claims (4)
1.一种基于DMD编码的哈达码变换高灵敏成像的方法,其特征在于包括以下步骤:
步骤一:确定编码长度l值,哈达码一维编码模板长度l=4m+3,m=0,1,2……,其中m为自然数;
步骤二:确定单一二维编码模板分辨率l+1=i×j,其中i,j为编码模板的行和列;步骤三:确定编码图像分辨率M×N=k2×(i×j),k=1,2,3……,其中M=k×i,N=k×j;
步骤四:采集编码图像光信号,左移步骤一中的一维编码模板,重复l次步骤二、三、四;
步骤五:对步骤四中采集的l个编码图像光信号进行哈达码逆变换处理,重建场景图像。
2.如权利要求1中所述的基于DMD编码的哈达码变换高灵敏成像的方法,其特征在于:步骤一中所述确定编码长度l值时,要求编码重建后的图像信噪比为高信噪比。
3.如权利1中所述的基于DMD编码的哈达码变换高灵敏成像的方法,其特征在于:步骤二中所述确定单一二维编码模板分辨率的具体的方法为:将一维编码模板长度加1,再将加1后的一维编码模板折叠成二维的矩阵。
4.如权利要求1中所述的基于DMD编码的哈达码变换高灵敏成像的方法,其特征在于:步骤三中所述确定编码图像分辨率的具体方法为:对步骤二中的单一二维编码模板进行拼接,将k2个单一二维编码模板拼接成一个k×k的编码模板,k值决定编码图像分辨率的大小。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410779177.6A CN105761212A (zh) | 2014-12-15 | 2014-12-15 | 基于dmd编码的哈达码变换高灵敏成像的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410779177.6A CN105761212A (zh) | 2014-12-15 | 2014-12-15 | 基于dmd编码的哈达码变换高灵敏成像的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105761212A true CN105761212A (zh) | 2016-07-13 |
Family
ID=56336866
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410779177.6A Pending CN105761212A (zh) | 2014-12-15 | 2014-12-15 | 基于dmd编码的哈达码变换高灵敏成像的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105761212A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108024037A (zh) * | 2017-11-28 | 2018-05-11 | 华中科技大学 | Hadamard矩阵感知成像系统及其成像方法 |
CN111130556A (zh) * | 2019-12-25 | 2020-05-08 | 北京航天控制仪器研究所 | 一种对Hadamard测量基进行排序的方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1999049412A1 (en) * | 1998-03-20 | 1999-09-30 | University Of Maryland | Method and apparatus for compressing and decompressing images |
CN101782430A (zh) * | 2010-04-12 | 2010-07-21 | 中国科学院西安光学精密机械研究所 | 一种基于哈达玛变换成像光谱仪的光谱复原方法 |
-
2014
- 2014-12-15 CN CN201410779177.6A patent/CN105761212A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1999049412A1 (en) * | 1998-03-20 | 1999-09-30 | University Of Maryland | Method and apparatus for compressing and decompressing images |
CN101782430A (zh) * | 2010-04-12 | 2010-07-21 | 中国科学院西安光学精密机械研究所 | 一种基于哈达玛变换成像光谱仪的光谱复原方法 |
Non-Patent Citations (1)
Title |
---|
贾辉: "星载太阳紫外光谱监视器辐射定标的研究", 《中国博士学位论文全文数据库工程科技Ⅱ辑》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108024037A (zh) * | 2017-11-28 | 2018-05-11 | 华中科技大学 | Hadamard矩阵感知成像系统及其成像方法 |
CN108024037B (zh) * | 2017-11-28 | 2020-12-18 | 华中科技大学 | Hadamard矩阵感知成像系统及其成像方法 |
CN111130556A (zh) * | 2019-12-25 | 2020-05-08 | 北京航天控制仪器研究所 | 一种对Hadamard测量基进行排序的方法 |
CN111130556B (zh) * | 2019-12-25 | 2023-11-10 | 北京航天控制仪器研究所 | 一种对Hadamard测量基进行排序的方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9900562B2 (en) | System and method for light-field imaging | |
CN104121990A (zh) | 基于随机光栅的压缩感知宽波段高光谱成像系统 | |
CN105589210B (zh) | 一种基于光瞳调制的数字化合成孔径成像方法 | |
JP2013535931A (ja) | 圧縮撮像装置の画像取得時間の減少 | |
CN102759408A (zh) | 一种单光子计数成像系统及其方法 | |
CN111344533A (zh) | 用于检测光的装置、系统和方法 | |
CN114721006B (zh) | 高精度单像素成像方法及系统 | |
US10446600B2 (en) | Imaging system and imaging device having a random optical filter array | |
CN116245726A (zh) | 基于深度学习框架的压缩感知偏振超分辨成像方法 | |
CN105120141B (zh) | 一种压缩感知光电成像方法及装置 | |
CN103743484B (zh) | 一种频谱编码成像装置 | |
CN105761212A (zh) | 基于dmd编码的哈达码变换高灵敏成像的方法 | |
CN116128931B (zh) | 基于掩模版空间调控的全谱段快速定位系统及方法 | |
CN109781260B (zh) | 超紧凑型快照式偏振光谱成像探测装置及探测方法 | |
CN110736539B (zh) | 一种基于压缩感知的凝视型光谱成像系统 | |
CN103558160A (zh) | 一种提高光谱成像空间分辨率的方法和系统 | |
CN105791643A (zh) | 基于dmd的多元探测成像系统 | |
KR101701138B1 (ko) | 해상도 및 초점 향상 | |
US12098954B2 (en) | Optical device and method | |
CN115278247A (zh) | 一种偏振光谱视频压缩采集系统 | |
Leblanc et al. | Interferometric lensless imaging: rank-one projections of image frequencies with speckle illuminations | |
CN110501069B (zh) | 时空联合调制光场光谱成像系统及方法 | |
CN103929577B (zh) | 基于压缩感知的紫外光与红外光成像系统和方法 | |
CN111854945A (zh) | 一种单像素紫外偏振成像方法及系统 | |
CN104992456B (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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20160713 |