ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/cv/bsdf2ttree.c
Revision: 2.6
Committed: Fri Nov 9 02:25:32 2012 UTC (11 years, 4 months ago) by greg
Content type: text/plain
Branch: MAIN
Changes since 2.5: +26 -45 lines
Log Message:
Made code a little more readable

File Contents

# User Rev Content
1 greg 2.1 #ifndef lint
2 greg 2.6 static const char RCSid[] = "$Id: bsdf2ttree.c,v 2.5 2012/11/09 02:16:29 greg Exp $";
3 greg 2.1 #endif
4     /*
5     * Load measured BSDF interpolant and write out as XML file with tensor tree.
6     *
7     * G. Ward
8     */
9    
10     #define _USE_MATH_DEFINES
11     #include <stdio.h>
12     #include <stdlib.h>
13     #include <math.h>
14     #include "platform.h"
15     #include "bsdfrep.h"
16     /* global argv[0] */
17     char *progname;
18     /* percentage to cull (<0 to turn off) */
19     int pctcull = 90;
20     /* sampling order */
21     int samp_order = 6;
22    
23     /* Interpolate and output isotropic BSDF data */
24     static void
25     interp_isotropic()
26     {
27     const int sqres = 1<<samp_order;
28     FILE *ofp = NULL;
29     char cmd[128];
30     int ix, ox, oy;
31     FVECT ivec, ovec;
32 greg 2.4 float bsdf;
33 greg 2.1 #if DEBUG
34     fprintf(stderr, "Writing isotropic order %d ", samp_order);
35     if (pctcull >= 0) fprintf(stderr, "data with %d%% culling\n", pctcull);
36     else fputs("raw data\n", stderr);
37     #endif
38     if (pctcull >= 0) { /* begin output */
39 greg 2.4 sprintf(cmd, "rttree_reduce -h -a -ff -r 3 -t %d -g %d",
40 greg 2.1 pctcull, samp_order);
41     fflush(stdout);
42     ofp = popen(cmd, "w");
43     if (ofp == NULL) {
44     fprintf(stderr, "%s: cannot create pipe to rttree_reduce\n",
45     progname);
46     exit(1);
47     }
48     SET_FILE_BINARY(ofp);
49     } else
50     fputs("{\n", stdout);
51     /* run through directions */
52     for (ix = 0; ix < sqres/2; ix++) {
53     RBFNODE *rbf;
54     SDsquare2disk(ivec, (ix+.5)/sqres, .5);
55     ivec[2] = input_orient *
56     sqrt(1. - ivec[0]*ivec[0] - ivec[1]*ivec[1]);
57     rbf = advect_rbf(ivec);
58     for (ox = 0; ox < sqres; ox++)
59     for (oy = 0; oy < sqres; oy++) {
60     SDsquare2disk(ovec, (ox+.5)/sqres, (oy+.5)/sqres);
61     ovec[2] = output_orient *
62     sqrt(1. - ovec[0]*ovec[0] - ovec[1]*ovec[1]);
63 greg 2.5 bsdf = eval_rbfrep(rbf, ovec) * output_orient/ovec[2];
64 greg 2.1 if (pctcull >= 0)
65     fwrite(&bsdf, sizeof(bsdf), 1, ofp);
66     else
67     printf("\t%.3e\n", bsdf);
68     }
69 greg 2.3 if (rbf != NULL)
70     free(rbf);
71 greg 2.1 }
72     if (pctcull >= 0) { /* finish output */
73     if (pclose(ofp)) {
74     fprintf(stderr, "%s: error running '%s'\n",
75     progname, cmd);
76     exit(1);
77     }
78     } else {
79     for (ix = sqres*sqres*sqres/2; ix--; )
80     fputs("\t0\n", stdout);
81     fputs("}\n", stdout);
82     }
83     }
84    
85     /* Interpolate and output anisotropic BSDF data */
86     static void
87     interp_anisotropic()
88     {
89     const int sqres = 1<<samp_order;
90     FILE *ofp = NULL;
91     char cmd[128];
92     int ix, iy, ox, oy;
93     FVECT ivec, ovec;
94 greg 2.4 float bsdf;
95 greg 2.1 #if DEBUG
96     fprintf(stderr, "Writing anisotropic order %d ", samp_order);
97     if (pctcull >= 0) fprintf(stderr, "data with %d%% culling\n", pctcull);
98     else fputs("raw data\n", stderr);
99     #endif
100     if (pctcull >= 0) { /* begin output */
101 greg 2.4 sprintf(cmd, "rttree_reduce -h -a -ff -r 4 -t %d -g %d",
102 greg 2.1 pctcull, samp_order);
103     fflush(stdout);
104     ofp = popen(cmd, "w");
105     if (ofp == NULL) {
106     fprintf(stderr, "%s: cannot create pipe to rttree_reduce\n",
107     progname);
108     exit(1);
109     }
110     } else
111     fputs("{\n", stdout);
112     /* run through directions */
113     for (ix = 0; ix < sqres; ix++)
114     for (iy = 0; iy < sqres; iy++) {
115     RBFNODE *rbf;
116     SDsquare2disk(ivec, (ix+.5)/sqres, (iy+.5)/sqres);
117     ivec[2] = input_orient *
118     sqrt(1. - ivec[0]*ivec[0] - ivec[1]*ivec[1]);
119     rbf = advect_rbf(ivec);
120     for (ox = 0; ox < sqres; ox++)
121     for (oy = 0; oy < sqres; oy++) {
122     SDsquare2disk(ovec, (ox+.5)/sqres, (oy+.5)/sqres);
123     ovec[2] = output_orient *
124     sqrt(1. - ovec[0]*ovec[0] - ovec[1]*ovec[1]);
125 greg 2.5 bsdf = eval_rbfrep(rbf, ovec) * output_orient/ovec[2];
126 greg 2.1 if (pctcull >= 0)
127     fwrite(&bsdf, sizeof(bsdf), 1, ofp);
128     else
129     printf("\t%.3e\n", bsdf);
130     }
131 greg 2.3 if (rbf != NULL)
132     free(rbf);
133 greg 2.1 }
134     if (pctcull >= 0) { /* finish output */
135     if (pclose(ofp)) {
136     fprintf(stderr, "%s: error running '%s'\n",
137     progname, cmd);
138     exit(1);
139     }
140     } else
141     fputs("}\n", stdout);
142     }
143    
144 greg 2.5 /* Output XML prologue to stdout */
145     static void
146     xml_prologue(int ac, char *av[])
147     {
148     static const char *bsdf_type[4] = {
149     "Reflection Back",
150     "Transmission Back",
151     "Transmission Front",
152     "Reflection Front"
153     };
154    
155 greg 2.6 puts("<?xml version=\"1.0\" encoding=\"UTF-8\"?>");
156     puts("<WindowElement xmlns=\"http://windows.lbl.gov\" xmlns:xsi=\"http://www.w3.org/2001/XMLSchema-instance\" xsi:schemaLocation=\"http://windows.lbl.gov/BSDF-v1.4.xsd\">");
157 greg 2.5 fputs("<!-- File produced by:", stdout);
158     while (ac-- > 0) {
159     fputc(' ', stdout);
160     fputs(*av++, stdout);
161     }
162     puts(" -->");
163 greg 2.6 puts("<WindowElementType>System</WindowElementType>");
164     puts("<FileType>BSDF</FileType>");
165     puts("<Optical>");
166     puts("<Layer>");
167     puts("\t<Material>");
168     puts("\t\t<Name>Name</Name>");
169     puts("\t\t<Manufacturer>Manufacturer</Manufacturer>");
170     puts("\t\t<DeviceType>Other</DeviceType>");
171     puts("\t</Material>");
172 greg 2.5 puts("\t<DataDefinition>");
173     printf("\t\t<IncidentDataStructure>TensorTree%c</IncidentDataStructure>\n",
174     single_plane_incident ? '3' : '4');
175     puts("\t</DataDefinition>");
176 greg 2.6 puts("\t<WavelengthData>");
177     puts("\t\t<LayerNumber>System</LayerNumber>");
178     puts("\t\t<Wavelength unit=\"Integral\">Visible</Wavelength>");
179     puts("\t\t<SourceSpectrum>CIE Illuminant D65 1nm.ssp</SourceSpectrum>");
180     puts("\t\t<DetectorSpectrum>ASTM E308 1931 Y.dsp</DetectorSpectrum>");
181     puts("\t\t<WavelengthDataBlock>");
182 greg 2.5 printf("\t\t\t<WavelengthDataDirection>%s</WavelengthDataDirection>\n",
183     bsdf_type[(input_orient>0)<<1 | (output_orient>0)]);
184 greg 2.6 puts("\t\t\t<AngleBasis>LBNL/Shirley-Chiu</AngleBasis>");
185     puts("\t\t\t<ScatteringDataType>BTDF</ScatteringDataType>");
186 greg 2.5 puts("\t\t\t<ScatteringData>");
187     }
188    
189     /* Output XML epilogue to stdout */
190     static void
191     xml_epilogue(void)
192     {
193 greg 2.6 puts("\t\t\t</ScatteringData>");
194     puts("\t\t</WavelengthDataBlock>");
195     puts("\t</WavelengthData>");
196     puts("</Layer>");
197     puts("</Optical>");
198     puts("</WindowElement>");
199 greg 2.5 }
200    
201 greg 2.1 /* Read in BSDF and interpolate as tensor tree representation */
202     int
203     main(int argc, char *argv[])
204     {
205     FILE *fpin = stdin;
206     int i;
207    
208 greg 2.5 progname = argv[0];
209     for (i = 1; i < argc-1 && argv[i][0] == '-'; i++)
210     switch (argv[i][1]) { /* get option */
211 greg 2.1 case 't':
212 greg 2.5 pctcull = atoi(argv[++i]);
213 greg 2.1 break;
214     case 'g':
215 greg 2.5 samp_order = atoi(argv[++i]);
216 greg 2.1 break;
217     default:
218     goto userr;
219     }
220 greg 2.5
221     if (i == argc-1) { /* open input if given */
222     fpin = fopen(argv[i], "r");
223 greg 2.2 if (fpin == NULL) {
224     fprintf(stderr, "%s: cannot open BSDF interpolant '%s'\n",
225     progname, argv[1]);
226     return(1);
227     }
228 greg 2.5 } else if (i < argc-1)
229 greg 2.1 goto userr;
230     SET_FILE_BINARY(fpin); /* load BSDF interpolant */
231     if (!load_bsdf_rep(fpin))
232     return(1);
233 greg 2.5 fclose(fpin);
234     xml_prologue(argc, argv); /* start XML output */
235 greg 2.1 if (single_plane_incident) /* resample dist. */
236     interp_isotropic();
237     else
238     interp_anisotropic();
239 greg 2.5 xml_epilogue(); /* finish XML output */
240 greg 2.1 return(0);
241     userr:
242     fprintf(stderr,
243     "Usage: %s [-t pctcull][-g log2grid] [bsdf.sir] > bsdf.xml\n",
244     progname);
245     return(1);
246     }