---

bingrid.c

Go to the documentation of this file.

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

/*----------------------------------------------------------------------------*/
/*-General Information--------------------------------------------------------*/
/*                                                                            */
/*   This function produces a grid on the input binary image.                 */
/*                                                                            */
/*----------------------------------------------------------------------------*/
/*-Interface Information------------------------------------------------------*/
void BINGRID (
IMAGE  *in,      /*  I   Pointer to the input image.                           */
IMAGE  **out     /*  O   Address of a pointer to the output image             */
/*----------------------------------------------------------------------------*/
) { register int i, j, k, m, n;
    char msg[SLEN];
    int pixsize = 5;

    /* if (TIMES) TIMING(T_BINGRID); */
    if (NAMES) {
        MESSAGE('I',"");
        MESSAGE('I',"BINGRID");
        MESSAGE('I',"");
        sprintf(msg," Input image:                       %s",in->text);
        MESSAGE('I',msg);
        MESSAGE('I'," ...............");
    }

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


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

    /* set background to black */
    for (j=0; j<(in->nlin * pixsize)+1; j++) {
        for (k=0; k<(in->npix * pixsize)+1; k++) {
            (*out)->data[0][j][k] = (PIXEL)0.0;
        }
    }

    /* insert binary data */
    for (j=0; j<in->nlin; j++) {
        for (k=0; k<in->npix; k++) {
            if (in->data[0][j][k] > in->gmin) {
                for (m=1; m<pixsize; m++) {
                    for (n=1; n<pixsize; n++) {
                        (*out)->data[0][(j*pixsize)+m][(k*pixsize)+n] = (PIXEL)1.0;
                    }
                }
            } else {
                for (m=1; m<pixsize; m++) {
                    for (n=1; n<pixsize; n++) {
                        (*out)->data[0][(j*pixsize)+m][(k*pixsize)+n] = (PIXEL)2.0;
                    }
                }
            }
        }
    }
        
    the_end:
    /* if (TIMES) TIMING(T_EXIT); */
}

---