---

nonmaxsuprx.c

Go to the documentation of this file.

#include "sadie.h"
#include "proto.h"

/*----------------------------------------------------------------------------*/
/*-General Information--------------------------------------------------------*/
/*                                                                            */
/* This function performs non-maxima suppression in the x-direction           */
/*                                                                            */
/*----------------------------------------------------------------------------*/
/*-Interface Information------------------------------------------------------*/
void NONMAXSUPRX (
IMAGE  *in,      /*  I   Pointer to the input image.                          */
IMAGE  **out     /*  O   Address of a pointer to the output image             */
/*----------------------------------------------------------------------------*/
) { register short i, j, k, n;
    char msg[SLEN];
        

    /* check input */
    if (!CHECKIMG(in)) {
        MESSAGE('E'," Can't identify image.");
        goto the_end;
    } else if (in->nbnd > 1) {
        MESSAGE('E'," Image must be single-band.");
        goto the_end;
    }


    /* create image of appropriate size */
    if (!CHECKIMG(*out)) GETMEM(in->nbnd,in->nlin,in->npix,out);
    if (!*out) goto the_end;



    /* perform non-maxima supression in the x-direction */
    for (j=0; j<in->nlin; j++) {
        for (k=0; k<in->npix; k++) {

            if ((k < 1)||(k > in->npix-2)) {
                /* ignore any data at the left and right edges of the window */
                (*out)->data[0][j][k] = 0.0;
            } else if (((k>0)&&(in->data[0][j][k-1] > in->data[0][j][k]))
                || ((k<in->npix-1)&&(in->data[0][j][k+1] > in->data[0][j][k]))) {
                    /* Supress this non-maxima edge pixel */
                    (*out)->data[0][j][k] = 0.0;
            } else {
                (*out)->data[0][j][k] = in->data[0][j][k];
            }
        }
    }
    
        
    the_end:
    /* if (TIMES) TIMING(T_EXIT); */
}

---