---

isolate.c

Go to the documentation of this file.

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

/*----------------------------------------------------------------------------*/
/*-General Information--------------------------------------------------------*/
/*                                                                            */
/*   This function isolates a single pixel value                 */
/*                                                                            */
/*----------------------------------------------------------------------------*/
/*-Interface Information------------------------------------------------------*/
void ISOLATE (
IMAGE  *in,      /*  I   Pointer to the input image.                           */
double  value,
IMAGE  **out     /*  O   Address of a pointer to the output image             */
/*----------------------------------------------------------------------------*/
) { register int i, j, k;
    char msg[SLEN];


    /* check input */
    if (!CHECKIMG(in)) {
        MESSAGE('E'," Can't identify image.");
        goto the_end;
    }

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


    for (i=0; i<in->nbnd; i++) {
        for (j=0; j<in->nlin; j++) {
            for (k=0; k<in->npix; k++) {
                if (in->data[i][j][k] > ((PIXEL)value + (PIXEL)0.01)) {
                    (*out)->data[i][j][k] = (PIXEL)0.0;
                } else if (in->data[i][j][k] < ((PIXEL)value - (PIXEL)0.01)) {
                    (*out)->data[i][j][k] = (PIXEL)0.0;
                } else {
                    (*out)->data[i][j][k] = (PIXEL)value;
                }
            }
        }
    }
        
    the_end:
    /* if (TIMES) TIMING(T_EXIT); */
}

---