数字图像实验报告 医学图像频域滤波.pdfVIP

数字图像实验报告医学图像频域滤波--第1页

实验8医学图像频域滤波与图像复原

实验目的：

1.熟悉医学图像离散傅里叶变化的原理和方法；

2.掌握医学图像频域滤波的原理；

3.掌握使用Matlab中的函数实现医学图像进行频域滤波的方法;

4.掌握使用Matlab中的图像退化与复原的方法;

实验内容：

一、医学图像频域滤波方法与实现

使用imnoise给图像BMRI1_24bit.bmp添加概率为0.2的椒盐噪声，对原图

像和加噪声后的图像进行傅氏变换并显示变换后的移中频谱图，然后分别使用

Butterworth低、高通滤波器对噪声图像进行低通和高通滤波，显示D为5，10，

0

20，40时的滤波效果图，并说明存两种滤波效果中所存在的差异及原因。

f=imread(BMRI1_24bit.bmp);

g=imnoise(f,saltpepper,0.2);

Tf=graythresh(f);

Tg=graythresh(g);

f1=im2bw(f,Tf);

g1=im2bw(g,Tg);

imshow(f1),title原图二值化;figure,imshow(g1),title噪声图像二值化;

数字图像实验报告医学图像频域滤波--第1页

数字图像实验报告医学图像频域滤波--第2页

fft2f=fft2(f1);

fft2g=fft2(g1);

imshow(abs(fft2f)),title原频谱图;figure,imshow(fft2g),title噪音频谱图;

shiftf=fftshift(fft2f);

shiftg=fftshift(fft2g);

imshow(abs(shiftf),[]),title原居中;

figure,imshow(abs(shiftg),[]),title噪音居中;

数字图像实验报告医学图像频域滤波--第2页

数字图像实验报告医学图像频域滤波--第3页

logf=log(1+abs(shiftf));

logg=log(1+abs(shiftg));

imshow(logf,[]),titlelog变换;figure,imshow(logg,[]),titlezaoyinlog;

f=imread(BMRI1_24bit.bmp);

g=imnoise(f,saltpepper,0.2);

BLPF(g);

D0=5和D0=10

D0=20和D0=40

数字图像实验报告医学图像频域滤波--第3页

数字图像实验报告医学图像频域滤波--第4页

GLPF(g);

D0=5和D0=10

D0=20和D0=40

数字图像实验报告医学图像频域滤波--第4页

数字图像实验报告医学图像频域滤波--第