--- ray/src/util/glare.h 1991/03/19 17:31:35 1.7 +++ ray/src/util/glare.h 2004/01/02 12:51:54 2.6 @@ -1,23 +1,34 @@ -/* Copyright (c) 1991 Regents of the University of California */ - -/* SCCSid "$SunId$ LBL" */ - +/* RCSid: $Id: glare.h,v 2.6 2004/01/02 12:51:54 schorsch Exp $ */ /* * Common data structures for glare source finding routines */ +#ifndef _RAD_GLARE_H_ +#define _RAD_GLARE_H_ #include "standard.h" #include "view.h" #include "color.h" #include "setscan.h" -#define GLAREBR 10.0 /* glare source is this * avg. lum. */ +#ifdef __cplusplus +extern "C" { +#endif -#define SAMPDENS 50 /* default samples per unit in image */ +#define GLAREBR 7.0 /* glare source is this * avg. lum. */ + +#define SAMPDENS 75 /* default samples per unit in image */ #define TSAMPSTEP 10 /* sample step to compute threshold */ -#define SEPS 2 /* sources this close ==> contig. */ +#define SEPS 1 /* sources this close ==> contig. */ +#define SAMIN .005 /* minimum solid angle for source */ +#define MAXBUDDY (4.*sqrt(SAMIN/PI)) /* max separation for pairing */ + +#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 */ @@ -32,18 +43,19 @@ 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 */ float lcos, lsin; /* cosine and sine to left view */ float rcos, rsin; /* cosine and sine to right view */ double sum; /* sum of indirect luminances */ - int n; /* number of values in sum */ + double n; /* number of values in sum */ } *indirect; /* array of indirect illuminances */ struct srcspan { @@ -61,4 +73,28 @@ extern struct source { struct source *next; /* next source in list */ } *donelist; /* finished sources */ -extern double getviewpix(); + +extern long npixinvw; /* number of samples in view */ +extern long npixmiss; /* number of samples missing */ + + /* defined in findglare.c */; +extern void memerr(char *s); +extern int compdir(FVECT vd, int x, int y); +extern double pixsize(int x, int y); + /* defined in glaresrc.c */ +extern void comp_thresh(void); +extern void analyze(void); +extern void absorb_specks(void); + /* defined in glareval.c */ +extern void open_pict(char *fn); +extern void fork_rtrace(char *av[]); +extern void close_pict(void); +extern void done_rtrace(void); +extern void getviewspan(int vv, float *vb); +extern double getviewpix(int vh, int vv); + +#ifdef __cplusplus +} +#endif +#endif /* _RAD_GLARE_H_ */ +