一类反问题—生物体发光断层成像—的理论分析和数值模拟

论文摘要

近年来,随着科学技术的发展,生物医学成像得到了长足的进步。分子成像是在细胞和分子水平上研究活体生理和病理过程的一个快速发展着的生物医学成像方法。它的目标是非入侵地准确描述活体细胞和分子过程。分子成像主要基于三种技术:核素成像技术,核磁共振成像技术和光学成像技术。与传统的成像方法,如计算机层析成像(CT)、正电子发射断层成像(PET)、核磁共振成像(MRI)以及它们的合成等,相比,作为光学成像模型,生物放光成像(BLI)有它独特的优点。对于活体人体研究,由于组织的吸收与发散特性,光学成像主要局限于表层位置,此时MRI或PET是比较合适的选择。然而,对于小动物来说,由于较短的路径长度,一大部分光子可以到达动物的外表,从而能以比MRI和PET低得多的成本探测到生物光和荧光信号。对于活体应用来说,用生物光源进行成像特别具有吸引力,因为它并不需要外部的刺激光源,从而背景噪音低而灵敏性高。生物发光断层成像(BLT)是一个非常有前途的BLI,因为他能实时地揭示生物体的分子和细胞活动。几年来,人们就BLT问题做了大量的理论分析和数值仿真。BLT的主要任务就是利用边界测量的光信号来决定小动物体内或一些大器官的表层的光子密度分布。利用它,我们能在三维空间中对小动物内部的光源进行定位并定量计算出生物光源的强度密度分布。进行BLT的第一步是确定组织的光学特性,这正是光学断层摄影(DOT)问题的的主要任务。因为光在任何实体中的传送都受到吸收和散射的影响,所以我们需要知道光子在生物组织中传送的精确表达式。一般来说,生物光光子在组织中的传播过程能被漫射转移方程(RTE)或蒙特卡罗模型(MCM)准确描述。然而,就当前而言,两者在计算上都是不可行的。因此,当光的波长在600 nm附近或更大时,人们通常用扩散近似方程来代替RTE。从数学上来讲,BLT是一个反问题,通常需要对它进行正则化处理。在这篇论文中,我们讨论了BLT和多谱BLT的理论分析和数值求解等问题.全文共分五章。在第一章里,我们首先简单介绍一下BLT问题的生物医学背景。然后,在1.2节中,我们给出了BLT问题的数学形式。在1.3节,我们给出了一些主要符号,它们将会在接下来的章节里被反复地用到.第二章主要是关于经典BLT的理论和数值方法的研究。2.1节陈述了经典BLT的形式。然后,在2.2节里,在一些假定下,利用共轭方法,我们将一个正则化的最小化问题转化为一个偏微分方程组。这为快速重构光源函数提供了可能。在第2.3节中,利用正则化方法,我们研究了BLT问题的一个更为广义框架。在这个框架下,我们建立了问题的适定性结果,考察了它与以前别人文章之间的关系。在第2.4节中,我们发展了BLT问题的一个新的变分形式。利用该变分形式,我们可以从理论上证明光源函数并不连续依赖于测量数据。在这一节里,我们还给出了正则化有限元解最优收敛阶的证明。另外,利用伴随方程,我们构造了一个简单而有效的迭代格式。第三章研究的是多谱BLT问题。首先,在3.2节,基于一个惩罚策略,我们提出了关于多谱BLT的全新的方法。这个新的数学框架的特点是其数值测量的使用涉及到两个相关但并不相同的边值问题(BVPs)。该数学框架包含了通常的多谱BLT模型。然后,在第3.2节,我们改进[1]中的误差阶。接着,在第3.3.2小节里,我们用线性而不是通常的常数有限元来离散光源函数空间.利用有效集策略,我们获得更高阶的有限元误差估计.第四章我们继续深入讨论BLT问题。首先,在4.2节,我们研究了具有随空间改变折射率的BLT问题。为此,利用Tikhonov正则化方法,我们介绍一个广义的框架,证明了它的适定性,并利用有限元方法(FEMs),建立了数值解的误差估计式。在第4.3节,我们提出了一个数学模型,它将BLT和DOT结合在一起,也就是说,同时进行两种类型的重构而不是一前一后。因为在实际中,光学参数通常被假定是常数或分片光滑函数,因此我们设定光学参数的求解范围是有界变差(BV)函数空间。我们证明了解的存在性和数值解的收敛性并给出了一个数值求解格式。最后,在第五章,我们作了总结并展望未来进一步要做的工作。

论文目录

• 致谢
• 摘要
• Abstract
• 目录
• 第一章 绪论
• 1.1 生物医学背景
• 1.2 BLT的数学模型
• 1.3 基本记号
• 第二章 BLT经典模型
• 2.1 引言
• 2.2 基于伴随方程的一个快速数值方法
• 2.2.1 问题的变形与正则化
• 2.2.2 微分方程组形式
• 2.2.3 有限元离散
• 2.2.4 数值例子
• 2.3 一个新的广义框架
• 2.3.1 广义格式的提出
• 2.3.2 极限情况
• 2.3.3 有限元近似
• 2.3.4 数值例子
• 2.4 基于变分不等式的BLT模型
• 2.4.1 BLT问题的不稳定性证明
• 2.4.2 一个新的正则化格式
• 2.4.3 有限元近似和误差估计
• 2.4.4 数值仿真
• 2.5 本章小节
• 第三章 多谱BLT模型
• 3.1 引言
• 3.2 基于罚方法的广义多谱BLT
• 3.2.1 新模型问题的提出
• 3.2.2 惩罚、正则解的适定性与极限行为
• 3.2.3 有限元近似
• 3.2.4 数值例子
• 3.3 多谱BLT的有限元误差分析
• 3.3.1 变分不等式形式
• 3.3.2 有限元离散及误差的改进
• 3.3.3 用线性有限元离散光源函数与误差估计
• 3.4 本章小节
• 第四章 BLT的进一步讨论
• 4.1 引言
• 4.2 折射率随空间改变的BLT
• 4.2.1 模型的描述及解的性质
• 4.2.2 有限元近似
• 4.2.3 数值实验
• 4.3 BLT与DOT的合成模型
• 4.3.1 基本知识
• 4.3.2 模型建立与解的存在性
• 4.3.3 数值近似
• 4.3.4 拟牛顿迭代算法
• 4.3.5 数值列子
• 4.4 本章小节
• 第五章 总结与展望
• 5.1 总结
• 5.2 展望
• 参考文献
• 发表文章目录
• 简历