在上一个代码基础上改写

基于上一个代码，将均值滤波器修改为4x4的方阵，可以按照以下方式修改代码：

clc; %清空命令窗口
clear all; %清空工作区变量
A = imread('pout.tif'); %读取图像
B = imnoise(A, 'salt & pepper', 0.2); %加入椒盐噪声
C = double(B); %将图像转换为double型
[m, n] = size(C);
D = C; %初始化去噪后的图像
%4x4均值滤波器
w = 1/16 * ones(4);
%对图像进行卷积操作
for i=3:m-2
    for j=3:n-2
        D(i,j) = sum(sum(w.*C(i-2:i+1,j-2:j+1)));
    end
end
D = uint8(D); %将图像转换为uint8型
%显示图像
subplot(1, 3, 1); imshow(A); title('原始图像');
subplot(1, 3, 2); imshow(B); title('加入椒盐噪声图像');
subplot(1, 3, 3); imshow(D); title('4x4均值滤波去噪图像');

其中，将均值滤波器修改为4x4的方阵`w = 1/16 * ones(4)`，遍历图像像素点时，需要将范围从原来的2x2调整为4x4，即`D(i,j) = sum(sum(w.*C(i-2:i+1,j-2:j+1)))`。其他部分与之前的代码类似。需要注意的是，在修改均值滤波器的大小时，需要根据具体情况调整，不能过度增加模板的大小，否则会导致图像细节信息的丢失。