% -*- texinfo -*-
% @deftypefn {Function File} {[@var{r}, @var{c}] =} immaximas (@var{im}, @var{radius})
% @deftypefnx{Function File} {[@var{r}, @var{c}] =} immaximas (@var{im}, @var{radius}, @var{thresh})
% @deftypefnx{Function File} {[@var{r}, @var{c}, ...] =} immaximas (...)
% @deftypefnx{Function File} {[..., @var{val}] =} immaximas (...)
% Finds local spatial maximas of the given image. A local spatial maxima is
% defined as an image point with a value that is larger than all neighbouring
% values in a square region of width 2*@var{radius}+1. By default @var{radius}
% is 1, such that a 3 by 3 neighbourhood is searched. If the @var{thresh} input
% argument is supplied, only local maximas with a value greater than @var{thresh}
% are retained.
%
% The output vectors @var{r} and @var{c} contain the row-column coordinates
% of the local maximas. The actual values are computed to sub-pixel precision
% by fitting a parabola to the data around the pixel. If @var{im} is
% @math{N}-dimensional, then @math{N} vectors will be returned.
%
% If @var{im} is @math{N}-dimensional, and @math{N}+1 outputs are requested,
% then the last output will contain the image values at the maximas. Currently
% this value is not interpolated.
%
% @seealso{ordfilt2, ordfiltn}
% @end deftypefn