WO2012146156A1 - 一种单光子计数成像系统及其方法 - Google Patents

Description

一种单光子计数成像系统及其方法

技术领域

本发明涉及极弱光探测的技术领域， 特别涉及一种单光子计数成像系统及其方法， 它采用压縮传感理论和 DLP技术， 利用点探测器可以实现极弱光对象的高质量二维成 像。 背景技术

普通成像设备通过记录观察对象上某点的光强和位置获得图像。当观察对象的光强 衰减到一定程度， 达到单光子水平， 就变成了离散脉冲信号。 单光子是一种极弱光， 被 认为是光不可分割的最小能量单位， 是可以探测的极限。单光子探测技术应用于生物自 发光、 医疗诊断、 非破坏性物质分析、 天文观测、 光谱测量、 量子光学等领域， 并在其 中扮演着重要角色。 研究极弱光成像探测技术在这些领域应用的发展非常有意义。

光子计数成像就是一种极弱光探测技术，通常它通过记录成像位置的光子计数以及 探测到光子的概率，在数据处理端进行累计和融合获得一幅图像，其核心是面元探测器， 面元探测器规模 （阵列大小）、 灵敏度范围、 以及响应波段直接影响能否获得单光子水 平的图象获取质量。但是， 用于单光子探测水平的面元探测器不但价格昂贵， 只能够在 少数波段可以实现， 且面元探测器灵敏度低，存在技术不甚成熟与极弱光对象二维成像 的强烈需求之间的矛盾。

压縮传感理论 （CS理论） 是由 E.J.Canc½_S等人提出的， 它打破了传统的线性采样 模式， 表明可压縮信号的少量线性随机投影中包含足够的信息来重建原信号。

CS理论包括两部分： 压縮采样和稀疏重建。

压縮采样是被测信号由高维向低维映射的过程。假设 是被测数据， ^是 观测数据， eR^ix "是随机投影矩阵（ <<")， e e ^是测量噪声， 那么， 压縮采样过 程可以描述为 （1) 式：

γ = χ + β (χ)

如果 是变换域稀疏的， 即 ·9 = Ψ Ψ是稀疏变换矩阵， 那么 （1) 式变化为 （2) 式：

y (2)

随机投影矩阵 Φ， 也叫测量矩阵， 需要满足 RIP (Restricted Isometry Property):

(l- ≤|K (i J χ 其中， 定义为使所有 s稀疏的向量； c都满足此不等式的最小常数， 并且 < 1。

另外， Φ与 Ψ越不相关， 采样所需的测量数 越小， 所以， 一般情况下 Φ设计为 随机矩阵。

稀疏重建实际上是在已知观测数据 y和测量矩阵 φ的条件下求解 a) 式中的 ， 这是一个 ill-po 一般用最优化方法求解， 可描述为 （3) 式：

如果 是 题可以描述为 （4) 式:

(3)式和 （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 ) 假设 是被测数据， y ^e ^是观测数据， e R^ix "是随机投影矩阵

( k « n -_{) >} 是测量噪声， 那么， 压縮采样的过程可以描述为 （1 ) 式： 如果 是变换域稀疏的， 即 >9 = Ψ ， Ψ是稀疏变换矩阵， 那么 （1) 式变化为 （2) 式：

y = ¹ ₁9 + e ₍₂₎

式中， Ψ为小波变换矩阵， Φ为 Gaussian随机矩阵；

12) 假设测量数为 ， 二维图像的像素个数为 n， 则 （1) 式中的测量矩阵则为

Φ = {Φ₁,···,Φ_!,···,Φ,}_; ^是0>的第 行， 把 ^x ^的二维图像的列首尾相连， 化成

"xl的一维列向量， 对应 （₁₎ 式中的 ， 其中的每一个元素代表相应位置处的光子密 度； DMD控制系统具有同样的分辨率， 它的列首尾相连， 化成 ^lx"的一维行向量， 对 应测量矩阵 Φ中的一行， 其中的每个元素代表相应位置处光子透射到第二透镜 2 的概 率；

13) 假设测量周期为: Γ， 在这段时间内， DMD控制系统保持不变， 单光子计数器 探测到的光子数为 W， 那么 就相当于光子密度图像与 DMD 系统上的随机数

阵列的内积值， 对应于（1)式中观察向量 y的一个元素 式中， "、 ^'分 别是 ^Φ '和 的第 '个元素； 根据测量矩阵改变 DMD控制系统， 重复 次测量， 就可以 得到整个观测数据 y ;

14)所述的稀疏重建是在已知观测数据 y和测量矩阵 Φ的条件下求解（ 1 )式中的 一般用最优化方法求解， 可描述为 （3) 式：

. ( 1

mini

2 - - (3)

如果 是变换域稀疏的， 对应于 （2) 式的重建问题可以描述为 （4) 式:

1

mini

(3)式和 （4)式中， 第一项是最小二乘约束， 记为 ^(^； 第二项是对 稀疏度的 一种约束； 两项之和是目标函数， 记为 ^^)。 作为上述技术方案的进一步的改进， 所述的最优化方法是采用 IWT (Iterative Wavelets Thresholding)算法， 把本次迭代的估计值做 DWT变换， 对变换系数做阈值处 理， 再做 DWT反变换得到下一次迭代的估计值； 如果阈值处理函数描述为¹ ^'^v)=^O^max^l^_v'^Q} ， 那么该算法可以描述为：

4

其中， 《= α^Φ 。

作为上述技术方案的另一种改进，所述的方法在所述的第二透镜 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控制系统上的测量矩阵经过最优 化算法重建光子密度图像。其中，光衰减器的设计是为了防止被测光子密度过大和单光 子计数器的门控时间过长引起的饱和。

为了便于数学上的理解， 假设测量数为 ， ( 1 ) 式中的测量矩阵写为

Φ = {Φ₁ , · · · , Φ_! , · · ·, Φ, } _; ^是 的第 行。 我们把

^的二维图像的列首尾相连， 化成" ^{x l}的一维列向量， 对应 （1 ) 式中的 其中的每一个元素代表相应位置处的光 子密度； DMD控制系统具有同样的分辨率， 它的列首尾相连， 化成 ^{l x} "的一维行向量， 对应测量矩阵 Φ中的一行，其中的每个元素代表相应位置处光子透射到第二透镜 2的概 率。 假设测量周期为 ， 在这段时间内， DMD控制系统保持不变， 单光子计数器探测 到的光子数为 W， 那么 就相当于光子密度图像与 DMD控制系统上的随机数阵列 的内积值， 对应于 （1 ) 式中观察向量 y的一个元素 ^ ' ( ^φ"、 ^'分别是 ^φ'' 和 的第 ^ '个元素）。 根据测量矩阵改变 DMD控制系统， 重复 次测量， 就可以得到整 个观测数据 ， 在物理上实现 （1 ) 式的过程。

根据光子学的知识， 在一个元面积 ^内， 任意时刻在 点观察到一个光子的概率 p(r> 正比于该处光强。 我们用生物芯片的灰度图像模拟光子密度图像。

生物芯片是典型的极弱光源， 目前主要通过荧光标记的方法使其便于观察。实际上 生物都有自发光的特性， 并且自发光光谱包含很多重要的信息。采用光子计数成像技术 就可以直接观测。 实验中选择分辨率为 64x64的生物芯片图像， 灰度为 256级， 最高灰 度级对应光子数为 O x ¹⁰² —¹。在假设不知道原图像的情况下，采用 Gaussian矩阵进行 压縮采样， IWT稀疏重建算法进行图像重建， 得到图 2所示的结果， 其中， （a)是原始 光子密度图像； （b) 是一次测量中 DMD控制系统上的随机矩阵， 黑点代表 0， 白点代 表 1， 灰点代表中间值； （c) 是 IWT算法的重建图像； （d) 是 IWT算法的残差图像。

图 2 (c) 所示的重建图像与原始图像图 2 (a) 的相关系数 Cov = 0.9783， 信噪比 PSNR=23.95dB, 满足对生物芯片图像的观测需求。 实验结果表明， 本发明硬件可以实 现压縮采样， 软件算法可以实现图像重建。

最后所应说明的是， 以上实施例仅用以说明本发明的技术方案而非限制。尽管参照 实施例对本发明进行了详细说明， 本领域的普通技术人员应当理解，对本发明的技术方 案进行修改或者等同替换， 都不脱离本发明技术方案的精神和范围，其均应涵盖在本发 明的权利要求范围当中。

Claims

权 利 要 求

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所述的单光子计数成像方法， 其特征在于， 所述的方法包括步 骤如下： 所述的压縮采样， 是被测信号由高维向低维映射的过程：

H ) 假设 e T是被测数据， 是观测数据， 0)eR^ix"是随机投影矩阵 (k«n), ee ^是测量噪声， 那么， 压縮采样的过程可以描述为 （1) 式： 如果 是变换域稀疏的， 即 >9 = Ψ ， Ψ是稀疏变换矩阵， 那么 （1) 式变化为 （2) 式：

y = ¹ ₁9 + e ₍₂₎

式中， Ψ为小波变换矩阵， Φ为 Gaussian随机矩阵；

12) 假设测量数为 ， 二维图像的像素个数为 n， 则 （1) 式中的测量矩阵则为

Φ = {Φ₁,···,Φ_!,···,Φ,}₎ ^是 的第 行， 把 ^x ^的二维图像的列首尾相连， 化成

"xl的一维列向量， 对应 （₁₎ 式中的 ， 其中的每一个元素代表相应位置处的光子密 度； DMD控制系统具有同样的分辨率， 它的列首尾相连， 化成 ^lx"的一维行向量， 对 应测量矩阵 Φ中的一行， 其中的每个元素代表相应位置处光子透射到第二透镜 （2) 的 概率；

13) 假设测量周期为: Γ， 在这段时间内， DMD控制系统保持不变， 单光子计数器 探测到的光子数为 W， 那么 就相当于光子密度图像与 DMD 系统上的随机数 阵列的内积值， 对应于（1)式中观察向量 y的一个元素

式中， "、 ^'分 别是 ^Φ '和 的第 '个元素； 根据测量矩阵改变 DMD控制系统， 重复 次测量， 就可以 得到整个观测数据 y;

14)所述的稀疏重建是在已知观测数据 y和测量矩阵 Φ的条件下求解（ 1 )式中的 一般用最优化方法求解， 可

如果 是变换域稀疏的， 2) 式的重建问题可以描述为 （4) 式:

(3)式和 （4)式中， 第一项是最小二乘约束， 记为 ^(^； 第二项是对 稀疏度的 一种约束； 两项之和是目标函数， 记为^ ^)。

5、 根据权利要求 4所述的单光子计数成像方法， 其特征在于， 所述的最优化方法 是采用 IWT算法， 把本次迭代的估计值做 DWT变换， 对变换系数做阈值处理， 再做 DWT反变换得到下一次迭代的估计值； 如果阈值处理函数描述为¹ ^'^v)⁼^"(")^max^l^_v'^Q} ，那么，该算法可以描述为:

= ψ- ¹. s _丄 ν/( ))，丄

a ) a

其中， α

6、 根据权利要求 3或 4所述的单光子计数成像方法， 其特征在于， 所述的方法在 所述的第二透镜 （2) 聚焦后至单光子计数器的光路上， 还通过光衰减器将光衰减到单 光子探测器的工作范围。

10