matlab双边滤波磨皮程序

Matlab中的双边滤波磨皮程序是一种图像处理算法，用于去除图像中的噪声和细节，同时保留图像的边缘信息。下面是一个简单的双边滤波磨皮程序的示例：

% 读取图像
image = imread('input.jpg');

% 将图像转换为灰度图像
gray_image = rgb2gray(image);

% 双边滤波参数设置
sigma_spatial = 5; % 空间域标准差
sigma_range = 0.1; % 强度域标准差

% 应用双边滤波
filtered_image = bilateralFilter(gray_image, sigma_spatial, sigma_range);

% 显示原始图像和处理后的图像
subplot(1, 2, 1);
imshow(gray_image);
title('原始图像');
subplot(1, 2, 2);
imshow(filtered_image);
title('双边滤波磨皮后的图像');

% 双边滤波函数定义
function filtered_image = bilateralFilter(image, sigma_spatial, sigma_range)
    [height, width] = size(image);
    filtered_image = zeros(height, width);
    
    % 对每个像素进行处理
    for i = 1:height
        for j = 1:width
            pixel = image(i, j);
            
            % 计算空间权重
            spatial_weights = exp(-((i - (1:height)).^2 + (j - (1:width)).^2) / (2 * sigma_spatial^2));
            
            % 计算强度权重
            range_weights = exp(-(pixel - image).^2 / (2 * sigma_range^2));
            
            % 计算归一化的双边滤波结果
            filtered_pixel = sum(sum(range_weights .* spatial_weights .* image)) / sum(sum(range_weights .* spatial_weights));
            
            filtered_image(i, j) = filtered_pixel;
        end
    end
end

这个程序首先读取输入图像，然后将其转换为灰度图像。接下来，通过调用`bilateralFilter`函数来应用双边滤波算法。最后，显示原始图像和处理后的图像。