---

dupl.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 enlarges an image by replicating bands, lines, and pixels. */
/*                                                                            */
/*----------------------------------------------------------------------------*/
/*-Interface Information------------------------------------------------------*/
void DUPL (
IMAGE *in,      /*  I   Pointer to the input image.                           */
short ibgn,     /*  I   First replicated band.                                */
short jbgn,     /*  I   First replicated line.                                */
short kbgn,     /*  I   First replicated pixel.                               */
short ifact,    /*  I   Band replication factor.                              */
short jfact,    /*  I   Line replication factor.                              */
short kfact,    /*  I   Pixel replication factor.                             */
short isize,    /*  I   Number of replicated bands.                           */
short jsize,    /*  I   Number of replicated lines.                           */
short ksize,    /*  I   Number of replicated pixels.                          */
IMAGE **out     /*  O   Address of a pointer to the output image.             */
/*----------------------------------------------------------------------------*/
) { register short i, j, k, l, m, n;
    char msg[SLEN];

    if (TIMES) TIMING(T_DUPL);
    if (NAMES) {
        MESSAGE('I',"");
        MESSAGE('I',"DUPL");
        MESSAGE('I',"");
        sprintf(msg," Input image:                                   %s",in->text);
        MESSAGE('I',msg);
        sprintf(msg," First replicated band, line, pixel:            %d, %d, %d",ibgn+1,jbgn+1,kbgn+1);
        MESSAGE('I',msg);  
        sprintf(msg," Replication factor for bands, lines, pixels:   %d, %d, %d",ifact,jfact,kfact);
        MESSAGE('I',msg);
        sprintf(msg," Number of replicated bands, lines, pixels:     %d, %d, %d",isize,jsize,ksize);
        MESSAGE('I',msg); 
        MESSAGE('I'," ...............");
    }

    /* check input */
    if (!CHECKIMG(in)) {
        MESSAGE('E'," Can't identify image.");
        goto the_end;
    } else if (0 > ibgn  ||  ibgn+isize > in->nbnd) {
        MESSAGE('E'," Block band range to be duplicated must be located completely inside input image.");
        goto the_end;
    } else if (0 > jbgn  ||  jbgn+jsize > in->nlin) {
        MESSAGE('E'," Block line range to be duplicated must be located completely inside input image.");
        goto the_end;
    } else if (0 > kbgn  ||  kbgn+ksize > in->npix) {
        MESSAGE('E'," Block pixel range to be duplicated must be located completely inside input image.");
        goto the_end;
    } else if (ifact <= 0) {
        MESSAGE('E'," Band duplication factor must be greater than zero.");
        goto the_end;
    } else if (jfact <= 0) {
        MESSAGE('E'," Line duplication factor must be greater than zero.");
        goto the_end;
    } else if (kfact <= 0) {
        MESSAGE('E'," Pixel duplication factor must be greater than zero.");
        goto the_end;
    } else if (isize <= 0) {
        MESSAGE('E'," Size of block in bands to be duplicated must be greater then zero.");
        goto the_end;
    } else if (jsize <= 0) {
        MESSAGE('E'," Size of block in lines to be duplicated must be greater then zero.");
        goto the_end;
    } else if (ksize <= 0) {
        MESSAGE('E'," Size of block in pixels/line to be duplicated must be greater then zero.");
        goto the_end;
    }

    /* create image of appropriate size */
    if (!CHECKIMG(*out)) GETMEM(isize*ifact,jsize*jfact,ksize*kfact,out);
    if (!*out) goto the_end;

    /* duplicate the image data */
    for (i=ibgn; i<ibgn+isize; i++) {
        for (j=jbgn; j<jbgn+jsize; j++) {
            for (k=kbgn; k<kbgn+ksize; k++) {
                for (l=0; l<ifact; l++) {
                    for (m=0; m<jfact; m++) {
                        for (n=0; n<kfact; n++) {
                            (*out)->data[i*ifact+l][j*jfact+m][k*kfact+n] = in->data[i][j][k];
                        }
                    }
                }
            }
        }
    }
    (*out)->gmin = in->gmin;
    (*out)->gmax = in->gmax;

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

---