---

gainadjmos.c

Go to the documentation of this file.

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

/*----------------------------------------------------------------------------*/
/*-General Information--------------------------------------------------------*/
/*                                                                            */
/*   This function adjusts the gain of, and assembles multiple images into    */
/*   a single composite mosaic.  It is based on MOSAIC and TFLINEAR from      */
/*   the SADIE library.                                                       */
/*                                                                            */
/*----------------------------------------------------------------------------*/
/*-Interface Information------------------------------------------------------*/
int GAINADJMOS (
IMAGE **in,     /*  I   Pointer to an array[nimg] containing pointers         */
                /*      to the input images.                                  */
int nimg,       /*  I   Number of input images.                               */
int nlin,       /*  I   Number of lines in the output image.                  */
int npix,       /*  I   Number of pixels/line in the output image.            */
int *jbgn,      /*  I   Address of an array[nimg], containing the starting    */
                /*      line coordinates of the input images in the output.   */
int *kbgn,      /*  I   Address of an array[nimg], containing the starting    */
                /*      pixel coordinates of the input images in the output.  */
PIXEL *bias,    /*  I   Address of an array[nimg], containing the bias        */
                /*      offsets of the input images relative to the output.   */
PIXEL *gain,    /*  I   Address of an array[nimg], containing the gain        */
                /*      offsets of the input images relative to the output.   */
PIXEL glev,     /*  I   Graylevel to fill otherwise empty pixels.             */
IMAGE **out     /*  O   Address of a pointer to the output image.             */
/*----------------------------------------------------------------------------*/
) { register int i, j, k, l, nbnd, elin, epix;
    char msg[SLEN];

    if (TIMES) TIMING(T_MOSAIC);
    if (NAMES) {
        MESSAGE('I',"");
        MESSAGE('I',"----------");
        MESSAGE('I',"");
        MESSAGE('I',"GAINADJMOS");
        MESSAGE('I',"");
        sprintf(msg," Number of input images: %d",nimg);
        MESSAGE('I',msg);
        sprintf(msg," Image size: lines:      %d",nlin);
        MESSAGE('I',msg);
        sprintf(msg,"             pixels:     %d",npix);
        MESSAGE('I',msg);
        for (i=0; i<nimg; i++) {
            sprintf(msg," Image: %s   Offset: lines: %-4d pixels: %-4d",in[i]->text,jbgn[i],kbgn[i]);
            MESSAGE('I',msg);
            sprintf(msg,"                   Contrast Adjustment: gain: %12.4e bias: %12.4e",gain[i],bias[i]);
            MESSAGE('I',msg);
        }
        sprintf(msg," Fill graylevel:       %12.4e",glev);
        MESSAGE('I',msg);
        MESSAGE('I'," ...............");
        MESSAGE('I',"");
    }

    /* check input */
    for (i=0; i<nimg; i++) {
        if (!CHECKIMG(in[i])) {
            MESSAGE('E'," Can't identify input image(s).");
            goto the_end;
        }
    }
    for (nbnd=in[0]->nbnd,i=1; i<nimg; i++) {
        if (nbnd != in[i]->nbnd) {
            MESSAGE('E'," Number of bands must be identical in all images.");
            goto the_end;
        }
    }
    for (i=0; i<nimg; i++) {
        if (0 > jbgn[i]  ||  jbgn[i]+in[i]->nlin > nlin) {
            MESSAGE('E'," Input image(s) must be located completely inside output image in the line dimension.");
            goto the_end;
        } else if (0 > kbgn[i]  ||  kbgn[i]+in[i]->npix > npix) {
            MESSAGE('E'," Input image(s) must be located completely inside output image in the pixel dimension.");
            goto the_end;
        }
    }


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

    /* fill the background */
    for (i=0; i<nbnd; i++) {
        for (j=0; j<nlin; j++) {
            for (k=0; k<npix; k++) {
                (*out)->data[i][j][k] = glev;
            }
        }
    }

    /* assemble the mosaic */
    for (i=0; i<nbnd; i++) {
        for (l=0; l<nimg; l++) {
            elin = min(jbgn[l]+in[l]->nlin,nlin);
            epix = min(kbgn[l]+in[l]->npix,npix);
            for (j=jbgn[l]; j<elin; j++) {
                for (k=kbgn[l]; k<epix; k++) {
                    (*out)->data[i][j][k] = (PIXEL)((double)gain[l]*(double)in[l]->data[i][j-jbgn[l]][k-kbgn[l]]+(double)bias[l]);
                }
            }
        }
    }

    the_end:
    if (TIMES) TIMING(T_EXIT);
}

---