---

createmosaic.c

Go to the documentation of this file.

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

/*----------------------------------------------------------------------------*/
/*-General Information--------------------------------------------------------*/
/*                                                                            */
/*   This function assembles multiple images into                             */
/*   a single composite mosaic.  Individual image are specified as                        */
/*                                                                            */
/*----------------------------------------------------------------------------*/
/*-Interface Information------------------------------------------------------*/
int CREATEMOSAIC (
char **imglist, /*  I   Pointer to an array[nimg] containing pointers         */
                /*      to the input image names.                             */
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 *xoff,      /*  I   Address of an array[nimg], containing the starting    */
                /*      line coordinates of the input images in the output.   */
int *yoff,      /*  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, elin, epix;
    char msg[SLEN];
    IMAGE* inimg=NULL;

    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 %d:   Offset: lines: %-4d pixels: %-4d",i,yoff[i],xoff[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',"");
    }

    for (i=0; i<nimg; i++) {
        printf(" Image %d:   Offset: lines: %-4d pixels: %-4d \n",i,yoff[i],xoff[i]);
        printf("                   Contrast Adjustment: gain: %12.4e bias: %12.4e \n",gain[i],bias[i]);
    }

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

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

    /* assemble the mosaic */
    for (l=0; l<nimg; l++) {
        DISK2IMG (imglist[l],-1,&inimg);
        
        elin = min(yoff[l]+inimg->nlin,nlin);
        epix = min(xoff[l]+inimg->npix,npix);
        for (j=yoff[l]; j<elin; j++) {
            for (k=xoff[l]; k<epix; k++) {
                (*out)->data[0][j][k] = (PIXEL)((double)gain[l]*(double)inimg->data[0][j-yoff[l]][k-xoff[l]]+(double)bias[l]);
            }
        }
        
        RELMEM(inimg);
        inimg = NULL;
    }

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

---