src/px/ra_rgbe.c

#ifndef lint
static const char       RCSid[] = "$Id: ra_rgbe.c,v 2.22 2024/06/11 17:24:01 greg Exp $";
#endif
/*
 *  program to convert from RADIANCE RLE to flat format
 */

#include  <math.h>

#include  "platform.h"
#include  "rtio.h"
#include  "paths.h"
#include  "color.h"
#include  "resolu.h"

#define dumpheader(fp)  putbinary(headlines, 1, headlen, fp)

int  bradj = 0;                         /* brightness adjustment */
int  doflat = 1;                        /* produce flat file */
int  force = 0;                         /* force file overwrite? */
int  findframe = 0;                     /* find a specific frame? */
int  frameno = 0;                       /* current frame number */
int  fmterr = 0;                        /* got input format error */
char  *headlines = NULL;                /* current header info. */
int  headlen1 = 0;                      /* length of initial frame header */
int  headlen = 0;                       /* current header length */
char  fmt[MAXFMTLEN];                   /* input format */

char  *progname;

static gethfunc addhline;
static int transfer(char *ospec);
static int loadheader(FILE *fp);


int
main(int  argc, char  *argv[])
{
        char    *ospec;
        int  i;

        progname = argv[0];

        for (i = 1; i < argc; i++)
                if (argv[i][0] == '-')
                        switch (argv[i][1]) {
                        case 'r':
                                doflat = !doflat;
                                break;
                        case 'e':
                                if (argv[i+1][0] != '+' && argv[i+1][0] != '-')
                                        goto userr;
                                bradj = atoi(argv[++i]);
                                break;
                        case 'n':
                                findframe = atoi(argv[++i]);
                                break;
                        case 'f':
                                force++;
                                break;
                        case '\0':
                                goto gotfile;
                        default:
                                goto userr;
                        }
                else
                        break;
gotfile:
        if (i < argc-2)
                goto userr;
        if (i <= argc-1 && strcmp(argv[i], "-") &&
                        freopen(argv[i], "r", stdin) == NULL) {
                fprintf(stderr, "%s: can't open input \"%s\"\n",
                                progname, argv[i]);
                exit(1);
        }
        SET_FILE_BINARY(stdin);
        ospec = i==argc-2 ? argv[i+1] : (char *)NULL;
        while (transfer(ospec))
                ;
        exit(0);
userr:
        fprintf(stderr,
        "Usage: %s [-r][-e +/-stops][-f][-n frame] [input [outspec]]\n",
                        progname);
        exit(1);
}


static int
transfer(                       /* transfer a Radiance picture */
        char    *ospec
)
{
        char    oname[PATH_MAX];
        FILE    *fp;
        int     order;
        int     xmax, ymax;
        COLR    *scanin;
        int     y;
                                        /* get header info. */
        if (!(y = loadheader(stdin)))
                return(0);
        if (y < 0 || (order = fgetresolu(&xmax, &ymax, stdin)) < 0) {
                fprintf(stderr, "%s: bad input format\n", progname);
                exit(1);
        }
                                        /* did we pass the target frame? */
        if (findframe && findframe < frameno)
                return(0);
                                        /* allocate scanline */
        scanin = (COLR *)malloc(xmax*sizeof(COLR));
        if (scanin == NULL) {
                perror(progname);
                exit(1);
        }
                                        /* skip frame? */
        if (findframe > frameno) {
                if (NCSAMP > 3) {
                        if (fseek(stdin, ymax*xmax*LSCOLR, SEEK_CUR) < 0) {
                                perror(progname);
                                exit(1);
                        }
                } else
                    for (y = ymax; y--; )
                        if (freadcolrs(scanin, xmax, stdin) < 0) {
                                fprintf(stderr,
                                        "%s: error reading input picture\n",
                                                progname);
                                exit(1);
                        }
                free(scanin);
                return(1);
        }
                                        /* open output file/command */
        if (ospec == NULL) {
                strcpy(oname, "<stdout>");
                fp = stdout;
        } else {
                sprintf(oname, ospec, frameno);
                if (oname[0] == '!') {
                        if ((fp = popen(oname+1, "w")) == NULL) {
                                fprintf(stderr, "%s: cannot start \"%s\"\n",
                                                progname, oname);
                                exit(1);
                        }
                } else {
                        if (!force && access(oname, 0) >= 0) {
                                fprintf(stderr,
                                        "%s: output file \"%s\" exists\n",
                                                progname, oname);
                                exit(1);
                        }
                        if ((fp = fopen(oname, "w")) == NULL) {
                                fprintf(stderr, "%s: ", progname);
                                perror(oname);
                                exit(1);
                        }
                }
        }
        SET_FILE_BINARY(fp);
        newheader("RADIANCE", fp);              /* put out header */
        dumpheader(fp);
        fputs(progname, fp);
        if (bradj)
                fprintf(fp, " -e %+d", bradj);
        if (!doflat)
                fputs(" -r", fp);
        fputc('\n', fp);
        if (bradj)
                fputexpos(pow(2.0, (double)bradj), fp);
        if (frameno)
                fprintf(fp, "FRAME=%d\n", frameno);
        if (fmt[0])
                fputformat(fmt, fp);
        fputc('\n', fp);
        fputresolu(order, xmax, ymax, fp);
                                        /* transfer picture */
        for (y = ymax; y--; ) {
                if (fread2colrs(scanin, xmax, stdin, NCSAMP, WLPART) < 0) {
                        fprintf(stderr, "%s: error reading input picture\n",
                                        progname);
                        exit(1);
                }
                if (bradj)
                        shiftcolrs(scanin, xmax, bradj);
                if (doflat ? (putbinary(scanin, sizeof(COLR), xmax, fp) != xmax) :
                                (fwritecolrs(scanin, xmax, fp) < 0))
                        goto writerr;
        }
        free(scanin);                   /* clean up */
        if (fflush(fp) == EOF)
                goto writerr;
        if (oname[0] == '!')
                pclose(fp);
        else if (ospec != NULL)
                fclose(fp);
        return(1);
writerr:
        fprintf(stderr, "%s: error writing output to \"%s\"\n",
                        progname, oname);
        exit(1);
}


static int
addhline(                       /* add a line to our info. header */
        char    *s,
        void    *p
)
{
        int     n;

        if (isheadid(s))
                return(0);
        if (!strncmp(s, "FRAME=", 6)) {
                frameno = atoi(s+6);
                return(0);
        }
        if (formatval(fmt, s)) {
                if (!strcmp(fmt, SPECFMT))
                        strcpy(fmt, COLRFMT);
                else
                        fmterr += !globmatch(PICFMT, fmt);
                return(0);
        }
        if (isncomp(s)) {
                NCSAMP = ncompval(s);
                return(NCSAMP - 3);
        }
        if (iswlsplit(s)) {
                wlsplitval(WLPART, s);
                return(0);
        }
        n = strlen(s);
        if (headlen)
                headlines = (char *)realloc((void *)headlines, headlen+n+1);
        else
                headlines = (char *)malloc(n+1);
        if (headlines == NULL) {
                perror(progname);
                exit(1);
        }
        strcpy(headlines+headlen, s);
        headlen += n;
        return(0);
}


static int
loadheader(                     /* load an info. header into memory */
        FILE    *fp
)
{
        fmterr = 0; frameno = 0;
                                /* revert to initial header length */
        if (!headlen1) headlen1 = headlen;
        else headlen = headlen1;

        if (getheader(fp, addhline, NULL) < 0)
                return(0);
        if (fmterr)
                return(-1);
        return(1);
}
Revision:	2.23
Committed:	Tue Sep 10 20:24:42 2024 UTC (7 months, 3 weeks ago) by greg
Content type:	text/plain
Branch:	MAIN
Changes since 2.22:	+25 -6 lines
Log Message:	feat(ra_bmp,ra_rgbe,ximage): Added quick conversion of spectral inputs to RGB
#	User	Rev	Content
1	greg	1.1	#ifndef lint
2	greg	2.23	static const char RCSid[] = "$Id: ra_rgbe.c,v 2.22 2024/06/11 17:24:01 greg Exp $";
3	greg	1.1	#endif
4			/*
5			* program to convert from RADIANCE RLE to flat format
6			*/
7
8	greg	2.4	#include <math.h>
9	schorsch	2.12
10			#include "platform.h"
11	greg	2.19	#include "rtio.h"
12	greg	2.18	#include "paths.h"
13	greg	2.4	#include "color.h"
14			#include "resolu.h"
15
16	greg	2.19	#define dumpheader(fp) putbinary(headlines, 1, headlen, fp)
17	gregl	2.8
18	greg	1.1	int bradj = 0; /* brightness adjustment */
19			int doflat = 1; /* produce flat file */
20	gregl	2.8	int force = 0; /* force file overwrite? */
21			int findframe = 0; /* find a specific frame? */
22			int frameno = 0; /* current frame number */
23			int fmterr = 0; /* got input format error */
24	greg	2.22	char headlines = NULL; / current header info. */
25			int headlen1 = 0; /* length of initial frame header */
26			int headlen = 0; /* current header length */
27			char fmt[MAXFMTLEN]; /* input format */
28	gregl	2.8
29	greg	1.1	char *progname;
30
31	schorsch	2.16	static gethfunc addhline;
32	schorsch	2.17	static int transfer(char *ospec);
33			static int loadheader(FILE *fp);
34	schorsch	2.16
35	greg	1.1
36	schorsch	2.17	int
37			main(int argc, char *argv[])
38	greg	1.1	{
39	gregl	2.8	char *ospec;
40	greg	1.1	int i;
41	gregl	2.7
42	greg	1.1	progname = argv[0];
43
44			for (i = 1; i < argc; i++)
45			if (argv[i][0] == '-')
46			switch (argv[i][1]) {
47			case 'r':
48			doflat = !doflat;
49			break;
50			case 'e':
51			if (argv[i+1][0] != '+' && argv[i+1][0] != '-')
52			goto userr;
53			bradj = atoi(argv[++i]);
54			break;
55	gregl	2.8	case 'n':
56			findframe = atoi(argv[++i]);
57			break;
58			case 'f':
59			force++;
60			break;
61			case '\0':
62			goto gotfile;
63	greg	1.1	default:
64			goto userr;
65			}
66			else
67			break;
68	gregl	2.8	gotfile:
69	greg	1.1	if (i < argc-2)
70			goto userr;
71	gregl	2.8	if (i <= argc-1 && strcmp(argv[i], "-") &&
72			freopen(argv[i], "r", stdin) == NULL) {
73	greg	1.1	fprintf(stderr, "%s: can't open input \"%s\"\n",
74			progname, argv[i]);
75			exit(1);
76			}
77	schorsch	2.12	SET_FILE_BINARY(stdin);
78	gregl	2.8	ospec = i==argc-2 ? argv[i+1] : (char *)NULL;
79			while (transfer(ospec))
80			;
81	greg	1.1	exit(0);
82			userr:
83	gregl	2.8	fprintf(stderr,
84			"Usage: %s [-r][-e +/-stops][-f][-n frame] [input [outspec]]\n",
85	greg	1.1	progname);
86			exit(1);
87			}
88
89
90	schorsch	2.17	static int
91			transfer( /* transfer a Radiance picture */
92			char *ospec
93			)
94	greg	1.1	{
95	schorsch	2.15	char oname[PATH_MAX];
96	gregl	2.8	FILE *fp;
97	greg	1.2	int order;
98	greg	2.3	int xmax, ymax;
99	greg	1.1	COLR *scanin;
100			int y;
101	gregl	2.8	/* get header info. */
102			if (!(y = loadheader(stdin)))
103			return(0);
104			if (y < 0 \|\| (order = fgetresolu(&xmax, &ymax, stdin)) < 0) {
105			fprintf(stderr, "%s: bad input format\n", progname);
106			exit(1);
107			}
108			/* did we pass the target frame? */
109			if (findframe && findframe < frameno)
110			return(0);
111			/* allocate scanline */
112	greg	2.23	scanin = (COLR )malloc(xmaxsizeof(COLR));
113	gregl	2.8	if (scanin == NULL) {
114			perror(progname);
115			exit(1);
116			}
117			/* skip frame? */
118			if (findframe > frameno) {
119	greg	2.23	if (NCSAMP > 3) {
120			if (fseek(stdin, ymaxxmaxLSCOLR, SEEK_CUR) < 0) {
121			perror(progname);
122			exit(1);
123			}
124			} else
125			for (y = ymax; y--; )
126	gregl	2.8	if (freadcolrs(scanin, xmax, stdin) < 0) {
127			fprintf(stderr,
128			"%s: error reading input picture\n",
129			progname);
130			exit(1);
131			}
132	greg	2.23	free(scanin);
133	gregl	2.8	return(1);
134			}
135			/* open output file/command */
136			if (ospec == NULL) {
137			strcpy(oname, "<stdout>");
138			fp = stdout;
139			} else {
140			sprintf(oname, ospec, frameno);
141			if (oname[0] == '!') {
142			if ((fp = popen(oname+1, "w")) == NULL) {
143			fprintf(stderr, "%s: cannot start \"%s\"\n",
144			progname, oname);
145			exit(1);
146			}
147			} else {
148			if (!force && access(oname, 0) >= 0) {
149			fprintf(stderr,
150			"%s: output file \"%s\" exists\n",
151			progname, oname);
152			exit(1);
153			}
154			if ((fp = fopen(oname, "w")) == NULL) {
155			fprintf(stderr, "%s: ", progname);
156			perror(oname);
157			exit(1);
158			}
159			}
160			}
161	schorsch	2.12	SET_FILE_BINARY(fp);
162	greg	2.22	newheader("RADIANCE", fp); /* put out header */
163			dumpheader(fp);
164	gregl	2.8	fputs(progname, fp);
165	greg	1.1	if (bradj)
166	gregl	2.8	fprintf(fp, " -e %+d", bradj);
167	gregl	2.7	if (!doflat)
168	gregl	2.8	fputs(" -r", fp);
169			fputc('\n', fp);
170	greg	2.2	if (bradj)
171	gregl	2.8	fputexpos(pow(2.0, (double)bradj), fp);
172	greg	2.22	if (frameno)
173			fprintf(fp, "FRAME=%d\n", frameno);
174			if (fmt[0])
175			fputformat(fmt, fp);
176	gregl	2.8	fputc('\n', fp);
177			fputresolu(order, xmax, ymax, fp);
178			/* transfer picture */
179	gregl	2.7	for (y = ymax; y--; ) {
180	greg	2.23	if (fread2colrs(scanin, xmax, stdin, NCSAMP, WLPART) < 0) {
181	gregl	2.8	fprintf(stderr, "%s: error reading input picture\n",
182			progname);
183			exit(1);
184			}
185	greg	1.1	if (bradj)
186			shiftcolrs(scanin, xmax, bradj);
187	greg	2.22	if (doflat ? (putbinary(scanin, sizeof(COLR), xmax, fp) != xmax) :
188			(fwritecolrs(scanin, xmax, fp) < 0))
189			goto writerr;
190	greg	1.1	}
191	greg	2.23	free(scanin); /* clean up */
192			if (fflush(fp) == EOF)
193	greg	2.22	goto writerr;
194	gregl	2.8	if (oname[0] == '!')
195			pclose(fp);
196			else if (ospec != NULL)
197			fclose(fp);
198			return(1);
199	greg	2.22	writerr:
200			fprintf(stderr, "%s: error writing output to \"%s\"\n",
201			progname, oname);
202			exit(1);
203	gregl	2.8	}
204
205
206	schorsch	2.16	static int
207			addhline( /* add a line to our info. header */
208			char *s,
209			void *p
210			)
211	gregl	2.8	{
212			int n;
213
214	greg	2.22	if (isheadid(s))
215			return(0);
216			if (!strncmp(s, "FRAME=", 6)) {
217			frameno = atoi(s+6);
218			return(0);
219			}
220			if (formatval(fmt, s)) {
221	greg	2.23	if (!strcmp(fmt, SPECFMT))
222			strcpy(fmt, COLRFMT);
223			else
224			fmterr += !globmatch(PICFMT, fmt);
225			return(0);
226			}
227			if (isncomp(s)) {
228			NCSAMP = ncompval(s);
229			return(NCSAMP - 3);
230			}
231			if (iswlsplit(s)) {
232			wlsplitval(WLPART, s);
233	greg	2.22	return(0);
234			}
235	gregl	2.8	n = strlen(s);
236			if (headlen)
237	greg	2.11	headlines = (char )realloc((void )headlines, headlen+n+1);
238	gregl	2.8	else
239			headlines = (char *)malloc(n+1);
240			if (headlines == NULL) {
241			perror(progname);
242			exit(1);
243			}
244			strcpy(headlines+headlen, s);
245			headlen += n;
246	gwlarson	2.9	return(0);
247	gregl	2.8	}
248
249
250	schorsch	2.17	static int
251			loadheader( /* load an info. header into memory */
252			FILE *fp
253			)
254	gregl	2.8	{
255			fmterr = 0; frameno = 0;
256	greg	2.22	/* revert to initial header length */
257			if (!headlen1) headlen1 = headlen;
258			else headlen = headlen1;
259
260	gregl	2.8	if (getheader(fp, addhline, NULL) < 0)
261			return(0);
262			if (fmterr)
263			return(-1);
264			return(1);
265	greg	1.1	}