标签：end 滤波 depth 图像处理 radius Matlab sigma result

网上的介绍太多了，但绝大多数Matlab代码应用了较多的循环，速度较慢。本代码充分应用矩阵运算，循环次数少，速度极快。

归一化到0到1之间以后，加入零均值，标准差为0.01的高斯白噪声，得到下图

设置 σ s = 5 ， σ r = 0.2 \sigma_s=5，\sigma_r=0.2 σs​=5，σr​=0.2，邻域边长 k = 9 ， k=9， k=9，得到滤波效果

笔记本CPU为i7-4710HQ，0.1秒，灰度图滤波速度很快，RGB彩色图即进行三次灰度图滤波，耗时乘3。

代码：

clear all
close all
f = imread('滤波原始图.jpg');
[row, col, depth] = size(f);
f = double(f);
f = imnoise(f / 255, 'gaussian', 0, 0.01);
figure, imshow(f);
%% 双边滤波
sigma_s = 5;
sigma_r = 0.2;
tic
result = zeros(row, col, depth);
A_sum = zeros(row, col, depth);
k = 9;
radius = floor(k / 2);
f_padded = padarray(f, [radius radius], 'symmetric');
for l = 1 : depth
    for i = -radius : radius
        for j = -radius : radius
            moving_f = f_padded(i + radius + 1 : end + i - radius, j + radius + 1 : end + j - radius, l);
            A = exp(-(i ^ 2 + j ^ 2) / (2 * sigma_s ^ 2) - (moving_f - f) .^ 2 / (2 * sigma_r ^ 2));
            result(:, :, l) = result(:, :, l) + moving_f .* A;
            A_sum = A_sum + A;
        end
    end
    result(:, :, l) = result(:, :, l) ./ A_sum;
end
result = result * 255;
result = uint8(result);
toc;
figure, imshow(result);

标签：end,滤波,depth,图像处理,radius,Matlab,sigma,result
来源： https://blog.csdn.net/u011861755/article/details/112146274

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。