Scuola di Dottorato di Ricerca in Parodontologia Sperimentale

% DOMAIN is equivalent to the structuring element used for
% binary image operations. It is a matrix containing only 1's
% and 0's; the 1's define the neighborhood for the filtering
% operation.
%
% For example, B=ORDFILT2(A,5,ONES(3,3)) implements a 3-by-3
% median filter; B=ORDFILT2(A,1,ONES(3,3)) implements a 3-by-3
% minimum filter; and B=ORDFILT2(A,9,ONES(3,3)) implements a
% 3-by-3 maximum filter. B=ORDFILT2(A,4,[0 1 0; 1 0 1; 0 1 0])
% replaces each element in A by the maximum of its north, east,
% south, and west neighbors.
%
% See also MEDFILT2.
% Copyright 1993-2004 The MathWorks, Inc.
% $Revision: 5.17.4.6 $ $Date: 2004/08/10 01:41:00 $
[A,order,domain,s,padopt,msg] = ParseInputs(varargin{:});
domainSize = size(domain);
center = floor((domainSize + 1) / 2);
[r,c] = find(domain);
r = r - center(1);
c = c - center(2);
padSize = max(max(abs(r)), max(abs(c)));
originalSize = size(A);
if (strcmp(padopt, 'zeros'))
A = padarray(A, padSize * [1 1], 0, 'both');
elseif (strcmp(padopt, 'ones'))
% padopt of 'ones' is for support of medfilt2; it is
% undocumented
A = padarray(A, padSize * [1 1], 1, 'both');
else
A = padarray(A, padSize * [1 1], 'symmetric', 'both');
end
Ma = size(A,1);
offsets = c*Ma + r;
% make sure that offsets are valid
if ~isreal(offsets) || any(floor(offsets) ~= offsets) || any(~isfinite(offsets))
98