引导图像滤波器 Matlab实现 绝对能用

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

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2010-11-02 14:06  引导滤波器
     文件         931  2010-07-08 19:30  引导滤波器oxfilter.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
     目录           0  2010-11-02 14:06  引导滤波器img_enhancement
     文件        5632  2010-11-02 14:06  引导滤波器img_enhancementThumbs.db
     文件     2880056  2010-03-05 16:09  引导滤波器img_enhancement ulips.bmp
     目录           0  2010-11-02 14:14  引导滤波器img_smoothing
     文件       77256  2010-02-27 16:16  引导滤波器img_smoothingcat.bmp
     文件        6656  2010-11-02 14:14  引导滤波器img_smoothingThumbs.db
     文件         354  2010-11-02 14:14  引导滤波器程序说明.txt