[发明专利]一种动态荧光分子断层图像的重建方法

权利要求书

1.一种动态荧光分子断层图像的重建方法，包括以下步骤：

1）设置一包括有计算机、小动物旋转平台装置和荧光成像激发与检测装置的小动物诱发荧光分子成像系统；

2）采用所述小动物诱发荧光分子成像系统，在不同时间点，全角度、等间隔采集反映荧光探针在成像对象体内的代谢分布的荧光投影图像，并将每一时间点所采集的不同角度的荧光投影图像存储为一矩阵序列，将所述矩阵序列标记为输入矩阵u；

3）根据输入矩阵u生成KL变换矩阵A；

4）基于求解输入矩阵u经KL变换后的数据式中，I_M为M×M的单位矩阵，为克罗内克积；

5）获取经KL变换后需要进行3D重建的L′个KL分量；

6）根据光场分布函数Φ_x和格林函数G生成权重矩阵W₁，其中，Φ_x描述了对应激发光波长的光子密度，G描述了在荧光谱段，光子在成像对象体内的传播；

7）根据公式采用3D FMT重建方法在KL域中对步骤5）中获取的每一个KL分量分别进行3D重建，其中，k＝1,2,...,L′，W^-1₁为权重矩阵W₁的逆矩阵；

8）对上述步骤7）中获取的KL分量的重建结果进行逆KL变换，得到完整的4D荧光断层图像序列。

2.如权利要求1所述的一种动态荧光分子断层图像的重建方法，其特征在于：所述步骤3）根据输入矩阵u生成KL变换矩阵A的方法为：

①构建输入矩阵u的时间协方差矩阵P^time；

②通过求解P^timeA^T=A^TD，生成KL变换矩阵A，式中d_l是协方差矩阵P^time的第l个特征值，L是动态成像中扫描的总帧数。

3.如权利要求1所述的一种动态荧光分子断层图像的重建方法，其特征在于：所述步骤3）根据输入矩阵u生成KL变换矩阵A的方法为：通过对输入矩阵u进行奇异值分解求得。

4.如权利要求1或2或3所述的一种动态荧光分子断层图像的重建方法，其特征在于：所述步骤5）中获取经KL变换后需要进行3D重建的L′个KL分量采用scree方法、Kaiser方法和基于特征值的累计方差方法中的一种。

5.如权利要求1或2或3所述的一种动态荧光分子断层图像的重建方法，其特征在于：所述步骤6）中在检测点r处获取的荧光投影图像采用下式表示：

Φ_m(r)＝∫_r′∈VG(r,r′)n(r′)Φ_x(r′)dr′

其中，W₁(r)＝G(r,r′)Φ_x(r′)ΔV(r)，ΔV是重建区域V被离散成的小网格；V是重建区域；Φ_m是发射荧光波长的光子密度；n(r′)为重建区域V中点r′处的待重建的荧光产额。

6.如权利要求4所述的一种动态荧光分子断层图像的重建方法，其特征在于：所述步骤6）中在检测点r处获取的荧光投影图像采用下式表示：

Φ_m(r)＝∫_r′∈VG(r,r′)n(r′)Φ_x(r′)dr′

其中，W₁(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方法中的一种。