背景技术
近红外光(600-1000nm)具备无损伤探测血红蛋白、脂肪、水等分子含量和流量的独特优势,目前利用近红外光进行成像的系统主要分为两种:第一种是血红蛋白含量成像仪,对含氧血红蛋白、无氧血红蛋白、水和脂肪的空间分布进行重建,这些信息可以用来对肿瘤或其他组织的新陈代谢进行估计,并可以用来追踪化疗的效果;第二种是血流成像仪,能够对体表或大脑表层血流分布进行测量,对局域血液循环不畅进行诊断,对大脑功能进行探测。血氧血流成像对于肿瘤治疗追踪有不可替代的优势,肿瘤区域一般都存在严重的缺氧和微血管增生现象,如肿瘤细胞被有效杀死,能够从血氧血流图像中看到明显改善。因为近红外光对人体无害,可以对放疗/化疗病人进行每日监察,分析治疗手段的有效性,及早终止无效治疗方案,尝试替代方案。
但是,无论是第一种还是第二种利用近红外光进行成像的系统,目前都不能同时进行深层血红蛋白分子含量和流动速率分布成像。
近红外血氧三维成像仪目前分为两种:一种利用CCD/APD/PMT等探测器,进行光波的强度测量。CCD的优势在于能够同时采集多个数据点,APD/PMT的优势在于能够进行光子计数,但是CCD的像素点很小,信号强度很弱,容易受到噪声影响,且由于近红外光容易扩散,本身对血氧含量的空间分辨能力很低,不需要探测器具有过于密集的数据采集点,因此CCD不是最优选择。传统APD/PMT价格昂贵,需要较高的工作电压(上千伏),体积较大,如需数百个探测点,需要昂贵的分光部件(Piezo-optical switch),牺牲时间分辨率换取空间分辨率。另一种血氧成像仪利用飞秒(fs)超短脉冲光和streak camera等仪器进行光脉冲展宽的测量,光脉冲展宽的程度和组织的吸收率,散射率有关,可以用来重建血氧的空间分布。但是,streak camera探测器非常昂贵(~$250,000),且利用PMT/APD作为单光子探测器依然需要分光部件和很长的扫描时间去获取数百个数据点,时间分辨率通常很低。
近红外血流探测仪目前主要集中在两个方向:第一是基于近红外光在组织表层反射时,击中运动粒子的多普勒效应,空间分辨率很高,例如激光光斑对比度血流成像仪(laser speckle contrast imaging),该种仪器仅能对表层组织的血流进行成像,因为近红外光在深层组织进行多次散射,多普勒模型失效。第二是基于近红外光的多次散射模型,对光子到达探测器时间进行自相关计算,根据不同探测器位置和脉冲到达时间的自相关程度,可以对不同深度的组织进行血流速率估计。但是,如果采用APD/PMT进行光子计数,则需要对探测点进行逐一扫描,每帧数据需要长达数秒的扫描时间,仅能对血流速率分布变化较慢的情况进行成像,比如,探测系统由10个APD组成,每个数据点采集积分0.5s,如需要测量400*(*数值取决于空间分辨率)个数据点,则至少需要约20秒时间来采集一帧数据。
随着探测器技术的进步,成本低廉的超快速单光子探测器阵列已经可以批量定制。譬如近期(2012/10)投入市场的SiPM,利用微加工制成每平方毫米数量上千的电子放大通道,每个通道相当于一个独立的APD。当光子击中探测器表面时,产生光电子,各个电子落入不同的放大通道,独立进行放大,能够达到106的增益。由于放大器尺寸微小,仅需30V左右的偏压。SiPM加工工艺具有均匀性,所以增益稳定,可以进行光子计数。通过连接快速泄流电阻,单光子电脉冲可以有非常窄的半高宽。
单光子探测器阵列与CCD相比的优势在于单个像素的独立性,互相之间干扰很小,不存在单个像素因为信号过高溢出影响周边像素的问题,如果有一个探测器损坏,可以很容易的替换。
单光子计数模式的另一个优势是不易受到背景噪声影响。一般的背景噪声会被计数的阈值过滤,在弱信号下依然有不错的信噪比。
由于采用低电压探测器,可以在探测模块端实现全部光电转换,采用柔性电缆FFC或者接插件进行供电和数据传输。现有技术获取多个数据点有两种方式:
(1)采用光纤进行信号传输,采用移动APD/PMT方式进行数据获取,但是仪器由于高压电缆的存在变得笨重和脆弱,且光纤的抖动会造成非常大的测量误差。
(2)采用APD/PMT和分光器件来获取多个数据点,取消了扫描,但是增加了光损失和仪器成本。
具体实施方式
下面详细描述本发明的实施例,所述实施例的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的,仅用于解释本发明,而不能理解为对本发明的限制。
在本发明的描述中,术语“内”、“外”、“纵向”、“横向”、“上”、“下”、“顶”、“底”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明而不是要求本发明必须以特定的方位构造和操作,因此不能理解为对本发明的限制。
请参阅图1,图1所示为本发明一实施方式中基于单光子计数器的近红外三维动态成像仪系统的结构图。
在本实施方式中,基于单光子计数器的近红外三维动态成像仪系统包括光源10、光切换器20、多模光纤30、单模光纤40、探测器50以及控制器60。
在本实施方式中,光源10用于产生不同波长的长相关激光,包括760nm、780nm、830nm、850nm等波长的长相关激光。
光切换器20一端与所述光源10光学连接,另一端与多模光纤30光学连接,本实施例中,所述光切换器20包括波长切换器与位置切换器。
所述生物体位于所述多模光纤30与所述单模光纤40之间。
所述探测器50和所述单模光纤41光学连接,用于对所述生物体进行检测。
控制器60,与所述光源10、所述光切换器20、所述探测器50电性连接,用于根据工作模式控制所述光源10产生相应的长相关激光以及所述探测器50采集相应的数据。
在本实施方式中,所述工作模式包括:血氧探测与血流探测。
在本实施方式中,所述控制器60在所述工作模式为血氧探测时控制所述光源10产生多个不同波长的激光并进行脉冲调制后经过所述波长切换器,由所述位置切换器经过所述多模光纤30,轮流输送相应波长的激光到入射光源点,所述探测器50记录单个脉冲展宽与到达时间。
在本实施方式中,所述探测器50利用单模光纤40采集光斑,输出光子到单光子计数器,并将光子转换为电信号,数据处理单元5022计算脉冲展宽及其相应的到达时间。
所述控制器60根据所述脉冲展宽与到达时间利用近红外光时域扩散方程进行三维重建以得到血氧分布。
在本实施方式中,近红外光时域扩散方程采用国际上比较认同的时域光扩散方程,如以下方式所示:
其中,v为光在组织中的传播速率,
为光量,
为光源强度在空间和时间的分布,这些参数都可以通过测量得出。可以看出,上述方程中的仅有未知量为散射率μ
s和吸收率μ
a(式中μ
t=μ
a+μ
s),探测点的采集光子数为光通量
可以根据探测器50的构型,给定的边界条件和边界的值,采取近似解析计算,也可以采取有限元分析的方法进行更为精密的重建。
比如被测组织为乳腺,光源10和探测器50可以在乳腺周围呈环形排布,单模光纤40埋在光吸收材料中,边界光反射量为零。通过不同位置入射光,探测器50的阵列给出的光通量值,可以形成一整套足以重建μs和μa的数据。
根据不同波长吸收率与血红蛋白含量的关系,可以通过简单的解方程组获取血氧含量信息
其中,∈i(λ)为分子在某个波长λ的光吸收率,c为分子的含量,在近红外某个波长λ的总光吸收率为含氧血红蛋白、无氧血红蛋白、脂肪及水等分子的综合。
在本实施方式中,所述控制器60在所述工作模式为血流探测时控制所述光源产生直流激光并进行脉冲调制后得到单光子脉冲,经过所述波长切换器,由所述为位置切换器经过所述多模光纤30,轮流输送相应波长的激光到入射光源点,所述探测器50记录单光子到达时间。
所述控制器60根据所述单光子的到达时间计算自相关函数,利用近红外光自相关方程进行三维重建以得到血流分布。
在本实施方式中,通过光扩散自相关方程计算自相关函数,光扩散自相关方程为:
在公认的光扩散自相关方程中,μs和μa已经由时域光扩散方程得出,可以认为是已知量,G1(r,τ)为每个探测点的未经归一化自相关函数,控制器60根据光子到达时间进行计算,k0为入射光的波数,S(r)为光源分布,<Δr2(τ)>是方程中唯一的未知量,也是反应血流速度的参数。
未经归一化的自相关函数Gn和经过归一化的自相关函数定义gn如下:
其中,n(t)为时间t测量到的光子数。为了节约计算时间,可以选取有限数目的τ值及对应的Gn,来计算组织的血流。
在本实施方式,单个光子脉冲到达必须进行阈值判选,最大限度消除暗电流的影响。
在本实施方式中,自相关函数计算可以采用FPGA,可以灵活选择光子计数时间段,对于不同血流流速采取不同的采样频率,最大限度的利用动态范围,增加采样率。
在本实施方式中,血流的三维重建需要测量每个数据点的自相关函数。
请参阅图2,为图1所示基于单光子计数器的近红外三维动态成像仪系统中光切换器20与多模光纤30的连接关系图。
在本实施方式中,光切换器20可以把输入的任意数据切换到多模光纤30中的其中一个或几个管线输出,控制器60利用编码,可以控制光切换器20自由选择输出组合。
在本实施方式中,多模光纤可最大限度地提高传输功率及性噪比。
请参阅图3,为图1所示基于单光子计数器的近红外三维动态成像仪系统中探测器50的结构图。
在本实施方式中,所述探测器50由多个单光子计数单元502组成,其中,所述多个单光子计数单元502一端与单模光纤40光学连接,所述单光子计数单元502通过所述单模光纤40对所述生物体进行探测。
在本实施方式中,单光子计数单元502包括数据处理单元5022、单模光纤5024以及遮光单元5026。
在本实施方式中,单模光纤5024用于探测光斑;遮光单元5026的中间设有一小孔,用于固定单模光纤5024。
在本实施方式中,数据处理单元5022用于处理单模光纤5024所探测到的光斑,并得到脉冲展宽与到达时间或单光子到达时间。
请参阅图4,图4所示为本发明一实施方式中自相关函数重建选取时间数据点的曲线图,从图4中可以看出,时间窗越小,自相关函数值对于血流计算的贡献越大,随着时间窗的增大,自相关函数值对血流计算的贡献会越来越小,因此,对自相关函数贡献最大的部分选取4-5个数据点就足够。
请参阅图5a~图5c,图5a~图5c所示分别为三种不同的分级光子计数时间窗的示意图。
在本实施方式中,由于计数时间窗不能无限短,且由于其指数衰减的特性,自相关时间很长的时候,g值计算不需要过于频繁的计数,τ值较大的情况下,密集计数工作量不会带来g函数计算的精确度增加,所以设定数据采集板对光子进行分级计数。计算τ值较小的g函数时,采用小的时间窗,计算τ值较大的g函数时,采用大的时间窗。对于血流较快的测量,可以减少时间窗个数,因为τ很大的时候,自相关函数为0,计算没有任何意义。
时间窗的选择也与个体差异有关,血流快的个体可以采用较小的时间窗(图5b),增加计算精确度,血流慢的个体采用较大的时间窗(图5a),减少计算负担,增加精确度(单个光子窗的积分时间增加)。
经过数据处理后,叠加(co-registered)的三维重建图可以同时显示如下信息:血氧分布、脂肪分布、水分布和血氧动态。由于采用同样的探测几何结构,和几乎同时的数据采集,参数之间相互关联极好。
本发明提供的实施方式中的基于单光子计数器的近红外三维动态成像仪系统,通过控制器60根据工作模式控制所述光源10产生相应的长相关激光以及控制所述探测器50采集相应的数据,有效的提高了数据采集的精度,并降低设备的成本。
虽然本发明参照当前的较佳实施方式进行了描述,但本领域的技术人员应能理解,上述较佳实施方式仅用来说明本发明,并非用来限定本发明的保护范围,任何在本发明的精神和原则范围之内,所做的任何修饰、等效替换、改进等,均应包含在本发明的权利保护范围之内。