一种基于自适应近邻正交最小二乘算法的光源重建方法

1.本发明属于分子影像领域，具体为一种基于自适应近邻正交最小二乘算法的光源重建方法。


背景技术：

2.荧光分子层析成像(以下简称fmt)或荧光分子断层成像，利用外界光源对生物体内的荧光分子探针进行激发使其发光，发出的荧光通过生物体外部的设备进行采集，利用光传输模型和数学方法可以反演出生物体内部荧光探针的三维分布。它是一种非侵入式的光学分子成像技术，具有很强的特异性和敏感性。因此，在肿瘤检测、药物开发和疗效评价的预临床研究中具有巨大的潜力。但是，光子的强散射和表面测量数据的不足，使得fmt重建结果不稳定，重建结果的质量还有待提高，而这些进一步阻碍了fmt在早期肿瘤检测的应用进程。
3.fmt三维重建存在严重的病态性，在数学上是典型的逆问题。因为近红外光在生物体组织中具有强散射特性，同时，获取的表面荧光数据有限，并且还会有噪声的干扰，这大大增加了问题的不适定性。为了获得稳定而准确的解，各种正则化方法应用到fmt重建。
4.现在比较流行的正则化方法主要有基于l2范数，l1范数和l0范数。其中，l2范数的方法设计简单，计算快捷，但其过平滑性往往导致重建分辨率差。fmt的一个典型应用就是肿瘤早期检测模型，该模型中肿瘤相对成像目标来说具有显著的稀疏特性。为了利用稀疏先验信息，l1范数和l0范数正则化方法开始应用到该模型中。它们有效克服了l2范数的过平滑弊端，解的稀疏性也较l2范数有了很大的提高。因此得到了研究者的广泛关注和重视。


技术实现要素：

5.为了提高荧光分子断层成像的重建结果质量，本发明提出了一种基于自适应近邻正交最小二乘算法的光源重建方法，求解基于l0范数的正则化问题。该算法是典型的贪婪算法，它不需要事先知道变量的稀疏度，而是能够自适应的获取稀疏度，同时采用近邻策略选取最优支撑集，最终提高了荧光分子断层成像的重建质量和空间分辨率。
6.具体流程包括：
7.一种基于自适应近邻正交最小二乘算法的荧光分子断层成像方法，包括以下步骤：
8.步骤一、获取测量数据：采集目标生物体的表面荧光数据和解剖结构信息，对重建目标生物体进行有限元网格化，然后基于光子传输模型的扩散方程来近似模型，建立表面荧光数据与重建目标内部荧光团分布的线性关系；
9.φ＝αx
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1)；
10.其中，φ为表面荧光数据，α表示系统矩阵，x是要求解的重建目标内部荧光团的分布；
11.步骤二、归一化测量数据：将步骤一中的线性关系转化为l0范数极小化问题：
[0012][0013]
||x||0是0范数，表示非零元素数量，ε为给定的阈值；
[0014]
步骤三、使用自适应近邻正交最小二乘算法进行重建，得到最优解x。
[0015]
可选的，所述的步骤三具体包括：
[0016]
x＝ki(si)；
[0017]ki
＝k
i-1
+s，i是迭代次数，s是步长，表示向上取整，s0＝k0＝6；
[0018]
si＝li(e
i-1
)∪s
i-1
，li(e
i-1
)表示从e
i-1
中选出li个节点，求出s
i-1
中所有节点的最近邻集合取并集，得到的集合记作近邻集合e
i-1
，将该过程用近邻算子n来描述，则e
i-1
＝n(s
i-1
)；
[0019]
i表示迭代次数，初始值为1；k表示稀疏度；si表示支撑集；
[0020]
每迭代一次i＝i+1，直到残差小于给定阈值或稀疏度k的变化步长为1时，停止迭代。
[0021]
优选的，所述的步骤三具体包括：
[0022]
3.1初始化参数：i表示迭代次数，初始值i＝1，稀疏度k0＝6，l0＝10，初始残差r0＝φ；
[0023]
元素选择标准qj，在迭代中，选择一个新的列索引jm：
[0024][0025][0026][0027]
其中i是列集合，aj是系统矩阵a的第j列，u
l
是正交基向量；根据上述公式(3)选出l0个元素构成初始支撑集s0；
[0028]
3.2自适应调整步长：采用非线性函数来调整步长，该函数是一个单调递减函数，步长表述如下：
[0029][0030]
其中s是步长，表示向上取整；设往支撑集si中添加的原子数为li，其迭代公式为li＝l
i-1-s，其中s0＝l0＝10，随着迭代次数的增加，向支撑集中添加的原子数逐渐变少，直到原子数变为1，则停止迭代；
[0031]
稀疏度迭代公式设为ki＝k
i-1
+s，其中s0＝k0＝6；随着迭代次数的增加，稀疏度逐渐增加，其增量开始阶段是较大步长，完成阶段是较小的步长，直到为1；
[0032]
3.3支撑集的更新：基于有限元理论，根据四面体单元结构构建最近邻节点集合。s
i-1
[0033]
是第i-1次迭代生成的支撑集，设任意节点gk∈s
i-1
，则该节点所在四面体的所有节点构成的集合就是最近邻集合，求出s
i-1
中所有节点的最近邻集合取并集，得到的集合记
作近邻集合e
i-1
，将该过程用近邻算子n来描述，则e
i-1
＝n(s
i-1
)，从e
i-1
中依据公式(3)，选出li个节点并入到支撑集s
i-1
，形成新的支撑集si，即si＝li(e
i-1
)∪s
i-1
，li(e
i-1
)表示从e
i-1
中选出li个节点；
[0034]
3.4残差更新：设定ri代表第i次迭代产生的残差向量，此外，下一次迭代所需的残差向量r
i+1
为：
[0035]ri+1
＝r
i-u
i+1
ꢀꢀꢀꢀ
(7)；
[0036]
其中，
[0037]
由于每次迭代需要选择li个节点，所以在第i次迭代需要构造li个线性无关向量
[0038]
3.5循环步骤3.1-3.4：每循环一次i＝i+1，直到残差ri《ε或稀疏度k的变化步长为1时，停止迭代；最终支撑集的前k列为解集，即
[0039]
x＝ki(si)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(8)；
[0040]
该解即为极小化问题的解。
[0041]
进一步的，利用tecplot成像软件对结果x进行展示。
[0042]
本发明具有以下优点：
[0043]
本发明适用于以重建荧光源为目标的光学分子断层成像，在重建过程中，利用提出的自适应近邻正交最小二乘算法进行重建，在保证精度的同时提高了空间分辨能力及鲁棒性，保证了重建结果的质量。
附图说明
[0044]
附图是用来提供对本公开的进一步理解，并且构成说明书的一部分，与下面的具体实施方式一起用于解释本公开，但并不构成对本公开的限制。在附图中：
[0045]
图1为基于本发明的基于自适应近邻正交最小二乘算法的光源重建方法流程图；
[0046]
图2为基于本发明的基于自适应近邻正交最小二乘算法的光源重建方法的实现流程图；
[0047]
图3为基于本发明的自适应近邻正交最小二乘算法的实现流程图；
[0048]
图4为基于本发明用于仿真实验的数字鼠模型；
[0049]
图5为基于本发明的fmt重建结果，其中深蓝色为重建的荧光团区域，为双目标的重建结果的三维显示；
[0050]
图6为本发明应用于fmt三维重建中，在使用不同数量的投影数据重建时获得的重建结果，双光源为z＝14mm截面的荧光目标重建分布图。
具体实施方式
[0051]
下面结合附图对本发明作进一步的详细描述，应指出的是，所描述的实施例仅旨在便于对本发明的理解，而对其不起任何限定作用。
[0052]
为了提高荧光分子断层成像的重建结果质量，本发明提出了一种基于自适应近邻的正交最小二乘算法来求解基于l0范数的正则化问题。该算法是典型的贪婪算法，它不需要事先知道变量的稀疏度，而是在迭代过程中国自适应的逼近稀疏度，同时采用近邻策略
选取最优支撑集，最终提高了荧光分子断层成像的重建质量和空间分辨率。具体流程包括：
[0053]
利用多模态光学分子影像系统，采集荧光数据和成像对象的解剖结构信息，对成像目标生物体进行有限元网格化，然后基于光传输模型的扩散近似方程，建立表面荧光数据与目标内部荧光团的线性关系。为了获得稳定稀疏的解，利用l0范数构建优化问题，并使用自适应近邻正交最小二乘法快速求解，最后展示重建结果。
[0054]
首先在多个位置激发获取多个角度的投影数据，以特定波长的激发光照射荧光探针标记过的生物体，用高灵敏度的光学探测器在体表采集荧光探针发出的荧光信号，构建基于有限元理论的扩散方程，建立荧光数据与重建目标内部荧光团分布的线性关系。
[0055]
φ＝αx
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1)；
[0056]
其中，φ为表面荧光数据，α表示系统矩阵，x是要求解的成像目标内部荧光团的分布。将上述线性关系转化为l0范数极小化问题：
[0057][0058]
||x||0是0范数，表示非零元素数量，ε为给定的阈值。本发明设计自适应近邻策略的正交最小二乘算法求解上述优化问题。
[0059]
本发明采用的技术方针如下：
[0060]
(a)初始化各参数：迭代次数i＝1，阈值ε，支撑集si,稀疏度k0＝6，以及l0＝10，初始残差r0＝φ，表面荧光光分布φ，
[0061]
元素选择标准qj。具体来说，在迭代中，它根据下式选择一个新的列索引jm：
[0062][0063][0064][0065]
其中i是列集合，aj是系统矩阵a的第j列，u
l
是正交基向量。
[0066]
根据上述公式(3)选出l0个元素构成初始支撑集s0。
[0067]
(b)自适应调整步长：采用非线性函数来调整步长，该函数是一个单调递减函数。步长表述如下：
[0068][0069]
其中i是迭代次数，s是步长，表示向上取整。设往支撑集si中添加的原子数为li，其迭代公式为li＝l
i-1-s，其中s0＝l0＝10，随着迭代次数的增加，向支撑集中添加的原子数逐渐变少，直到原子数变为1，则停止迭代。稀疏度迭代公式设为ki＝k
i-1
+s，其中s0＝k0＝6。随着迭代次数的增加，稀疏度逐渐增加，其增量开始阶段是较大步长，完成阶段是较小的步长，直到为1。
[0070]
(c)支撑集的更新：基于有限元理论，根据四面体单元结构构建最近邻节点集合。s
i-1
是第i-1次迭代生成的支撑集，设任意节点gk∈s
i-1
，则该节点所在四面体的所有节点构
成的集合就是最近邻集合，求出s
i-1
中所有节点的最近邻集合取并集，得到的集合记作近邻集合e
i-1
，将该过程用近邻算子n来描述，则e
i-1
＝n(s
i-1
)，从e
i-1
中依据公式(3)，选出li个节点并入到支撑集s
i-1
，形成新的支撑集si，即si＝li(e
i-1
)∪s
i-1
，li(e
i-1
)表示从e
i-1
中选出li个节点。
[0071]
(d)残差更新：设定ri代表第i次迭代产生的残差向量，此外，下一次迭代所需的残差向量r
i+1
为：
[0072]ri+1
＝r
i-u
i+1
ꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(7)；
[0073]
其中，
[0074]
由于每次迭代需要选择li个节点，所以在第i次迭代需要构造li个线性无关向量
[0075]
(e)循环步骤(b)-(d)，每循环一次i＝i+1，直到残差ri《ε或稀疏度k的变化步长为1时，停止迭代。最终支撑集的前k列为解集，即
[0076]
x＝ki(si)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(8)；
[0077]
该解即为极小化问题(2)的解。
[0078]
(f)利用tecplot成像软件对结果进行展示。
[0079]
实施例一：
[0080]
参见图1，本发明荧光断层成像重建方法的流程如下：
[0081]
(1)获取表面荧光测量数据，建立荧光数据与内部荧光目标分布的线性关系：
[0082]
(2)使用多模态成像系统，采集荧光数据和成像对象的解剖结构信息，
[0083]
(3)利用多模态光学分子影像系统，采集荧光数据和成像对象的解剖结构信息，对成像目标生物体进行有限元网格化，然后基于光子传输模型的扩散方程来近似模型，
[0084]
φ＝αx
[0085]
其中，φ为表面荧光数据，α表示系统矩阵，x是要求解的成像目标内部荧光团的分布。将上述线性关系转化为l0范数极小化问题：
[0086][0087]
||x||0是0范数，表达非零元素数量，ε为给定的阈值。
[0088]
(4)利用自适应近邻正交最小二乘算法得到重建结果，使用tecplot软件进行显示。
[0089]
参见图2为基于本发明的荧光分子断层成像重建实现的流程图。
[0090]
参见图3为基于本发明的自适应近邻正交最小二乘算法的实现流程图。
[0091]
结合图4-图6对基于本发明的重建结果做进一步的描述。
[0092]
图4是用于仿真实验的数字鼠模型。我们只选取了小鼠的躯干部分，同时分割出了小鼠的主要器官，包括心脏，肺，肝脏，胃，肾脏和肌肉组织。
[0093]
图5是基于本发明的fmt双光源不同边到边距离重建的结果，图中绿色的区域就是各阶段的重建结果。可以看到，双光源重建结果随着边边距的增大，重建效果越好。
[0094]
图6是本发明应用于fmt三维重建，双光源为6个投影点下的在z＝14mm截面的荧光目标重建分布图。双光源重建可在边到边距离为2mm至4mm时完成重建，本发明在贪婪策略
的基础上引入了非线性函数自适应策略来逼近真实稀疏度，并且使用近邻策略完成支撑集的精确选取，可以准确有效的进行重建，具有精度高和空间分辨率高的特点。
[0095]
以上结合附图选择优先详细论述最佳实施例，并不用于限制本发明。在上述描述的各个具体技术特征，在不矛盾的情况下，可以通过任意合适的形式进行组合，本发明不在一一赘述。任何本领域技术人员在不脱离技术方案范围内的前提下采取对技术方案进行任意组合或同等替换等简单修改或修饰的手段，并不影响其技术方案的本质仍属于本发明的各实施例代表的技术方案的保护范围之内。