--- ray/src/util/glare.h 1991/04/12 10:28:39 1.14 +++ ray/src/util/glare.h 2003/02/22 02:07:30 2.3 @@ -1,7 +1,4 @@ -/* Copyright (c) 1991 Regents of the University of California */ - -/* SCCSid "$SunId$ LBL" */ - +/* RCSid: $Id: glare.h,v 2.3 2003/02/22 02:07:30 greg Exp $ */ /* * Common data structures for glare source finding routines */ @@ -13,7 +10,7 @@ #define GLAREBR 7.0 /* glare source is this * avg. lum. */ -#define SAMPDENS 50 /* default samples per unit in image */ +#define SAMPDENS 75 /* default samples per unit in image */ #define TSAMPSTEP 10 /* sample step to compute threshold */ #define SEPS 1 /* sources this close ==> contig. */ @@ -23,6 +20,9 @@ #define TOOSMALL(s) ((s)->brt*(s)->dom < threshold*SAMIN) +#define SABIG .025 /* solid angle of splittable source */ +#define LCORR .12 /* linearity of splittable source */ + extern VIEW ourview; /* our view */ extern VIEW pictview; /* picture view */ extern VIEW leftview, rightview; /* leftmost and rightmost views */ @@ -37,11 +37,12 @@ extern ANGLE glarang[]; /* glare calculation angles * extern int nglarangs; extern double maxtheta; /* maximum glare angle (in radians) */ extern int hsize; /* horizontal size */ -extern int hlim; /* horizontal limit of central view */ #define nglardirs (2*nglarangs+1) #define vsize (sampdens-1) -#define h_theta(h) (-(double)(h)/(double)sampdens) +#define hscale(v) sqrt((double)(sampdens*sampdens - (v)*(v))) +#define hlim(v) (int)(maxtheta*hscale(v)) +#define h_theta(h,v) (-(h)/hscale(v)) extern struct illum { float theta; /* glare direction */ @@ -66,12 +67,8 @@ extern struct source { struct source *next; /* next source in list */ } *donelist; /* finished sources */ -typedef struct { - double err; /* cumulative error */ - double prob; /* target probability */ -} SPANERR; /* probability record for computing spans */ - extern double getviewpix(); +extern double pixsize(); extern long npixinvw; /* number of samples in view */ extern long npixmiss; /* number of samples missing */