引导图像滤波器 Matlab实现

Matlab demo code for Guided Image Filtering 引导滤波器

Guided Image Filtering 需要两幅图像， 原图和引导图

function imDst = boxfilter（imSrc r）

%   BOXFILTER   O（1） time box filtering using cumulative sum
%
%   - Definition imDst（x y）=sum（sum（imSrc（x-r:x+ry-r:y+r）））;
%   - Running time independent of r; 
%   - Equivalent to the function: colfilt（imSrc [2*r+1 2*r+1] ‘sliding‘ @sum）;
%   - But much faster.

[hei wid] = size（imSrc）;
imDst = zeros（size（imSrc））;

%cumulative sum over Y axis
imCum = cumsum（imSrc 1）;
%difference over Y axis
imDst（1:r+1 :） = imCum（1+r:2*r+1 :）;
imDst（r+2:hei-r :） = imCum（2*r+2:hei :） - imCum（1:hei-2*r-1 :）;
imDst（hei-r+1:hei :） = repmat（imCum（hei :） [r 1]） - imCum（hei-2*r:hei-r-1 :）;

%cumulative sum over X axis
imCum = cumsum（imDst 2）;
%difference over Y axis
imDst（: 1:r+1） = imCum（: 1+r:2*r+1）;
imDst（: r+2:wid-r） = imCum（: 2*r+2:wid） - imCum（: 1:wid-2*r-1）;
imDst（: wid-r+1:wid） = repmat（imCum（: wid） [1 r]） - imCum（: wid-2*r:wid-r-1）;
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        931  2010-07-08 19:30  引导滤波器\boxfilter.m

     文件        452  2010-07-08 20:22  引导滤波器\example_enhancement.m

     文件        285  2010-07-08 19:46  引导滤波器\example_smoothing.m

     文件        957  2010-07-08 21:57  引导滤波器\guidedfilter.m

     文件       2461  2010-07-08 22:11  引导滤波器\guidedfilter_color.m

    ..A.SH.      5632  2010-11-02 14:06  引导滤波器\img_enhancement\Thumbs.db

     文件    2880056  2010-03-05 16:09  引导滤波器\img_enhancement\tulips.bmp

     文件      77256  2010-02-27 16:16  引导滤波器\img_smoothing\cat.bmp

    ..A.SH.      6656  2010-11-02 14:14  引导滤波器\img_smoothing\Thumbs.db

     文件        354  2010-11-02 14:14  引导滤波器\程序说明.txt

     目录          0  2010-11-02 14:06  引导滤波器\img_enhancement

     目录          0  2010-11-02 14:14  引导滤波器\img_smoothing

     目录          0  2010-11-02 14:06  引导滤波器

----------- ---------  ---------- -----  ----

              2975040                    13