超声层析成像的理论与实现.ppt

单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,超 声 层 析 成 像 的 理 论 与 实 现,答辩人:刘 超,指导教师:汪元美 教 授,浙江大学生物医学工程系,二三年九月二日,浙江大学博士论文答辩,英国从事超声成像的专家P.N.T Wells在2000年的文章,超声成像技术的现状与未来,一文中指出：“在最近的十几年里，有关超声成像技术的研究在医学成像领域至少占百分之二十五以上的份额，并且这种趋势还在继续增长Wells还指出：“目前成功地应用于医学领域的超声成像设备大都是基于反射波，且其成像也只是定性的，根据超声散射波的信息，定量地生成人体内部的结构图，是超声应用技术的研究者追求的新目标未来的超声成像技术应该是制造出不需成像专家或医学专家才能识别的反映客观现实真实图像的超声成像设备，即使是这种设备是不完美的主要内容,一.超声层析成像技术的发展历史,二.超声层析成像技术的基本模型及方法,三.问题的不适定性及其正则化,四.模型噪声的判断方法Picard准则,五.静态正则化技术在超声层析技术中的应用,六.迭代正则化技术在超声层析技术中的应用,七.总结与展望,一.,超声层析成像的发展历史,1.折射系数层析成像方法,2.衰减系数层析成像方法,3.射线跟踪方法,4.透射式衍射层析成像及反射式衍射,层析成像方法,5.基于精确场描述的层析成像方法,1.折射系数层析成像方法,Refractive-index tomography,2.超声衰减系数层析成像,Attenuation tomography,衰减系数,综合衰减系数,3.射线跟踪方法,Ray Tracing Method,4.透射式衍射层析成像及反射式衍 射层析成像方法,物体,傅里叶变换,频域,空域,入射波,前向散射场,从不同方向照射物体时,前向散射场数据的傅里叶变换,5.基于精确场描述的层析成像方法,二.超声层析成像技术的基本模型 及方法,非齐次亥姆霍兹方程,(Helmholtz Equation),1.波动方程及其解,全场方程,(Total Field Equation)(第二类Fredholm积分方程),散射场方程(Scattering Field Equation),探测器方程(Detector Equation),2.积分方程的离散化,矩量法,向量形式:,3.波动方程的近似,Born近似,Born逆解,O,应满足的条件:,Rytov近似,应满足的条件:,4.基本方法,Born迭代算法(BI),Levenberg-Marquardt和Newton-Kantorovich方法,变形Born迭代方法(DBI),Born迭代算法(BI),求Born逆解,O,由全场方程,确定全场,由散射场方程,求散射场,并计算,由方程,求改变量,求,求Born逆解,O,由全场方程,确定全场,由散射场方程,求散射场,并计算,变形Born迭代算法(DBI),根据最新求得的,O,k,改变散射方程的系数矩阵,D,求,由方程,求改变量,Levenberg-Marquardt和Newton-Kantorovich方法,代入,三.问题的不适定性及其正则化,适定性,问题是指：,对于连续算子方程,Kx,=,y,，如果解,x,满足：,(1).存在；,(2).唯一；,(3).连续地依赖于数据,y,。

否则，即上述三个条件有一个不满足，则称其为,不适定的（Ill-posed）,离散不适定问题,(Discrete Ill-Posed Problem),若:,(1).矩阵,A,的条件数非常大，或者说矩,阵,A,的最大奇异值和最小奇异值之比,非常大；,(2).矩阵,A,的奇异值逐渐下降趋于零对于线性方程组,Ax,=,b,或最小二乘问题:,Tikhonov正则化,L=I,n,，,x,0,=0,时，称为,Tikhonov,正则化的标准形式,其解可表示为:,四.模型噪声的判断方法,:Picard,准则,离散,Picard,准则：,若方程组,Ax,=,b,的傅里叶系数 趋于零的速度在平均意义下快于矩阵,A,的奇异值趋于零的速度的话，则称该方程组满足离散,Picard,准则(条件)最小二乘解:,Tikhonov正则化解:,受噪声污染和无噪声污染的,Picard,图,污染严重,污染较轻,A,对比度为30时,对比度为20时,对比度为10时,五.静态正则化技术,1.截断奇异值分解正则化方法,Truncated Singular Value Decomposition(TSVD),2.截断完全最小二乘正则化方法,Truncated Total Least Squares,(TTLS),1.截断奇异值分解正则化方法,(TSVD),对于线性方程组,Ax,=,b,或最小二乘问题,最小二乘解:,Tikhonov,正则化解:,TSVD,正则化解:,正则化参数的选取方法,离差原理(Discrepancy Principle)方法,广义交叉验证(GCV)方法,L曲线(L-Curve)方法,减小时,增加,时,由L曲线方法确定,k,采用一维搜索的方法确定更精确的,k,TSVD方法的数值仿真结果,B,A,C,D,E,对比度为10时,对比度为20时,对比度为30时,原始图像,迭代过程的相对误差和相对残差曲线,迭代过程的相对误差和相对残差曲线,2.截断完全最小二乘正则化方法,满足：,最小二乘问题：,完全最小二乘问题：,满足：,截断完全最小二乘的步骤,1.首先，计算增广矩阵,(,A,，,b,),的奇异值分解：,2确定截断参数,k,min,(,n,，rank(A，b),使得：,3.记,q,=,n-k,+1,将矩阵分块,4.则完全最小二乘问题的解为：,TTLS方法的数值仿真结果,原始图像,对比度为10时,对比度为20时,对比度为30时,迭代过程的相对误差和相对残差曲线,六.迭代正则化技术,1.求解最小二乘问题的共轭梯度方法,(cgls),2.LSQR方法,1.求解最小二乘问题的共轭梯度 方 法,(cgls),将共轭梯度法应用于法方程,相当于在,Krylov,子空间：,产生的序列,x,k,，使得：,cgls方法的解可表示为:,的,k,-1次多项式，其系数的确定,是,其中：,依赖于：(1).方程的右侧项,b,的特征；(2).矩阵,A,的奇,异值的分布；,(3).迭代的次数,迭代次数增加，残差变化不大，但解的范数受影响较大,正则化参数对迭代的影响,cgls方法的数值仿真结果,B,A,C,D,E,对比度为10时,原始图像,对比度为20时,对比度为30时,迭代过程的相对误差和相对残差曲线,图5.15采用clgs方法，五种不同图像在对比度为30时的,相对残差(RRE)曲线，cgls迭代次数为10,LSQR迭代方法,Lanczos,三对角过程,Lanczos,应用于,将矩阵,A,双对角化,Golub和Kahan(1965),Paige和Saunders(1982),线性方程组,Ax,=,b,和,应用于,LSQR方法的优点:,1.速度快,2.对不适定性问题数值稳定,3.从迭代过程很容易求得数值分析的数值,B,A,C,D,E,原始图像,LSQR方法的数值仿真结果,对比度为10时,对比度为20时,对比度为30时,七.总结与展望,首先利用Picard理论，分析了超声层析成像问题的中的模型噪声问题，给出了入射波的确定方法、以及正则化方法的适用范围的判断方法。

采用了两类四种正则化方法对超声层析成像问题中的不适定性问题进行了研究，通过对正则化参数选择的修正，完成了较大对比度物体的成像问题结论：静态正则化方法数值稳定，但速度慢；迭代正则化方法速度快，但数值稳定性不如静态方法今后需要进一步研究的工作,1.前向散射问题的研究（波动方程的精确程度）2.离散化方法 有限元法、边界元法 矩量法中基函数的确定3.正则化问题 基于非对成方程的,Krylov,子空间方法,谢 谢！,。