基于小波变换的医学超声图像去噪方法研究程序Word格式文档下载.docx
- 文档编号:19915153
- 上传时间:2023-01-12
- 格式:DOCX
- 页数:20
- 大小:724.21KB
基于小波变换的医学超声图像去噪方法研究程序Word格式文档下载.docx
《基于小波变换的医学超声图像去噪方法研究程序Word格式文档下载.docx》由会员分享,可在线阅读,更多相关《基于小波变换的医学超声图像去噪方法研究程序Word格式文档下载.docx(20页珍藏版)》请在冰豆网上搜索。
A=log(1+J);
subplot(2,2,4);
imshow(A,[]);
对数变换后的图像'
imshow(y);
J=imnoise(y,'
subplot(2,2,2);
imshow(J);
[m,n]=size(J);
J=double(J);
k1=filter2(fspecial('
average'
3),J)/255;
%3*3均值滤波
k2=filter2(fspecial('
7),J)/255;
%7*7均值滤波
subplot(2,2,3);
imshow(k1,[]);
3*3均值滤波'
subplot(2,2,4);
imshow(k2,[]);
7*7均值滤波'
[PSNR1,MSE1]=psnr(y,k1);
[PSNR2,MSE2]=psnr(y,k2);
display('
MSE1:
'
MSE1
PSNR1:
PSNR1
MSE2:
MSE2
PSNR2:
PSNR2
MSE1=
1.0948e+004
PSNR1=
7.7376
MSE2=
PSNR2=
7.7373
x1=medfilt2(J,[33]);
%3*3中值滤波
x2=medfilt2(J,[77]);
%7*7中值滤波
imshow(x1,[]);
3*3中值滤波'
imshow(x2,[]);
7*7中值滤波'
[PSNR1,MSE1]=psnr(y,x1);
[PSNR2,MSE2]=psnr(y,x2);
156.7310
26.1793
99.1610
28.1674
closeall;
clc
c1=wiener2(J,[33]);
%3*3维纳滤波
c2=wiener2(J,[77]);
%7*7维纳滤波
figure;
imshow(y,[]);
imshow(J,[]);
imshow(c1,[]);
3*3维纳滤波'
imshow(c2,[]);
7*7维纳滤波'
[PSNR1,MSE1]=psnr(y,c1);
[PSNR2,MSE2]=psnr(y,c2);
118.0365
27.4106
88.1850
28.6769
clc;
[c,l]=wavedec2(J,2,'
sym4'
a1=uint8(wrcoef2('
a'
c,l,'
coif2'
2));
[gc,gl]=wavedec2(a1,2,'
n=[1,2];
%设置尺度向量
p=[10.28,24.08];
%设置阈值向量
nc=wthcoef2('
h'
gc,gl,n,p,'
s'
v'
d'
mc=wthcoef2('
d1=waverec2(nc,gl,'
imshow(x1,[]);
第一次消噪图像'
d2=waverec2(mc,gl,'
imshow(x2,[]);
第二次消噪图像'
[m,n]=size(J);
A=log(1+J);
THR=thselect(A,'
heursure'
[c,l]=wavedec2(A,2,'
a1=wrcoef2('
2);
[gc,gl]=wavedec2(a1,2,'
n=[1,2];
p=[THR,THR];
nc=wthcoef2('
mc=wthcoef2('
nc,gl,n,p,'
e1=waverec2(nc,gl,'
E1=exp(e1)-1;
e2=waverec2(mc,gl,'
E2=exp(e2)-1;
imshow(E1,[]);
改进后第一次消噪图像'
imshow(E2,[]);
改进后第二次消噪图像'
[PSNR1,MSE1]=psnr(y,E1);
[PSNR2,MSE2]=psnr(y,E2);
84.9459
28.8394
b1=medfilt2(J,[33]);
b2=medfilt2(J,[77]);
imshow(b1,[]);
imshow(b2,[]);
THR=thselect(J,'
[c,l]=wavedec2(J,2,'
a1=uint8(wrcoef2('
d1=waverec2(nc,gl,'
imshow(d1,[]);
imshow(d2,[]);
[PSNR3,MSE3]=psnr(y,b1);
[PSNR4,MSE4]=psnr(y,b2);
MSE3:
MSE3
PSNR3:
PSNR3
MSE4:
MSE4
PSNR4:
PSNR4
[PSNR5,MSE5]=psnr(y,c1);
[PSNR6,MSE6]=psnr(y,c2);
MSE5:
MSE5
PSNR5:
PSNR5
MSE6:
MSE6
PSNR6:
PSNR6
[PSNR7,MSE7]=psnr(y,d1);
[PSNR8,MSE8]=psnr(y,d2);
MSE7:
MSE7
PSNR7:
PSNR7
MSE8:
MSE8
PSNR8:
PSNR8
[PSNR9,MSE9]=psnr(y,E1);
[PSNR0,MSE0]=psnr(y,E2);
MSE9:
MSE9
PSNR9:
PSNR9
MSE0:
MSE0
PSNR0:
PSNR0
MSE3=
156.5938
PSNR3=
26.1831
MSE4=
99.0738
PSNR4=
28.1712
MSE5=
PSNR5=
MSE6=
PSNR6=
MSE7=
76.0333
PSNR7=
29.3208
MSE8=
PSNR8=
MSE9=
PSNR9=
MSE0=
PSNR0=
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 变换 医学 超声 图像 方法 研究 程序
![提示](https://static.bdocx.com/images/bang_tan.gif)