一种荧光分子断层成像可行域优化方法技术

本发明专利技术属于光学分子成像技术领域，公开了一种荧光分子断层成像可行域优化方法；使用有限元计算软件对仿体进行四面体网格划分，建立表面测量值与内部荧光目标分布的线性关系；利用不完全变量截断共轭梯度法、共轭梯度最小二乘法求解系统方程，得到两组重建结果；根据三支决策理论划分出两组重建结果的正域，合并构成目标可行区，作为下一级重建的区域；再次使用不完全变量截断共轭梯度法、共轭梯度最小二乘法进行重建，得到改进后的最终结果。本发明专利技术采用两种不同的重建算法进行重建；提取出两组结果中的正域；在此区域内完成最终重建。本发明专利技术能够获取到准确的目标范围区域，有效缓解了问题的病态性，上述两种重建方法的图像质量均得到明显提升。

A feasible region optimization method for fluorescence molecular tomography

The invention belongs to the field of optical molecular imaging technology, and discloses an optimization method for the feasible region of fluorescent molecular tomography; uses finite element calculation software to divide tetrahedral meshes of the imitated body, establishes the linear relationship between the surface measurement value and the distribution of internal fluorescent targets; uses incomplete variable truncated conjugate gradient method and conjugate gradient least square method to solve the system equation, and obtains two groups. According to the three-branch decision theory, two groups of positive domains of reconstruction results are divided, and the target feasible area is merged to form the next reconstruction area. Then the incomplete variable truncated conjugate gradient method and conjugate gradient least squares method are used to reconstruct, and the improved final results are obtained. The invention adopts two different reconstruction algorithms to reconstruct, extracts the positive domain of two sets of results, and completes the final reconstruction in this region. The invention can acquire an accurate target range area, effectively alleviate the ill-conditioned nature of the problem, and the image quality of the two reconstruction methods is obviously improved.

【技术实现步骤摘要】
一种荧光分子断层成像可行域优化方法
本专利技术属于光学分子成像
，尤其涉及一种荧光分子断层成像可行域优化方法。
技术介绍
目前，业内常用的现有技术是这样的：作为光学分子成像技术的重要分支，荧光分子断层重建(FMT)借助特定波长的激发光照射被荧光标记物的生物体，然后基于全角度非接触式成像系统，应用高性能CCD采集生物体表面的光强分布，结合合适的数学模型反演出荧光标记物在生物体内的三维分布和荧光浓度，该技术能够在分子水平上实现对生物体生理过程和病理变化的在体动态监测。荧光分子断层重建的成像灵敏度高，成本低以及荧光探针标记技术发展较为成熟，从而使其在肿瘤早期检测、药物研制和治疗评估等方面具有很大优势。荧光分子断层重建主要涉及两大问题：正向问题和逆向问题。由于光在生物组织中传播时具有高散射、低吸收的特点，因此正向问题的主要任务为建立数学模型来描述光在生物组织中的传输过程，通常使用辐射传输方程(RTE)来描述这一过程。但是在实际使用中，辐射传输方程的求解十分复杂，在大多数情况下无法求解出解析解，因此一般使用漫射方程来近似逼近辐射传输方程，然后使用有限元这一数值方法求解漫射方程，获得每个节点的能量分布值。逆向问题则是指基于正向问题建立的光传输模型，选择一定的重建算法，求解出荧光标记物在生物体内的位置和浓度分布信息。由于在逆问题求解过程中，未知数量远远大于方程数量，因此逆问题具有严重的不适定性。目前主要围绕降低逆问题病态性，更加精确有效的定位荧光目标的位置及分布。但是基于传统方法，无法一次性精确定位到荧光目标的准确位置，为了提高图像质量，减少计算过程中的数据量，通常可以借助可行域，基于经典重建算法，可以达到改进重建方法的目的。采用传统重建算法得到的结果中，虽然能够逼近荧光目标，但通常无法给出荧光标记物的精确分布，这使得重建结果在后续研究中的参考意义大打折扣，降低了荧光分子断层重建技术的应用价值。通常在使用目标可行区提高重建质量时需要参考初次重建结果，因此初次重建结果的精确度决定着可行域的选取范围，如果未能准确提供荧光目标存在的范围，将造成下一级重建结果失准，继而无法得到更加准确的目标定位。综上所述，现有技术存在的问题是：采用传统重建算法得到的结果无法给出荧光标记物的精确分布，使得重建结果在后续研究中的参考意义大打折扣，降低了荧光分子断层重建技术的应用价值；在小动物全身成像、人体乳腺器官成像的临床及预临床研究中无法提供有效的重建结果。解决上述技术问题的难度和意义：借助三支决策理论的思想可以将重建区域划分为三部分：正域、负域及边界域，设置初始阈值，根据各节点的能量值，筛选出处于正域中的节点构成可行域，作为下一次重建的区域。在这个过程中，去掉荧光能量值较小的节点，保留能量值较大的节点进行二次重建，有效地去除了噪声数据，能够进一步提升可行域的准确度。最后使用上述两种方法进行二次重建时，本专利技术能够有效提升两种方法的重建精度。在本专利技术中，如何合理地设置阈值参数使决策风险损失最小一直是三支决策理论的挑战性难题，选择正域能够尽可能地保证重建范围不会存在较大偏差，接近目标所在区域。在使用其指导下一级重建时，能够提供准确的目标范围依据，有利于提高重建质量。
技术实现思路
针对现有技术存在的问题，本专利技术提供了一种荧光分子断层成像可行域优化方法。本专利技术是这样实现的，一种荧光分子断层成像可行域优化方法，所述荧光分子断层成像可行域优化方法包括：第一步，使用有限元计算软件对仿体进行四面体网格划分，建立表面测量值与内部荧光目标分布的线性关系；第二步，利用不完全变量截断共轭梯度法、共轭梯度最小二乘法求解系统方程，得到两组重建结果；第三步，根据三支决策理论划分出两组重建结果的正域，合并构成目标可行区，作为下一级重建的区域；第四步，再次使用不完全变量截断共轭梯度法、共轭梯度最小二乘法进行重建，得到改进后的最终结果。进一步，所述第一步具体包括：(1)基于全角度非接触型成像系统，应用高性能CCD相机采集所有角度的荧光数据；(2)建立圆柱仿体模型，使用有限元网格剖分软件COMSOLMultiphysics对仿体进行离散化，将成像区域Ω划分为m个互不重叠的四面体单元(Ω1,Ω2,...,Ωm)和n个网格节点(p1,p2,...,pn)；第i个结点表示为pi，1≤i≤n，第i个节点的荧光产额值表示为Yi；(3)建立表面测量值与荧光目标分布的线性关系，构建荧光分子断层重建问题的系统方程；表示边界测量值，A为系统矩阵，X表示所要求解的荧光目标的分布。进一步，所述第二步具体包括：基于同一网格，分别使用不完全变量截断共轭梯度法、共轭梯度最小二乘法进行重建，得到两组全域重建结果，各节点的荧光产额值，分别记做：进一步，所述第三步具体包括：(1)设置初始阈值α(0＜α＜1)，根据三支决策理论，筛选出两组重建结果output1，output2中属于正域的节点集合，分别记为：PS_R1＝{pi|Y1i≥max(Y1i)*α,1≤i≤n}；PS_R2＝{pi|Y2i≥max(Y2i)*α,1≤i≤n}；(2)合并两组重建结果的正域，对PS_R1，PS_R2取并集，得到荧光目标的可行域PR＝PS_R1∪PS_R2，根据可行区PR定义一个列矩阵D＝[P1,P2,...,Pn]T：矩阵D表示若第i个结点若属于可行域，则标记为i，否则置为0，去除D中的零值，得到矩阵D'。进一步，所述第四步具体包括：(1)基于形成的目标可行区D'，产生新的系统矩阵A'，利用正则化方法，将重建问题转化为拉普拉斯正则化的l1范数极小化问题：其中τ是正则化参数，Φ'为表面的荧光测量值，X表示要求解的荧光目标三维分布与浓度；(2)再次使用不完全变量截断共轭梯度法、共轭梯度最小二乘法进行重建，得到改进后的最终结果。进一步，所述荧光分子断层成像可行域优化方法的定量分析指标包括定位误差、归一化方差误差；定位误差表示重建出的荧光目标位置与真实荧光目标位置之间的误差，归一化方差误差代表重建的荧光产额值与真实荧光产额值之间的误差；其中(x,y,z)为重建光源的中心位置，(x0,y0,z0)为真实光源的中心位置；xrec表示重建结果的荧光产额，xorg为真实的荧光产额；若定位误差和归一化均方差误差越接近于零值，则表示重建结果更加精确。本专利技术的另一目的在于提供一种应用所述荧光分子断层成像可行域优化方法的光学分子成像平台。本专利技术的另一目的在于提供一种应用所述荧光分子断层成像可行域优化方法的荧光分子断层重建平台。本专利技术的另一目的在于提供一种应用所述荧光分子断层成像可行域优化方法的切伦科夫发光断层成像平台。综上所述，本专利技术的优点及积极效果为：本专利技术的主要应用领域包括：荧光分子断层重建、生物发光断层重建及契伦科夫断层成像等。受三支决策理论启发，本专利技术采用两种不同的重建算法进行重建，获得初始重建结果，然后提取出两组结果中的正域，合并构成目标可能存在的范围，最后在此区域内完成最终重建，能够获取到更加准确的目标范围区域，有效缓解了问题的病态性，提升了下一级重建图像的质量。最后的圆柱仿真结果也表明通过本专利技术提取的可行域有助于上述两种方法的重建，有效改进了两组实验的重建精度。圆柱仿真实验结果中：IVTCG方法的位置误差降为0.67mm，本文档来自技高网...

【技术保护点】
1.一种荧光分子断层成像可行域优化方法，其特征在于，所述荧光分子断层成像可行域优化方法包括：第一步，使用有限元计算软件对仿体进行四面体网格划分，建立表面测量值与内部荧光目标分布的线性关系；第二步，利用不完全变量截断共轭梯度法、共轭梯度最小二乘法求解系统方程，得到两组重建结果；第三步，根据三支决策理论划分出两组重建结果的正域，合并构成目标可行区，作为下一级重建的区域；第四步，再次使用不完全变量截断共轭梯度法、共轭梯度最小二乘法进行重建，得到改进后的最终结果。

【技术特征摘要】
1.一种荧光分子断层成像可行域优化方法，其特征在于，所述荧光分子断层成像可行域优化方法包括：第一步，使用有限元计算软件对仿体进行四面体网格划分，建立表面测量值与内部荧光目标分布的线性关系；第二步，利用不完全变量截断共轭梯度法、共轭梯度最小二乘法求解系统方程，得到两组重建结果；第三步，根据三支决策理论划分出两组重建结果的正域，合并构成目标可行区，作为下一级重建的区域；第四步，再次使用不完全变量截断共轭梯度法、共轭梯度最小二乘法进行重建，得到改进后的最终结果。2.如权利要求1所述的荧光分子断层成像可行域优化方法，其特征在于，所述第一步具体包括：(1)基于全角度非接触型成像系统，应用高性能CCD相机采集所有角度的荧光数据；(2)建立圆柱仿体模型，使用有限元网格剖分软件COMSOLMultiphysics对仿体进行离散化，将成像区域Ω划分为m个互不重叠的四面体单元(Ω1,Ω2,...,Ωm)和n个网格节点(p1,p2,...,pn)；第i个结点表示为pi，1≤i≤n，第i个节点的荧光产额值表示为Yi；(3)建立表面测量值与荧光目标分布的线性关系，构建荧光分子断层重建问题的系统方程；表示边界测量值，A为系统矩阵，X表示所要求解的荧光目标的分布。3.如权利要求1所述的荧光分子断层成像可行域优化方法，其特征在于，所述第二步具体包括：基于同一网格，分别使用不完全变量截断共轭梯度法、共轭梯度最小二乘法进行重建，得到两组全域重建结果，各节点的荧光产额值，分别记做：4.如权利要求1所述的荧光分子断层成像可行域优化方法，其特征在于，所述第三步具体包括：(1)设置初始阈值α(0＜α＜1)，根据三支决策理论，筛选出两组重建结果output1，o...

【专利技术属性】
技术研发人员：宋小磊，焦璞，易黄建，梁嘉翔，赵凤军，曹欣，贺小伟，侯榆青，
申请(专利权)人：西北大学，
类型：发明
国别省市：陕西,61