图像中值处理MATLAB实现

I=imread('cameraman.tif');
[M N]=size(I);%求出图片大小
I1=I;
%3x3模型第一行一列，最后一行一列不处理保留原值
for i=2:M-1
    for j=2:N-1
        summary=0;
%中值处理
        for p=-1:1 
            for q=-1:1
                summary=summary+I(i+p,j+q)/9;
            end
        end
        I1(i,j)=summary;
    end
end
%画图
subplot(121),imshow(I)
subplot(122),imshow(I1)