图像频域增强技巧研究wpy.docVIP

设计题目：图像频域增强方法研究 学院： 班级： 学号： 设计人： 指导老师： 设计时间：2014年12月14日 图像频域增强方法研究 设计目的 掌握图像频域增强的概念及其计算方法； 熟悉掌握傅里叶变换和卷积的计算过程； 熟练掌握频域滤波中常用的Butterworth低通滤波器； 利用MATLAB程序进行图像增强； 掌握常用频域高通滤波器的设计； 加深理解和掌握图像频谱的特点和频域高通滤波的原理。 设计要求 熟练掌握MATLAB的运行环境，并能熟练应用； 掌握并理解MATLAB的各种编程方式及函数定义； 加深对数字图像处理研究。 设计方案 图像频域增强原理 图像增强是按特定的需要突出一幅图像中的某些信息，同时，消弱或去除某些不需要的信息的处理方法。其主要目的是处理后的图像对某些特定的应用比原理的图像更加有效。 2.实现步骤 1 计算需增强图的傅里叶变换； 2 将其与1个（根据需要设计的）转移函数相乘； 3 再将结果傅里叶反变换以得到增强的图。 3.实现方法 常用的频域增强方法有：①低通滤波；②高通滤波；③同态滤波。 4.具体设计内容及其结果 1.频域理想低通滤波器处理图像 I=imread(13.jpg); J=rgb2gray(I); subplot(331); imshow(J); I=double(J); f=fft2(J); g=fftshift(f); subplot(332); imshow(log(abs(g)), [ ]), color(jet(64)); [M, N]=size(f); n1=floor(M/2); n2= floor(N/2); d0=35; for i=1:M for j=1:N d=sqrt((i-n1)^2+(j-n2)^2); if d=d0 h=1; else h=0; end g(i,j)=h*g(i,j); end end g=ifftshift(g); g=uint8(real(ifft2(g))); subplot(333); imshow(g); 2.用巴特沃斯低通滤波器处理图像 I=imread(13.jpg); I1=rgb2gray(I); J=imnoise(I1,salt pepper,0.02); subplot(121);imshow(J); title(含有椒盐噪声的图像) J=double(J); f=fft2(J); g=fftshift(f); [M,N]=size(f); n=3; d0=40 n1=floor(M/2) n2=floor(N/2) for i=1:M for j=1:N d=sqrt((i-n1)^2+(j-n2)^2) h=1/(1+(d/d0)^(2*n)); g(i,j)=h*g(i,j); end end g=ifftshift(g); g=uint8(real(ifft2(g))); subplot(122); imshow(g); 由实验结果知，经频域低通滤波后，图中大部分能量集中在低频分量里。低通滤波器能够让低频信息通过而过滤掉高频分量，以此实现图像平滑去除噪声。 3.频域理想高通滤波器增强图像 I1=imread(jh.jpg); I=rgb2gray(I1);subplot(221);imshow(I); J=imnoise(I, gaussian); subplot(222);imshow(J); title(模糊图像); J=double(J); f=fft2(J); g=fftshift(f); [M, N]=size(f); n1=floor(M/2); n2= floor(N/2); d0=20; for i=1:M for j=1:N d=sqrt((i-n1)^2+(j-n2)^2); if d=d0 h1=1; h2=1+0.5; else h1=0; h2=0.5; end g1(i,j)=h1*g(i,j); g2(i,j)=h2*g(i,j); end end g1=if