| 16 |
|
|
| 17 |
|
#include "color.h" |
| 18 |
|
|
| 19 |
– |
|
| 19 |
|
extern char *malloc(); |
| 20 |
|
|
| 21 |
|
#define CHECKRAD 1.5 /* radius to check for filtering */ |
| 45 |
|
char *tfname = NULL; |
| 46 |
|
|
| 47 |
|
int xres, yres; /* resolution of input */ |
| 48 |
+ |
double inpaspect = 1.0; /* pixel aspect ratio of input */ |
| 49 |
|
|
| 50 |
|
int xrad; /* x window size */ |
| 51 |
|
int yrad; /* y window size */ |
| 64 |
|
extern char *mktemp(); |
| 65 |
|
extern double atof(), pow(); |
| 66 |
|
extern long ftell(); |
| 67 |
< |
extern int quit(); |
| 67 |
> |
extern int quit(), headline(); |
| 68 |
|
FILE *fin; |
| 69 |
|
long fpos; |
| 70 |
+ |
double outaspect = 0.0; |
| 71 |
|
double d; |
| 72 |
|
int i; |
| 73 |
|
|
| 103 |
|
} else |
| 104 |
|
nrows = atoi(argv[i]); |
| 105 |
|
break; |
| 106 |
+ |
case 'p': |
| 107 |
+ |
i++; |
| 108 |
+ |
outaspect = atof(argv[i]); |
| 109 |
+ |
break; |
| 110 |
|
case 'e': |
| 111 |
|
if (argv[i+1][0] == '+' || argv[i+1][0] == '-') |
| 112 |
|
d = pow(2.0, atof(argv[i+1])); |
| 136 |
|
case '2': |
| 137 |
|
singlepass = 0; |
| 138 |
|
break; |
| 139 |
< |
case 'p': |
| 139 |
> |
case 'n': |
| 140 |
|
npts = atoi(argv[++i]) / 2; |
| 141 |
|
break; |
| 142 |
|
case 's': |
| 190 |
|
fprintf(stderr, "%s: bad # file arguments\n", progname); |
| 191 |
|
quit(1); |
| 192 |
|
} |
| 193 |
< |
/* copy header */ |
| 194 |
< |
copyheader(fin, stdout); |
| 193 |
> |
/* get header */ |
| 194 |
> |
getheader(fin, headline); |
| 195 |
|
/* add new header info. */ |
| 196 |
|
printargs(i, argv, stdout); |
| 197 |
|
/* get picture size */ |
| 199 |
|
fprintf(stderr, "%s: bad picture size\n", progname); |
| 200 |
|
quit(1); |
| 201 |
|
} |
| 202 |
< |
if (ncols > 0) |
| 203 |
< |
x_c = (double)ncols/xres; |
| 199 |
< |
else |
| 202 |
> |
/* compute output resolution */ |
| 203 |
> |
if (ncols <= 0) |
| 204 |
|
ncols = x_c*xres + .5; |
| 205 |
< |
if (nrows > 0) |
| 202 |
< |
y_r = (double)nrows/yres; |
| 203 |
< |
else |
| 205 |
> |
if (nrows <= 0) |
| 206 |
|
nrows = y_r*yres + .5; |
| 207 |
+ |
if (outaspect > .01) { |
| 208 |
+ |
d = inpaspect * yres/xres / outaspect; |
| 209 |
+ |
if (d * ncols > nrows) |
| 210 |
+ |
ncols = nrows / d; |
| 211 |
+ |
else |
| 212 |
+ |
nrows = ncols * d; |
| 213 |
+ |
} |
| 214 |
+ |
x_c = (double)ncols/xres; |
| 215 |
+ |
y_r = (double)nrows/yres; |
| 216 |
|
|
| 217 |
< |
if (singlepass) { |
| 207 |
< |
/* skip exposure, etc. */ |
| 217 |
> |
if (singlepass) { /* skip exposure, etc. */ |
| 218 |
|
pass1default(); |
| 219 |
|
pass2(fin); |
| 220 |
|
quit(0); |
| 234 |
|
} |
| 235 |
|
|
| 236 |
|
|
| 237 |
+ |
headline(s) /* process line from header */ |
| 238 |
+ |
char *s; |
| 239 |
+ |
{ |
| 240 |
+ |
fputs(s, stdout); /* copy to output */ |
| 241 |
+ |
if (isaspect(s)) /* get aspect ratio */ |
| 242 |
+ |
inpaspect *= aspectval(s); |
| 243 |
+ |
} |
| 244 |
+ |
|
| 245 |
+ |
|
| 246 |
|
copyfile(in, out) /* copy a file */ |
| 247 |
|
register FILE *in, *out; |
| 248 |
|
{ |
| 331 |
|
|
| 332 |
|
scan2init() /* prepare scanline arrays */ |
| 333 |
|
{ |
| 334 |
< |
double e; |
| 334 |
> |
double d; |
| 335 |
|
register int i; |
| 336 |
|
|
| 337 |
|
if (rad <= 0.0) { |
| 360 |
|
fprintf(stderr, "%s: out of memory\n", progname); |
| 361 |
|
quit(1); |
| 362 |
|
} |
| 363 |
< |
e = bright(exposure); |
| 364 |
< |
if (e < 1-1e-7 || e > 1+1e-7) /* record exposure */ |
| 365 |
< |
fputexpos(e, stdout); |
| 363 |
> |
/* record pixel aspect and exposure */ |
| 364 |
> |
d = x_c / y_r; |
| 365 |
> |
if (d < .99 || d > 1.01) |
| 366 |
> |
fputaspect(d, stdout); |
| 367 |
> |
d = bright(exposure); |
| 368 |
> |
if (d < .995 || d > 1.005) |
| 369 |
> |
fputexpos(d, stdout); |
| 370 |
|
printf("\n"); |
| 371 |
|
fputresolu(YMAJOR|YDECR, ncols, nrows, stdout); /* resolution */ |
| 372 |
|
} |