基于变分法的医学超声图像去噪研究

1、基于变分法的医学超声图像去噪研究【摘 要】目的：医学超声图像中斑点噪声的存在， 降低了图像质量，本文着重讨论医学超声图像的去噪问题。 方法：针对超声图像的斑点噪声，本文基于全变分正则化模 型，首先对超声图像进行对数变换，将乘性斑点噪声转化为 加性噪声，再对对数变换后的图像进行全变分正则化处理， 最后通过指数变换重构超声图像。结果：子宫超声图像去噪 实验中，将全变分法与常用的中值滤波和小波变换去噪方法 进行对比，结果显示全变分法的去噪性能指标明显优于其余 方法。结论：采用基于全变分正则化的方法，不仅很大程度 上抑制了医学超声图像的斑点噪声，而且保留了清晰的边缘 细节信息，具有重要的学术价值和现实

2、意义。【关键词】超声图像处理；斑点噪声去噪；全变分正则 化模型0引言与X光透视、CT、MRI等医学成像方法相比，医学超 声成像因非侵入无创伤性、成像速度快、成本低、操作简便 等优点，成为目前普遍应用的医学成像技术。超声诊断作为 一种理想的无损检查方法，有着广阔的发展前景。据报道， 近十年，世界医学超声仪器的数量以15%左右的速度增长1o然而，由于成像机制的限制，超声图像存在固有的斑点 噪声，极大地降低了超声图像质量，增加了图像特征分析的 难度，影响了疾病诊断的准确度。因此，超声图像中对斑点 噪声的抑制具有重要的学术价值和现实意义。目前，超声图像去噪方法常见的有中值滤波2、直方图3和小波变换4等

3、。基于小波域的去噪方法，以其良好的时频特性，广泛运用于超声图像去噪。该方法主要基于图像 中有用信息和噪声之间的频率特性存在差异的假设，进行频 域分析去噪。但实际上假设条件并不总是成立，图像中的有 用信息部分和噪声往往在频带上存在重迭。基于小波域的去 噪方法容易丢失部分高频分量一一图像中的细节和边缘等 有用信息，限制了图像质量的提高。最近的研究表明，全变 分正则化法对稀疏或梯度稀疏图像的重构效果显著，很好地 保留了图像的边缘信息。Rudin等人首次将全变分去噪法引入到图像处理中5 o由于全变分正则化在边缘检测中的巨大 应用价值，近几年研究人员提由很多基于全变分的去噪算法 6-7 o针对超声图像的

4、去噪特殊性，本文提由基于全变分正 则化的超声图像去噪方法。首先对超声图像进行对数变换， 将乘性噪声转化为加性噪声，再对对数变换后的图像进行全 变分正则化处理，最后通过指数变换重构超声图像。本方法 去除斑点噪声的同时，能够很好地保留图像的边缘信息。1变分法原理变分法本质上为泛函的极值问题。设 F (x, y (x), y'(x)是独立变量 x, y (x), y' (x)在区间x0 , x1上的已知函数，且二阶连续可微。变分法的基本问题为：在边界 条件为y (x0) =y0, y (x1) =y1时，求使得泛函这样求泛函极值的问题转化为求解在满足边界条件y(x0) =y0, y

5、(x1) =y1下的定解问题。2基于变分法的超声图像去噪超声图像一般由原始图像和乘性噪声构成，但乘性噪声 不便于处理。对超声图像进行对数变换，可将其中混有的乘 性噪声转变为加性的高斯噪声 n8,即u0=u+n (3)其中U0和u分别为混噪图像和原始图像的对数变换。基于全变分的图像去噪，目的是要还原图像U,使得它的全变分值小于U0的全变分，但又尽量接近 U0O该方法基 于含噪声图像的变分值显著高于无噪声图像的变分值，使用 最优化的方法使含噪声图像的变分值最小化就可以降低图 像的噪声。基于全变分正则化的图像去噪模型为上面的能量泛函方程由两项组成：第一部分即为正则 项，是对图像能量的估计，起到平滑的

6、作用；第二部分是保 真项，表示被噪声污染的图像和原图像图像之间的距离，使 含噪声图像与去噪后图像距离不被缩小，即是对高频突变部 分图像的补偿。参数 人起到平衡正则项与保真项的作用，随着人的增大，求生的u越接近U0O该能量泛函对应的欧拉- 拉格朗日方程为：其中，1/?萃U为扩散系数，表示朝图像中各个方向的 扩散能力。从该方程可看由，在边缘处，？萃 U较大，扩散 系数较小，因此沿边缘方向的扩散较弱，从而保留了边缘； 在平滑区，？萃U较小，扩散系数较大，因此在平滑区的扩 散能力较强，从而去除了噪声。全变分正则化法实际上是对 图像的高频边缘部分进行了补偿，因此该去噪声方法具有良 好的边缘细节保留效果。

7、3实验与分析本节将全变分法应用于超声图像去噪实验分析，并与中 值滤波、小波变换方法去噪方法进行效果对比。图1 (a)为原始含噪的子宫超声图像，图像大小为256X 256像素，取自剑桥大学网页9。图1 (b)为中值滤波去噪结果，图1 (c) 为小波去噪结果。可以看由，中值滤波和小波变换虽在一定 程度上达到去噪效果，使噪声部分平滑，但导致图像高频部 分丢失，图形边缘特征不明显，对目标区域的可分辨性较低。 图1 (d)为利用本文的全变分正则化方法进行去噪处理，可 以看由，采用本文的方法不仅有效地去除了医学超声图像的 斑点噪声，而且增强了目标区域内轮廓边缘的清晰度，在分 辨效果上体现了明显的优越性。(

8、a)原始含噪图像(b)中值滤波(c)小波变换(d)全变分法图1去噪前后的子宫超声图像在该实验中，以超声图像中目标区域的均值与标准差之 比(the mean-to-standard-deviation ratio ,简称 MSR)来作为 衡量去噪图像质量的性能指标。MSR值越大，去噪效果越好，图像质量越高。从图1中，截取感兴趣的目标区域：方形区 域1和矩形区域2。表1所示为不同去噪方法所得图像目标 区域内的MSR值，通过比较可知全变分法的MSR值最高。表1不同去噪方法所得目标区域的MSR比较注：最佳结果由黑体表示.4结论医学超声图像中人体组织结构明确的边界是医生诊断 的重要依据，因此在去除斑点噪

9、声的同时须保留边界信息。 相对于中值滤波、小波变换等图像去噪方法，全变分正则化 去噪方法更适合用于医学图像处理，它能够在去除斑点噪声 的同时很好地保留高频的图像边缘信息。本文从理论推导和 实验分析两个角度，验证了全变分法应用于超声图像去噪的 可行性。但医学超声图像中斑点噪声的形成受多因素影响， 为了更加深入地分析全变分正则化方法在医学图像中的应 用效果，还需要大量的临床实验来进一步完善该方法，这是 今后需要研究的方向。【参考文献】1邢立伟.心脏B超的医学图像处理D.西安：西北工业大学，2002: 12-132关新平，赵立兴&唐英干.图像去噪混合滤波方法J. 中国图象图形学报，2005,

10、 10: 332-337.3刘培锋，刘红秀.超声图像预处理方法的改进研究J.检验医学与临床，2010, 7: 686-688.4郭敏，马远良，朱建.基于小波变换的医学超声图像去 噪及增强方法J.中国医学影像技术，2006, 22: 1435-1437.5Rudin L I , Osher S, Fatemi E. Nonlinear total variation based noise removal algorithmsJ. Physica D , 1992, 60 (1-4): 259-2686吴亚东；孙世新；张红英，等.一种基于图割的全变 差图像去噪算法J.电子学报， 2007, 35, 265-268.7马少贤&江成顺.基于四阶偏微分方程的盲图像恢 复模型J.中国图象图形学报，2010, 15, 26-30.8H. Rabbani , M. Vafadust, P. Abolmaesumi, and S. Gazor.Speckle noise reduction of medica