如何用matlab对三维图像进行中值滤波，写一个程序

以下是用MATLAB编写的三维图像中值滤波程序：

% 读取三维图像
img = imread('example.img');

% 设定滤波器大小
filter_size = 3;

% 计算滤波器边缘大小
filter_edge = floor(filter_size/2);

% 获取图像大小
[n, m, p] = size(img);

% 创建新的图像矩阵
new_img = zeros(n, m, p);

% 对每个像素进行滤波处理
for i = 1+filter_edge:n-filter_edge
    for j = 1+filter_edge:m-filter_edge
        for k = 1+filter_edge:p-filter_edge
            % 获取滤波器覆盖范围内的像素值
            neighborhood = img(i-filter_edge:i+filter_edge, j-filter_edge:j+filter_edge, k-filter_edge:k+filter_edge);
            % 对像素值进行排序
            sorted = sort(neighborhood(:));
            % 取中值作为新的像素值
            new_img(i,j,k) = sorted(floor(numel(neighborhood)/2)+1);
        end
    end
end

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

其中，`example.img` 是一个三维图像文件。程序首先读取该文件，并设定滤波器大小。然后，程序遍历图像中的每个像素，在滤波器覆盖范围内获取像素值，并对这些像素值进行排序。最后，程序取排序后的中值作为新的像素值，并将其存储在新的图像矩阵中。程序最后显示原始图像和处理后的图像。