CN102759408A - 一种单光子计数成像系统及其方法 - Google Patents
一种单光子计数成像系统及其方法 Download PDFInfo
- Publication number
- CN102759408A CN102759408A CN2011101035593A CN201110103559A CN102759408A CN 102759408 A CN102759408 A CN 102759408A CN 2011101035593 A CN2011101035593 A CN 2011101035593A CN 201110103559 A CN201110103559 A CN 201110103559A CN 102759408 A CN102759408 A CN 102759408A
- Authority
- CN
- China
- Prior art keywords
- mrow
- lens
- control system
- single photon
- photons
- 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
- 238000003384 imaging method Methods 0.000 title claims abstract description 33
- 238000000034 method Methods 0.000 title claims abstract description 26
- 239000011159 matrix material Substances 0.000 claims abstract description 35
- 238000005259 measurement Methods 0.000 claims abstract description 34
- 238000005070 sampling Methods 0.000 claims abstract description 19
- 238000012545 processing Methods 0.000 claims abstract description 14
- 238000005457 optimization Methods 0.000 claims abstract description 11
- 230000003287 optical effect Effects 0.000 claims description 15
- 238000001514 detection method Methods 0.000 claims description 14
- 230000006835 compression Effects 0.000 claims description 9
- 238000007906 compression Methods 0.000 claims description 9
- 230000008569 process Effects 0.000 claims description 9
- 238000013507 mapping Methods 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 2
- 230000003071 parasitic effect Effects 0.000 abstract 1
- 230000035945 sensitivity Effects 0.000 description 8
- 238000000018 DNA microarray Methods 0.000 description 4
- 230000006872 improvement Effects 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 2
- 230000002238 attenuated effect Effects 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 238000004020 luminiscence type Methods 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000007123 defense Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000007850 fluorescent dye Substances 0.000 description 1
- 238000001215 fluorescent labelling Methods 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 238000001748 luminescence spectrum Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000002269 spontaneous effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J1/00—Photometry, e.g. photographic exposure meter
- G01J1/42—Photometry, e.g. photographic exposure meter using electric radiation detectors
- G01J1/44—Electric circuits
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J1/00—Photometry, e.g. photographic exposure meter
- G01J1/02—Details
- G01J1/04—Optical or mechanical part supplementary adjustable parts
- G01J1/0407—Optical elements not provided otherwise, e.g. manifolds, windows, holograms, gratings
- G01J1/0414—Optical elements not provided otherwise, e.g. manifolds, windows, holograms, gratings using plane or convex mirrors, parallel phase plates, or plane beam-splitters
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J1/00—Photometry, e.g. photographic exposure meter
- G01J1/42—Photometry, e.g. photographic exposure meter using electric radiation detectors
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J1/00—Photometry, e.g. photographic exposure meter
- G01J1/42—Photometry, e.g. photographic exposure meter using electric radiation detectors
- G01J1/44—Electric circuits
- G01J2001/4413—Type
- G01J2001/442—Single-photon detection or photon counting
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)
- Photometry And Measurement Of Optical Pulse Characteristics (AREA)
Abstract
本发明涉及一种单光子计数成像系统及其方法,该系统包括:滤光片、第一透镜(1)、DMD控制系统、第二透镜(2)、单光子计数器和数据处理单元;其中,DMD结合第一透镜(1)和第二透镜(2),用于将二维图像数据转化为一维数据序列完成被测信号的压缩采样,极弱光通过滤光片滤除杂光,经第一透镜(1)在DMD控制系统处成像,并由DMD控制系统控制光子被反射到第二透镜(2)的概率,经过第二透镜(2)控制光子聚焦;数据处理单元结合单光子计数器完成稀疏重建,数据处理单元根据单光子计数器在一定时间内对光子进行计数折算成探测到光子数的概率作为测量值,和DMD控制系统上的测量矩阵经过最优化算法重建光子密度图像,解算出二维图像。
Description
技术领域
本发明涉及极弱光探测的技术领域,特别涉及一种单光子计数成像系统及其方法,它采用压缩传感理论和DLP技术,利用点探测器可以实现极弱光对象的高质量二维成像。
背景技术
普通成像设备通过记录观察对象上某点的光强和位置获得图像。当观察对象的光强衰减到一定程度,达到单光子水平,就变成了离散脉冲信号。单光子是一种极弱光,被认为是光不可分割的最小能量单位,是可以探测的极限。单光子探测技术应用于生物自发光、医疗诊断、非破坏性物质分析、天文观测、光谱测量、量子光学等领域,并在其中扮演着重要角色。研究极弱光成像探测技术在这些领域应用的发展非常有意义。
光子计数成像就是一种极弱光探测技术,通常它通过记录成像位置的光子计数以及探测到光子的概率,在数据处理端进行累计和融合获得一幅图像,其核心是面元探测器,面元探测器规模(阵列大小)、灵敏度范围、以及响应波段直接影响能否获得单光子水平的图象获取质量。但是,用于单光子探测水平的面元探测器不但价格昂贵,只能够在少数波段可以实现,且面元探测器灵敏度低,存在技术不甚成熟与极弱光对象二维成像的强烈需求之间的矛盾。
压缩传感理论(CS理论)是由E.J.Candes等人提出的,它打破了传统的线性采样模式,表明可压缩信号的少量线性随机投影中包含足够的信息来重建原信号。
CS理论包括两部分:压缩采样和稀疏重建。
压缩采样是被测信号由高维向低维映射的过程。假设x∈Rn是被测数据,y∈Rk是观测数据,Φ∈Rk×n是随机投影矩阵(k<<n),e∈Rk是测量噪声,那么,压缩采样过程可以描述为(1)式:
y=Φx+e (1)
如果x是变换域稀疏的,即θ=Ψx,Ψ是稀疏变换矩阵,那么(1)式变化为(2)式:
y=ΦΨ-1θ+e (2)
随机投影矩阵Φ,也叫测量矩阵,需要满足RIP(Restricted Isometry Property):
其中,δs定义为使所有s稀疏的向量x都满足此不等式的最小常数,并且δs<1。
另外,Φ与Ψ越不相关,采样所需的测量数k越小,所以,一般情况下Φ设计为随机矩阵。
稀疏重建实际上是在已知观测数据y和测量矩阵Φ的条件下求解(1)式中的x,这是一个ill-posed问题,一般用最优化方法求解,可描述为(3)式:
如果x是变换域稀疏的,对应于(2)式的重建问题可以描述为(4)式:
DLP技术是美国德州仪器公司(TI)提出的一项技术,它与数字视频或图形信号结合起来,其微镜和透镜系统可以将数字图像反射到屏幕或其它表面,其核心是DLP芯片——数字微镜器件(Digital Micro-mirror Device,简称DMD控制系统),这目前是世界上最精密的光开关。它包含一个多达200万个安装在铰链上的微镜的矩阵,每个微镜的大小小于人的头发丝的五分之一,每一个微镜都可以在一定角度范围内(通常-12°和+12°)摆动,如果把这两种状态记为0和1,那么,用脉宽调制波(PWM)来驱动微镜,使其在0和1之间高速抖动,就可以实现中间状态。DMD控制系统及其相关的精密的电子元件就是所谓的DLP技术,该技术有成熟产品,广泛用于投影仪器等产品中。
这种“先采样,后重建”的思想使得将二维信号转换为随时间分布的一维信号,并且由单个探测器采样成为可能。而点探测器无论在探测灵敏度、波长范围具有更宽的选择范围,成本优势明显,利用点探测器实现单光子记数成像成为未来单光子水平成像的重要发展趋势。
发明内容
本发明的目的在于,为解决目前面元探测器灵敏度低、技术不甚成熟与极弱光对象二维成像的强烈需求之间的矛盾,从而提供一种单光子计数成像系统及其方法,以压缩传感(Compressive Sensing,CS)理论为基础,采用DLP技术将图象信号随机变化变成随机的光强信号,再以单光子计数器为探测元件获得记数信号,实现了用点探测器对极弱光对象进行二维成像。
为实现上述目的,本发明提供了一种单光子计数成像系统,其特征在于,该单光子计数成像系统采用压缩传感理论和DLP技术,并以单光子计数器为探测元件,实现了单光子级别的极弱光对象的二维成像,所述的单光子计数成像系统包括:滤光片、第一透镜1、DMD控制系统、第二透镜2、单光子计数器和数据处理单元;其中,
DMD结合第一透镜1和第二透镜2,用于将二维图像数据转化为一维数据序列完成被测信号的压缩采样,极弱光通过滤光片滤除杂光,经第一透镜1在DMD控制系统处成像,并由DMD控制系统控制光子被反射到第二透镜2的概率,经过第二透镜2控制光子聚焦;
数据处理单元结合单光子计数器完成稀疏重建,数据处理单元根据单光子计数器在一定时间内对光子进行计数折算成探测到光子数的概率作为测量值,和DMD控制系统上的测量矩阵经过最优化算法重建光子密度图像,解算出二维图像。
作为上述技术方案的一种改进,所述的第二透镜2聚焦后至单光子计数器的光路上还设有光衰减器,用于将光衰减到单光子探测器的工作范围。该光衰减器的设计是为了防止被测光子密度过大和单光子计数器的门控时间过长引起的饱和。
为实现上述的另一发明目的,本发明还提供了一种单光子计数成像方法,该方法采用了压缩传感理论和DLP技术,并以单光子计数器为探测元件,实现了单光子级别的极弱光对象的二维成像,其步骤包括:
1)压缩采样的步骤;
所述的压缩采样是由DMD控制系统结合第一透镜1和第二透镜2,用于将二维图像数据转化为一维数据序列完成被测信号的压缩采样,
极弱光通过滤光片滤除杂光,经第一透镜1在DMD控制系统处成像,并由DMD控制系统控制光子被反射到第二透镜2的概率,经过第二透镜2控制光子聚焦;
2)稀疏重建的步骤;
所述的稀疏重建是由数据处理单元根据单光子计数器在一定时间内对光子进行计数折算成探测到光子数的概率作为测量值,和DMD控制系统上的测量矩阵经过最优化算法重建光子密度图像,解算出二维图像。
作为上述技术方案的一种改进,所述的方法包括步骤如下:
所述的压缩采样,是被测信号由高维向低维映射的过程:
11)假设x∈Rn是被测数据,y∈Rk是观测数据,Φ∈Rk×n是随机投影矩阵(k<<n),e∈Rk是测量噪声,那么,压缩采样的过程可以描述为(1)式:
y=Φx+e (1)
如果x是变换域稀疏的,即θ=Ψx,Ψ是稀疏变换矩阵,那么(1)式变化为(2)式:
y=ΦΨ-1θ+e (2)
式中,Ψ为小波变换矩阵,Φ为Gaussian随机矩阵;
12)假设测量数为k,二维图像的像素个数为n,则(1)式中的测量矩阵则为Φ={Φ1,…,Φi,…,Φk},Φi是Φ的第i行,把的二维图像的列首尾相连,化成n×1的一维列向量,对应(1)式中的x,其中的每一个元素代表相应位置处的光子密度;DMD控制系统具有同样的分辨率,它的列首尾相连,化成1×n的一维行向量,对应测量矩阵Φ中的一行,其中的每个元素代表相应位置处光子透射到第二透镜2的概率;
13)假设测量周期为T,在这段时间内,DMD控制系统保持不变,单光子计数器探测到的光子数为N,那么N/T就相当于光子密度图像与DMD控制系统上的随机数阵列的内积值,对应于(1)式中观察向量y的一个元素式中,Φi,j、xj分别是Φi和x的第j个元素;根据测量矩阵改变DMD控制系统,重复k次测量,就可以得到整个观测数据y;
14)所述的稀疏重建是在已知观测数据y和测量矩阵Φ的条件下求解(1)式中的x,一般用最优化方法求解,可描述为(3)式:
如果x是变换域稀疏的,对应于(2)式的重建问题可以描述为(4)式:
作为上述技术方案的进一步的改进,所述的最优化方法是采用IWT(IterativeWavelets Thresholding)算法,把本次迭代的估计值做DWT变换,对变换系数做阈值处理,再做DWT反变换得到下一次迭代的估计值;
如果阈值处理函数描述为S(u,v)=sign(u)max{|u|-v,0},那么该算法可以描述为:
其中,α=||ΦTΦ||2。
作为上述技术方案的另一种改进,所述的方法在所述的第二透镜2聚焦后至单光子计数器的光路上,还通过光衰减器将光衰减到单光子探测器的工作范围。
本发明的优点在于:本发明采用压缩传感理论(Compressive Sensing,简称CS)结合数字光处理(Digital Light Processing,简称DLP)技术的方案来解决利用点探测器实现高探测灵敏度的成像问题,其灵敏度可以达到单光子水平,分辨率与DMD控制系统直接相关,而DMD控制系统目前可以达到很高的分辨率。本发明可广泛应用于生物自发光检测、医疗诊断、非破坏性物质分析、天文观测、国防军事、光谱测量、量子电子学等领域。
本发明以压缩传感(Compressive Sensing,CS)理论为基础,以单光子点探测器为探测元件,用单光子计数器实现极弱光的二维成像,结构简单,灵敏度可以达到单光子水平,分辨率与DMD控制系统直接相关,而DMD控制系统目前可以达到很高的分辨率,解决了目前该领域中焦平面传感器灵敏度低、阵列规模小、探测波长范围相对单一与极弱光对象二维成像的强烈需求之间的矛盾。
附图说明
图1是本发明的单光子计数成像系统的结构示意图。
图2是对本发明的模拟实验结果,图2(a)是原始光子密度图像;图2(b)是一次测量中DMD控制系统上的随机矩阵,黑点代表0,白点代表1,灰点代表中间值;图2(c)是IWT算法的重建图像,与原始图像的相关系数Cov=0.9783,信噪比PSNR=23.95dB;图2(d)是IWT算法的残差图像。
具体实施方式
以下结合附图对本发明作进一步的详细说明。
如图1所示,观察对象发出的极弱光被滤光片滤除杂光,经第一透镜1在DMD控制系统处成像,由DMD控制系统控制光子被反射到第二透镜2的概率,经过第二透镜2的控制光子会聚点,光衰减器作用是当光过强时将光衰减到单光子探测器的工作范围,由单光子计数器在一定时间内对光子进行计数,该数值可以折算成探测到光子数的概率作为测量值。最后由数据处理单元根据测量值和DMD控制系统上的测量矩阵经过最优化算法重建光子密度图像。其中,光衰减器的设计是为了防止被测光子密度过大和单光子计数器的门控时间过长引起的饱和。
为了便于数学上的理解,假设测量数为k,(1)式中的测量矩阵写为Φ={Φ1,…,Φi,…,Φk},Φi是Φ的第i行。我们把的二维图像的列首尾相连,化成n×1的一维列向量,对应(1)式中的x,其中的每一个元素代表相应位置处的光子密度;DMD控制系统具有同样的分辨率,它的列首尾相连,化成1×n的一维行向量,对应测量矩阵Φ中的一行,其中的每个元素代表相应位置处光子透射到第二透镜2的概率。假设测量周期为T,在这段时间内,DMD控制系统保持不变,单光子计数器探测到的光子数为N,那么N/T就相当于光子密度图像与DMD控制系统上的随机数阵列的内积值,对应于(1)式中观察向量y的一个元素(Φi,j、xj分别是Φi和x的第j个元素)。根据测量矩阵改变DMD控制系统,重复k次测量,就可以得到整个观测数据y,在物理上实现(1)式的过程。
根据光子学的知识,在一个元面积dA内,任意时刻在r点观察到一个光子的概率p(r)dA正比于该处光强。我们用生物芯片的灰度图像模拟光子密度图像。
生物芯片是典型的极弱光源,目前主要通过荧光标记的方法使其便于观察。实际上生物都有自发光的特性,并且自发光光谱包含很多重要的信息。采用光子计数成像技术就可以直接观测。实验中选择分辨率为64×64的生物芯片图像,灰度为256级,最高灰度级对应光子数为4.0×102s-1。在假设不知道原图像的情况下,采用Gaussian矩阵进行压缩采样,IWT稀疏重建算法进行图像重建,得到图2所示的结果,其中,(a)是原始光子密度图像;(b)是一次测量中DMD控制系统上的随机矩阵,黑点代表0,白点代表1,灰点代表中间值;(c)是IWT算法的重建图像;(d)是IWT算法的残差图像。
图2(c)所示的重建图像与原始图像图2(a)的相关系数Cov=0.9783,信噪比PSNR=23.95dB,满足对生物芯片图像的观测需求。实验结果表明,本发明硬件可以实现压缩采样,软件算法可以实现图像重建。
最后所应说明的是,以上实施例仅用以说明本发明的技术方案而非限制。尽管参照实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,对本发明的技术方案进行修改或者等同替换,都不脱离本发明技术方案的精神和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (6)
1.一种单光子计数成像系统,其特征在于,该单光子计数成像系统采用压缩传感理论和DLP技术,并以单光子计数器为探测元件,实现了单光子级别的极弱光对象的二维成像,所述的单光子计数成像系统包括:滤光片、第一透镜(1)、DMD控制系统、第二透镜(2)、单光子计数器和数据处理单元;其中,
DMD结合第一透镜(1)和第二透镜(2),用于将二维图像数据转化为一维数据序列完成被测信号的压缩采样,极弱光通过滤光片滤除杂光,经第一透镜(1)在DMD控制系统处成像,并由DMD控制系统控制光子被反射到第二透镜(2)的概率,经过第二透镜(2)控制光子聚焦;
数据处理单元结合单光子计数器完成稀疏重建,数据处理单元根据单光子计数器在一定时间内对光子进行计数折算成探测到光子数的概率作为测量值,和DMD控制系统上的测量矩阵经过最优化算法重建光子密度图像,解算出二维图像。
2.根据权利要求1所述的单光子计数成像系统,其特征在于,所述的第二透镜(2)聚焦后至单光子计数器的光路上还设有光衰减器,用于将光衰减到单光子探测器的工作范围。
3.一种单光子计数成像方法,该方法采用了压缩传感理论和DLP技术,并以单光子计数器为探测元件,实现了单光子级别的极弱光对象的二维成像,其步骤包括:
1)压缩采样的步骤;
所述的压缩采样是由DMD控制系统结合第一透镜(1)和第二透镜(2),用于将二维图像数据转化为一维数据序列完成被测信号的压缩采样,
极弱光通过滤光片滤除杂光,经第一透镜(1)在DMD控制系统处成像,并由DMD控制系统控制光子被反射到第二透镜(2)的概率,经过第二透镜(2)控制光子聚焦;
2)稀疏重建的步骤;
所述的稀疏重建是由数据处理单元根据单光子计数器在一定时间内对光子进行计数折算成探测到光子数的概率作为测量值,和DMD控制系统上的测量矩阵经过最优化算法重建光子密度图像,解算出二维图像。
4.根据权利要求3所述的单光子计数成像方法,其特征在于,所述的方法包括步骤如下:
所述的压缩采样,是被测信号由高维向低维映射的过程:
11)假设x∈Rn是被测数据,y∈Rk是观测数据,Φ∈Rk×n是随机投影矩阵(k<<n),e∈Rk是测量噪声,那么,压缩采样的过程可以描述为(1)式:
y=Φx+e (1)
如果x是变换域稀疏的,即θ=Ψx,Ψ是稀疏变换矩阵,那么(1)式变化为(2)式:
y=ΦΨ-1θ+e (2)
式中,Ψ为小波变换矩阵,Φ为Gaussian随机矩阵;
12)假设测量数为k,二维图像的像素个数为n,则(1)式中的测量矩阵则为Φ={Φ1,…,Φi,…,Φk},Φi是Φ的第i行,把的二维图像的列首尾相连,化成n×1的一维列向量,对应(1)式中的x,其中的每一个元素代表相应位置处的光子密度;DMD控制系统具有同样的分辨率,它的列首尾相连,化成1×n的一维行向量,对应测量矩阵Φ中的一行,其中的每个元素代表相应位置处光子透射到第二透镜(2)的概率;
13)假设测量周期为T,在这段时间内,DMD控制系统保持不变,单光子计数器探测到的光子数为N,那么N/T就相当于光子密度图像与DMD控制系统上的随机数阵列的内积值,对应于(1)式中观察向量y的一个元素式中,Φi,j、xj分别是Φi和x的第j个元素;根据测量矩阵改变DMD控制系统,重复k次测量,就可以得到整个观测数据y;
14)所述的稀疏重建是在已知观测数据y和测量矩阵Φ的条件下求解(1)式中的x,一般用最优化方法求解,可描述为(3)式:
如果x是变换域稀疏的,对应于(2)式的重建问题可以描述为(4)式:
5.根据权利要求4所述的单光子计数成像方法,其特征在于,所述的最优化方法是采用1WT算法,把本次迭代的估计值做DWT变换,对变换系数做阈值处理,再做DWT反变换得到下一次迭代的估计值;
如果阈值处理函数描述为S(u,v)=sign(u)max{|u|-v,0},那么,该算法可以描述为:
其中,α=||ΦTΦ||2。
6.根据权利要求3或4所述的单光子计数成像方法,其特征在于,所述的方法在所述的第二透镜(2)聚焦后至单光子计数器的光路上,还通过光衰减器将光衰减到单光子探测器的工作范围。
Priority Applications (5)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110103559.3A CN102759408B (zh) | 2011-04-25 | 2011-04-25 | 一种单光子计数成像系统及其方法 |
JP2014506741A JP6211512B2 (ja) | 2011-04-25 | 2012-04-23 | 単一光子計数イメージングシステム及び方法 |
EP12776896.8A EP2685227B1 (en) | 2011-04-25 | 2012-04-23 | Single photon-counting imaging system and method thereof |
US14/000,421 US8723130B2 (en) | 2011-04-25 | 2012-04-23 | Single photon-counting imaging system and method thereof |
PCT/CN2012/074533 WO2012146156A1 (zh) | 2011-04-25 | 2012-04-23 | 一种单光子计数成像系统及其方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110103559.3A CN102759408B (zh) | 2011-04-25 | 2011-04-25 | 一种单光子计数成像系统及其方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102759408A true CN102759408A (zh) | 2012-10-31 |
CN102759408B CN102759408B (zh) | 2015-04-15 |
Family
ID=47053946
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110103559.3A Expired - Fee Related CN102759408B (zh) | 2011-04-25 | 2011-04-25 | 一种单光子计数成像系统及其方法 |
Country Status (5)
Country | Link |
---|---|
US (1) | US8723130B2 (zh) |
EP (1) | EP2685227B1 (zh) |
JP (1) | JP6211512B2 (zh) |
CN (1) | CN102759408B (zh) |
WO (1) | WO2012146156A1 (zh) |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103207015A (zh) * | 2013-04-16 | 2013-07-17 | 华东师范大学 | 一种光谱重构方法及其光谱仪装置 |
CN103645491A (zh) * | 2013-11-25 | 2014-03-19 | 中国科学院高能物理研究所 | 放射源定位方法、装置及系统 |
CN103969693A (zh) * | 2014-04-30 | 2014-08-06 | 中国科学院长春光学精密机械与物理研究所 | 光电探测成像系统及其成像方法 |
CN104796674A (zh) * | 2015-04-17 | 2015-07-22 | 南京理工大学 | 基于压缩感知的彩色成像装置及方法 |
CN105915869A (zh) * | 2016-04-22 | 2016-08-31 | 南京理工大学 | 一种彩色自适应压缩计算鬼成像系统及方法 |
CN105915868A (zh) * | 2016-04-22 | 2016-08-31 | 南京理工大学 | 基于扩展小波树的彩色成像系统及方法 |
CN106054038A (zh) * | 2016-06-29 | 2016-10-26 | 南京理工大学 | 一种连通域标记的紫外单光子计数的方法 |
CN106772430A (zh) * | 2016-12-30 | 2017-05-31 | 南京理工大学 | 基于多分辨率小波逼近的单像素光子计数三维成像系统及方法 |
CN106911893A (zh) * | 2017-02-23 | 2017-06-30 | 北京建筑大学 | 一种单像素计算成像方法 |
CN107078011A (zh) * | 2014-10-28 | 2017-08-18 | 株式会社日立高新技术 | 带电粒子束装置和信息处理装置 |
CN109361833A (zh) * | 2018-10-08 | 2019-02-19 | 南昌大学 | 一种单光子压缩视频的传输装置及方法 |
CN109477796A (zh) * | 2016-05-25 | 2019-03-15 | 徕卡显微系统复合显微镜有限公司 | 允许更高光强度的利用时间相关单光子计数的荧光寿命成像显微学方法 |
CN111103062A (zh) * | 2019-12-06 | 2020-05-05 | 太原理工大学 | 一种基于单光子计数二维成像的装置及方法 |
CN111208531A (zh) * | 2020-01-19 | 2020-05-29 | 中国科学技术大学 | 一种基于宽谱光源的单光子成像激光雷达系统 |
CN111896125A (zh) * | 2020-07-09 | 2020-11-06 | 武汉大学 | 一种单光子计数成像的偏振去噪方法 |
Families Citing this family (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2755327A4 (en) | 2011-10-25 | 2015-11-04 | Ct Space Sci & Applied Res Cas | MULTI-DIMENSIONAL MULTI-DIMENSIONAL IMAGING SYSTEM AND METHOD WITH TEMPORARY RESOLUTION OR ULTRA-LIGHT LIGHT |
WO2015143378A1 (en) * | 2014-03-20 | 2015-09-24 | Kla-Tencor Corporation | Compressive sensing with illumination patterning |
CN104238996B (zh) * | 2014-09-04 | 2017-08-11 | 清华大学 | 源无关量子随机数的产生方法及装置 |
CN104375158B (zh) * | 2014-09-10 | 2017-09-05 | 北京交通大学 | 一种单原子的探测方法及装置 |
TWI546533B (zh) * | 2015-07-14 | 2016-08-21 | 龍彩科技股份有限公司 | 即時空間與時間光譜量測系統及其量測模組 |
CN109816603B (zh) * | 2018-12-30 | 2023-04-11 | 天津大学 | 单光子计数成像的图像传感器图像还原方法 |
CN110187498B (zh) * | 2019-05-27 | 2021-08-17 | 中国科学院国家空间科学中心 | 一种真热光关联成像系统 |
CN110443883B (zh) * | 2019-07-08 | 2023-04-07 | 杭州电子科技大学 | 一种基于dropblock的单张彩色图片平面三维重建方法 |
CN110852981B (zh) * | 2019-11-19 | 2023-08-22 | 中智科仪(北京)科技有限公司 | 一种单光子成像图像处理方法 |
CN112957011B (zh) * | 2021-02-01 | 2022-06-24 | 西安电子科技大学 | 高灵敏度微弱荧光信号探测系统、方法、存储介质及应用 |
CN113534094B (zh) * | 2021-05-28 | 2023-04-18 | 西安电子科技大学 | 一种基于纠缠态的量子探测恒虚警检测系统及检测方法 |
CN114486746B (zh) * | 2021-11-25 | 2023-12-08 | 中国科学院西安光学精密机械研究所 | 基于压缩感知的高分辨率光子集成成像系统及成像方法 |
CN115442505B (zh) * | 2022-08-30 | 2023-07-21 | 山西大学 | 一种单光子压缩感知成像系统及其方法 |
CN117705728B (zh) * | 2023-12-12 | 2024-08-02 | 成都信息工程大学 | 基于光自旋的单光子计数手性分子检测装置及检测方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7197523B2 (en) * | 2001-05-09 | 2007-03-27 | Magiq Technologies, Inc. | Efficient use of detectors for random number generation |
US6700108B2 (en) * | 2001-11-01 | 2004-03-02 | Northrop Grumman Corporation | Means for protecting optical focal plane sensor arrays against excessive irradiation |
JP2005250714A (ja) * | 2004-03-03 | 2005-09-15 | Univ Nihon | 光子乱数発生器 |
GB2414881A (en) * | 2004-06-01 | 2005-12-07 | Imp College Innovations Ltd | Imaging system capable of reproducing a wide range of intensities |
WO2006116134A2 (en) * | 2005-04-21 | 2006-11-02 | William Marsh Rice University | Method and apparatus for compressive imaging device |
JP2011075491A (ja) * | 2009-10-01 | 2011-04-14 | Panasonic Corp | 光量測定装置 |
-
2011
- 2011-04-25 CN CN201110103559.3A patent/CN102759408B/zh not_active Expired - Fee Related
-
2012
- 2012-04-23 WO PCT/CN2012/074533 patent/WO2012146156A1/zh active Application Filing
- 2012-04-23 EP EP12776896.8A patent/EP2685227B1/en not_active Not-in-force
- 2012-04-23 US US14/000,421 patent/US8723130B2/en not_active Expired - Fee Related
- 2012-04-23 JP JP2014506741A patent/JP6211512B2/ja not_active Expired - Fee Related
Non-Patent Citations (3)
Title |
---|
MACRO F.DUARTE,ET AL: "SINGLE-PIXEL IMAGING VIA COMPRESSIVE SAMPLING", 《IEEE SIGNAL PROCESSING MAGZINE》 * |
杜克铭: "基于压缩传感的光子计数成像系统研究", 《2010国防空天信息技术前沿论坛论文集》 * |
高睿: "基于压缩传感的匹配追踪重建算法研究", 《北京交通大学硕士学位论文》 * |
Cited By (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103207015A (zh) * | 2013-04-16 | 2013-07-17 | 华东师范大学 | 一种光谱重构方法及其光谱仪装置 |
CN103645491B (zh) * | 2013-11-25 | 2016-12-07 | 中国科学院高能物理研究所 | 放射源定位方法、装置及系统 |
CN103645491A (zh) * | 2013-11-25 | 2014-03-19 | 中国科学院高能物理研究所 | 放射源定位方法、装置及系统 |
CN103969693A (zh) * | 2014-04-30 | 2014-08-06 | 中国科学院长春光学精密机械与物理研究所 | 光电探测成像系统及其成像方法 |
CN107078011B (zh) * | 2014-10-28 | 2018-11-02 | 株式会社日立高新技术 | 带电粒子束装置 |
CN107078011A (zh) * | 2014-10-28 | 2017-08-18 | 株式会社日立高新技术 | 带电粒子束装置和信息处理装置 |
CN104796674A (zh) * | 2015-04-17 | 2015-07-22 | 南京理工大学 | 基于压缩感知的彩色成像装置及方法 |
CN105915868A (zh) * | 2016-04-22 | 2016-08-31 | 南京理工大学 | 基于扩展小波树的彩色成像系统及方法 |
CN105915869A (zh) * | 2016-04-22 | 2016-08-31 | 南京理工大学 | 一种彩色自适应压缩计算鬼成像系统及方法 |
CN109477796A (zh) * | 2016-05-25 | 2019-03-15 | 徕卡显微系统复合显微镜有限公司 | 允许更高光强度的利用时间相关单光子计数的荧光寿命成像显微学方法 |
CN109477796B (zh) * | 2016-05-25 | 2022-05-31 | 徕卡显微系统复合显微镜有限公司 | 允许更高光强度的利用时间相关单光子计数的荧光寿命成像显微学方法 |
CN106054038A (zh) * | 2016-06-29 | 2016-10-26 | 南京理工大学 | 一种连通域标记的紫外单光子计数的方法 |
CN106054038B (zh) * | 2016-06-29 | 2019-01-15 | 南京理工大学 | 一种连通域标记的紫外单光子计数的方法 |
CN106772430A (zh) * | 2016-12-30 | 2017-05-31 | 南京理工大学 | 基于多分辨率小波逼近的单像素光子计数三维成像系统及方法 |
CN106911893A (zh) * | 2017-02-23 | 2017-06-30 | 北京建筑大学 | 一种单像素计算成像方法 |
CN109361833A (zh) * | 2018-10-08 | 2019-02-19 | 南昌大学 | 一种单光子压缩视频的传输装置及方法 |
CN109361833B (zh) * | 2018-10-08 | 2020-08-11 | 南昌大学 | 一种单光子压缩视频传输装置的传输方法 |
CN111103062A (zh) * | 2019-12-06 | 2020-05-05 | 太原理工大学 | 一种基于单光子计数二维成像的装置及方法 |
CN111208531A (zh) * | 2020-01-19 | 2020-05-29 | 中国科学技术大学 | 一种基于宽谱光源的单光子成像激光雷达系统 |
CN111896125A (zh) * | 2020-07-09 | 2020-11-06 | 武汉大学 | 一种单光子计数成像的偏振去噪方法 |
Also Published As
Publication number | Publication date |
---|---|
JP2014512540A (ja) | 2014-05-22 |
WO2012146156A1 (zh) | 2012-11-01 |
EP2685227B1 (en) | 2018-10-10 |
CN102759408B (zh) | 2015-04-15 |
US8723130B2 (en) | 2014-05-13 |
JP6211512B2 (ja) | 2017-10-11 |
EP2685227A1 (en) | 2014-01-15 |
EP2685227A4 (en) | 2014-10-01 |
US20130341487A1 (en) | 2013-12-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102759408B (zh) | 一种单光子计数成像系统及其方法 | |
CN102353449A (zh) | 一种极弱光多光谱成像方法及其系统 | |
Studer et al. | Compressive fluorescence microscopy for biological and hyperspectral imaging | |
CN109297925B (zh) | 一种基于分块压缩感知的太赫兹高分辨率快速成像装置 | |
Hojman et al. | Photoacoustic imaging beyond the acoustic diffraction-limit with dynamic speckle illumination and sparse joint support recovery | |
US11561134B2 (en) | Compressed-sensing ultrafast spectral photography systems and methods | |
US20120330162A1 (en) | Modulated aperture imaging for automatic moving target detection | |
CN110793632A (zh) | 一种用于火焰拍摄的高速高精度光谱视频系统及方法 | |
Yu et al. | Single-photon compressive imaging with some performance benefits over raster scanning | |
Edgar et al. | Real-time computational photon-counting LiDAR | |
CN103954357A (zh) | 压缩光谱成像系统测量矩阵的获取方法 | |
Ouyang et al. | Compressive line sensing imaging system in a controlled hybrid scattering environment | |
CN110632002A (zh) | 一种基于压缩感知的孔径编码光谱探测装置 | |
CN111123284B (zh) | 探测方法及探测设备 | |
Thapa et al. | Less is more: compressive sensing in optics and image science | |
Du Bosq et al. | An overview of joint activities on computational imaging and compressive sensing systems by NATO SET-232 | |
Leblanc et al. | Interferometric lensless imaging: rank-one projections of image frequencies with speckle illuminations | |
KR101766328B1 (ko) | 현미경 | |
CN115442505A (zh) | 一种单光子压缩感知成像系统及其方法 | |
Zhao et al. | Passive ghost imaging using caustics modeling | |
CN111640063A (zh) | 基于空间频域多尺度调制与重建的压缩成像系统及方法 | |
Cao et al. | Quantum projection ghost imaging: a photon-number-selection method | |
Wu | Compressive optical imaging systems | |
Huang et al. | Enhancing Compressive Single-Pixel Imaging with Zig-Zag-Ordered Walsh-Hadamard Light Modulation | |
Leblanc et al. | Interferometric Lensless Imaging |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CP01 | Change in the name or title of a patent holder |
Address after: 100190 No. two south of Zhongguancun, Haidian District, Beijing 1 Patentee after: NATIONAL SPACE SCIENCE CENTER, CAS Address before: 100190 No. two south of Zhongguancun, Haidian District, Beijing 1 Patentee before: NATIONAL SPACE SCIENCE CENTER, CHINESE ACADEMY OF SCIENCES |
|
CP01 | Change in the name or title of a patent holder | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150415 |
|
CF01 | Termination of patent right due to non-payment of annual fee |