src/cv/pkgBSDF.c

#ifndef lint
static const char RCSid[] = "$Id: pkgBSDF.c,v 2.8 2018/12/01 19:45:43 greg Exp $";
#endif
/*
 * Take BSDF XML file and generate a referencing Radiance object
 */

#include "rtio.h"
#include "paths.h"
#include "bsdf.h"

int     do_stdout = 0;                  /* send Radiance object to stdout? */
int     do_instance = 0;                /* produce instance/octree pairs? */

char    *progname;                      /* global argv[0] */

/* Return appropriate suffix for Z offset */
static const char *
get_suffix(double zval)
{
        return((zval > FTINY) ? "_f" : (zval < -FTINY) ? "_b" : "");
}


/* Produce a single rectangle corresponding to the given BSDF dimensions */
static void
faceBSDF(const SDData *bsp, double zval)
{
        const char      *sfx = get_suffix(zval);
                                                        /* output material */
        printf("void BSDF m_%s%s\n", bsp->name, sfx);
        printf("6 %g %s.xml 0 1 0 .\n", zval*1.00004, bsp->name);
        printf("0\n0\n\n");
                                                        /* output surface */
        zval *= (zval < 0) + 0.00002;
        printf("m_%s%s polygon %s%s\n", bsp->name, sfx, bsp->name, sfx);
        printf("0\n0\n12\n");
        printf("\t%.6e\t%.6e\t%.6e\n",
                                -.5*bsp->dim[0], -.5*bsp->dim[1], zval);
        printf("\t%.6e\t%.6e\t%.6e\n",
                                .5*bsp->dim[0], -.5*bsp->dim[1], zval);
        printf("\t%.6e\t%.6e\t%.6e\n",
                                .5*bsp->dim[0], .5*bsp->dim[1], zval);
        printf("\t%.6e\t%.6e\t%.6e\n\n",
                                -.5*bsp->dim[0], .5*bsp->dim[1], zval);
}


/* Convert BSDF (MGF) geometry to Radiance description */
static int
geomBSDF(const SDData *bsp)
{
        char    tmpfile[64];
        char    command[SDnameLn+64];
        int     fd;
                                        /* write MGF to temp file */
        fd = open(mktemp(strcpy(tmpfile,TEMPLATE)), O_WRONLY|O_CREAT|O_EXCL, 0600);
        if (fd < 0) {
                fprintf(stderr, "Cannot open temp file '%s'\n", tmpfile);
                return(0);
        }
        (void)write(fd, bsp->mgf, strlen(bsp->mgf));
        close(fd);
                                        /* set up command */
        if (do_instance) {
                sprintf(command, "mgf2rad -s %s | oconv -f - > %s.oct",
                                        tmpfile, bsp->name);
        } else {
                fflush(stdout);
                sprintf(command, "mgf2rad -s %s", tmpfile);
        }
        if (system(command)) {
                fprintf(stderr, "Error running: %s\n", command);
                return(0);
        }
        unlink(tmpfile);                /* remove temp file */
        if (do_instance) {              /* need to output instance? */
                printf("void instance %s_g\n", bsp->name);
                printf("1 %s.oct\n0\n0\n\n", bsp->name);
        }
        return(1);
}


/* Load a BSDF XML file and produce a corresponding Radiance object */
static int
cvtBSDF(char *fname)
{
        int     retOK;
        SDData  myBSDF;
        char    *pname, *fnbeg;
                                        /* find and load the XML file */
        retOK = strlen(fname);
        if (retOK < 5 || strcasecmp(fname+retOK-4, ".xml")) {
                fprintf(stderr, "%s: input does not end in '.xml'\n", fname);
                return(0);
        }
        pname = getpath(fname, getrlibpath(), R_OK);
        if (pname == NULL) {
                fprintf(stderr, "%s: cannot find BSDF file\n", fname);
                return(0);
        }
        fnbeg = strrchr(fname, DIRSEP);
        if (fnbeg != NULL)              /* eliminate directory */
                fname = fnbeg+1;
        SDclearBSDF(&myBSDF, fname);
        if (SDreportError(SDloadFile(&myBSDF, pname), stderr))
                return(0);
        retOK = (myBSDF.dim[0] > FTINY) & (myBSDF.dim[1] > FTINY);
        if (!retOK) {
                fprintf(stderr, "%s: zero width or height\n", fname);
        } else {
                if (!do_stdout) {
                        char    rname[SDnameLn+4];
                        strcpy(rname, myBSDF.name);
                        strcat(rname, ".rad");
                        retOK = (freopen(rname, "w", stdout) != NULL);
                }
                if (retOK) {
                        printf("# Produced by: %s %s\n", progname, fname);
                        if (myBSDF.matn[0] && myBSDF.makr[0])
                                printf("# Material '%s' by '%s'\n\n",
                                                myBSDF.matn, myBSDF.makr);
                        if (myBSDF.mgf == NULL) {
                                faceBSDF(&myBSDF, .0);
                        } else {
                                faceBSDF(&myBSDF, myBSDF.dim[2]);
                                if (myBSDF.rb != NULL)
                                        faceBSDF(&myBSDF, -myBSDF.dim[2]);
                                retOK = geomBSDF(&myBSDF);
                        }
                }
        }
        SDfreeBSDF(&myBSDF);            /* clean up and return */
        return(retOK);
}


/* Generate BSDF-referencing Radiance scenes, one per XML input */
int
main(int argc, char *argv[])
{
        int     status = 0;
        int     i;

        progname = argv[0];

        for (i = 1; i < argc && argv[i][0] == '-'; i++)
                switch (argv[i][1]) {
                case 'i':
                        do_instance = 1;
                        break;
                case 's':
                        do_stdout = 1;
                        break;
                default:
                        goto userr;
                }
        if (i >= argc) {
                fprintf(stderr, "%s: missing XML input\n", argv[0]);
                goto userr;
        }
        if (i < argc-1 && do_stdout) {
                fprintf(stderr, "%s: cannot send multiple BSDFs to stdout\n",
                                argv[0]);
                return(1);
        }
        for ( ; i < argc; i++)
                if (!cvtBSDF(argv[i])) {
                        fprintf(stderr, "%s: conversion of '%s' failed\n",
                                        argv[0], argv[i]);
                        status = 1;
                }
        return(status);
userr:
        fprintf(stderr, "Usage: %s [-i][-s] input.xml ..\n", argv[0]);
        return(1);
}
Revision:	2.9
Committed:	Sun Dec 8 18:33:23 2024 UTC (4 months, 3 weeks ago) by greg
Content type:	text/plain
Branch:	MAIN
CVS Tags:	HEAD
Changes since 2.8:	+3 -3 lines
Log Message:	feat(pkgBSDF): Added spectral conversion for MGF materials
#	User	Rev	Content
1	greg	2.1	#ifndef lint
2	greg	2.9	static const char RCSid[] = "$Id: pkgBSDF.c,v 2.8 2018/12/01 19:45:43 greg Exp $";
3	greg	2.1	#endif
4			/*
5			* Take BSDF XML file and generate a referencing Radiance object
6			*/
7
8	greg	2.2	#include "rtio.h"
9	greg	2.1	#include "paths.h"
10			#include "bsdf.h"
11
12			int do_stdout = 0; /* send Radiance object to stdout? */
13			int do_instance = 0; /* produce instance/octree pairs? */
14
15	greg	2.8	char progname; / global argv[0] */
16	greg	2.1
17			/* Return appropriate suffix for Z offset */
18			static const char *
19			get_suffix(double zval)
20			{
21			return((zval > FTINY) ? "_f" : (zval < -FTINY) ? "_b" : "");
22			}
23
24
25			/* Produce a single rectangle corresponding to the given BSDF dimensions */
26			static void
27			faceBSDF(const SDData *bsp, double zval)
28			{
29			const char *sfx = get_suffix(zval);
30			/* output material */
31			printf("void BSDF m_%s%s\n", bsp->name, sfx);
32			printf("6 %g %s.xml 0 1 0 .\n", zval*1.00004, bsp->name);
33			printf("0\n0\n\n");
34			/* output surface */
35			zval *= (zval < 0) + 0.00002;
36			printf("m_%s%s polygon %s%s\n", bsp->name, sfx, bsp->name, sfx);
37			printf("0\n0\n12\n");
38			printf("\t%.6e\t%.6e\t%.6e\n",
39			-.5bsp->dim[0], -.5bsp->dim[1], zval);
40			printf("\t%.6e\t%.6e\t%.6e\n",
41			.5bsp->dim[0], -.5bsp->dim[1], zval);
42			printf("\t%.6e\t%.6e\t%.6e\n",
43			.5bsp->dim[0], .5bsp->dim[1], zval);
44			printf("\t%.6e\t%.6e\t%.6e\n\n",
45			-.5bsp->dim[0], .5bsp->dim[1], zval);
46			}
47
48
49			/* Convert BSDF (MGF) geometry to Radiance description */
50			static int
51			geomBSDF(const SDData *bsp)
52			{
53			char tmpfile[64];
54			char command[SDnameLn+64];
55			int fd;
56			/* write MGF to temp file */
57	greg	2.3	fd = open(mktemp(strcpy(tmpfile,TEMPLATE)), O_WRONLY\|O_CREAT\|O_EXCL, 0600);
58	greg	2.1	if (fd < 0) {
59			fprintf(stderr, "Cannot open temp file '%s'\n", tmpfile);
60			return(0);
61			}
62	greg	2.3	(void)write(fd, bsp->mgf, strlen(bsp->mgf));
63	greg	2.1	close(fd);
64			/* set up command */
65			if (do_instance) {
66	greg	2.9	sprintf(command, "mgf2rad -s %s \| oconv -f - > %s.oct",
67	greg	2.1	tmpfile, bsp->name);
68			} else {
69			fflush(stdout);
70	greg	2.9	sprintf(command, "mgf2rad -s %s", tmpfile);
71	greg	2.1	}
72			if (system(command)) {
73			fprintf(stderr, "Error running: %s\n", command);
74			return(0);
75			}
76			unlink(tmpfile); /* remove temp file */
77			if (do_instance) { /* need to output instance? */
78			printf("void instance %s_g\n", bsp->name);
79			printf("1 %s.oct\n0\n0\n\n", bsp->name);
80			}
81			return(1);
82			}
83
84
85			/* Load a BSDF XML file and produce a corresponding Radiance object */
86			static int
87	greg	2.2	cvtBSDF(char *fname)
88	greg	2.1	{
89			int retOK;
90			SDData myBSDF;
91	greg	2.6	char pname, fnbeg;
92	greg	2.2	/* find and load the XML file */
93	greg	2.1	retOK = strlen(fname);
94	greg	2.7	if (retOK < 5 \|\| strcasecmp(fname+retOK-4, ".xml")) {
95	greg	2.1	fprintf(stderr, "%s: input does not end in '.xml'\n", fname);
96			return(0);
97			}
98	greg	2.2	pname = getpath(fname, getrlibpath(), R_OK);
99			if (pname == NULL) {
100			fprintf(stderr, "%s: cannot find BSDF file\n", fname);
101			return(0);
102			}
103	greg	2.6	fnbeg = strrchr(fname, DIRSEP);
104			if (fnbeg != NULL) /* eliminate directory */
105			fname = fnbeg+1;
106	greg	2.1	SDclearBSDF(&myBSDF, fname);
107	greg	2.5	if (SDreportError(SDloadFile(&myBSDF, pname), stderr))
108	greg	2.1	return(0);
109			retOK = (myBSDF.dim[0] > FTINY) & (myBSDF.dim[1] > FTINY);
110			if (!retOK) {
111			fprintf(stderr, "%s: zero width or height\n", fname);
112			} else {
113			if (!do_stdout) {
114			char rname[SDnameLn+4];
115			strcpy(rname, myBSDF.name);
116			strcat(rname, ".rad");
117			retOK = (freopen(rname, "w", stdout) != NULL);
118			}
119			if (retOK) {
120	greg	2.8	printf("# Produced by: %s %s\n", progname, fname);
121	greg	2.4	if (myBSDF.matn[0] && myBSDF.makr[0])
122			printf("# Material '%s' by '%s'\n\n",
123			myBSDF.matn, myBSDF.makr);
124	greg	2.1	if (myBSDF.mgf == NULL) {
125			faceBSDF(&myBSDF, .0);
126			} else {
127			faceBSDF(&myBSDF, myBSDF.dim[2]);
128			if (myBSDF.rb != NULL)
129			faceBSDF(&myBSDF, -myBSDF.dim[2]);
130			retOK = geomBSDF(&myBSDF);
131			}
132			}
133			}
134			SDfreeBSDF(&myBSDF); /* clean up and return */
135			return(retOK);
136			}
137
138
139			/* Generate BSDF-referencing Radiance scenes, one per XML input */
140			int
141			main(int argc, char *argv[])
142			{
143			int status = 0;
144			int i;
145
146	greg	2.8	progname = argv[0];
147
148	greg	2.1	for (i = 1; i < argc && argv[i][0] == '-'; i++)
149			switch (argv[i][1]) {
150			case 'i':
151			do_instance = 1;
152			break;
153			case 's':
154			do_stdout = 1;
155			break;
156			default:
157			goto userr;
158			}
159			if (i >= argc) {
160			fprintf(stderr, "%s: missing XML input\n", argv[0]);
161			goto userr;
162			}
163			if (i < argc-1 && do_stdout) {
164			fprintf(stderr, "%s: cannot send multiple BSDFs to stdout\n",
165			argv[0]);
166			return(1);
167			}
168			for ( ; i < argc; i++)
169			if (!cvtBSDF(argv[i])) {
170			fprintf(stderr, "%s: conversion of '%s' failed\n",
171			argv[0], argv[i]);
172			status = 1;
173			}
174			return(status);
175			userr:
176			fprintf(stderr, "Usage: %s [-i][-s] input.xml ..\n", argv[0]);
177			return(1);
178			}