---

treering_detect.c

Go to the documentation of this file.

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

extern  PIXEL   thresh;


/*-General Information--------------------------------------------------------*/
/*                                                                            */
/*  Description:  This procedure detects tree ring boundaries and             */
/*      links fragments.                                                      */
/*                                                                            */
/*----------------------------------------------------------------------------*/
/*-Interface Information------------------------------------------------------*/
void TREERING_DETECT (
IMAGE  *in,     /*  I   Pointer to the input image.                           */
IMAGE  **out    /*  O   Address of a pointer to the output image.             */
/*----------------------------------------------------------------------------*/
) { register short i, j, k;
    char   msg[SLEN];
    IMAGE *ringmag=NULL;         /* Tree ring magnitudes prior to thresholding */
    IMAGE *thresh1img=NULL;      /* High threshold, to remove most noise */
    IMAGE *thresh2img=NULL;      /* Low threshold, including noise and tree rings */
    IMAGE *sizethresh1img=NULL;  /* High threshold, also with a size threshold */
    IMAGE *double_thresh_rings=NULL;
    IMAGE *double_thresh_links=NULL;
    IMAGE *linkrings_rings=NULL;
    IMAGE *linkrings_links=NULL;
    IMAGE *combineoutput=NULL;
    short (*tf)(IMAGE*, IMAGE **);
    tf = TFTHRESH;

    /* if (TIMES) TIMING(T_TREERING_DETECT); */
    if (NAMES) {
        MESSAGE('I',"");
        MESSAGE('I',"TREERING_DETECT");
        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->nbnd > 1) {
        MESSAGE('E'," Only first band of input image will be used.");
    }
    
    /* Detect rings and return gradient magnitudes of possible ring boundaries */
    /* Use a gaussian of size 11 and std. deviation 3.0 */
    TREERINGMAG(in,3.0,11,&ringmag);
    if (!CHECKIMG(ringmag)) {
        MESSAGE('E'," Problem determining tree ring boundary magnitudes.");
        goto the_end;
    } else {
        sprintf(ringmag->text, "%s - Ring Magnitude", in->text);
    }

    /* For some reason this is necessary to make the following threshold work */
    /* Without it, the following threshold returns a all-zero image always
        under linux, and after the third run under Solaris ???                */
    RANGE(ringmag);


    /* Create high threshold image */
    //thresh = (PIXEL)9.0;   /* Need to set this automatically */
    thresh = (PIXEL)25.0;   /* Need to set this automatically */
    STRETCH(ringmag, tf, &thresh1img);
    if (!CHECKIMG(thresh1img)) {
        MESSAGE('E'," Problem creating threshold 1.");
        goto the_end;
    } else {
        sprintf(thresh1img->text, "%s - Threshold 1", in->text);
    }

    
    /* Perform size threshold */
    /* Don't pass any fragment smaller than 25 pixels */
    SIZEFILTER(thresh1img, 25, &sizethresh1img);
    if (!CHECKIMG(sizethresh1img)) {
        MESSAGE('E'," Problem creating threshold 1.");
        goto the_end;
    } else {
        sprintf(sizethresh1img->text, "%s - Size Filtered Threshold 1", in->text);
    }


    /* Create low threshold image */
    thresh = (PIXEL)0.001;
    STRETCH(ringmag, tf, &thresh2img);
    if (!CHECKIMG(thresh2img)) {
        MESSAGE('E'," Problem creating threshold 2.");
        goto the_end;
    } else {
        sprintf(thresh2img->text, "%s - Threshold 2", in->text);
    }
    
    
    /* Perform double threshold linking */
    LINK_DOUBLETHRESH(sizethresh1img,thresh2img,5,&double_thresh_links,&double_thresh_rings);
    if (!CHECKIMG(double_thresh_links) || !CHECKIMG(double_thresh_rings)  ||  !CHECKIMG(in)) {
        MESSAGE('E'," Problem with double threshold linking.");
        goto the_end;
    } else {
        sprintf(double_thresh_links->text, "%s - Double Threshold Links", in->text);
        sprintf(double_thresh_rings->text, "%s - Double Threshold Rings", in->text);
    }

    
    /* Combine the output into a single edge map */
    COMPARE (double_thresh_links, double_thresh_rings, OR, out);

goto the_end;
    COMPARE (double_thresh_links, double_thresh_rings, OR, &combineoutput);
    if (!CHECKIMG(combineoutput)) {
        MESSAGE('E'," Problem with combining double threshold linking output.");
        goto the_end;
    } else {
        sprintf(combineoutput->text, "%s - Double Threshold Links + Rings", in->text);
    }


    
    /* Perform tree ring fragment linking */
    LINK_TREERINGS(combineoutput,1,&linkrings_links, &linkrings_rings);
    if (!CHECKIMG(linkrings_links) || !CHECKIMG(linkrings_rings)) {
        MESSAGE('E'," Problem with treering linking.");
        goto the_end;
    } else {
        sprintf(linkrings_links->text, "%s - TreeRing Links", in->text);
        sprintf(linkrings_rings->text, "%s - TreeRing Rings", in->text);
    }
    
    /* Combine the output into a single edge map */
    COMPARE (linkrings_links, linkrings_rings, OR, out);
    if (!CHECKIMG(*out)) {
        MESSAGE('E'," Problem with combining treering linking output.");
        goto the_end;
    }

    


    the_end:

    /* if (TIMES)   TIMING(T_EXIT); */
    if (CHECKIMG(ringmag))  RELIMG(&ringmag);
    if (CHECKIMG(thresh1img))  RELIMG(&thresh1img);
    if (CHECKIMG(thresh2img))  RELIMG(&thresh2img);
    if (CHECKIMG(sizethresh1img))  RELIMG(&sizethresh1img);
    if (CHECKIMG(double_thresh_rings))  RELIMG(&double_thresh_rings);
    if (CHECKIMG(double_thresh_links))  RELIMG(&double_thresh_links);
    if (CHECKIMG(linkrings_rings))  RELIMG(&linkrings_rings);
    if (CHECKIMG(linkrings_links))  RELIMG(&linkrings_links);
    if (CHECKIMG(combineoutput))  RELIMG(&combineoutput);
}

---