Superiorization算法的理论研究和改进及其在CT和MRI中的应用

来源 :河南大学 | 被引量 : 0次 | 上传用户:usaend
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
图像重建技术是一种可以获得物体内部结构或者功能的图像而不对物体造成损伤的技术。目前,不完全数据重建问题和快速成像问题是图像重建技术的研究热点,这些问题都可以归结为不适定问题,研究这类问题的求解算法具有十分重要的应用价值,正则化方法是解决这类不适定问题的常用的方法,已经在医学、工业等领域得到了广泛的应用。本文主要基于不完全数据重建问题和快速成像问题,研究了一种求解大规模约束正则化模型的superiorization算法,提出了两种新的改进算法,并证明了第一种改进算法的收敛性。本文首先简单介绍了医学图像重建、正则化重建技术以及superiorization算法的研究背景、意义和国内外研究现状;其次介绍了计算机断层成像(computed tomography,CT)、核磁共振成像(magnetic resonance imaging,MRI)和并行核磁共振成像(partially parallel imaging,PPI)的基本原理、重建模型、重建算法及其相关的正则化重建模型。然后,本文对superiorization算法进行了研究。Superiorization算法的思想是在迭代过程中加扰动参数来求解约束最小化模型的近似最优解。对于该算法,本文从理论分析和大量的数值实验发现了其存在以下一些问题:原算法需要使迭代结果满足目标函数下降和在约束条件范围内,如果满足这两个迭代条件,则迭代进行一次,否则调整参数的大小,直到满足条件,这不仅耗时而且影响迭代结果;另外,由于在迭代后期扰动参数变得很小起不到正则化的作用,这使得重建时间增加,且影响最终的迭代结果。为了解决这些问题,本文提出了两种改进算法。两种改进算法为:(1)最优参数superiorization算法—–改变参数的选择方法,其主要思想是用无约束最小化模型的解代替原算法扰动后的迭代结果,那么替换后的结果一定可以满足目标函数下降,我们就不需要再验证这一条件,从而节省大量时间、提高重建质量;并证明了改进算法的收敛性。(2)阈值约束的superiorization算法,其思想主要是在迭代过程中对扰动参数进行阈值限制,如果扰动参数小于这一阈值,令扰动参数等于这一阈值,这就使得正则化在整个迭代过程中都起到了作用,从而节省了成像时间并提高了重建质量。此外,本文还证明了基于PPI成像的投影迭代算法满足BPR(Bounded perturbation resilience,BPR)条件,说明superiorization算法同样适用于PPI成像。最后,将改进算法应用到基于全变分(total variation,TV)的正则化重建模型的CT和PPI图像重建中,分析、比较模拟仿真实验结果,实验结果表明改进后的superiorization算法优于原算法,不仅提高了重建速度,而且提高了重建图像质量。
其他文献
本文主要研究了有理函数的双曲维数和Poincaré指数之间的关系以及TCE条件和大导数条件的性质。  第一章简要介绍了复解析动力系统的发展历史和研究的主要内容,并且介绍了本
本文主要讨论双曲型积分微分方程的最小二乘混合有限元方法,同时利用最小二乘混合有限元方法与特征线方法相结合的技巧研究Sobolev方程。最小二乘混合元方法可以避免标准混合
基因芯片技术是一种高通量的基因表达分析平台,因此基因表达数据的分析被广泛应用于生物医学研究、疾病诊断和药物筛选等多个领域。但是在基因表达数据的分析过程中,经常存在着
本文主要运用了拟线性化方法和广义拟线性化方法讨论了不同类型的积分微分方程的解的收敛性,达到了提高解的收敛速度的目的。全文其分五章。   第一章简述了积分微分方程系
图像处理广泛地运用于现代生活的各个领域.图像复原是图像处理的基础步骤,图像分割是从图像处理到图像分析的关键步骤.因此,改进图像复原和分割的方法,提出相应的模型和算法
本文考虑高维空间的区域中的具有两个特征矩阵的Beltrami方程组Dtf(x)H(x)Df(x)= J(x,f)2/nG(x)这里H(x)= diag(21 (x),22(x), ···,2n(x))为对角阵,detH(x)= 1,and G(x)∈ S(n)为正
本文主要研究的是仿射Weyl群E8a-值等于6的双边胞腔W2⑹中的左胞腔,找出了双边胞腔W2⑹中的三组左胞腔代表元系,画出了它们的左胞腔图,并给出了W2⑹中所含的特异对合元。    
本论文主要研究结构动力模型修正中一类二次特征值反问题的中心对称解及其最佳逼近.设,如果,则称A为n阶中心对称矩阵.本文研究的一类二次特征值反问题具体是指构造n阶中心对称矩