---

radialavg.c

Go to the documentation of this file.

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


PIXEL globarray_original[50][100];
PIXEL globarray_percent[50][100];
PIXEL globarray_normalized[50][100];
int globframecount;
int globfreqcount;


/*----------------------------------------------------------------------------*/
/*-General Information--------------------------------------------------------*/
/*                                                                            */
/*   This function computes the radial average intensities around the center  */
/*   of the input image.  It is useful for representing a 2-D power spectrum  */
/*   in 1-D.                                                                  */
/*----------------------------------------------------------------------------*/
/*-Interface Information------------------------------------------------------*/
void RADIALAVG (
IMAGE  *in     /*  I   Pointer to the input image.                           */
/*----------------------------------------------------------------------------*/
) { register short i, j, k;
    char msg[SLEN];
    int length, sublength;
    PIXEL *value_array=NULL, *value_subarray=NULL;
    int *count_array=NULL;
    int radius;
    double delta_y, delta_x;
    PIXEL maxval, tempval;

    /* if (TIMES) TIMING(T_RADIALAVG); */
    if (NAMES) {
        MESSAGE('I',"");
        MESSAGE('I',"RADIALAVG");
        MESSAGE('I',"");
        sprintf(msg," Input image:                       %s",in->text);
        MESSAGE('I',msg);
        MESSAGE('I'," ...............");
    }

    /* check input */
    if (!CHECKIMG(in)) {
        MESSAGE('E'," Can't identify image.");
        goto the_end;
    } else if (in->npix != in->nlin) {
        MESSAGE('E'," Input image must have the same width and height.");
        goto the_end;
    } else if (in->nbnd > 1) {
        MESSAGE('E'," Input image must be single band!");
        goto the_end;
    }
    
    length = (int) rnd((double)in->nlin / 2.0);
    
    value_array = malloc((length+1) * sizeof(PIXEL));
    count_array = malloc((length+1) * sizeof(int));
    
    /* Initialize arrays */
    for (i=0; i<=length; i++) {
        value_array[i] = (PIXEL)0.0;
        count_array[i] = 0;
    }


    for (j=0; j<in->nlin; j++) {
        delta_y = (double)(j-length);
        for (k=0; k<in->npix; k++) {

            delta_x = (double)(k-length);

            radius = (int) rnd(sqrt((delta_x)*(delta_x) + (delta_y)*(delta_y)));
            
            if (radius <= length) {
                value_array[radius] += in->data[0][j][k];
                count_array[radius] += 1;
            }
        }
    }

        

    for (i=0; i<=length; i++) {
        if (count_array[i] > 0) {
            value_array[i] /= (PIXEL) count_array[i];
        }
        globarray_original[globframecount][i] = value_array[i];
    }



    maxval = value_array[0];
    for (i=0; (i<=length)&&(i<100); i++) {
        globarray_normalized[globframecount][i] = value_array[i]/maxval;
        globarray_percent[globframecount][i] = value_array[i]/maxval;
        //printf("%2.8f ",globarray_normalized[globframecount][i]);
    }
    //printf("\n");

    
    //sublength = (int) rnd((double)length / 10.0);
    //value_subarray = malloc((sublength) * sizeof(PIXEL));
    //for (i=0; i<sublength; i++) {
    //    value_subarray[i] = (PIXEL)0.0;
    //    
    //    for (j=0; j<10; j++) {
    //        value_subarray[i] += value_array[j+(i*10)]/maxval;
    //    }
    //    
    //    printf("%2.8f ",value_subarray[i]/10.0);
    //}
    //printf("\n");

//
//Divide each value by the DC value to normalize --> convert to percent ??
//
//Print out array with each row containing the results for an image and each column containing a radius?
//


//    PLOT(value_array,length,LOG);


    globframecount++;
    if (globfreqcount > length) {
        globfreqcount = (int)min(length,100);
    }

    the_end:
    /* if (TIMES) TIMING(T_EXIT); */
    if (count_array) free(count_array);
    if (value_array) free(value_array);
    if (value_subarray) free(value_subarray);
}

---