图像处理中的卷积与插值技术实现

16、编写一个名为greduce的函数，该函数接受一个大小为M × N的方形图像img作为输入，将img的行和列与滤波器内核1/20[1 5 8 5 1]ᵀ进行卷积，然后对卷积结果进行下采样，以生成大小为M/2 × N/2的输出图像。在你选择的图像上演示该函数。

以下是实现 greduce 函数的 MATLAB 代码示例：

function output_img = greduce(img)
    % 定义滤波器内核
    kernel = (1/20) * [1 5 8 5 1]';

    % 对图像的行进行卷积
    img_conv_rows = conv2(img, kernel', 'same');

    % 对图像的列进行卷积
    img_conv = conv2(img_conv_rows, kernel, 'same');

    % 下采样，生成大小为 M/2 × N/2 的图像
    output_img = img_conv(1:2:end, 1:2:end);
end

你可以使用以下代码在一个示例图像上演示该函数：

% 读取图像，这里以 MATLAB 自带的示例图像为例
img = imread('cameraman.tif');

% 调用 greduce 函数
output = greduce(img);

% 显示原始图像和处理后的图像
subplot(1,2,1);
imshow(img);
title('原始图像');

subplot(1,2,2);
imshow(output);
title('处理后的图像');

上述代码首先定义了 greduce 函数，该函数实现了对输入图像的卷积和下采样操作。然后，使用 MATLAB 自带的示例图像 cameraman.tif 进行演示，将原始图像和处理后的图像显示在同一窗口中。你可以根据需要替换为自己选择的图像。