一种基于迭代收缩的双模态融合断层成像方法技术

本发明专利技术属于医学分子影像领域，涉及自发荧光断层成像与计算机断层成像两种模态融合的方法，尤其是一种基于迭代收缩的双模态融合断层成像方法。这种技术用于定量并定位重建目标体内部的光源强度与位置，解决由目标体表面有限的光强分布反推得到内部全部的光强分布的逆向问题。本技术方案的要点在于，利用自发荧光断层成像中获得的表面光强信息，融合计算机断层成像中获得的内部几何结构信息，采用迭代收缩法将重建中复杂的多维求优过程转化为一维并行求优的高效循环过程，从整体上得到对正则参数、lp范数、噪声、初值鲁棒的精确重建结果。该技术能有效地应用于目标体全身生理代谢的研究，重建效率高且适用于成像系统性能较低的情况。

【技术实现步骤摘要】

本专利技术属于医学分子影像领域，涉及自发荧光断层成像与计算机断层成像两种模态融合的方法，尤其是。
技术介绍
自发荧光断层成像是近年来新兴的一种分子影像技术，是一种可以在细胞分子水平上实现生物组织内部生物物理过程无创实时动态连续成像的一种技术。通过获取到的目标体表面出射光强分布，探测在目标体内部的光源分布则能够对目标体进行定量化研究， 为荧光光源提供三维的定位信息。目前普遍采用的断层成像技术大多是基于梯度的传统方法，由于自发荧光断层成像的信号十分微弱，导致了采集的数据集噪声较大，加之该问题的固有病态性，造成在待探测过程中对正则参数的选择、Ip范数的选择、初始值的选择、可行域的选择以及噪声的容忍度都有着很强的要求。这限制了自发荧光断层成像在实际中的应用。寻找一种高效、鲁棒、 全域的探测方法一直是该领域的热点、难点问题，到目前为止仍然没有一种通用的技术能够很好地解决这个问题。迭代收缩法是一种新兴的数值计算方法，与传统的迭代优化方法相比，它更具高效性，且近年来已被陆续应用于图像降噪以及逆问题求解中，推动了这些领域的发展。
技术实现思路
本专利技术的目的是提供一种基于迭代收缩的双模态融合断层成像技术，利用高性能探测器捕捉目标体在自发荧光断层成像中出射的表面光强分布，融合计算机断层成像中获得的目标体内部几何结构信息，使用扩散方程建立光子传输模型，采用迭代收缩法将重建中复杂的多维求优过程转化为对每个变量逐次进行一维并行求优的高效循环过程，从整体上得到对正则参数、Ip范数、噪声、初值鲁棒的精确重建结果。为了实现上述目的，本专利技术涉及的基于迭代收缩的双模态融合断层成像方法具体步骤如下步骤a，采用波尔兹曼传输方程的一阶球谐近似得到的扩散方程建立光子在目标体内部的传输模型，即扩散方程传输模型；步骤b，利用有限元方法将扩散方程传输模型转换成线性矩阵方程的形式；步骤c，基于线性矩阵方程形式的扩散方程传输模型，利用正则理论，建立优化目标函数f (S)，其中，S为待求解的目标体内部光强分布；步骤d，对f (S)进行降维处理，得到一维优化目标函数；步骤e，基于所述一维优化目标函数，以迭代的方式求解得到最优目标体内部光强分布S。本专利技术的有益效果是通过使用基于迭代收缩的双模态融合断层成像技术，融合自发荧光断层成像与计算机断层成像两种模态的信息降低逆问题求解的病态性，由并行计算迭代收缩算子的复杂的多维求优过程转化为对每个变量逐次进行一维并行求优的高效循环过程，最终重建得到目标体内部光源分布的定量和定位信息。使用该方法的优势为，在不同正则参数、不同Ip范数、不同初值、较强噪声以及全部目标区域下均能计算得到精确的重建结果；同时，通过迭代算子的降维并行处理，使这种双模态融合的断层成像技术具有高于传统方法的重建效率。附图说明图1为基于迭代收缩的双模态融合断层成像技术框架流程图。图2为在体非勻质小鼠的自发荧光断层成像实验的表面光强图。图3为在体非勻质小鼠的计算机断层成像实验的解剖结构图。图4为平滑算子θ (w)的参数w与其拟合对象的参数ρ的关系图。图5为使用本技术在不同正则参数下的重建结果图。图6为使用本技术在不同Ip范数下的重建结果图。图7为目标优化函数值与迭代次数的演化曲线图。具体实施例方式为使本专利技术的目的、技术方案和优点更加清楚明白，以下结合具体实施例，并参照附图，对本专利技术进一步详细说明。下面结合附图详细说明本专利技术技术方案中所涉及的各个细节问题。应指出的是， 所描述的实例仅旨在便于对本专利技术的理解，而对其不起任何限定作用。本专利技术是一种基于迭代收缩的双模态融合断层成像技术。在自发荧光断层成像中，例如通过对小动物转基因，使得小动物在病灶处产生荧光蛋白酶，进而通过静脉注射荧光底物，从而在该蛋白酶的催化作用下发射出可见或近红外光，光在生物组织中的传输将历经多次散射。最终有一部分光子到达小动物皮肤并出射到自由空间中被探测器接收，作为自发荧光断层成像中获取到的表面光强分布信息。再融合通过计算机断层成像获取到的小动物内部解剖结构信息，借助不同组织器官的光学特性参数和分割、配准、网格剖分技术建立非勻质小动物模型，最终重建出小动物体内的光源强度与位置信息。本专利技术能在不同正则参数、不同Ip范数、不同初值、较强的噪声的情况下，精确重建出小动物全身的内部光源{曰息。下面结合附图，详细描述本专利技术的断层图像重建方法。作为一种具体的实现技术方案，其框架流程图参阅图1，步骤如下步骤1 采用波尔兹曼传输方程的一阶球谐近似得到的扩散方程建立光子在目标体内部的传输模型，即扩散方程传输模型- V. +//a(r)0(r)= ^(r)式中，r是位置向量，D是散射系数，Pa是吸收系数，Φ是目标体表面光强分布，S 表示待求解的目标体内部光强分布。步骤2 利用自发荧光断层成像获取目标体表面出射的光强分布测量值Φ* ；利用计算机断层成像获取目标体内部几何结构信息，通过分割出目标体内部不同的器官结构， 并为不同器官设置相应的光学参数，包括吸收系数和散射系数，从而建立非勻质目标体模型。步骤3 对非勻质目标体模型进行几何结构的网格剖分，并利用配准技术将表面光强信息与几何结构信息进行位置融合与能量融合。请参阅图2和图3为非勻质小鼠断层成像实验获取的信息。图2给出在体非勻质小鼠表面光强信息，该数据采集于自发荧光断层成像实验，并经过配准等后处理操作。图3给出在体非勻质小鼠的解剖结构信息，该数据采集于计算机断层成像实验，并经过分割等后处理操作。融合图2和图3中的数据，并作网格剖分处理，建立起非勻质小鼠的模型，包含了躯干、心脏、肺、肝脏、肌肉、骨骼等组织器官。步骤4 基于步骤3获得的信息，比如离散化几何信息和能量信息，利用有限元方法将扩散方程传输模型转换成线性矩阵方程MS= Φ的形式，其中，M是描述模型性质的系统矩阵。步骤5 为了以优化的方式迭代求解目标体内部光强分布S，要基于线性矩阵方程形式的扩散方程传输模型，利用正则理论，建立优化目标函数权利要求1.，其特征在于，该方法包括以下步骤步骤a，采用波尔兹曼传输方程的一阶球谐近似得到的扩散方程建立光子在目标体内部的传输模型，即扩散方程传输模型；步骤b，利用有限元方法将扩散方程传输模型转换成线性矩阵方程的形式； 步骤c，基于线性矩阵方程形式的扩散方程传输模型，利用正则理论，建立优化目标函数f (S)，其中，S为待求解的目标体内部光强分布；步骤d，对f (S)进行降维处理，得到一维优化目标函数；步骤e，基于所述一维优化目标函数，以迭代的方式求解得到最优目标体内部光强分布S。2.根据权利要求1所述的基于迭代收缩的双模态融合断层成像方法，其特征在于，所述扩散方程传输模型表示为-V · +//a(r)0(r)=对r)，其中，r是位置向量，D是散射系数，是吸收系数，Φ是目标体表面光强分布，S表示待求解的目标体内部光强分布。3.根据权利要求1所述的基于迭代收缩的双模态融合断层成像方法，其特征在于，步骤b之前还包括建立非勻质目标体模型的步骤，并且所述建立非勻质目标体模型的步骤具体为利用自发荧光断层成像获取目标体表面出射的光强分布测量值Φ* ；利用计算机断层成像获取目标体内部几何结构信息，通过分割出目标体内部不同的器官结构，并为不同器官设置相应的光学参数，包括吸收系数和散射系数本文档来自技高网...

【技术保护点】
１．一种基于迭代收缩的双模态融合断层成像方法，其特征在于，该方法包括以下步骤：步骤ａ，采用波尔兹曼传输方程的一阶球谐近似得到的扩散方程建立光子在目标体内部的传输模型，即扩散方程传输模型；步骤ｂ，利用有限元方法将扩散方程传输模型转换成线性矩阵方程的形式；步骤ｃ，基于线性矩阵方程形式的扩散方程传输模型，利用正则理论，建立优化目标函数ｆ（Ｓ），其中，Ｓ为待求解的目标体内部光强分布；步骤ｄ，对ｆ（Ｓ）进行降维处理，得到一维优化目标函数；步骤ｅ，基于所述一维优化目标函数，以迭代的方式求解得到最优目标体内部光强分布Ｓ。

【技术特征摘要】

【专利技术属性】
技术研发人员：田捷，秦承虎，吴萍，杨鑫，刘凯，
申请(专利权)人：中国科学院自动化研究所，
类型：发明
国别省市：11