面向小动物分子成像的时域荧光扩散层析系统

OPT OpenIR > 半导体激光器专利数据库

	面向小动物分子成像的时域荧光扩散层析系统
其他题名	面向小动物分子成像的时域荧光扩散层析系统
	高峰; 马艺闻; 赵会娟; 张丽敏; 和慧园
	2011-06-01
专利权人	天津大学
公开日期	2011-06-01
授权国家	中国
专利类型	授权发明
摘要	本发明涉及小动物分子成像领域，特别涉及面向小动物分子成像的时域荧光扩散层析系统。为提供能获得高定量精度和空间分辨率的三维层析、以及稳健的双参数和多组分成像能力的方法和装置，本发明采用的技术方案是：提供所需波长超短激光的两个皮秒半导体激光器及荧光探针；用于选择前述两个皮秒半导体激光器出射的两路激光的一个2∶1光开关；用于将2∶1光开关出射激光强度衰减到所需测量范围的可变衰减器；用于将可变衰减器出射的所选波长光源切换成多个光源点的光开关；用于将光源点出射激光投射到目标体的入射光纤；用于接收来自目标体的反射或透射激光的接收光纤；本发明主要用于分子成像的时域荧光扩散层析。
其他摘要	本发明涉及小动物分子成像领域，特别涉及面向小动物分子成像的时域荧光扩散层析系统。为提供能获得高定量精度和空间分辨率的三维层析、以及稳健的双参数和多组分成像能力的方法和装置，本发明采用的技术方案是：提供所需波长超短激光的两个皮秒半导体激光器及荧光探针；用于选择前述两个皮秒半导体激光器出射的两路激光的一个2∶1光开关；用于将2∶1光开关出射激光强度衰减到所需测量范围的可变衰减器；用于将可变衰减器出射的所选波长光源切换成多个光源点的光开关；用于将光源点出射激光投射到目标体的入射光纤；用于接收来自目标体的反射或透射激光的接收光纤；本发明主要用于分子成像的时域荧光扩散层析。
主权项	一种面向小动物分子成像的时域荧光扩散层析系统，其特征在于，包括：提供所需波长超短激光的两个皮秒半导体激光器及荧光探针；用于选择前述两个皮秒半导体激光器出射的两路激光的一个2∶1光开关；用于将2∶1光开关出射激光强度衰减到所需测量范围的可变衰减器；用于将可变衰减器出射的所选波长光源切换成多个光源点的光开关；用于将光源点出射激光投射到目标体的入射光纤；用于接收来自目标体的反射或透射激光的接收光纤；用于接收来自接收光纤的检验模块，检验模块进一步依次包括：滤除相应激发光的带通滤波器、PMT光电倍增管检测器、TCSPC单光子计数模块；用于由检验模块的输出产生测量平面或三维荧光发射率和寿命的图像的荧光分子层析时域FMT图像重建模块，FMT图像重建模块进一步包括下列模块： ①对检测量Γx(ζd，ξs，t)和Γm(ζd，ξs，t)进行拉普拉斯变换，得到Γx(ζd，ξs，p)和Γm(ζd，ξs，p)的模块； ②求出荧光波长检测量相对于激发波长检测量的玻恩比(Bom ratio)：的模块； ③在Robin边界条件下，对激发光波长的拉普拉斯变换时域扩散方程求解，获得ξs处激励下r和ζd处之光密度拉普拉斯变换值Φx(r，ξs，p)和Φx(ζd，ξs，p)的模块，即求解 [▿·Dx(r)▿-μax(r)c-p]Φx(r,ξs,p)=-δ(r-ξs)cΦx(r,ξs,p)+2KDx(r)n·▿Φx(r,ξs,p)\|r∈∂Ω=0Φx(ζd,ξs,p)=Φx(r,ξs,p)\|r=ζd 其中μax和Dx为激发波长下的吸收系数和扩散系数；为成像域Ω之边界；n为边界检测点之外法向；K为内反射控制因子；c为组织中的光速； ④在Robin边界条件下，获得荧光波长下的拉普拉斯变换时域扩散方程之格林函数Gm(ζd，r，p)的模块；即求解 [▿·Dm(r′)▿-μam(r′)c-p]Gm′(r′,r,p)=-δ(r′-r)cGm′(r′,r,p)+2KDm(r′)n·▿Gm′(r′,r,p)\|r′∈∂Ω=0Gm(ζd,r,p)=(cK/2)Gm′(r′,r,p)\|r′=ζd 其中r和r′为两个独立的空间矢量，μam和Dm为荧光波长下的吸收系数和扩散系数，为r处激励，r′处出射的光密度格林函数拉普拉斯变换值； ⑤由荧光扩散方程 [▿·Dm(r)▿-μam(r)c-p]Φm(r,rs,p)=-cΦx(r,rs,p)x(r,p)x(r,p)=ημaf(r)/[1+pτ(r)] 得出荧光扩散成像积分方程并进行空间离散化获得成像矩阵方程Γm(p)＝W(p)x(p)的模块，其中x(p)＝[x1(p)，x2(p)，…，xN(p)]T；Γm(p)＝[Γm(ζ1，ξ1，p)，Γm(ζ2，ξ1，p)，…，Γm(ζD，ξS，p)]T；W(p)为SD×N维矩阵；N为离散后之体元数； ⑥由任务②和⑤得Born比形式成像矩阵方程之模块，其中中元素的计算式为F(Γ)(ζd，ξs，p)为模型计算值； ⑦基于代数重建技术(Algebraic Reconstruction Technique，ART)的成像矩阵方程求解模块，即 xk+1(p)=xk(p)+λ[Γ^m(k+1)(p)-W(k+1)(p)xk(p)][W(k+1)(p)]·[W(k+1)(p)]T[W(k+1)(p)]Tk=0,1,2,...,(S×D-1) 其中为的第k个元素；W(k)(p)为W(p)的第k行；λ为迭代松弛因子； ⑦将拉普拉斯变换因子p分别取p1和p2，利用上述ART方法的荧光产率ημaf(r)和荧光寿命τ(r)成像计算模块，即 ημaf(r)=(p1-p2)x(r,p1)x(r,p2)/[p1x(r,p1)-p2x(r,p2)]τ(r)=-[x(r,p1)-x(r,p2)]/[p1x(r,p1)-p2x(r,p2)] 上述公式中r和r′为两个独立的空间矢量，ξs代表光源的位置，ζd为探测器的位置，Γx(ζd，ξs，t)为光源ξs激励，ζd点探测时激发光的检测光流量，Γm(ζd，ξs，t)为光源ξs激励，ζd点探测时出射光的检测光流量；Φx(r，ξs，p)为在外推边界条件下光源ξs激励，r位置处的激发光的光密度；Φx(ζd，ξs，p)为在外推边界条件下光源ξs激励，ζd位置处探测到的激发光的光密度；Gm(ζd，r，p)为r位置激励时ζd位置处探测到的出射光光流量的格林函数解；ημaf(r)表示荧光产率；η为量子效率；μaf(r)为荧光吸收系数；τ(r)为荧光寿命，δ(r，rs)为rs处的弥向点光源，Φx(r，rs，p)为rs处光源激励，r处激发光的光密度，Φm(r，rs，p)为rs处光源激励，r处荧光的光密度。
申请日期	2008-03-04
专利号	CN101236160B
专利状态	授权
申请号	CN200810052363.4
公开（公告）号	CN101236160B
IPC 分类号	G01N21/64 \| G06F17/10
专利代理人	刘国威
代理机构	天津市北洋有限责任专利代理事务所
文献类型	专利
条目标识符	http://ir.opt.ac.cn/handle/181661/45259
专题	半导体激光器专利数据库
作者单位	天津大学
推荐引用方式 GB/T 7714	高峰,马艺闻,赵会娟,等. 面向小动物分子成像的时域荧光扩散层析系统. CN101236160B[P]. 2011-06-01.

条目包含的文件
文件名称/大小	文献类型	版本类型	开放类型	使用许可
CN101236160B.PDF（1007KB）	专利		开放获取	CC BY-NC-SA	请求全文