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

来源 :浙江大学理学院 浙江大学 | 被引量 : 2次 | 上传用户:erikwg
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
近年来,随着科学技术的发展,生物医学成像得到了长足的进步.分子成像是在细胞和分子水平上研究活体生理和病理过程的一个快速发展着的生物医学成像方法.它的目标是非入侵地准确描述活体细胞和分子过程.分子成像主要基于三种技术:核素成像技术,核磁共振成像技术和光学成像技术.与传统的成像方法,如计算机层析成像(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)函数空间.我们证明了解的存在性和数值解的收敛性并给出了一个数值求解格式.   最后,在第五章,我们作了总结并展望未来进一步要做的工作.
其他文献
本论文主要包括以下三个方面的内容.  首先,利用非协调EQrot1元以及零阶Raviart-Thomas(R-T)元,我们研究了双相滞热传导方程的两类不同的全离散有限元格式.其一是构造它的一
多面体样条是B样条的推广,有着诸如局部支集性,光滑性,灵活性等非常好的性质.鉴于规则网格数据的曲面设计方法无法满足实物的复杂性需求,而以细分法为代表的任意网格数据的曲
背景误差协方差矩阵(B)是变分同化系统中一个很重要的组成部分。它反映了不同变量之间的关系,并对同化系统输出的分析场有极大地影响。如何及时吸收观测信息并有效更新当前的
本文主要研究半环的结构和同余,探讨了半环强分配格的性质及同余,给出了双半环的拟强分配格的定义,并得到了双半环的拟强分配格的结构分解,找到了一些半环的半格同余及分配格同余
将库存控制与价格确定联合起来加以讨论的问题已经被很多专家学者进行了研究,同时,允许顾客提前订货的库存管理策略也在被人们深入的研究,并且已经被很多企业进行了实践。对
在Lotka—Volterra竞争模型和相似的竞争模型中,如果种群间的竞争比种群内部的竞争弱,那么种群之间才能共存。先验模型和实地研究表明:通过空间或时间的生态变异来减弱种群间的