CN103188989B - 一种动态荧光分子断层图像的重建方法 - Google Patents
一种动态荧光分子断层图像的重建方法 Download PDFInfo
- Publication number
- CN103188989B CN103188989B CN201280002807.8A CN201280002807A CN103188989B CN 103188989 B CN103188989 B CN 103188989B CN 201280002807 A CN201280002807 A CN 201280002807A CN 103188989 B CN103188989 B CN 103188989B
- Authority
- CN
- China
- Prior art keywords
- reconstructing
- fluorescence
- reconstruction
- matrix
- fmt
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 90
- 238000003325 tomography Methods 0.000 title abstract description 6
- 239000011159 matrix material Substances 0.000 claims abstract description 67
- 238000003384 imaging method Methods 0.000 claims abstract description 32
- 230000009466 transformation Effects 0.000 claims abstract description 19
- 230000003287 optical effect Effects 0.000 claims abstract description 7
- 238000005315 distribution function Methods 0.000 claims abstract description 5
- 238000006243 chemical reaction Methods 0.000 claims description 29
- 239000007850 fluorescent dye Substances 0.000 claims description 14
- 238000013398 bayesian method Methods 0.000 claims description 5
- 238000012360 testing method Methods 0.000 claims description 5
- 238000000354 decomposition reaction Methods 0.000 claims description 4
- 238000000799 fluorescence microscopy Methods 0.000 claims description 4
- 230000004060 metabolic process Effects 0.000 claims description 4
- 230000005284 excitation Effects 0.000 claims description 3
- 230000003595 spectral effect Effects 0.000 claims description 2
- 230000008569 process Effects 0.000 abstract description 8
- 230000002503 metabolic effect Effects 0.000 abstract description 5
- 239000000523 sample Substances 0.000 abstract 1
- 238000004613 tight binding model Methods 0.000 abstract 1
- 238000002474 experimental method Methods 0.000 description 12
- 238000010586 diagram Methods 0.000 description 8
- 230000000694 effects Effects 0.000 description 7
- 241000699670 Mus sp. Species 0.000 description 5
- 230000008859 change Effects 0.000 description 5
- 238000001506 fluorescence spectroscopy Methods 0.000 description 5
- 210000004185 liver Anatomy 0.000 description 5
- 230000002123 temporal effect Effects 0.000 description 4
- 230000006870 function Effects 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 2
- 241000699666 Mus <mouse, genus> Species 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- MOFVSTNWEDAEEK-UHFFFAOYSA-M indocyanine green Chemical compound [Na+].[O-]S(=O)(=O)CCCCN1C2=CC=C3C=CC=CC3=C2C(C)(C)C1=CC=CC=CC=CC1=[N+](CCCCS([O-])(=O)=O)C2=CC=C(C=CC=C3)C3=C2C1(C)C MOFVSTNWEDAEEK-UHFFFAOYSA-M 0.000 description 2
- 229960004657 indocyanine green Drugs 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000004140 cleaning Methods 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 238000013144 data compression Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 230000036267 drug metabolism Effects 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000001727 in vivo Methods 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
- 230000031068 symbiosis, encompassing mutualism through parasitism Effects 0.000 description 1
- 230000036962 time dependent Effects 0.000 description 1
- 210000003462 vein Anatomy 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/006—Inverse problem, transformation from projection-space into object-space, e.g. transform methods, back-projection, algebraic methods
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/62—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
- G01N21/63—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited
- G01N21/64—Fluorescence; Phosphorescence
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/412—Dynamic
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Algebra (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)
Abstract
本发明涉及一种动态荧光分子断层图像的重建方法,包括以下步骤:1)设置一小动物诱发荧光分子成像***;2)采用所述***,在不同时间点全角度、等间隔采集反映荧光探针在成像对象体内代谢分布的荧光投影图像,并将采集图像存储为输入矩阵u;3)根据u生成KL变换矩阵A;4)求解输入矩阵u经KL变换后的数据5)获取经KL变换后需要进行3D重建的L′个KL分量;6)根据光场分布函数Φx和格林函数G生成权重矩阵W1;7)根据公式采用3D FMT重建方法对每一个KL分量分别进行3D重建;8)对KL分量的重建结果进行逆KL变换,得到4D荧光断层序列。本发明可以广泛应用于动态荧光分子图像的重建过程中。
Description
技术领域
本发明涉及一种荧光分子断层图像的重建方法,特别是关于一种动态荧光分子断层图像的重建方法。
背景技术
目前动态(4D)荧光分子断层图像的重建大都是利用逐帧(frame-by-frame)重建方法完成,即对采集到的每一帧(圈)的荧光投影图像进行逐个重建。此方法的优点是可以直接利用已有的各种诱发荧光分子断层成像(Fluorescence MolecularTomography,FMT)重建方法完成图像重建。然而,在动态FMT成像研究中,例如在药物代谢的研究中,药物(荧光探针)的浓度是随时间变化的,因此所采集到的每一帧测量数据(荧光投影图像)在时间上是高度相关的,如果采用上述逐帧重建方法对荧光投影图像进行4D FMT重建,则忽略了重建过程中荧光探针浓度沿时间轴的相关性,但是如果直接将采集的所有荧光投影图像序列作为一个整体直接进行4D重建,其计算量通常已超出普通计算机的计算能力,例如对一个35帧的荧光投影图像序列(假设在每帧FMT成像过程中,等间隔采集24张不同角度的荧光投影图像)直接进行4D FMT重建,其用于重建的空间-时间权重矩阵将占用大约3234GB的内存,这大大超出了普通计算机的计算能力。
发明内容
针对上述问题,本发明的目的是提供一种充分考虑荧光投影图像序列的时间相关性,且能够有效提高动态荧光分子断层图像重建速度的动态荧光分子断层图像的重建方法。
为实现上述目的,本发明采取以下技术方案:一种动态荧光分子断层图像的重建方法,包括以下步骤:1)设置一包括有计算机、小动物旋转平台装置和荧光成像激发与检测装置的小动物诱发荧光分子成像***;2)采用所述小动物诱发荧光分子成像***,在不同时间点,全角度、等间隔采集反映荧光探针在成像对象体内的代谢分布的荧光投影图像,并将每一时间点所采集的不同角度的荧光投影图像存储为一矩阵序列,将所述矩阵序列标记为输入矩阵u;3)根据输入矩阵u生成KL变换矩阵A;4)基于求解输入矩阵u经KL变换后的数据式中,IM为M×M的单位矩阵,为克罗内克积;5)获取经KL变换后需要进行3D重建的L′个KL分量;6)根据光场分布函数Φx和格林函数G生成权重矩阵W1,其中,Φx描述了对应激发光波长的光子密度,G描述了在荧光谱段,光子在成像对象体内的传播;7)根据公式采用3D FMT重建方法在KL域中对步骤5)中获取的每一个KL分量分别进行3D重建,其中,k=1,2,...,L′,W-1 1为权重矩阵W1的逆矩阵;8)对上述步骤7)中获取的KL分量的重建结果进行逆KL变换,得到完整的4D荧光断层图像序列。
所述步骤3)根据输入矩阵u生成KL变换矩阵A的方法为:①构建输入矩阵u的时间协方差矩阵Ptime;②通过求解PtimeAT=ATD,生成KL变换矩阵A,式中dl是协方差矩阵Ptime的第l个特征值,L是动态成像中扫描的总帧数。
所述步骤3)根据输入矩阵u生成KL变换矩阵A的方法为:通过对输入矩阵u进行奇异值分解求得。
所述步骤5)中获取经KL变换后需要进行3D重建的L′个KL分量采用scree方法、Kaiser方法和基于特征值的累计方差方法中的一种。
所述步骤6)中在检测点r处获取的荧光投影图像采用下式表示:
Φm(r)=∫r′∈VG(r,r′)n(r′)Φx(r′)dr′
其中,W1(r)=G(r,r′)Φx(r′)ΔV(r),ΔV是重建区域V被离散成的小网格;V是重建区域;Φm是发射荧光波长的光子密度;n(r′)为重建区域V中点r′处的待重建的荧光产额。
所述步骤7)采用3D FMT重建方法在KL域中对步骤5)中获取的每一个KL分量分别进行3D重建中,3D FMT重建方法采用ART方法、LSQR方法和Bayesian方法中的一种。
本发明由于采取以上技术方案,其具有以下优点:1、本发明在不同时间点,全角度、连续采集荧光探针在成像对象体内的代谢分布,且将采集的所有荧光投影图像作为一个整体进行重建,因此在重建过程中充分考虑到了荧光投影图像在时间上的相关性,弥补了现有重建方法对时间相关性考虑的不足。2、本发明通过对测量得到的荧光投影数据进行KL变换,去掉了荧光投影数据间的时间相关性,将复杂的4D FMT重建问题简化为一系列KL分量在KL域中的3D重建问题,因此与直接对测量数据进行4D重建的方法相比降低了4D FMT重建的计算规模。3、本发明采用的KL变换具有很高的数据压缩性,将荧光投影数据进行KL变换后,有用的信号被保存在前几个KL分量中,采用现有的scree方法或其它类似方法提取需要进行FMT重建的KL分量的个数,并对提取的每一KL分量分别采用现有的3D FMT重建方法进行重建,然后将重建后的数据进行逆KL变换则可以恢复出完整的4D断层图像序列,因此有效提高4D FMT重建方法的速度,经本发明的仿真试验表明,采用本发明的重建方法,4D FMT的重建时间大约为传统重建方法的1.5%。本发明可以广泛应用于动态荧光分子图像的重建过程中。
附图说明
图1是本发明基于KL变换的4D FMT重建方法结构框图;
图2是本发明的仿真实验结构示意图,其中,图2(a)是主视示意图,图2(b)是图2(a)的俯视示意图;
图3是吲哚菁绿ICG(Indocyanine Green,ICG)在小鼠肝部的代谢曲线示意图,也是本发明仿真实验中所使用的ICG浓度时间曲线示意图,横坐标为时间,纵坐标为ICG在小鼠内的浓度值;
图4(a)~(d)是采用本发明方法计算所得仿真实验的FMT重建结果的效果示意图,图4(a)是第2帧(2min)荧光投影图像的重建结果效果示意图,图4(b)是第20帧(20min)荧光投影图像重建结果效果示意图;图4(c)是第50帧(50min)荧光投影图像重建结果效果示意图;图4(d)是第120帧(120min)荧光投影图像重建结果效果示意图;
图5是采用本发明的重建方法所获得的ICG在仿真实验中的浓度时间变化曲线示意图,其中,圆形为图3中的实际的ICG浓度,圆点为采用本发明方法计算所得到的ICG浓度,横坐标为采集时间点。
具体实施方式
下面结合附图和实施例对本发明进行详细的描述。
对于FMT成像,基于一阶波恩近似,在检测点r处获取的荧光投影图像可以采用下式表示:
Φm(r)=∫r′∈VG(r,r′)n(r′)Φx(r′)dr′ (1)
式中,V是重建区域;Φx为光场分布函数,描述了对应激发光波长的光子密度,Φm为发射(荧光)波长的光子密度;n(r′)为重建区域V中点r′处的待重建的荧光产额,通常n(r′)正比于荧光探针的浓度;格林函数G(r,r′)描述在荧光波段,光子在成像对象体内的传播。
对于动态的FMT,重建区域V被离散成小的网格ΔV后,公式(1)可被离散为下述线性方程:
uk=W1nk (2)
式中,uk是一个M×1的向量,用以描述第k帧的荧光投影图像;nk是一个N×1的向量,用以描述待重建的荧光断层图像;W1(r)=G(r,r′)Φx(r′)ΔV(r)是基于漫射方程生成的M×N权重矩阵,用以描述光子在成像对象体内的传播。
对于4D荧光断层图像重建,现有的方法是采用公式(2)逐帧进行重建,得到每帧的断层图像nk,然后将所有的断层图像进行组合生成动态的FMT图像序列,上述重建方法忽略了荧光投影图像uk在时间轴上的相关性。为了解决上述问题,可以将所有测量得到的荧光投影图像序列作为整体进行4D FMT重建,重建公式为:
u=Wn (3)
式中,是所采集到的完整荧光投影图像序列,其中描述了第i帧的荧光投影图像;L是动态成像中扫描的总帧数(总次数);一系列待重建的荧光断层图像,其中描述了第i时刻荧光探针在成像体内的浓度分布;W=diag[W1,W1,…,W1]是动态FMT成像的空间-时间***矩阵,其大小为ML×NL;但是此种直接对荧光投影图像序列进行4D重建的方法,W的内存占用及随后的计算量都超出了普通计算机能力。
为了解决上述问题,本发明提出了基于KL变换的4D FMT重建方法,即通过对公式(3)进行KL(Karhunen–Loève)变换。在动态FMT的研究中,本发明令λ=[λ1,λ2,...,λL]T表示像素的时间强度曲线(Time Activity Curve,TAC),KL变换和KL逆变换的公式分别为:
KL变换:
KL逆变换:
式中,是经KL变换后的数据;A是KL变换矩阵,AT是KL变换的逆变换矩阵(在KL变换中KL的逆变换矩阵近似等于KL变换矩阵的转置矩阵);KL变换矩阵A可以通过以下公式进行求解:
PtimeAT=ATD(6)
式中,Ptime为所有TAC的协方差矩阵,矩阵大小为L×L:
式中, dl是协方差矩阵Ptime的第l个特征值。在动态FMT重建中,Ptime的维数通常较小(L<200),因此在公式(6)中,特征矢量能够被快速计算。
为了描述荧光投影图像u在时间方向上的KL变换,本发明定义如下KL变换矩阵AM:
式中,IM是一个M×M的单位矩阵;是克罗内克积(kronecker product)。
将公式(2)的两侧分别乘以AM,基于Wernick(韦尼克)推导,得到以下公式:
令:
将公式(10)和公式(11)代入公式(9),公式(9)可以写为:
由于KL变换能够去除荧光投影图像u的时间相关性。因此可以对KL变换后的每个KL分量基于公式(12)在KL域中单独进行重建,即利用KL变换的去相关性,将复杂的4D荧光投影数据的重建转化为一系列KL分量在KL域的3D重建问题,降低了4D FMT重建的计算规模。
此外,由于公式(12)与公式(3)具有相同的成像模型,因此可以采用现有各种的FMT重建方法完成各个KL分量的重建,而不需对其进行修改,对重建后的KL分量再进行逆KL变换,就能够获取完整的4D断层图像序列n:
式中,IN是一个N×N的单位矩阵,由于KL变化具有压缩属性,经过KL变换后,信号被保存在前几个KL分量,更高阶的KL分量几乎没有包含有用的信号,通过抛弃高阶KL分量,即仅通过重建前几个KL分量,然后对重建后的KL分量进行逆KL变就可以恢复出完整的4D断层序列,在选择需保留的KL分量个数时可以采用现有的“scree”方法、“Kaiser”方法或基于特征值的累计方差方法。
在实验过程中获得荧光投影图像可以采用现有的非接触、全角度小动物诱发荧光分子成像***,小动物诱发荧光分子成像***包括有计算机、小动物旋转平台装置和荧光成像激发与检测装置。
如图1所示,下面通过具体实施例进一步说明对动态荧光投影图像序列进行4D荧光分子断层图像重建的方法流程,具体包括以下步骤:
1)在不同时间点t1、t2…tL′…tL,使用全角度、非接触FMT成像***,等时间间隔采集反映荧光探针在成像对象体内的代谢分布的荧光投影图像(Proj.1、Proj.2…Proj.K)。对于每个成像时刻(帧)所采集的不同角度的荧光投影图像,将其转化为一个列向量,通过变换所有帧,得到KL变换的输入矩阵u(图1中纵向虚线框内的方框分别表示各采集时刻所采集到不同角度的荧光投影图像序列1、2…L′…L)。
如图2(a)~(b)所示,本发明以模拟ICG在小鼠肝中的代谢情况(如图3所示)的仿真实验为实施例进一步说明动态荧光投影图像序列的采集和重建过程。
本发明参照小鼠形体大小和组织的平均光学参数,仿真实验设置一直径为3.0cm、高为4.5cm的圆柱作为仿体1,仿体1的光学约化散射系数为10cm-1,光学吸收系数为0.3cm-1,将一个直径为0.43cm、高为0.5cm的小圆柱作为荧光棒2放置在仿体1内后,将仿体1放置在小动物旋转平台装置上,开启激光光源照射荧光棒2,利用荧光成像激发与检测装置对成像对象进行全角度的荧光投影图像采集,成像过程中,为了模拟ICG在小鼠肝部的时间代谢曲线,在成像的不同时刻,根据小鼠肝部代谢曲线将不同时刻的ICG浓度赋值给荧光棒(如图3所示),对于每一时间点(帧)成像,等间隔采集24个角度的荧光投影图像,即对每一时间点共生成24个不同角度的荧光投影图像。考虑到由于在在体实验中,一帧的成像时间近似为1min,因此在本发明的仿真实验中,设置每帧图像的采集间隔为1min。具体的动态仿真过程为,首先,根据给定的ICG的浓度时间曲线(如图3所示),在不同的时间点(从1min到120min)为荧光棒2设置不同的ICG浓度。然后,根据公式(2)生成120帧(每一min为一帧)的荧光投影图像序列,用以描述模拟ICG在小鼠肝内不同时刻的浓度分布信息。最后,将所有采集到的荧光投影图像序列(总计120帧,每帧24个不同角度的荧光投影数据)组合成一个整体,形成动态的荧光投影图像序列,作为KL变换的输入矩阵u。
2)根据公式(7),构建输入矩阵u的时间协方差矩阵Ptime。
3)通过求解PtimeAT=ATD,生成L×L的KL变换矩阵A,在动态FMT成像中,由于总的扫描次数L通常较小,因此变换矩阵A的计算是快速的。
4)基于公式(10),获取经KL变换后的数据(图1中细线框内的 表示经KL变换后得到的KL分量)。
5)采用“scree”方法,保留KL分量中的前L′个KL分量(L′=2表示取前两个KL分量)。
6)根据光场分布函数Φx和格林函数G生成权重矩阵W1。
7)基于公式k=1,2,...,L′,使用3D FMT重建方法,对步骤5)中获取的每一个KL分量分别进行3D重建得到重建结果(如图1中粗实线框内得到的L′个KL分量的重建结果)。其中,3D FMT重建方法为现有的方法,可以采用代数重建方法(Algebraic Reconstruction Technique,ART)方法、最小二乘QR分解迭代法(LeastSquares QR,LSQR)方法或Bayesian方法。
8)如图4(a)~(d)所示,基于公式(13),将步骤7)的重建结果进行逆KL变换,得到完整的4D荧光断层图像序列n1、n2…nL(如图1所示)。
上述实施例中,步骤2)和3)根据输入矩阵u生成KL变换矩阵A的方法为还可以采用通过对输入矩阵u进行奇异值分解求得。
上述实施例中,本发明的实施例采用仿真实验,在仿真实验中,改变荧光探针的浓度只需改变浓度设定值即可;在仿体实验中,改变荧光探针的浓度需要将放置荧光探针的玻璃棒取出,清洗,注入新的不同浓度的荧光探针,然后在放回仿体内;在活体实验中,在尾静脉注射荧光探针后,荧光探针浓度随时间自行改变。
如图5所示,将各个时刻重建的ICG的浓度取平均值,生成采用本方法计算得到的不同时刻所对应的ICG浓度的分布图,图中圆点表示120个采集时间点所对应的使用本方法计算得到ICG的浓度,图中的圆形为图3的给定的实际ICG的浓度时间曲线,可以看出两者是基本接近的,充分说明了本发明方法的有效性。
上述各实施例仅用于说明本发明,其中本发明的方法的实施步骤等都是可以有所变化的,凡是在本发明技术方案的基础上进行的等同变换和改进,均不应排除在本发明的保护范围之外。
Claims (9)
1.一种动态荧光分子断层图像的重建方法,包括以下步骤:
1)设置一包括有计算机、小动物旋转平台装置和荧光成像激发与检测装置的小动物诱发荧光分子成像***;
2)采用所述小动物诱发荧光分子成像***,在不同时间点,全角度、等时间间隔采集反映荧光探针在成像对象体内的代谢分布的荧光投影图像,并将每一时间点所采集的不同角度的荧光投影图像存储为一矩阵序列,将所述矩阵序列标记为输入矩阵u;
3)根据输入矩阵u生成KL变换矩阵A;
4)基于求解输入矩阵u经KL变换后的数据式中,IM为M×M的单位矩阵,为克罗内克积;
5)获取经KL变换后需要进行3D重建的L'个KL分量;
6)根据光场分布函数Φx和格林函数G生成权重矩阵W1,其中,Φx描述了对应激发光波长的光子密度,G描述了在荧光谱段,光子在成像对象体内的传播;
7)根据公式采用3D FMT重建方法在KL域中对步骤5)中获取的每一个KL分量分别进行3D重建,其中,k=1,2,...,L',W-1 1为权重矩阵W1的逆矩阵;
8)对上述步骤7)中获取的KL分量的重建结果进行逆KL变换,得到完整的4D荧光断层图像序列。
2.如权利要求1所述的一种动态荧光分子断层图像的重建方法,其特征在于:所述步骤3)根据输入矩阵u生成KL变换矩阵A的方法为:
①构建输入矩阵u的时间协方差矩阵Ptime;
②通过求解PtimeAT=ATD,生成KL变换矩阵A,式中dl是协方差矩阵Ptime的第l个特征值,L是动态成像中扫描的总帧数。
3.如权利要求1所述的一种动态荧光分子断层图像的重建方法,其特征在于:所述步骤3)根据输入矩阵u生成KL变换矩阵A的方法为:通过对输入矩阵u进行奇异值分解求得。
4.如权利要求1或2或3所述的一种动态荧光分子断层图像的重建方法,其特征在于:所述步骤5)中获取经KL变换后需要进行3D重建的L'个KL分量采用基于特征值的累计方差方法。
5.如权利要求1或2或3所述的一种动态荧光分子断层图像的重建方法,其特征在于:所述步骤6)中在检测点r处获取的荧光投影图像采用下式表示:
其中,W1(r)=G(r,r')Φx(r')ΔV(r),ΔV是重建区域V被离散成的小网格;V是重建区域;Φm是发射荧光波长的光子密度;n(r')为重建区域V中点r′处的待重建的荧光产额。
6.如权利要求4所述的一种动态荧光分子断层图像的重建方法,其特征在于:所述步骤6)中在检测点r处获取的荧光投影图像采用下式表示:
其中,W1(r)=G(r,r')Φx(r')ΔV(r),ΔV是重建区域V被离散成的小网格;V是重建区域;Φm是发射荧光波长的光子密度;n(r')为重建区域V中点r′处的待重建的荧光产额。
7.如权利要求1或2或3或6所述的一种动态荧光分子断层图像的重建方法,其特征在于:所述步骤7)采用3D FMT重建方法在KL域中对步骤5)中获取的每一个KL分量分别进行3D重建中,3D FMT重建方法采用ART方法、LSQR方法和Bayesian方法中的一种。
8.如权利要求4所述的一种动态荧光分子断层图像的重建方法,其特征在于:所述步骤7)采用3D FMT重建方法在KL域中对步骤5)中获取的每一个KL分量分别进行3D重建中,3D FMT重建方法采用ART方法、LSQR方法和Bayesian方法中的一种。
9.如权利要求5所述的一种动态荧光分子断层图像的重建方法,其特征在于:所述步骤7)采用3D FMT重建方法在KL域中对步骤5)中获取的每一个KL分量分别进行3D重建中,3D FMT重建方法采用ART方法、LSQR方法和Bayesian方法中的一种。
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/CN2012/000576 WO2013159250A1 (zh) | 2012-04-28 | 2012-04-28 | 一种动态荧光分子断层图像的重建方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103188989A CN103188989A (zh) | 2013-07-03 |
CN103188989B true CN103188989B (zh) | 2015-05-27 |
Family
ID=48679799
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201280002807.8A Expired - Fee Related CN103188989B (zh) | 2012-04-28 | 2012-04-28 | 一种动态荧光分子断层图像的重建方法 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN103188989B (zh) |
WO (1) | WO2013159250A1 (zh) |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104434035B (zh) * | 2014-12-04 | 2017-08-29 | 中国人民解放军第四军医大学 | 基于光学影像的非侵入式测定裸鼠肿瘤受体密度的方法 |
DE102015112628A1 (de) * | 2015-07-31 | 2017-02-02 | Carl Zeiss Microscopy Gmbh | Verfahren zur Erzeugung eines digitalen Fluoreszenzbildes |
CN106683180B (zh) * | 2017-01-03 | 2019-08-27 | 清华大学 | 图像处理方法及*** |
CN107184181A (zh) * | 2017-05-15 | 2017-09-22 | 清华大学 | 动态荧光分子断层成像的处理方法和*** |
CN113327305B (zh) * | 2021-05-31 | 2022-06-24 | 华中科技大学 | 一种模型驱动的深度学习荧光分子断层成像方法及*** |
CN113724351B (zh) * | 2021-08-24 | 2023-12-01 | 南方医科大学 | 一种光声图像衰减校正方法 |
CN113951831B (zh) * | 2021-11-19 | 2023-08-29 | 复旦大学 | 基于随机重构策略的超高分辨荧光分子断层成像方法 |
CN113951835B (zh) * | 2021-12-01 | 2024-02-13 | 复旦大学 | 一种基于光学断层重建策略的三维荧光显微成像方法 |
CN114463457B (zh) * | 2022-01-07 | 2023-03-21 | 苏州大学 | 基于小波和舒尔分解的荧光分子断层重建方法及*** |
CN115607108A (zh) * | 2022-09-30 | 2023-01-17 | 西北大学 | 一种基于表面测量信号盲源分离的多目标重建方法 |
CN116124752B (zh) * | 2023-04-04 | 2023-07-18 | 之江实验室 | 一种基于多光谱调控的组织仿生模体及其生成方法 |
CN116823836B (zh) * | 2023-08-30 | 2023-11-07 | 南京诺源医疗器械有限公司 | 适用于细胞核靶向上转换荧光探针的数据处理方法 |
CN117946516B (zh) * | 2024-03-26 | 2024-07-09 | 之江实验室 | 聚氨酯荧光光学仿体 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1994229A (zh) * | 2006-12-01 | 2007-07-11 | 清华大学 | 旋转平台式小动物在体多模成像检测*** |
CN101057788A (zh) * | 2006-09-29 | 2007-10-24 | 清华大学 | 核素与荧光双模一体的小动物分子成像的数据获取*** |
CN101396262A (zh) * | 2008-10-31 | 2009-04-01 | 清华大学 | 一种基于线性关系的荧光分子断层成像重建方法 |
CN101984928A (zh) * | 2010-09-29 | 2011-03-16 | 北京大学 | 一种多模态分子层析成像*** |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8918163B2 (en) * | 2008-03-25 | 2014-12-23 | Visen Medical, Inc. | Animal holder for in vivo tomographic imaging with multiple modalities |
US20110176715A1 (en) * | 2010-01-21 | 2011-07-21 | Foos David H | Four-dimensional volume imaging system |
-
2012
- 2012-04-28 CN CN201280002807.8A patent/CN103188989B/zh not_active Expired - Fee Related
- 2012-04-28 WO PCT/CN2012/000576 patent/WO2013159250A1/zh active Application Filing
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101057788A (zh) * | 2006-09-29 | 2007-10-24 | 清华大学 | 核素与荧光双模一体的小动物分子成像的数据获取*** |
CN1994229A (zh) * | 2006-12-01 | 2007-07-11 | 清华大学 | 旋转平台式小动物在体多模成像检测*** |
CN101396262A (zh) * | 2008-10-31 | 2009-04-01 | 清华大学 | 一种基于线性关系的荧光分子断层成像重建方法 |
CN101984928A (zh) * | 2010-09-29 | 2011-03-16 | 北京大学 | 一种多模态分子层析成像*** |
Non-Patent Citations (4)
Title |
---|
Can T. Xu,et al.Fluorescence diffuse optical tomography using upconverting nanoparticles.《APPLIED PHYSICS LETTERS》.2009,第94卷 * |
Ralf B. Schulz,et al.Hybrid System for Simultaneous Fluorescence and X-Ray Computed Tomography.《 IEEE TRANSACTIONS ON MEDICAL IMAGING》.2010,第29卷(第2期), * |
姚俊杰,等.用于荧光分子断层成像的小动物躯干部分三维表面轮廓重建研究.《中国生物医学工程学报》.2008,第27卷(第3期),360-365. * |
邹玮,等.荧光分子断层成像正向问题的并行计算.《光学学报》.2007,第27卷(第3期),443-450. * |
Also Published As
Publication number | Publication date |
---|---|
WO2013159250A1 (zh) | 2013-10-31 |
CN103188989A (zh) | 2013-07-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103188989B (zh) | 一种动态荧光分子断层图像的重建方法 | |
Dehghani et al. | Near infrared optical tomography using NIRFAST: Algorithm for numerical model and image reconstruction | |
Tong et al. | Domain transform network for photoacoustic tomography from limited-view and sparsely sampled data | |
US9495516B2 (en) | Systems, methods, and devices for image reconstruction using combined PDE-constrained and simplified spherical harmonics algorithm | |
Abdoulaev et al. | Optical tomography as a PDE-constrained optimization problem | |
CN102334979B (zh) | 一种基于迭代收缩的双模态融合断层成像方法 | |
CN103027705B (zh) | 产生运动补偿的ct 图像数据组的方法和*** | |
US9672639B1 (en) | Bioluminescence tomography reconstruction based on multitasking bayesian compressed sensing | |
CN104665770B (zh) | 一种用于近红外光脑功能研究的自引导扩散光层析成像方法 | |
CN104903933A (zh) | 投影数据降噪 | |
Mohajerani et al. | Optical and optoacoustic model-based tomography: theory and current challenges for deep tissue imaging of optical contrast | |
CN103300829A (zh) | 一种基于迭代重加权的生物自发荧光断层成像方法 | |
Liu et al. | 4-D reconstruction for dynamic fluorescence diffuse optical tomography | |
Li et al. | Multispectral interlaced sparse sampling photoacoustic tomography based on directional total variation | |
Zheng et al. | Quantitative photoacoustic tomography with light fluence compensation based on radiance Monte Carlo model | |
CN106023124A (zh) | 一种基于先验图像约束的磁共振稀疏重建图像方法 | |
CN108830915A (zh) | 口腔影像3d模拟动画实现方法及装置 | |
Guan et al. | Optical tomography reconstruction algorithm based on the radiative transfer equation considering refractive index: Part 2. Inverse model | |
Imms et al. | A high performance biometric signal and image processing method to reveal blood perfusion towards 3D oxygen saturation mapping | |
CN204072055U (zh) | 头戴式分子影像导航*** | |
US20200329974A1 (en) | Method and system for fluence matching in photoacoustic imaging | |
CN110473145A (zh) | 被动式太赫兹图像融合超高分辨率重建的检测方法及装置 | |
CN113724351B (zh) | 一种光声图像衰减校正方法 | |
Liu et al. | Parametric diffuse optical imaging in reflectance geometry | |
CN102599910B (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 | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150527 Termination date: 20190428 |