空间滤波

1.1 创建二维滤波器：fspecial函数
语法格式：h=fspecial(type)
h=fspecial(type,parameters)

I=imread('dragon.jpg');
subplot(2,2,1);
imshow(I);
title('sharpened image');
H=fspecial('motion',20,50);
MotionBlur=imfilter(I,H,'replicate');
subplot(2,2,2);
imshow(MotionBlur);
title('motion blurred image');
H=fspecial('disk',10);
blurred=imfilter(I,H,'replicate');
subplot(2,2,3);
imshow(blurred);
title('blurred image');
H=fspecial('unsharp');
sharpened=imfilter(I,H,'replicate');
subplot(2,2,4);
imshow(sharpened);
title('sharpened image');

1.2 进行二维线性滤波操作: filter2函数
语法格式：Y=filter2(h,X)
Y=filter2(h,X,shape) 使用二维FIR滤波器h对矩阵X进行滤波操作字符串参量shape指定返回值Y的形式，取值为full（Y的维数大于 X）,same（Y的维数等于X）和valid（Y的维数小于X），默认值为same

>> x=magic(3);%创建一个3*3的魔幻矩阵
>> h=fspecial('motion',20,45);%生成长度20、角度45度的近似线性移动滤波
>> y=filter2(h,x) %二维线性滤波

y =

    0.4094    0.4793    0.9311
    0.4793    0.9311    0.9365
    0.9311    0.9365    0.2985

1.3 对图像进行滤波：imfilter函数
语法格式;g=imfilter(f,w,filtering_mode,boundary_options,size_options)

f=imread('dragon.jpg');
g=double(f);
w=ones(31) %31*31的简单滤波器
gd=imfilter(g,w)%设置默认的边界选项，默认零填充
imshow(gd,[]);
gr=imfilter(g,w,'replicate');%复制外边界来填充
figure,imshow(gr,[]);
gs=imfilter(g,w,'symmetric');%镜像反射其边界来填充
figure,imshow(gs,[]);
gc=imfilter(g,w,'circular');%二值周期函数的周期填充
figure,imshow(gc,[]);
g8=im2uint8(g);
g8r=imfilter(g8,w,'replicate');
figure,imshow(g8r,[])

1.4 二维中值滤波函数：medfilt2函数
语法格式：g=medfilt2(f,[m n],padopt)
g=medfilt2(f)