| 32 |
|
|
| 33 |
|
int rev_orient = 0; /* shall we reverse surface orientation? */ |
| 34 |
|
|
| 35 |
+ |
double lim_graze = 0; /* limit scattering near grazing (deg above) */ |
| 36 |
+ |
|
| 37 |
|
/* Compare incident angles */ |
| 38 |
|
static int |
| 39 |
|
cmp_indir(const void *p1, const void *p2) |
| 180 |
|
fclose(fp); |
| 181 |
|
return(1); |
| 182 |
|
} |
| 183 |
< |
if (rev_orient) { /* reverse Z-axis to face outside */ |
| 183 |
> |
/* check if scatter angle is too low */ |
| 184 |
> |
if (fabs(theta_out - 90.) < lim_graze-FTINY) |
| 185 |
> |
continue; |
| 186 |
> |
if (rev_orient) { /* reverse Z-axis to face outside? */ |
| 187 |
|
theta_out = 180. - theta_out; |
| 188 |
|
phi_out = 360. - phi_out; |
| 189 |
|
} |
| 229 |
|
main(int argc, char *argv[]) |
| 230 |
|
{ |
| 231 |
|
extern int nprocs; |
| 232 |
+ |
int auto_grazing = 0; |
| 233 |
|
const char *symmetry = "0"; |
| 234 |
|
int ninpfiles, totinc; |
| 235 |
|
int a, i; |
| 236 |
< |
/* start header */ |
| 231 |
< |
SET_FILE_BINARY(stdout); |
| 232 |
< |
newheader("RADIANCE", stdout); |
| 233 |
< |
printargs(argc, argv, stdout); |
| 234 |
< |
fputnow(stdout); |
| 236 |
> |
|
| 237 |
|
progname = argv[0]; /* get options */ |
| 238 |
|
for (a = 1; a < argc && argv[a][0] == '-'; a++) |
| 239 |
|
switch (argv[a][1]) { |
| 246 |
|
case 's': |
| 247 |
|
symmetry = argv[++a]; |
| 248 |
|
break; |
| 249 |
+ |
case 'g': |
| 250 |
+ |
if (toupper(argv[a+1][0]) == 'A') |
| 251 |
+ |
auto_grazing = 1; |
| 252 |
+ |
else |
| 253 |
+ |
lim_graze = atof(argv[a+1]); |
| 254 |
+ |
++a; |
| 255 |
+ |
break; |
| 256 |
|
default: |
| 257 |
|
goto userr; |
| 258 |
|
} |
| 264 |
|
inpfile = (PGINPUT *)malloc(sizeof(PGINPUT)*totinc); |
| 265 |
|
if (inpfile == NULL) |
| 266 |
|
return(1); |
| 267 |
< |
for (i = 0; i < ninpfiles; i++) |
| 267 |
> |
if (auto_grazing) |
| 268 |
> |
lim_graze = 90.; |
| 269 |
> |
for (i = 0; i < ninpfiles; i++) { |
| 270 |
|
if (!init_pabopto_inp(i, argv[a+i])) |
| 271 |
|
return(1); |
| 272 |
< |
|
| 272 |
> |
if (auto_grazing && fabs(inpfile[i].theta - 90.) < lim_graze) |
| 273 |
> |
lim_graze = fabs(inpfile[i].theta - 90.); |
| 274 |
> |
} |
| 275 |
|
for (i = ninpfiles; i < totinc; i++) { /* copy for "up" symmetry */ |
| 276 |
|
inpfile[i] = inpfile[i-ninpfiles]; |
| 277 |
|
inpfile[i].phi += 180.; /* invert duplicate data */ |
| 318 |
|
quadrant_sym[inp_coverage]); |
| 319 |
|
#endif |
| 320 |
|
build_mesh(); /* create interpolation */ |
| 321 |
< |
save_bsdf_rep(stdout); /* write it out */ |
| 321 |
> |
SET_FILE_BINARY(stdout); /* start header */ |
| 322 |
> |
newheader("RADIANCE", stdout); |
| 323 |
> |
printargs(argc, argv, stdout); |
| 324 |
> |
fputnow(stdout); |
| 325 |
> |
save_bsdf_rep(stdout); /* complete header + data */ |
| 326 |
|
return(0); |
| 327 |
|
userr: |
| 328 |
< |
fprintf(stderr, "Usage: %s [-t][-n nproc][-s symmetry] meas1.dat meas2.dat .. > bsdf.sir\n", |
| 328 |
> |
fprintf(stderr, "Usage: %s [-t][-n nproc][-s symmetry][-g angle|'A'] meas1.dat meas2.dat .. > bsdf.sir\n", |
| 329 |
|
progname); |
| 330 |
|
return(1); |
| 331 |
|
} |
