% INTEGAVERAGE - performs averaging filtering  using an integral image
%
% Usage:  avim = integaverage(intim,rad)
%
% Arguments:  initim  - Integral image
%             rad     - 'Radius' of square region over which averaging is
%                       performed  (rad = 1 implies a 3x3 average) 
% Returns:    avim    - Averaged image
%
% See also:  INTEGRALIMAGE, INTEGRALFILTER, INTFILTTRANSPOSE

% Reference: Paul Viola and Michael Jones, "Robust Real-Time Face Detection",
% IJCV 57(2). pp 137-154. 2004.

% Copyright (c) 2007 Peter Kovesi
% School of Computer Science & Software Engineering
% The University of Western Australia
% pk at csse uwa edu au
% http://www.csse.uwa.edu.au/
% 
% Permission is hereby granted, free of charge, to any person obtaining a copy
% of this software and associated documentation files (the "Software"), to deal
% in the Software without restriction, subject to the following conditions:
% 
% The above copyright notice and this permission notice shall be included in 
% all copies or substantial portions of the Software.
%
% The Software is provided "as is", without warranty of any kind.

% October 2007

function avim = integaverage(intim, rad)
    
    [rows, cols] = size(intim);
    avim = zeros(rows, cols);

    % Fiddle with indices to ensure we calculate the average over a square
    % region that has an odd No of pixels on each side (ie has a centre pixel
    % located over each pixel of interest)
    
    down = rad;   % offset to 'lower' indices
    up = down+1;  % offset to 'upper' indices
    
    for r = 1+up:rows-down
        for c = 1+up:cols-down
        
            avim(r,c) = intim(r+down,c+down) - intim(r-up,c+down) ...
                       - intim(r+down,c-up) + intim(r-up,c-up);
            
        end
    end
    
    avim = avim/(down+up)^2;