基于小波變換的醫(yī)學(xué)超聲圖像去噪方法研究程序

《基于小波變換的醫(yī)學(xué)超聲圖像去噪方法研究程序》由會(huì)員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在教育資源-天天文庫(kù)。

1、I=imread('xian.bmp');subplot(2,2,1);imshow(I);title('原始圖像');>>if~isempty(size(I,3))y=rgb2gray(I);end>>subplot(2,2,2);>>imshow(y);>>title('原始灰度圖像');>>J=imnoise(y,'speckle');>>subplot(2,2,3);>>imshow(J);>>title('加噪后的圖像');>>J=double(J);>>A=log(1+J);>>subplot(2,2,4);>>imshow(A,[]);>>title('對(duì)數(shù)變換后

2、的圖像');I=imread('xian.bmp');if~isempty(size(I,3))y=rgb2gray(I);endsubplot(2,2,1);imshow(y);title('原始灰度圖像');J=imnoise(y,'speckle');subplot(2,2,2);imshow(J);title('加噪后的圖像');[m,n]=size(J);J=double(J);k1=filter2(fspecial('average',3),J)/255;%3*3均值濾波k2=filter2(fspecial('average',7),J)/255;%7*7均值濾波

3、subplot(2,2,3);imshow(k1,[]);title('3*3均值濾波');subplot(2,2,4);imshow(k2,[]);title('7*7均值濾波');[PSNR1,MSE1]=psnr(y,k1);[PSNR2,MSE2]=psnr(y,k2);display('MSE1:');MSE1display('PSNR1:');PSNR1display('MSE2:');MSE2display('PSNR2:');PSNR2MSE1:MSE1=1.0948e+004PSNR1:PSNR1=7.7376MSE2:MSE2=1.0948e+004PSNR

4、2:PSNR2=7.7373I=imread('xian.bmp');if~isempty(size(I,3))y=rgb2gray(I);endsubplot(2,2,1);imshow(y);title('原始灰度圖像');J=imnoise(y,'speckle');subplot(2,2,2);imshow(J);title('加噪后的圖像');[m,n]=size(J);J=double(J);x1=medfilt2(J,[33]);%3*3中值濾波x2=medfilt2(J,[77]);%7*7中值濾波subplot(2,2,3);imshow(x1,[]);tit

5、le('3*3中值濾波');subplot(2,2,4);imshow(x2,[]);title('7*7中值濾波');[PSNR1,MSE1]=psnr(y,x1);[PSNR2,MSE2]=psnr(y,x2);display('MSE1:');MSE1display('PSNR1:');PSNR1display('MSE2:');MSE2display('PSNR2:');PSNR2MSE1:MSE1=156.7310PSNR1:PSNR1=26.1793MSE2:MSE2=99.1610PSNR2:PSNR2=28.1674closeall;clcI=imread('x

6、ian.bmp');if~isempty(size(I,3))y=rgb2gray(I);endJ=imnoise(y,'speckle');[m,n]=size(J);J=double(J);c1=wiener2(J,[33]);%3*3維納濾波c2=wiener2(J,[77]);%7*7維納濾波figure;subplot(2,2,1);imshow(y,[]);title('原始灰度圖像');subplot(2,2,2);imshow(J,[]);title('加噪后的圖像');subplot(2,2,3);imshow(c1,[]);title('3*3維納濾波');

7、subplot(2,2,4);imshow(c2,[]);title('7*7維納濾波');[PSNR1,MSE1]=psnr(y,c1);[PSNR2,MSE2]=psnr(y,c2);display('MSE1:');MSE1display('PSNR1:');PSNR1display('MSE2:');MSE2display('PSNR2:');PSNR2MSE1:MSE1=118.0365PSNR1:PSNR1=27.4106MSE2:MSE2=88.1850PSNR2:PSNR2=