matlab 中值滤波

1、中值滤波的原理：对于一串连续输入的信号（量化后是一组数据）。如下图所示，是输入的原信号。中值滤波的原理为，重新计算每一个x的输出值(y)，新的输出值。

相当于y=new(x)，new的 *** 作是，从在以x为中心，长度为2k的原信号中（区间为[x-k+1,x+k]），提取出这段区间内中间的那个值，作为y=new(x)的结果。

2、举例来说，输入：Y[1-10]:1,2,3,4,5,6,7,8,9,10.取区间2k=4,所以k=2执行中值滤波K=中值滤波(Y)、由x-k+1>=1，所以当k=2时，x>=2、滤波时：

K[1]=Y[1]

K[2]=(Y[1]、Y[2]、Y[3]、Y[4])的中间值，即为2或3。

3、matlab的中值滤波实现方式：调用函数：A=medfilt1(B,n)、B为输入信号，A为滤波后的信号，即结果。

4、对于输入信号（最开头的图），以下分别为设置区间n=8和n=16得到的滤波图像。中值滤波可以过滤尖峰脉冲。目的在于我们对于滤波后的数据更感兴趣。滤波后的数据保留的原图像的变化趋势，同时去除了尖峰脉冲对分析造成的影响。

// 中值滤波对椒盐噪声

RGB=imread('peppers','PNG')

I=rgb2gray(RGB)

J=imnoise(I,'salt &pepper',0.02)

subplot(121),imshow(J)

L=medfilt2(J,[3 3])

subplot(122),imshow(L)

I=imread('lena.bmp')

inf=imfinfo('lena.bmp')

figure,imshow(I)

X=grayslice(I,64)

imshow(X,pink(64))

load trees

figure,image(10,10,X)

imwrite(X,map,'trees.bmp')

imfinfo('trees.bmp')

figure,imshow(X,map)

BW=im2bw(X,map,0.6)

figure,imshow(BW)

I=imread('lena.bmp')

inf=imfinfo('lena.bmp')

figure,imshow(I)

X=grayslice(I,64)

figure,imshow(X,pink(64))

A=imread('lena.bmp')

imshow(A)

B=fftshift(fft2(A))

figure

imshow(log(abs(B)),[8,10])

clc

clear all

I=imread('lena.tif')

% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% %用中值滤波，多维滤波，使用中心为-4,-8的拉普

% %拉斯滤波器，高斯低通滤波，拉普拉斯滤波器进行滤波处理

% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

figure%figure1

subplot(2,2,1)

imshow(I)

title('原始图像')

J=imnoise(I,'salt &pepper',0.04)%加椒盐噪声

title('加椒盐噪声')

subplot(2,2,2)

imshow(J)

K=medfilt2(J,[4,4])%进行中值滤波

subplot(2,2,3)

imshow(K)

title('进行中值滤波')

h=ones(3,3)/9%多维滤波

w=h

g=imfilter(I,w,'conv','replicate')

subplot(2,2,4)

imshow(g)

title('多维滤波')

%使用中心为-4,-8的拉普拉斯滤波器,

w4=[1 1 1

1 -4 1

1 1 1]

w8=[1 1 1

1 -8 1

1 1 1]

f=im2double(I)

g4=f-imfilter(f,w4,'replicate')

g8=f-imfilter(f,w8,'replicate')

imshow(f)

figure%figure2

subplot(1,2,1)

imshow(g4)

title('中心为-4的拉普拉斯滤波')

subplot(1,2,2)

imshow(g8)

title('中心为-8的拉普拉斯滤波')

h3=fspecial('gaussian',[3,3],0.5)%高斯低通滤波

figure%figure3

B4=filter2(h3,I)

subplot(1,2,1)

imshow(B4,[ ])

title('高斯低通滤波')

h4=fspecial('laplacian',0)%使用拉普拉斯滤波器

B5=filter2(h4,I)

subplot(1,2,2)

imshow(B5,[ ])

title('拉普拉斯滤波器')

% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% %从空域的角度进行亮度变换

% %把灰度等级是10-100的变化到10-255

% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

figure%figure4

subplot(2,2,1)

imshow(I)

title('原始图像')

J2=imadjust(I,[],[],0.5)% 增强低灰度级

subplot(2,2,2)

imshow(J2)

title('增强低灰度级')

J3=imadjust(I,[ ],[ ],2)%增强高灰度级

subplot(2,2,3)

imshow(J3)

title('增强高灰度级')

a1=100/255%把灰度等级是10-100的变化到10-255

a2=255/255

a3=10/255

J2=imadjust(I,[a3,a1],[a3,a2],[])

subplot(2,2,4)

imshow(J2)

title('把灰度等级是10-100的变化到10-255')

% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% %从频域的角度进行亮度变换

% %fft2

% %由于能量主要集中在低频部分

% %所以对低频进行处理可以得到理想的效果

% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

I=imread('lena.tif')

up=0.5%设置处理频率上限

down=0.09%%设置处理频率下限

figure%figure5

subplot(421)

imshow(I)

title('原始图像')

f=double(I)

subplot(4,2,3)

imshow(log(abs(f)),[])

title('unit8转化为double')

g=fft2(f)

g=fftshift(g)

subplot(4,2,5)

imshow(log(abs(g)),[])

title('FFT2变化后的图像')

[M,N]=size(g)% 转换数据矩阵

y1=max(max(abs(g)))%求出最大频率

y2=min(min(abs(g)))%%求出最小频率

y3=(y1-y2)*up+y2%设置滤波上限

y4=(y1-y2)*down+y2%%设置滤波下限

for i=1:M

for j=1:N

if (abs(g(i,j))<y4)

g(i,j)=g(i,j)^1.1%对低频部分进行灰度增强

end

end

end

result=ifftshift(g)

J2=ifft2(result)

J3=uint8(abs(J2))

subplot(427)

imshow(J3,[ ])

title('频域处理后的图像')

subplot(422)

imhist(I,64)

subplot(424)

imhist(f,64)

subplot(426)

imhist(g,64)

subplot(428)

imhist(J3,64)

---