小波分形插值医学图像处理.pdfVIP

维普资讯 国耗技 2001年第6期 裂ESEARC秘＆DEVEI研OP究M与E开N发I 文章编号 ：1001—893X(2001)06—0067一O4 小 波 分 形插 值 医学 图像 处 理 陈德元，涂国防 (中国科学技术大学研究生院．北京 10co39) 摘【要】为了得到高分辩率图像，本文拳j用小液变换后各方 向上的高频图像具有相似性的特点，提出 了小波分彤插值处理医学图像的方法。仿真实验表明，这种处理方法比常用的方法(如双线性插值 等)具有更好的效果 关键词 ：图像处理 ；分形布 朗运动 ；小波变换 ；分形插值 中图分类号 ：TN919．81 文葡标识码 ：A 的多分辨分析(MRA)V~，使得f(x，Y)可分解为如下 一 、 引 言 形式 在医学图像成像时，由于技术等原因，图像分辨 5f，(n，m)： ，Y)0 (一 ，一y))(2～ ，2～m) 率较低 。从低分辨率数据得到高分辨率数据，最直 ，(n，m)： ，)0 (一 ，一y))(2～n，2m) 接的方法就是插值。但医学图像 的高度复杂及丰富 称为二维 图像 f(x，Y)的离散小渡变换 (DWT)。其 的纹理促使人们寻求新的插值方法代替传统插值方 中m，i,i∈Z2，i=1，2，3，符号0表线性卷积 ，一1≤J 法(如双线性插值 、样条插值等)。J、Grevera等 提 ≤一J， (x，Y)为相应的 MRA的尺 度函数 中(x， 出基于形状的多维灰度级图像插值，先将 rt维 图像 v)的膨胀．即 “提升”(1ifting)为 11+1维 图像，对 D-+1维图像进行 ( ， )：2 (2 ，2 ) 插值 ，然后再将 rt+1维插值图像坍塌 (collapsing)为 (x，v)为相应于尺度函数 中(x，v)的小波 函 rt维图像。这种处理方法所放大的图像会丢失图像 数， (x，Y)为小波 函数的膨胀 ，即 的纹理。有 的文献 提出 了用分形插值 处理医学 图 (，y)：2J (2 ，2Y) 像，使得图像的纹理得以保持 ，但没有考虑对图像噪 为计算方便 ，一般取可分离的小渡 函数 ，即 声的抑制。 雪(，Y)= c) () 本文考虑到小波变换 (WT)具有过滤噪声的作 用 ，且图像的小渡变换具有时频特性 ，对小渡变换后 (，Y)= () (Y) 的图像进行插值处理 ，再作逆小渡变换 (IWT)得处 (，Y)= () (Y) 理后的图像。这样 ，在插值处理时可避免同时放大 (，Y)= () (Y) 噪声。最后，对脑 CT的模拟实验表 明，本文提出的 在上述分解 中，Sjf(D-，m)为 f图像 (x，Y)的低 方法具有 良好效果