%简单地说:

%A为给定图像,归一化到[,]的矩阵

%W为双边滤波器(核)的边长/

%定义域方差σd记为SIGMA(),值域方差σr记为SIGMA()

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% Pre-process input and select appropriate filter.

function B = bfilter2(A,w,sigma)

% Verify that the input image exists and is valid.

if ~exist('A','var') || isempty(A)

   error('Input image A is undefined or invalid.');

end

if ~isfloat(A) || ~sum([,] == size(A,)) || ...

      min(A(:)) <  || max(A(:)) > 

   error(['Input image A must be a double precision ',...

          'matrix of size NxMx1 or NxMx3 on the closed ',...

          'interval [0,1].']);      

end

% Verify bilateral filter window size.

if ~exist('w','var') || isempty(w) || ...

      numel(w) ~=  || w < 

   w = ;

end

w = ceil(w);

% Verify bilateral filter standard deviations.

if ~exist('sigma','var') || isempty(sigma) || ...

      numel(sigma) ~=  || sigma() <=  || sigma() <= 

   sigma = [ 0.1];

end

% Apply either grayscale or color bilateral filtering.

if size(A,) == 

   B = bfltGray(A,w,sigma(),sigma());

else

   B = bfltColor(A,w,sigma(),sigma());

end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% Implements bilateral filtering for grayscale images.

function B = bfltGray(A,w,sigma_d,sigma_r)

% Pre-compute Gaussian distance weights.

[X,Y] = meshgrid(-w:w,-w:w);

%创建核距离矩阵,e.g.

%  [x,y]=meshgrid(-:,-:)

% 

% x =

% 

%     -          

%     -          

%     -          

% 

% 

% y =

% 

%     -    -    -

%                

%                

%计算定义域核

G = exp(-(X.^+Y.^)/(*sigma_d^));

% Create waitbar.

h = waitbar(,'Applying bilateral filter...');

set(h,'Name','Bilateral Filter Progress');

% Apply bilateral filter.

%计算值域核H 并与定义域核G 乘积得到双边权重函数F

dim = size(A);

B = zeros(dim);

for i = :dim()

   for j = :dim()

         % Extract local region.

         iMin = max(i-w,);

         iMax = min(i+w,dim());

         jMin = max(j-w,);

         jMax = min(j+w,dim());

         %定义当前核所作用的区域为(iMin:iMax,jMin:jMax)

         I = A(iMin:iMax,jMin:jMax);%提取该区域的源图像值赋给I

         % Compute Gaussian intensity weights.

         H = exp(-(I-A(i,j)).^/(*sigma_r^));

         % Calculate bilateral filter response.

         F = H.*G((iMin:iMax)-i+w+,(jMin:jMax)-j+w+);

         B(i,j) = sum(F(:).*I(:))/sum(F(:));

   end

   waitbar(i/dim());

end

% Close waitbar.

close(h);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% Implements bilateral filter for color images.

function B = bfltColor(A,w,sigma_d,sigma_r)

% Convert input sRGB image to CIELab color space.

if exist('applycform','file')

   A = applycform(A,makecform('srgb2lab'));

else

   A = colorspace('Lab<-RGB',A);

end

% Pre-compute Gaussian domain weights.

[X,Y] = meshgrid(-w:w,-w:w);

G = exp(-(X.^+Y.^)/(*sigma_d^));

% Rescale range variance (using maximum luminance).

sigma_r = *sigma_r;

% Create waitbar.

h = waitbar(,'Applying bilateral filter...');

set(h,'Name','Bilateral Filter Progress');

% Apply bilateral filter.

dim = size(A);

B = zeros(dim);

for i = :dim()

   for j = :dim()

         % Extract local region.

         iMin = max(i-w,);

         iMax = min(i+w,dim());

         jMin = max(j-w,);

         jMax = min(j+w,dim());

         I = A(iMin:iMax,jMin:jMax,:);

         % Compute Gaussian range weights.

         dL = I(:,:,)-A(i,j,);

         da = I(:,:,)-A(i,j,);

         db = I(:,:,)-A(i,j,);

         H = exp(-(dL.^+da.^+db.^)/(*sigma_r^));

         % Calculate bilateral filter response.

         F = H.*G((iMin:iMax)-i+w+,(jMin:jMax)-j+w+);

         norm_F = sum(F(:));

         B(i,j,) = sum(sum(F.*I(:,:,)))/norm_F;

         B(i,j,) = sum(sum(F.*I(:,:,)))/norm_F;

         B(i,j,) = sum(sum(F.*I(:,:,)))/norm_F;

   end

   waitbar(i/dim());

end

% Convert filtered image back to sRGB color space.

if exist('applycform','file')

   B = applycform(B,makecform('lab2srgb'));

else  

   B = colorspace('RGB<-Lab',B);

end

% Close waitbar.

close(h);

双边滤波Matlab代码的更多相关文章

  1. 快速双边滤波 附完整C代码

    很早之前写过<双边滤波算法的简易实现bilateralFilter>. 当时学习参考的代码来自cuda的样例. 相关代码可以参阅: https://github.com/johng12/c ...

  2. opencv-12-高斯滤波-双边滤波(附C++代码实现)

    开始之前 这几天由于自己的原因没有写, 一个是因为自己懒了, 一个是感觉这里遇到点问题不想往下写了, 我们先努力结束这个章节吧, 之前介绍了比较常用而且比较好理解的均值和中值滤波, 但是呢,在例程Sm ...

  3. o(1)复杂度之双边滤波算法的原理、流程、实现及效果。

    一.引言     双边滤波在图像处理领域中有着广泛的应用,比如去噪.去马赛克.光流估计等等,最近,比较流行的Non-Local算法也可以看成是双边滤波的一种扩展.自从Tomasi et al等人提出该 ...

  4. Bilateral Filtering(双边滤波) for SSAO(转)

    原文链接:http://blog.csdn.net/bugrunner/article/details/7170471 另外一篇相似的英文资料:http://homepages.inf.ed.ac.u ...

  5. Bilateral Filtering(双边滤波) for SSAO

    原网址:http://blog.csdn.net/bugrunner/article/details/7170471 1. 简介 图像平滑是一个重要的操作,而且有多种成熟的算法.这里主要简单介绍一下B ...

  6. OpenCV计算机视觉学习(4)——图像平滑处理(均值滤波,高斯滤波,中值滤波,双边滤波)

    如果需要处理的原图及代码,请移步小编的GitHub地址 传送门:请点击我 如果点击有误:https://github.com/LeBron-Jian/ComputerVisionPractice &q ...

  7. 利用联合双边滤波或引导滤波进行升采样(Upsampling)技术提高一些耗时算法的速度。

    这十年来,在图像处理领域提出了很多新的图像分析和处理方法,包括是自动的以及一些需要有人工参与的,典型的比如stereo depth computations.image colorization.to ...

  8. 学习 opencv---(8)非线性滤波:中值滤波,双边滤波

    正如我们上一篇文章中讲到的,线性滤波可以实现很多种不同的图像变换.然而非线性滤波,如中值滤波器和双边滤波器,有时可以达到更好的实现效果. 邻域算子的其他一些例子还有对 二值图像进行操作的形态学算子,用 ...

  9. 双边滤波算法的简易实现bilateralFilter

    没怎么看过双边滤波的具体思路,动手写一写,看看能不能突破一下. 最后,感觉算法还是要分开 水平 与 垂直 方向进行分别处理,才能把速度提上去. 没耐性写下去了,发上来,给大伙做个参考好了. 先上几张效 ...

随机推荐

  1. Kubernetes组件及网络基础

    在前面的部分了解了Pod的创建删除 ,查看信息等.那么我们怎么去管理Pod呢?我们可以通过 ReplicationController 去管理维护 Pod. Replication Controlle ...

  2. java中时间与时间戳的相互转换

    package com.test.one; import java.text.ParseException; import java.text.SimpleDateFormat; import jav ...

  3. chrome浏览器截长图的方法

    1.首先打开一个你想要截图的长页面 2.然后按下F12 3.按Ctrl+Shift+P打开console菜单 4.在有个红对勾的位置输入Capture full size screenshot,找到C ...

  4. Springboot+Jedis+Ehcache整合

    项目结构概览: 1. 导包 <parent> <groupId>org.springframework.boot</groupId> <artifactId& ...

  5. oracle数据库中的游标

    oracle中的游标,游标的概念与作用,游标的分类,游标的使用. 一,游标的概念与作用 摘自百度百科:游标(Cursor)是处理数据的一种方法,为了查看或者处理结果集中的数据,游标提供了在结果集中一次 ...

  6. HTML弹性布局

    1.弹性布局的使用 (1)  display:flex:给父容器添加这个属性: (2)  display:flex; 容器添加弹性布局后,显示为块级元素: display:inline-flex; 容 ...

  7. 九、LaTex中的浮动体

    未利用浮动体:

  8. MySQL之表查询

    语法执行顺序 from >>>从那张表 where >>> 全局的筛选条件 group by>>> 必须是用在where 之后一组就是为了接下来我 ...

  9. Ubuntu 14.04 下的MAC OS X 主题安装

    Ubuntu 14.04 下的MAC OS X 主题安装 安装 MAC OS X 主题会帮助你的 Ubuntu 14.04 看起来更像MAC OS X.在这里我们介绍的Macbuntu安装包包含了GT ...

  10. oracle比较两个查询结果的差异

    可以使用minus select * from A minus select * from B; select * from B minus select * from A;