空间滤波（模板运算）（线性相关与卷积）

I=imread(‘lena.png’);
[h,w] = size(I);
L=3;%定义尺寸
N=(2L+1)^2;%定义模板的尺寸
mask=(1/N)ones(2L+1,2L+1);%均值算式
II=MirrorImage(double(I),L);%边界扩展
IInew=zeros(h+2L,w+2L);%
for i= (L+1):(h+L)
for j =(L+1):(w+L)
region = II(i-L:i+L,j-L:j+L);
tempI=region.mask;
temp = sum(tempI(😃);%（笑脸改成’:’)
IInew(i,j)=temp;
end
end
Inew=zeros(h,w);
Inew=uint8(IInew(((L+1):(h+L)),((L+1):(L+w))));
figure;
subplot(1,2,1);imshow(I)
title(‘原始图像’);
subplot(1,2,2);imshow(Inew)
title(‘均值滤波后的图像’);

库函数（镜像）
function Ioutput=MirrorImage(Iinput,wmax)
% to create a mirror image, wmax is the numb