1 |
< |
/* Copyright (c) 1991 Regents of the University of California */ |
1 |
> |
/* Copyright (c) 1994 Regents of the University of California */ |
2 |
|
|
3 |
|
#ifndef lint |
4 |
|
static char SCCSid[] = "$SunId$ LBL"; |
20 |
|
|
21 |
|
#include "resolu.h" |
22 |
|
|
23 |
– |
#ifndef BSD |
24 |
– |
#define vfork fork |
25 |
– |
#endif |
26 |
– |
|
23 |
|
#define pscan(y) (ourpict+(y)*hresolu) |
24 |
|
#define zscan(y) (ourzbuf+(y)*hresolu) |
25 |
|
|
165 |
|
goto badopt; |
166 |
|
check(3,"s"); |
167 |
|
gotvfile = viewfile(argv[++i], &ourview, 0, 0); |
168 |
< |
if (gotvfile < 0) { |
169 |
< |
perror(argv[i]); |
170 |
< |
exit(1); |
175 |
< |
} else if (gotvfile == 0) { |
168 |
> |
if (gotvfile < 0) |
169 |
> |
syserror(argv[i]); |
170 |
> |
else if (gotvfile == 0) { |
171 |
|
fprintf(stderr, "%s: bad view file\n", |
172 |
|
argv[i]); |
173 |
|
exit(1); |
186 |
|
if (fillsamp == 1) |
187 |
|
fillo &= ~F_BACK; |
188 |
|
/* set view */ |
189 |
< |
if (err = setview(&ourview)) { |
189 |
> |
if ((err = setview(&ourview)) != NULL) { |
190 |
|
fprintf(stderr, "%s: %s\n", progname, err); |
191 |
|
exit(1); |
192 |
|
} |
195 |
|
ourpict = (COLR *)bmalloc(hresolu*vresolu*sizeof(COLR)); |
196 |
|
ourzbuf = (float *)bmalloc(hresolu*vresolu*sizeof(float)); |
197 |
|
if (ourpict == NULL || ourzbuf == NULL) |
198 |
< |
syserror(); |
198 |
> |
syserror(progname); |
199 |
|
bzero((char *)ourzbuf, hresolu*vresolu*sizeof(float)); |
200 |
+ |
/* new header */ |
201 |
+ |
newheader("RADIANCE", stdout); |
202 |
|
/* get input */ |
203 |
|
for ( ; i < argc; i += 2) |
204 |
|
addpicture(argv[i], argv[i+1]); |
209 |
|
fillpicture(fillfunc); |
210 |
|
/* close calculation */ |
211 |
|
caldone(); |
212 |
+ |
/* aft clipping */ |
213 |
+ |
clipaft(); |
214 |
|
/* add to header */ |
215 |
|
printargs(argc, argv, stdout); |
216 |
|
if (gotvfile) { |
245 |
|
{ |
246 |
|
char fmt[32]; |
247 |
|
|
248 |
< |
if (isformat(s)) { |
249 |
< |
formatval(fmt, s); |
248 |
> |
if (isheadid(s)) |
249 |
> |
return; |
250 |
> |
if (formatval(fmt, s)) { |
251 |
|
wrongformat = strcmp(fmt, COLRFMT); |
252 |
|
return; |
253 |
|
} |
274 |
|
struct position *plast; |
275 |
|
int y; |
276 |
|
/* open picture file */ |
277 |
< |
if ((pfp = fopen(pfile, "r")) == NULL) { |
278 |
< |
perror(pfile); |
279 |
< |
exit(1); |
280 |
< |
} |
277 |
> |
if ((pfp = fopen(pfile, "r")) == NULL) |
278 |
> |
syserror(pfile); |
279 |
|
/* get header with exposure and view */ |
280 |
|
theirexp = 1.0; |
281 |
|
gotview = 0; |
301 |
|
plast = (struct position *)calloc(scanlen(&tresolu), |
302 |
|
sizeof(struct position)); |
303 |
|
if (scanin == NULL || zin == NULL || plast == NULL) |
304 |
< |
syserror(); |
304 |
> |
syserror(progname); |
305 |
|
/* get z specification or file */ |
306 |
|
if ((zfd = open(zspec, O_RDONLY)) == -1) { |
307 |
|
double zvalue; |
308 |
|
register int x; |
309 |
< |
if (!isfloat(zspec) || (zvalue = atof(zspec)) <= 0.0) { |
310 |
< |
perror(zspec); |
313 |
< |
exit(1); |
314 |
< |
} |
309 |
> |
if (!isfloat(zspec) || (zvalue = atof(zspec)) <= 0.0) |
310 |
> |
syserror(zspec); |
311 |
|
for (x = scanlen(&tresolu); x-- > 0; ) |
312 |
|
zin[x] = zvalue; |
313 |
|
} |
382 |
|
float *zline; |
383 |
|
struct position *lasty; /* input/output */ |
384 |
|
{ |
389 |
– |
extern double sqrt(); |
385 |
|
FVECT pos; |
386 |
|
struct position lastx, newpos; |
387 |
|
register int x; |
470 |
|
movepixel(pos) /* reposition image point */ |
471 |
|
FVECT pos; |
472 |
|
{ |
473 |
+ |
double d0, d1; |
474 |
|
FVECT pt, direc; |
475 |
|
|
476 |
|
if (pos[2] <= 0) /* empty pixel */ |
477 |
|
return(-1); |
478 |
+ |
if (normdist && theirview.type == VT_PER) { /* adjust distance */ |
479 |
+ |
d0 = pos[0] + theirview.hoff - .5; |
480 |
+ |
d1 = pos[1] + theirview.voff - .5; |
481 |
+ |
pos[2] /= sqrt(1. + d0*d0*theirview.hn2 + d1*d1*theirview.vn2); |
482 |
+ |
} |
483 |
|
if (hasmatrix) { |
484 |
|
pos[0] += theirview.hoff - .5; |
485 |
|
pos[1] += theirview.voff - .5; |
486 |
|
if (theirview.type == VT_PER) { |
486 |
– |
if (normdist) /* adjust for eye-ray distance */ |
487 |
– |
pos[2] /= sqrt( 1. |
488 |
– |
+ pos[0]*pos[0]*theirview.hn2 |
489 |
– |
+ pos[1]*pos[1]*theirview.vn2 ); |
487 |
|
pos[0] *= pos[2]; |
488 |
|
pos[1] *= pos[2]; |
489 |
|
} |
498 |
|
pos[1] += .5 - ourview.voff; |
499 |
|
return(0); |
500 |
|
} |
501 |
< |
if (viewray(pt, direc, &theirview, pos[0], pos[1]) < 0) |
501 |
> |
if (viewray(pt, direc, &theirview, pos[0], pos[1]) < -FTINY) |
502 |
|
return(-1); |
503 |
|
pt[0] += direc[0]*pos[2]; |
504 |
|
pt[1] += direc[1]*pos[2]; |
520 |
|
/* get back buffer */ |
521 |
|
yback = (int *)malloc(hresolu*sizeof(int)); |
522 |
|
if (yback == NULL) |
523 |
< |
syserror(); |
523 |
> |
syserror(progname); |
524 |
|
for (x = 0; x < hresolu; x++) |
525 |
|
yback[x] = -2; |
526 |
|
/* |
617 |
|
} |
618 |
|
|
619 |
|
|
620 |
+ |
clipaft() /* perform aft clipping as indicated */ |
621 |
+ |
{ |
622 |
+ |
register int x, y; |
623 |
+ |
double tstdist; |
624 |
+ |
double yzn2, vx; |
625 |
+ |
|
626 |
+ |
if (ourview.vaft <= FTINY) |
627 |
+ |
return; |
628 |
+ |
tstdist = ourview.vaft; |
629 |
+ |
for (y = 0; y < vresolu; y++) { |
630 |
+ |
if (ourview.type == VT_PER) { /* adjust distance */ |
631 |
+ |
yzn2 = (y+.5)/vresolu + ourview.voff - .5; |
632 |
+ |
yzn2 = 1. + yzn2*yzn2*ourview.vn2; |
633 |
+ |
tstdist = ourview.vaft * sqrt(yzn2); |
634 |
+ |
} |
635 |
+ |
for (x = 0; x < hresolu; x++) |
636 |
+ |
if (zscan(y)[x] > tstdist) { |
637 |
+ |
if (ourview.type == VT_PER) { |
638 |
+ |
vx = (x+.5)/hresolu + ourview.hoff - .5; |
639 |
+ |
if (zscan(y)[x] <= ourview.vaft * |
640 |
+ |
sqrt(vx*vx*ourview.hn2 + yzn2)) |
641 |
+ |
continue; |
642 |
+ |
} |
643 |
+ |
bzero(pscan(y)[x], sizeof(COLR)); |
644 |
+ |
zscan(y)[x] = 0.0; |
645 |
+ |
} |
646 |
+ |
} |
647 |
+ |
} |
648 |
+ |
|
649 |
+ |
|
650 |
|
writepicture() /* write out picture */ |
651 |
|
{ |
652 |
|
int y; |
654 |
|
fprtresolu(hresolu, vresolu, stdout); |
655 |
|
for (y = vresolu-1; y >= 0; y--) |
656 |
|
if (fwritecolrs(pscan(y), hresolu, stdout) < 0) |
657 |
< |
syserror(); |
657 |
> |
syserror(progname); |
658 |
|
} |
659 |
|
|
660 |
|
|
661 |
|
writedistance(fname) /* write out z file */ |
662 |
|
char *fname; |
663 |
|
{ |
637 |
– |
extern double sqrt(); |
664 |
|
int donorm = normdist && ourview.type == VT_PER; |
665 |
|
int fd; |
666 |
|
int y; |
667 |
|
float *zout; |
668 |
|
|
669 |
< |
if ((fd = open(fname, O_WRONLY|O_CREAT|O_TRUNC, 0666)) == -1) { |
670 |
< |
perror(fname); |
645 |
< |
exit(1); |
646 |
< |
} |
669 |
> |
if ((fd = open(fname, O_WRONLY|O_CREAT|O_TRUNC, 0666)) == -1) |
670 |
> |
syserror(fname); |
671 |
|
if (donorm |
672 |
|
&& (zout = (float *)malloc(hresolu*sizeof(float))) == NULL) |
673 |
< |
syserror(); |
673 |
> |
syserror(progname); |
674 |
|
for (y = vresolu-1; y >= 0; y--) { |
675 |
|
if (donorm) { |
676 |
|
double vx, yzn2; |
685 |
|
} else |
686 |
|
zout = zscan(y); |
687 |
|
if (write(fd, (char *)zout, hresolu*sizeof(float)) |
688 |
< |
< hresolu*sizeof(float)) { |
689 |
< |
perror(fname); |
666 |
< |
exit(1); |
667 |
< |
} |
688 |
> |
< hresolu*sizeof(float)) |
689 |
> |
syserror(fname); |
690 |
|
} |
691 |
|
if (donorm) |
692 |
|
free((char *)zout); |
732 |
|
cp = combuf; |
733 |
|
wp = argv; |
734 |
|
for ( ; ; ) { |
735 |
< |
while (isspace(*cp)) cp++; |
736 |
< |
if (!*cp) break; |
737 |
< |
*wp++ = cp; |
738 |
< |
while (!isspace(*cp)) |
739 |
< |
if (!*cp++) goto done; |
740 |
< |
*cp++ = '\0'; |
735 |
> |
while (isspace(*cp)) /* nullify spaces */ |
736 |
> |
*cp++ = '\0'; |
737 |
> |
if (!*cp) /* all done? */ |
738 |
> |
break; |
739 |
> |
*wp++ = cp; /* add argument to list */ |
740 |
> |
while (*++cp && !isspace(*cp)) |
741 |
> |
; |
742 |
|
} |
720 |
– |
done: |
743 |
|
*wp = NULL; |
744 |
|
/* start process */ |
745 |
|
if ((rval = open_process(PDesc, argv)) < 0) |
746 |
< |
syserror(); |
746 |
> |
syserror(progname); |
747 |
|
if (rval == 0) { |
748 |
|
fprintf(stderr, "%s: command not found\n", argv[0]); |
749 |
|
exit(1); |
784 |
|
register float *fbp; |
785 |
|
register int i; |
786 |
|
|
787 |
+ |
if (queuesiz == 0) |
788 |
+ |
return; |
789 |
|
fbp = fbuf; |
790 |
|
for (i = 0; i < queuesiz; i++) { |
791 |
|
viewray(orig, dir, &ourview, |
819 |
|
} |
820 |
|
|
821 |
|
|
822 |
< |
syserror() /* report error and exit */ |
822 |
> |
syserror(s) /* report error and exit */ |
823 |
> |
char *s; |
824 |
|
{ |
825 |
< |
perror(progname); |
825 |
> |
perror(s); |
826 |
|
exit(1); |
827 |
|
} |