用matlab编写程序实现对一幅真彩色图像的双边滤波
时间: 2023-10-09 07:11:14 浏览: 133
图像双边滤波的matlab实现与仿真.zip
好的,以下是用 MATLAB 编写程序实现对一幅真彩色图像的双边滤波的示例代码:
```matlab
% 读取原始图像
img = imread('image.jpg');
% 定义滤波器的参数
window_size = 15;
sigma_color = 75;
sigma_space = 75;
% 对图像进行双边滤波
filtered_img = bilateralFilter(img, window_size, sigma_color, sigma_space);
% 显示原始图像和滤波后的图像
subplot(1, 2, 1);
imshow(img);
title('Original Image');
subplot(1, 2, 2);
imshow(filtered_img);
title('Filtered Image');
% 定义双边滤波函数
function [filtered_img] = bilateralFilter(img, window_size, sigma_color, sigma_space)
% 将图像转换为双精度类型
img = double(img);
% 获取图像的大小
[height, width, channels] = size(img);
% 计算滤波器的半径
radius = floor(window_size / 2);
% 初始化输出图像
filtered_img = zeros(height, width, channels);
% 对每个像素进行滤波
for i = 1:height
for j = 1:width
% 定义滤波器的范围
i_min = max(i - radius, 1);
i_max = min(i + radius, height);
j_min = max(j - radius, 1);
j_max = min(j + radius, width);
% 获取滤波器的窗口
window = img(i_min:i_max, j_min:j_max, :);
% 计算颜色和空间距离权重
color_dist = sum((window - img(i, j, :)) .^ 2, 3);
space_dist = (i_min:i_max)' * ones(1, length(j_min:j_max)) - i;
space_dist = space_dist .^ 2 + (ones(length(i_min:i_max), 1) * (j_min:j_max) - j) .^ 2;
space_dist = exp(-space_dist / (2 * sigma_space ^ 2));
weight = space_dist .* exp(-color_dist / (2 * sigma_color ^ 2));
% 对像素进行加权平均
filtered_img(i, j, :) = sum(sum(weight .* window, 1), 2) / sum(sum(weight, 1), 2);
end
end
% 将输出图像转换为无符号 8 位整型
filtered_img = uint8(filtered_img);
end
```
说明:
- `imread()` 函数用于读取图像文件。
- `bilateralFilter()` 函数用于对图像进行双边滤波处理,其中 `window_size` 表示滤波器的窗口大小,`sigma_color` 表示颜色空间的标准方差,`sigma_space` 表示坐标空间的标准方差。
- `subplot()` 函数用于在同一窗口显示多张图像,`imshow()` 函数用于显示图像,`title()` 函数用于设置图像标题。
- 在代码中,我定义了一个名为 `bilateralFilter()` 的函数来实现双边滤波算法。该函数输入原始图像和滤波器的参数,输出滤波后的图像。函数中使用了两个嵌套的循环来对每个像素进行滤波,并使用了矩阵运算来提高计算效率。
阅读全文