---

insert.c

Go to the documentation of this file.

#include        "sadie.h"

/*-Copyright Information------------------------------------------------------*/
/* Copyright (c) 1988 by the University of Arizona Digital Image Analysis Lab */
/*----------------------------------------------------------------------------*/
/*-General Information--------------------------------------------------------*/
/*                                                                            */
/*   This function inserts one image into another by replacing overlapping    */
/*   pixels in the "background" image with those from the "foreground" image. */
/*                                                                            */
/*----------------------------------------------------------------------------*/
/*-Interface Information------------------------------------------------------*/
void INSERT (
IMAGE *in1,     /*  I   Pointer to the background image.                      */
IMAGE *in2,     /*  I   Pointer to the foreground image.                      */
short jbgn,     /*  I   Line offset in the background image.                  */
short kbgn,     /*  I   Pixel offset in the background image.                 */
IMAGE **out     /*  O   Address of a pointer to the output image.             */
                /*      The background image may be overwritten (out = &in1). */
/*----------------------------------------------------------------------------*/
) { register short i, j, k, elin, epix;
    char msg[SLEN];

    if (TIMES) TIMING(T_INSERT);
    if (NAMES) {
        MESSAGE('I',"");
        MESSAGE('I',"INSERT");
        MESSAGE('I',"");
        sprintf(msg," Background input image:   %s",in1->text);
        MESSAGE('I',msg);
        sprintf(msg," Foreground input image:   %s",in2->text);
        MESSAGE('I',msg);
        sprintf(msg," Offset: lines:            %d",jbgn);
        MESSAGE('I',msg);
        sprintf(msg,"         pixels:           %d",kbgn);
        MESSAGE('I',msg);
        MESSAGE('I'," ...............");
    }

    /* check input */
    if (!CHECKIMG(in1)) {
        MESSAGE('E'," Can't identify background image.");
        goto the_end;
    } else if (!CHECKIMG(in2)) {
        MESSAGE('E'," Can't identify foreground image.");
        goto the_end;
    } else if (in1->nbnd != in2->nbnd) {
        MESSAGE('E'," Number of bands must be identical in both images.");
        goto the_end;
    } else if (0 > jbgn  ||  jbgn >= in1->nlin  ||  jbgn+in2->nlin > in1->nlin) {
        MESSAGE('E'," Foreground image must be located completely inside the background image in the line dimension.");
        goto the_end;
    } else if (0 > kbgn  ||  kbgn >= in1->npix  ||  kbgn+in2->npix > in1->npix) {
        MESSAGE('E'," Foreground image must be located completely inside the background image in the pixel dimension.");
        goto the_end;
    }

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

    /* copy image */
    for (i=0; i<in1->nbnd; i++) {
        for (j=0; j<in1->nlin; j++) {
            for (k=0; k<in1->npix; k++) {
                (*out)->data[i][j][k] = in1->data[i][j][k];
            }
        }
    }

    /* replace overlapping pixels */
    elin = min(in1->nlin,jbgn+in2->nlin);
    epix = min(in1->npix,kbgn+in2->npix);
    for (i=0; i<in1->nbnd; i++) {
        for (j=jbgn; j<elin; j++) {
            for (k=kbgn; k<epix; k++) {
                (*out)->data[i][j][k] = in2->data[i][j-jbgn][k-kbgn];
            }
        }
    }

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

---