[发明专利]组织结构引导的复合正则化生物发光断层成像重建方法

摘要

组织结构引导的复合正则化生物发光断层成像重建方法，本方法属于医学图像处理领域；针对现有生物发光断层成像重建算法中存在的上述问题，本方法可以在已知肿瘤所在的生物组织而不知具体位置时很好的实现光源的重建；另外在不知肿瘤所在生物组织时，通过均分权重系数来实现重建。首先采用有限元方法仿真光在生物组织中的传输，并基于生物发光荧光断层成像中光源的稀疏分布特性，采用结构成像引导的复合正则化方法融入更多先验信息以降低BLT重建的病态性，最后在此基础上实现荧光光源的准确重建。在两个波段设置了双光源进行仿真实验，以验证算法的有效性和准确性。本方法不但可以对荧光光源进行准确重建定位，而且可以极大地提高计算效率。

主权项

1.组织结构引导的复合正则化生物发光断层成像重建方法，其特征在于：首先采用有限元方法仿真光在生物组织中的传输，并基于生物发光荧光断层成像中光源的稀疏分布特性，采用结构成像引导的复合正则化方法融入更多先验信息以降低BLT重建的病态性，最后在此基础上实现荧光光源的准确重建；具体包括下述步骤：步骤一，光在生物组织中的传输规律用扩散方程加以表示；根据变分原理，结合Robin边界条件，建立与扩散近似方程和边值问题等价的弱形式： ∫ Ω ( ( D ( r ) ( ▿ Φ ( r ) ) · ▿ Ψ ( r ) ) + μ a Φ ( r ) Ψ ( r ) ) d r + ∫ ∂ Ω Φ ( r ) Ψ ( r ) / ( 2 I n ( r ) ) d r = ∫ Ω S ( r ) Ψ ( r ) d r - - - ( 1 ) ]]>其中D(r)为扩散系数分布，Φ(r)为光子密度分布，μa为光吸收系数，Ψ(r)为任意测试函数，S(r)为光源分布，In(r)为使用Robin边界条件时引入的一个量，Rn≈-1.4399n-2+0.7099n-1+0.6681+0.636n表示扩算传输内反射系数，n是一个常数，与边界内外的光学折射系数偏差有关；步骤二，有限元离散；对成像体Ω进行网格剖分，实现离散化；本方法采用二维仿体进行仿真，由于三角网格单元能够实现使用较少数量来近似描述不规则区域，所以，本方法使用三角网格剖分成像物体；根据对近似解精度的要求、单元的几何形状、节点个数和节点的自由度信息来选择单元基函数；由于函数均表示为基函数的线性组合，假设成像物体所划分的网格包含Nn个节点和Nt个三角单元，光子密度分布Φ(r)表示为：其中，φk表示节点i上的函数值，表示节点i处的基函数；同样，光源函数S(r)表示为：其中，sk表示节点i上的光源值；步骤三，合成总体刚度矩阵；将(2)、(3)式代入(1)式中，得到：该矩阵方程简单表示为：([K]+[C]+[B])＝MΦ＝FS (5)其中，M为刚度矩阵，其他各矩阵元素通过下面公式计算得到：整理(6)式，得到：M-1FS＝Φ (7)进一步地，上式写成如下形式：AS＝Φ (8)其中，A＝M-1F；对于生物发光断层成像，探测器所能够采集到的测量数据仅仅是成像物体表面的荧光信息，因此在(7)式中，需要将光子密度向量Φ中的非表面元素去掉，并将系统