int
viewfile( /* get view from file */
char *fname,
VIEW *vp,
RESOLU *rp
)
{
struct myview mvs;
FILE *fp;
if (fname == NULL || !strcmp(fname, "-"))
fp = stdin;
else if ((fp = fopen(fname, "r")) == NULL)
return(-1);
mvs.hv = vp;
mvs.ok = 0;
getheader(fp, gethview, &mvs);
if (rp != NULL && !fgetsresolu(rp, fp))
mvs.ok = 0;
fclose(fp);
return(mvs.ok);
}
static PICTURE *
openpicture( /* open/check Radiance picture file */
char *fname
)
{
FILE *fp;
register PICTURE *pp;
register char *cp;
/* check filename suffix */
if (fname == NULL) return(NULL);
for (cp = fname; *cp; cp++)
;
while (cp > fname && cp[-1] != '.')
if (*--cp == '/') {
cp = fname;
break;
}
if (cp > fname && !strncmp(cp, "tif", 3))
return(NULL); /* assume it's a TIFF */
/* else try opening it */
if ((fp = fopen(fname, "r")) == NULL)
return(NULL);
SET_FILE_BINARY(fp);
/* allocate struct */
if ((pp = (PICTURE *)malloc(sizeof(PICTURE))) == NULL)
return(NULL); /* serious error -- should exit? */
pp->fp = fp; pp->fmt[0] = '\0'; pp->pa = 1.;
/* load header */
if (getheader(fp, headline, pp) < 0) {
closepicture(pp);
return(NULL);
}
if (!pp->fmt[0]) /* assume RGBE if unspecified */
strcpy(pp->fmt, COLRFMT);
if (!globmatch(PICFMT, pp->fmt) || !fgetsresolu(&pp->rs, fp)) {
closepicture(pp); /* failed test -- close file */
return(NULL);
}
rewind(fp); /* passed test -- rewind file */
return(pp);
}
static void
checkfile(void) /* ready a file */
{
register int i;
/* process header */
gotview = 0;
if (echoheader) {
fputs(input[nfiles].name, stdout);
fputs(":\n", stdout);
}
getheader(input[nfiles].fp, headline, NULL);
if (wrongformat < 0) {
eputs(input[nfiles].name);
eputs(": not a Radiance picture\n");
quit(1);
}
if (wrongformat > 0) {
wputs(input[nfiles].name);
wputs(": warning -- incompatible picture format\n");
}
if (!gotview || setview(&input[nfiles].vw) != NULL)
input[nfiles].vw.type = 0;
if (!fgetsresolu(&input[nfiles].rs, input[nfiles].fp)) {
eputs(input[nfiles].name);
eputs(": bad picture size\n");
quit(1);
}
if (xmax == 0 && ymax == 0) {
xmax = scanlen(&input[nfiles].rs);
ymax = numscans(&input[nfiles].rs);
} else if (scanlen(&input[nfiles].rs) != xmax ||
numscans(&input[nfiles].rs) != ymax) {
eputs(input[nfiles].name);
eputs(": resolution mismatch\n");
quit(1);
}
/* allocate scanlines */
for (i = 0; i < WINSIZ; i++)
input[nfiles].scan[i] = (COLOR *)emalloc(xmax*sizeof(COLOR));
}
static void
getahead(void) /* load picture header */
{
char *err;
getheader(infp, headline, NULL);
if (lumf == NULL || !fgetsresolu(&inpres, infp)) {
fprintf(stderr, "%s: %s: not a Radiance picture\n",
progname, infn);
exit(1);
}
if (lumf == rgblum)
comprgb2xyzWBmat(inrgb2xyz, inprims);
else if (mbcalfile != NULL) {
fprintf(stderr, "%s: macbethcal only works with RGB pictures\n",
progname);
exit(1);
}
if (!gotview || ourview.type == VT_PAR ||
(ourview.horiz <= 5.) | (ourview.vert <= 5.)) {
ourview = stdview;
ourview.type = VT_PER;
if (pixaspect*inpres.yr < inpres.xr) {
ourview.horiz = 40.0;
ourview.vert = 2.*180./PI *
atan(.3639702*pixaspect*inpres.yr/inpres.xr);
} else {
ourview.vert = 40.0;
ourview.horiz = 2.*180./PI *
atan(.3639702*inpres.xr/pixaspect/inpres.yr);
}
}
if ((err = setview(&ourview)) != NULL) {
fprintf(stderr, "%s: view error in picture \"%s\": %s\n",
progname, infn, err);
exit(1);
}
}
int
main(
int argc,
char **argv
)
{
double d, expval = 1.0;
int i;
progname = argv[0];
mybright = &rgb_bright; /* default */
for (i = 1; i < argc; i++)
if (argv[i][0] == '-' || argv[i][0] == '+')
switch (argv[i][1]) {
case 'h': /* header */
header = argv[i][0] == '+';
break;
case 'H': /* resolution string */
resolution = argv[i][0] == '+';
break;
case 's': /* skip bytes in header */
skipbytes = atol(argv[++i]);
break;
case 'u': /* unique values */
uniq = argv[i][0] == '-';
break;
case 'o': /* original values */
original = argv[i][0] == '-';
break;
case 'g': /* gamma correction */
gamcor = atof(argv[i+1]);
if (argv[i][0] == '+')
gamcor = 1.0/gamcor;
i++;
break;
case 'e': /* exposure correction */
d = atof(argv[i+1]);
if (argv[i+1][0] == '-' || argv[i+1][0] == '+')
d = pow(2.0, d);
if (argv[i][0] == '-')
expval *= d;
scalecolor(exposure, d);
doexposure++;
i++;
break;
case 'R': /* reverse byte sequence */
if (argv[i][0] == '-') {
ord[0]=BLU; ord[1]=GRN; ord[2]=RED;
} else {
ord[0]=RED; ord[1]=GRN; ord[2]=BLU;
}
break;
case 'r': /* reverse conversion */
reverse = argv[i][0] == '-';
break;
case 'n': /* non-interleaved RGB */
interleave = argv[i][0] == '+';
break;
case 'b': /* brightness values */
putprim = argv[i][0] == '-' ? BRIGHT : ALL;
break;
case 'p': /* primary controls */
switch (argv[i][2]) {
/* these two options affect -r conversion */
case '\0':
myprims[RED][CIEX] = atof(argv[++i]);
myprims[RED][CIEY] = atof(argv[++i]);
myprims[GRN][CIEX] = atof(argv[++i]);
myprims[GRN][CIEY] = atof(argv[++i]);
myprims[BLU][CIEX] = atof(argv[++i]);
myprims[BLU][CIEY] = atof(argv[++i]);
myprims[WHT][CIEX] = atof(argv[++i]);
myprims[WHT][CIEY] = atof(argv[++i]);
outprims = myprims;
break;
case 'x': case 'X': outprims = NULL; break;
/* the following options affect +r only */
case 'r': case 'R': putprim = RED; break;
case 'g': case 'G': putprim = GRN; break;
case 'b': case 'B': putprim = BLU; break;
default: goto unkopt;
}
break;
case 'd': /* data only (no indices) */
dataonly = argv[i][0] == '-';
switch (argv[i][2]) {
case '\0':
case 'a': /* ascii */
format = 'a';
fmtid = "ascii";
break;
case 'i': /* integer */
format = 'i';
fmtid = "ascii";
break;
case 'b': /* byte */
dataonly = 1;
format = 'b';
fmtid = "byte";
break;
case 'W': /* 16-bit swapped */
swapbytes = 1;
case 'w': /* 16-bit */
dataonly = 1;
format = 'w';
fmtid = "16-bit";
break;
case 'F': /* swapped floats */
swapbytes = 1;
case 'f': /* float */
dataonly = 1;
format = 'f';
fmtid = "float";
break;
case 'D': /* swapped doubles */
swapbytes = 1;
case 'd': /* double */
dataonly = 1;
format = 'd';
fmtid = "double";
break;
default:
goto unkopt;
}
break;
case 'x': /* x resolution */
case 'X': /* x resolution */
resolution = 0;
if (argv[i][0] == '-')
picres.rt |= XDECR;
picres.xr = atoi(argv[++i]);
break;
case 'y': /* y resolution */
case 'Y': /* y resolution */
resolution = 0;
if (argv[i][0] == '-')
picres.rt |= YDECR;
if (picres.xr == 0)
picres.rt |= YMAJOR;
picres.yr = atoi(argv[++i]);
break;
default:
unkopt:
fprintf(stderr, "%s: unknown option: %s\n",
progname, argv[i]);
quit(1);
break;
}
else
break;
/* recognize special formats */
if (dataonly && format == 'b') {
if (putprim == ALL)
fmtid = "24-bit_rgb";
else
fmtid = "8-bit_grey";
}
if (dataonly && format == 'w') {
if (putprim == ALL)
fmtid = "48-bit_rgb";
else
fmtid = "16-bit_grey";
}
/* assign reverse ordering */
rord[ord[0]] = 0;
rord[ord[1]] = 1;
rord[ord[2]] = 2;
/* get input */
if (i == argc) {
fin = stdin;
} else if (i < argc) {
if ((fin = fopen(argv[i], "r")) == NULL) {
fprintf(stderr, "%s: can't open file \"%s\"\n",
progname, argv[i]);
quit(1);
}
if (reverse && putprim != BRIGHT && i == argc-3) {
if ((fin2 = fopen(argv[i+1], "r")) == NULL) {
fprintf(stderr, "%s: can't open file \"%s\"\n",
progname, argv[i+1]);
quit(1);
}
if ((fin3 = fopen(argv[i+2], "r")) == NULL) {
fprintf(stderr, "%s: can't open file \"%s\"\n",
progname, argv[i+2]);
quit(1);
}
interleave = -1;
} else if (i != argc-1)
fin = NULL;
if (reverse && putprim != BRIGHT && !interleave) {
fin2 = fopen(argv[i], "r");
fin3 = fopen(argv[i], "r");
}
if (skipbytes && (fseek(fin, skipbytes, 0) || (fin2 != NULL &&
(fseek(fin2, skipbytes, 0) ||
fseek(fin3, skipbytes, 0))))) {
fprintf(stderr, "%s: cannot skip %ld bytes on input\n",
progname, skipbytes);
quit(1);
}
}
if (fin == NULL) {
fprintf(stderr, "%s: bad # file arguments\n", progname);
quit(1);
}
if (reverse) {
#ifdef _WIN32
SET_FILE_BINARY(stdout);
if (format != 'a' && format != 'i')
SET_FILE_BINARY(fin);
#endif
/* get header */
if (header) {
if (checkheader(fin, fmtid, stdout) < 0) {
fprintf(stderr, "%s: wrong input format\n",
progname);
quit(1);
}
if (fin2 != NULL) {
getheader(fin2, NULL, NULL);
getheader(fin3, NULL, NULL);
}
} else
newheader("RADIANCE", stdout);
/* get resolution */
if ((resolution && !fgetsresolu(&picres, fin)) ||
picres.xr <= 0 || picres.yr <= 0) {
fprintf(stderr, "%s: missing resolution\n", progname);
quit(1);
}
if (resolution && fin2 != NULL) {
RESOLU pres2;
if (!fgetsresolu(&pres2, fin2) ||
pres2.rt != picres.rt ||
pres2.xr != picres.xr ||
pres2.yr != picres.yr ||
!fgetsresolu(&pres2, fin3) ||
pres2.rt != picres.rt ||
pres2.xr != picres.xr ||
pres2.yr != picres.yr) {
fprintf(stderr, "%s: resolution mismatch\n",
progname);
quit(1);
}
}
/* add to header */
printargs(i, argv, stdout);
if (expval < .99 || expval > 1.01)
fputexpos(expval, stdout);
if (outprims != NULL) {
if (outprims != stdprims)
fputprims(outprims, stdout);
fputformat(COLRFMT, stdout);
} else /* XYZ data */
fputformat(CIEFMT, stdout);
putchar('\n');
fputsresolu(&picres, stdout); /* always put resolution */
valtopix();
} else {
#ifdef _WIN32
SET_FILE_BINARY(fin);
if (format != 'a' && format != 'i')
SET_FILE_BINARY(stdout);
#endif
/* get header */
getheader(fin, checkhead, NULL);
if (wrongformat) {
fprintf(stderr,
"%s: input not a Radiance RGBE picture\n",
progname);
quit(1);
}
if (!fgetsresolu(&picres, fin)) {
fprintf(stderr, "%s: missing resolution\n", progname);
quit(1);
}
if (header) {
printargs(i, argv, stdout);
if (expval < .99 || expval > 1.01)
fputexpos(expval, stdout);
fputformat(fmtid, stdout);
putchar('\n');
}
if (resolution) /* put resolution */
fputsresolu(&picres, stdout);
pixtoval();
}
quit(0);
return 0; /* pro forma return */
}
void
addpicz( /* add a picture + depth-buffer */
char *pcf,
char *zbf
)
{
FILE *pfp;
int zfd;
COLR *cscn;
float *zscn;
struct phead phd;
int eshft;
double emult;
RESOLU prs;
RREAL vl[2];
FVECT ro, rd;
double aftd;
COLOR ctmp;
int j;
register int i;
/* open files */
if ((pfp = fopen(pcf, "r")) == NULL) {
sprintf(errmsg, "cannot open picture file \"%s\"", pcf);
error(SYSTEM, pcf);
}
if ((zfd = open(zbf, O_RDONLY)) < 0) {
sprintf(errmsg, "cannot open depth file \"%s\"", zbf);
error(SYSTEM, pcf);
}
/* load picture header */
phd.vw = stdview;
phd.expos = 1.0;
phd.badfmt = phd.gotview = phd.altprims = 0;
if (getheader(pfp, picheadline, &phd) < 0 ||
phd.badfmt || !fgetsresolu(&prs, pfp)) {
sprintf(errmsg, "bad format for picture file \"%s\"", pcf);
error(USER, errmsg);
}
if (!phd.gotview || setview(&phd.vw) != NULL) {
sprintf(errmsg, "missing/illegal view in picture \"%s\"",
pcf);
error(USER, errmsg);
}
if (phd.altprims) {
sprintf(errmsg, "ignoring primary values in picture \"%s\"",
pcf);
error(WARNING, errmsg);
}
/* figure out what to do about exposure */
if ((phd.expos < 0.99) | (phd.expos > 1.01)) {
emult = -log(phd.expos)/log(2.);
eshft = emult >= 0. ? emult+.5 : emult-.5;
emult -= (double)eshft;
if ((emult <= 0.01) & (emult >= -0.01))
emult = -1.;
else {
emult = 1./phd.expos;
eshft = 0;
}
} else {
emult = -1.;
eshft = 0;
}
/* allocate buffers */
cscn = (COLR *)malloc(scanlen(&prs)*sizeof(COLR));
zscn = (float *)malloc(scanlen(&prs)*sizeof(float));
if ((cscn == NULL) | (zscn == NULL))
error(SYSTEM, "out of memory in addpicz");
/* read and process each scanline */
for (j = 0; j < numscans(&prs); j++) {
i = scanlen(&prs); /* read colrs */
if (freadcolrs(cscn, i, pfp) < 0) {
sprintf(errmsg, "error reading picture \"%s\"", pcf);
error(USER, errmsg);
}
if (eshft) /* shift exposure */
shiftcolrs(cscn, i, eshft);
i *= sizeof(float); /* read depth */
if (read(zfd, (char *)zscn, i) != i) {
sprintf(errmsg, "error reading depth file \"%s\"", zbf);
error(USER, errmsg);
}
for (i = scanlen(&prs); i--; ) { /* do each pixel */
pix2loc(vl, &prs, i, j);
aftd = viewray(ro, rd, &phd.vw, vl[0], vl[1]);
if (aftd < -FTINY)
continue; /* off view */
if (aftd > FTINY && zscn[i] > aftd)
continue; /* aft clipped */
if (emult > 0.) { /* whatta pain */
colr_color(ctmp, cscn[i]);
scalecolor(ctmp, emult);
setcolr(cscn[i], colval(ctmp,RED),
colval(ctmp,GRN), colval(ctmp,BLU));
}
addray(ro, rd, (double)zscn[i], cscn[i]);
}
}
/* write output and free beams */
hdflush(NULL);
/* clean up */
free((void *)cscn);
free((void *)zscn);
fclose(pfp);
close(zfd);
}
