---

real2com.c

Go to the documentation of this file.

#include        "sadie.h"

/*-Copyright Information------------------------------------------------------*/
/* Copyright (c) 1990 by the University of Arizona Digital Image Analysis Lab */
/*----------------------------------------------------------------------------*/
/*-General Information--------------------------------------------------------*/
/*                                                                            */
/*   This function combines two "real" images into one "complex" image.       */
/*                                                                            */
/*----------------------------------------------------------------------------*/
/*-Interface Information------------------------------------------------------*/
void REAL2COM (
IMAGE *in1,     /*  I   Pointer to the first input image.                     */
IMAGE *in2,     /*  I   Pointer to the second input image.                    */
short option,   /*  I   Switch, indicating the type of real images:           */
                /*      REIMAG   -   real and imaginary images are given.     */
                /*      MAGPHA   -   magnitude and phase images are given.    */
IMAGE **out     /*  O   Address of a pointer to the output image.             */
/*----------------------------------------------------------------------------*/
) { register short i, j, k, n;
    char msg[SLEN];

    if (TIMES) TIMING(T_REAL2COM);
    if (NAMES) {
        MESSAGE('I',"");
        MESSAGE('I',"REAL2COM");
        MESSAGE('I',"");
        sprintf(msg," First input image:    %s",in1->text);
        MESSAGE('I',msg);
        sprintf(msg," Second input image:   %s",in2->text);
        MESSAGE('I',msg);
        if (option == REIMAG) {
            MESSAGE('I'," Real and imaginary images are given");
        } else if (option == MAGPHA) {
            MESSAGE('I'," Magnitude and phase images are given");
        }
        MESSAGE('I'," ...............");
    }

    /* check input */
    if (!CHECKIMG(in1)) {
        MESSAGE('E'," Can't identify first image.");
        goto the_end;
    } else if (!CHECKIMG(in2)) {
        MESSAGE('E'," Can't identify second 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 (in1->nlin != in2->nlin) {
        MESSAGE('E'," Number of lines/band must be identical in both images.");
        goto the_end;
    } else if (in1->npix != in2->npix) {
        MESSAGE('E'," Number of pixels/line must be identical in both images.");
        goto the_end;
    } else if (option != REIMAG  &&  option != MAGPHA) {
        MESSAGE('E'," Option must be either REIMAG or MAGPHA.");
        goto the_end;
    }

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

    if (option == REIMAG) {
        /* combine real and imaginary parts */
        for (i=0; i<in1->nbnd; i++) {
            for (j=0; j<in1->nlin; j++) {
                for (k=0,n=0; k<in1->npix; k+=1,n+=2) {
                    (*out)->data[i][j][n]   = in1->data[i][j][k];
                    (*out)->data[i][j][n+1] = in2->data[i][j][k];
                }
            }
        }
    } else /* (option == MAGPHA) */ {
        /* combine magnitude and phase parts */
        for (i=0; i<in1->nbnd; i++) {
            for (j=0; j<in1->nlin; j++) {
                for (k=0,n=0; k<in1->npix; k+=1,n+=2) {
                    (*out)->data[i][j][n]   = (PIXEL)(double)in1->data[i][j][k]*cos((double)in2->data[i][j][k]);
                    (*out)->data[i][j][n+1] = (PIXEL)(double)in1->data[i][j][k]*sin((double)in2->data[i][j][k]);
                }
            }
        }
    }       

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

---