HELP CORNERS_HS                                     David Young
                                                    October 2003

LIB * CORNERS_HS implements the Harris-Stephens (also known as the
"Harris" or "Plessey") corner detector. Harris, C. and Stephens, M.
(1988) A combined corner and edge detector, ACV88 Proceedings of the
Fourth Alvey Vision Conference, pp 147-151.

corner_response_hs(_I, _____sigma, _k, _R) -> _R
corner_response_hs(_I, ______region, _____sigma, _k, _R) -> _R
        Computes the corner response _R for the image and parameters
        given. This is defined by the following equations:

                              2
                _R = Det - _k Tr

        where
                            2
                Det = AB - C
                Tr = A + B

                     2
                A = X  * G(_____sigma)
                     2
                B = Y  * G(_____sigma)
                C = (XY) * G(_____sigma)

                X = _I * Dx
                Y = _I * Dy

        and where "*" means the convolution operation, Dx means the
        [-1 0 1] mask and Dy its vertical equivalent, and G(_____sigma) means
        the Gaussian smoothing mask with width parameter _____sigma.

        The arguments are:

        _I: the input image, an array containing numbers

        ______region: if present, a 4-element list specifying the region of
        the image from which to take data, in *boundslist style. If
        absent, the whole of the array is used.

        _____sigma: a positive number, the smoothing parameter

        _k: a positive number, which controls sensitivity to edges or
        corners. The best value for an application should be found
        empirically - a reasonable starting point appears to be 0.04 for
        images with grey levels in the range 0 to 1, but Harris and
        Stephens fail to discuss this.

        _R: on entry to the procedure this may be:

            * an array, which will be filled with the results and
            returned. Its boundslist must be the same as ______region, if
            given; otherwise its boundslist must be the same as that of
            _I;

            * <false> - a new array will be created and returned;

            * anything else - an array will be created using
            oldsfloatarray (see *newsfloatarray) using the input
            argument as the tag. This will be returned.


corners_hs(_I, _____sigma, _k, ___thr) -> ____list
        Returns a list of corner positions. _I, _____sigma and _k are as for
        corner_response_hs above. ___thr is a number giving a threshold in
        _R for local maximum detection. ____list is a list of possible corner
        positions, taken to be those local maxima of the response array
        that exceed the threshold, in the format used by *array_peaks.


--- __________________________$popvision/help/corners_hs
--- _________Copyright __________University __of ______Sussex _____2003. ___All ______rights _________reserved.