---

hmosaic.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 mosaics two images horizontally.                           */
/*                                                                            */
/*----------------------------------------------------------------------------*/
/*-Interface Information------------------------------------------------------*/
void HMOSAIC (
IMAGE *in1,     /*  I   Pointer to the left input image.                      */
IMAGE *in2,     /*  I   Pointer to the right input image.                     */
short jbgn,     /*  I   Number of lines (positive or negative) to offset      */
                /*      right image with respect to left image.               */
PIXEL glev,     /*  I   Graylevel to fill otherwise empty pixels.             */
IMAGE **out     /*  O   Address of a pointer to the output image.             */
/*----------------------------------------------------------------------------*/
) { register short i, j, k, nbnd, nlin, npix;
    char msg[SLEN];

    if (TIMES) TIMING(T_HMOSAIC);
    if (NAMES) {
        MESSAGE('I',"");
        MESSAGE('I',"HMOSAIC");
        MESSAGE('I',"");
        sprintf(msg," Left input image:     %s",in1->text);
        MESSAGE('I',msg);
        sprintf(msg," Right input image:    %s",in2->text);
        MESSAGE('I',msg);
        sprintf(msg," Vertical offset:      %d",jbgn);
        MESSAGE('I',msg);
        sprintf(msg," Fill graylevel:     %12.4e",glev);
        MESSAGE('I',msg);
        MESSAGE('I'," ...............");
    }

    /* check input */
    if (!CHECKIMG(in1)) {
        MESSAGE('E'," Can't identify left image.");
        goto the_end;
    } else if (!CHECKIMG(in2)) {
        MESSAGE('E'," Can't identify right image.");
        goto the_end;
    } else if ((nbnd=in1->nbnd) != in2->nbnd) {
        MESSAGE('E'," Number of bands must be identical in both images.");
        goto the_end;
    }

    /* create image of appropriate size */
    nlin = jbgn > 0  ?  max(in1->nlin,in2->nlin+jbgn) : max(in1->nlin-jbgn,in2->nlin);
    npix = in1->npix + in2->npix;
    if (!CHECKIMG(*out)) GETMEM(nbnd,nlin,npix,out);
    if (!*out) goto the_end;

    /* concatenate images */
    for (i=0; i<nbnd; i++) {
        for (j=0; j<nlin; j++) {
            for (k=0; k<npix; k++) {
                (*out)->data[i][j][k] = glev;
            }
        }    
        if (jbgn > 0) {
            for (j=0; j<in1->nlin; j++) {
                for (k=0; k<in1->npix; k++) {
                    (*out)->data[i][j][k]=in1->data[i][j][k];
                }
            }
            for (j=0; j<in2->nlin; j++) {
                for (k=0; k<in2->npix; k++) {
                    (*out)->data[i][j+jbgn][k+in1->npix]=in2->data[i][j][k];
                }
            }
        } else /* (jbgn <= 0) */ {
            for (j=0; j<in1->nlin; j++) {
                for (k=0; k<in1->npix; k++) {
                    (*out)->data[i][j-jbgn][k]=in1->data[i][j][k];
                }
            }
            for (j=0; j<in2->nlin; j++) {
                for (k=0; k<in2->npix; k++) {
                    (*out)->data[i][j][k+in1->npix]=in2->data[i][j][k];
                }
            }
        }
    }

    /* compute new minimum and maximum */
    if (in1->nlin == in2->nlin  &&  jbgn == 0) {
        (*out)->gmin = min(in1->gmin,in2->gmin);
        (*out)->gmax = max(in1->gmax,in2->gmax);
    } else {
        (*out)->gmin = min(glev,min(in1->gmin,in2->gmin));
        (*out)->gmax = max(glev,max(in1->gmax,in2->gmax));
    }

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

---