---

findcenter.c

Go to the documentation of this file.

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

void FINDCENTER (
IMAGE  *in     /*  I   Pointer to the input image.                           */
) { register short i, j, k;
    int xsum,ysum,count;
    float xcenter,ycenter;


    /* check input */
    if (!CHECKIMG(in)) {
        MESSAGE('E'," Can't identify image.");
        goto the_end;
    }
    RANGE(in);
    
    
    xsum = 0;
    ysum = 0;
    count = 0;
    
    for (i=0; i<in->nlin; i++) {
        for (j=0; j<in->npix; j++) {
            if (in->data[0][i][j] > in->gmin) {
                xsum = xsum + (j+1);
                ysum = ysum + (i+1);
                count = count + 1;
            }
        }
    }
    
    xcenter = ((float)xsum / (float)count);
    ycenter = ((float)ysum / (float)count);
    
    
    printf("xsum = %d, ysum = %d\n",xsum,ysum);
    printf("Center = (%.2f,%.2f)\n",xcenter,ycenter);
    
    
        
    the_end:
    /* if (TIMES) TIMING(T_EXIT); */
}

---