| 14 |
|
#include <math.h> |
| 15 |
|
#include "rtio.h" |
| 16 |
|
#include "platform.h" |
| 17 |
+ |
#include "paths.h" |
| 18 |
|
#include "bsdfrep.h" |
| 18 |
– |
/* global argv[0] */ |
| 19 |
– |
char *progname; |
| 19 |
|
|
| 20 |
|
typedef struct { |
| 21 |
|
const char *fname; /* input file path */ |
| 31 |
|
|
| 32 |
|
int rev_orient = 0; /* shall we reverse surface orientation? */ |
| 33 |
|
|
| 34 |
+ |
double lim_graze = 0; /* limit scattering near grazing (deg above) */ |
| 35 |
+ |
|
| 36 |
|
/* Compare incident angles */ |
| 37 |
|
static int |
| 38 |
|
cmp_indir(const void *p1, const void *p2) |
| 179 |
|
fclose(fp); |
| 180 |
|
return(1); |
| 181 |
|
} |
| 182 |
< |
if (rev_orient) { /* reverse Z-axis to face outside */ |
| 182 |
> |
/* check if scatter angle is too low */ |
| 183 |
> |
if (fabs(theta_out - 90.) < lim_graze-FTINY) |
| 184 |
> |
continue; |
| 185 |
> |
if (rev_orient) { /* reverse Z-axis to face outside? */ |
| 186 |
|
theta_out = 180. - theta_out; |
| 187 |
|
phi_out = 360. - phi_out; |
| 188 |
|
} |
| 228 |
|
main(int argc, char *argv[]) |
| 229 |
|
{ |
| 230 |
|
extern int nprocs; |
| 231 |
+ |
static char gval_buf[16]; |
| 232 |
+ |
char * auto_grazing = NULL; |
| 233 |
|
const char *symmetry = "0"; |
| 234 |
|
int ninpfiles, totinc; |
| 235 |
|
int a, i; |
| 236 |
< |
/* start header */ |
| 237 |
< |
SET_FILE_BINARY(stdout); |
| 238 |
< |
newheader("RADIANCE", stdout); |
| 233 |
< |
printargs(argc, argv, stdout); |
| 234 |
< |
fputnow(stdout); |
| 235 |
< |
progname = argv[0]; /* get options */ |
| 236 |
> |
/* set global progname */ |
| 237 |
> |
fixargv0(argv[0]); |
| 238 |
> |
/* get options */ |
| 239 |
|
for (a = 1; a < argc && argv[a][0] == '-'; a++) |
| 240 |
|
switch (argv[a][1]) { |
| 241 |
|
case 't': |
| 247 |
|
case 's': |
| 248 |
|
symmetry = argv[++a]; |
| 249 |
|
break; |
| 250 |
+ |
case 'g': |
| 251 |
+ |
if (toupper(argv[a+1][0]) == 'A') |
| 252 |
+ |
auto_grazing = argv[a+1] = gval_buf; |
| 253 |
+ |
else |
| 254 |
+ |
lim_graze = atof(argv[a+1]); |
| 255 |
+ |
++a; |
| 256 |
+ |
break; |
| 257 |
|
default: |
| 258 |
|
goto userr; |
| 259 |
|
} |
| 265 |
|
inpfile = (PGINPUT *)malloc(sizeof(PGINPUT)*totinc); |
| 266 |
|
if (inpfile == NULL) |
| 267 |
|
return(1); |
| 268 |
< |
for (i = 0; i < ninpfiles; i++) |
| 268 |
> |
if (auto_grazing) |
| 269 |
> |
lim_graze = 90.; |
| 270 |
> |
for (i = 0; i < ninpfiles; i++) { |
| 271 |
|
if (!init_pabopto_inp(i, argv[a+i])) |
| 272 |
|
return(1); |
| 273 |
< |
|
| 273 |
> |
if (auto_grazing && fabs(inpfile[i].theta - 90.) < lim_graze) |
| 274 |
> |
lim_graze = fabs(inpfile[i].theta - 90.); |
| 275 |
> |
} |
| 276 |
> |
if (auto_grazing) |
| 277 |
> |
sprintf(auto_grazing, "%.2f", lim_graze); |
| 278 |
|
for (i = ninpfiles; i < totinc; i++) { /* copy for "up" symmetry */ |
| 279 |
|
inpfile[i] = inpfile[i-ninpfiles]; |
| 280 |
|
inpfile[i].phi += 180.; /* invert duplicate data */ |
| 281 |
|
inpfile[i].up_phi -= 180.; |
| 282 |
|
set_grid_pos(&inpfile[i]); /* grid location for sorting */ |
| 283 |
|
} |
| 284 |
< |
qsort(inpfile, ninpfiles, sizeof(PGINPUT), cmp_indir); |
| 284 |
> |
qsort(inpfile, totinc, sizeof(PGINPUT), cmp_indir); |
| 285 |
|
/* compile measurements */ |
| 286 |
|
for (i = 0; i < totinc; i++) |
| 287 |
|
if (!add_pabopto_inp(i)) |
| 321 |
|
quadrant_sym[inp_coverage]); |
| 322 |
|
#endif |
| 323 |
|
build_mesh(); /* create interpolation */ |
| 324 |
< |
save_bsdf_rep(stdout); /* write it out */ |
| 324 |
> |
SET_FILE_BINARY(stdout); /* start header */ |
| 325 |
> |
newheader("RADIANCE", stdout); |
| 326 |
> |
printargs(argc, argv, stdout); |
| 327 |
> |
fputnow(stdout); |
| 328 |
> |
save_bsdf_rep(stdout); /* complete header + data */ |
| 329 |
|
return(0); |
| 330 |
|
userr: |
| 331 |
< |
fprintf(stderr, "Usage: %s [-t][-n nproc][-s symmetry] meas1.dat meas2.dat .. > bsdf.sir\n", |
| 331 |
> |
fprintf(stderr, "Usage: %s [-t][-n nproc][-s symmetry][-g angle|'A'] meas1.dat meas2.dat .. > bsdf.sir\n", |
| 332 |
|
progname); |
| 333 |
|
return(1); |
| 334 |
|
} |
