| 1 | < | /* Copyright (c) 1993 Regents of the University of California */ | 
| 1 | > | /* Copyright (c) 1997 Regents of the University of California */ | 
| 2 |  |  | 
| 3 |  | #ifndef lint | 
| 4 |  | static char SCCSid[] = "$SunId$ LBL"; | 
| 27 |  | #include  <X11/Xutil.h> | 
| 28 |  | #include  <X11/Xatom.h> | 
| 29 |  |  | 
| 30 | < | #include  "color.h" | 
| 30 | > | #undef  NOPROTO | 
| 31 | > | #define NOPROTO 1 | 
| 32 | > | #include  "tonemap.h" | 
| 33 |  | #include  "view.h" | 
| 34 |  | #include  "x11raster.h" | 
| 35 |  | #include  "random.h" | 
| 36 |  | #include  "resolu.h" | 
| 37 |  |  | 
| 36 | – | #ifdef  __alpha | 
| 37 | – | #define  int4           int | 
| 38 | – | #endif | 
| 38 |  | #ifndef  int4 | 
| 39 | < | #define  int4           long | 
| 39 | > | #define  int4           int             /* most int's are 32-bit */ | 
| 40 |  | #endif | 
| 41 |  |  | 
| 42 |  | #define  FONTNAME       "8x13"          /* text font we'll use */ | 
| 55 |  | #define  redraw(x,y,w,h) patch_raster(wind,(x)-xoff,(y)-yoff,x,y,w,h,ourras) | 
| 56 |  |  | 
| 57 |  | double  gamcor = 2.2;                   /* gamma correction */ | 
| 58 | + | char  *gamstr = NULL;                   /* gamma value override */ | 
| 59 |  |  | 
| 60 |  | int  dither = 1;                        /* dither colors? */ | 
| 61 |  | int  fast = 0;                          /* keep picture in Pixmap? */ | 
| 73 |  | int  yoff = 0;                          /* y image offset */ | 
| 74 |  |  | 
| 75 |  | int  parent = 0;                        /* number of children, -1 if child */ | 
| 76 | + | int  sequential = 0;                    /* display images in sequence */ | 
| 77 |  |  | 
| 78 | + | char  *tout = "od";                     /* output of 't' command */ | 
| 79 | + | int  tinterv = 0;                       /* interval between mouse reports */ | 
| 80 | + |  | 
| 81 | + | int  tmflags = -1;                      /* tone mapping flags (-1 for none) */ | 
| 82 | + |  | 
| 83 |  | VIEW  ourview = STDVIEW;                /* image view parameters */ | 
| 84 |  | int  gotview = 0;                       /* got parameters from file */ | 
| 85 |  |  | 
| 86 |  | COLR  *scanline;                        /* scan line buffer */ | 
| 87 | + | TMbright  *lscan;                       /* encoded luminance scanline */ | 
| 88 | + | BYTE  *cscan;                           /* encoded chroma scanline */ | 
| 89 | + | BYTE  *pscan;                           /* compute pixel scanline */ | 
| 90 |  |  | 
| 91 |  | RESOLU  inpres;                         /* input resolution and ordering */ | 
| 92 |  | int  xmax, ymax;                        /* picture dimensions */ | 
| 121 |  |  | 
| 122 |  | char  errmsg[128]; | 
| 123 |  |  | 
| 124 | < | extern BYTE  clrtab[256][3];            /* global color map */ | 
| 124 | > | BYTE  clrtab[256][3];                   /* global color map */ | 
| 125 |  |  | 
| 126 |  | extern long  ftell(); | 
| 127 |  |  | 
| 128 | + | extern char  *getenv(); | 
| 129 | + |  | 
| 130 |  | Display  *thedisplay; | 
| 131 |  | Atom  closedownAtom, wmProtocolsAtom; | 
| 132 |  |  | 
| 133 | < | int  noop() {} | 
| 133 | > | int  sigrecv; | 
| 134 |  |  | 
| 135 | + | int  onsig() { sigrecv++; } | 
| 136 |  |  | 
| 137 | + |  | 
| 138 |  | main(argc, argv) | 
| 139 |  | int  argc; | 
| 140 |  | char  *argv[]; | 
| 141 |  | { | 
| 129 | – | extern char  *getenv(); | 
| 130 | – | char  *gv; | 
| 142 |  | int  headline(); | 
| 143 |  | int  i; | 
| 144 |  | int  pid; | 
| 145 |  |  | 
| 146 |  | progname = argv[0]; | 
| 136 | – | if ((gv = getenv("GAMMA")) != NULL) | 
| 137 | – | gamcor = atof(gv); | 
| 147 |  |  | 
| 148 |  | for (i = 1; i < argc; i++) | 
| 149 |  | if (argv[i][0] == '-') | 
| 150 |  | switch (argv[i][1]) { | 
| 151 | < | case 'c': | 
| 151 | > | case 'c':                       /* number of colors */ | 
| 152 |  | maxcolors = atoi(argv[++i]); | 
| 153 |  | break; | 
| 154 | < | case 'b': | 
| 154 | > | case 'b':                       /* greyscale only */ | 
| 155 |  | greyscale = !greyscale; | 
| 156 |  | break; | 
| 157 | < | case 'm': | 
| 157 | > | case 'm':                       /* monochrome */ | 
| 158 | > | greyscale = 1; | 
| 159 |  | maxcolors = 2; | 
| 160 |  | break; | 
| 161 | < | case 'd': | 
| 161 | > | case 'd':                       /* display or dither */ | 
| 162 |  | if (argv[i][2] == 'i') | 
| 163 |  | dispname = argv[++i]; | 
| 164 |  | else | 
| 165 |  | dither = !dither; | 
| 166 |  | break; | 
| 167 | < | case 'f': | 
| 167 | > | case 'f':                       /* save pixmap */ | 
| 168 |  | fast = !fast; | 
| 169 |  | break; | 
| 170 | < | case 'e': | 
| 171 | < | if (argv[i+1][0] != '+' && argv[i+1][0] != '-') | 
| 170 | > | case 's':                       /* one at a time */ | 
| 171 | > | sequential = !sequential; | 
| 172 | > | break; | 
| 173 | > | case 'o':                       /* 't' output */ | 
| 174 | > | tout = argv[i]+2; | 
| 175 | > | break; | 
| 176 | > | case 't':                       /* msec interval */ | 
| 177 | > | tinterv = atoi(argv[++i]); | 
| 178 | > | break; | 
| 179 | > | case 'e':                       /* exposure comp. */ | 
| 180 | > | i++; | 
| 181 | > | if (argv[i][0] == 'a') { | 
| 182 | > | tmflags = TM_F_CAMERA; | 
| 183 | > | break; | 
| 184 | > | } | 
| 185 | > | if (argv[i][0] == 'h') { | 
| 186 | > | tmflags = TM_F_HUMAN; | 
| 187 | > | break; | 
| 188 | > | } | 
| 189 | > | if (argv[i][0] != '+' && argv[i][0] != '-') | 
| 190 |  | goto userr; | 
| 191 | < | scale = atoi(argv[++i]); | 
| 191 | > | scale = atoi(argv[i]); | 
| 192 |  | break; | 
| 193 | < | case 'g': | 
| 193 | > | case 'g':                       /* gamma comp. */ | 
| 194 |  | if (argv[i][2] == 'e') | 
| 195 |  | geometry = argv[++i]; | 
| 196 |  | else | 
| 197 | < | gamcor = atof(argv[++i]); | 
| 197 | > | gamstr = argv[++i]; | 
| 198 |  | break; | 
| 199 |  | default: | 
| 200 |  | goto userr; | 
| 207 |  | if (i > argc) | 
| 208 |  | goto userr; | 
| 209 |  | while (i < argc-1) { | 
| 210 | + | sigrecv = 0; | 
| 211 | + | signal(SIGCONT, onsig); | 
| 212 |  | if ((pid=fork()) == 0) {        /* a child for each picture */ | 
| 213 |  | parent = -1; | 
| 214 |  | break; | 
| 216 |  | if (pid < 0) | 
| 217 |  | quiterr("fork failed"); | 
| 218 |  | parent++; | 
| 219 | < | signal(SIGCONT, noop); | 
| 220 | < | pause();                /* wait for wake-up call */ | 
| 219 | > | while (!sigrecv) | 
| 220 | > | pause();        /* wait for wake-up call */ | 
| 221 |  | i++; | 
| 222 |  | } | 
| 223 |  | if (i < argc) {                 /* open picture file */ | 
| 224 |  | fname = argv[i]; | 
| 225 |  | fin = fopen(fname, "r"); | 
| 226 | < | if (fin == NULL) { | 
| 227 | < | sprintf(errmsg, "cannot open file \"%s\"", fname); | 
| 198 | < | quiterr(errmsg); | 
| 199 | < | } | 
| 226 | > | if (fin == NULL) | 
| 227 | > | quiterr("cannot open picture file"); | 
| 228 |  | } | 
| 229 |  | /* get header */ | 
| 230 |  | getheader(fin, headline, NULL); | 
| 241 |  |  | 
| 242 |  | init(argc, argv);                       /* get file and open window */ | 
| 243 |  |  | 
| 216 | – | if (parent < 0) | 
| 217 | – | kill(getppid(), SIGCONT);       /* signal parent if child */ | 
| 218 | – |  | 
| 244 |  | for ( ; ; ) | 
| 245 |  | getevent();             /* main loop */ | 
| 246 |  | userr: | 
| 247 |  | fprintf(stderr, | 
| 248 | < | "Usage: %s [-di disp][[-ge] spec][-b][-m][-d][-f][-c nclrs][-e +/-stops] pic ..\n", | 
| 248 | > | "Usage: %s [-di disp][[-ge] spec][-b][-m][-d][-f][-c nclrs][-e spec][-g gamcor][-s][-ospec][-t intvl] pic ..\n", | 
| 249 |  | progname); | 
| 250 |  | exit(1); | 
| 251 |  | } | 
| 258 |  |  | 
| 259 |  | if (isexpos(s)) | 
| 260 |  | exposure *= exposval(s); | 
| 261 | < | else if (isformat(s)) { | 
| 237 | < | formatval(fmt, s); | 
| 261 | > | else if (formatval(fmt, s)) | 
| 262 |  | wrongformat = strcmp(fmt, COLRFMT); | 
| 263 | < | } else if (isview(s) && sscanview(&ourview, s) > 0) | 
| 263 | > | else if (isview(s) && sscanview(&ourview, s) > 0) | 
| 264 |  | gotview++; | 
| 265 |  | } | 
| 266 |  |  | 
| 294 |  | name += i+1; | 
| 295 |  | if ((thedisplay = XOpenDisplay(dispname)) == NULL) | 
| 296 |  | quiterr("cannot open display"); | 
| 297 | + | /* set gamma value */ | 
| 298 | + | if (gamstr == NULL)             /* get it from the X server */ | 
| 299 | + | gamstr = XGetDefault(thedisplay, "radiance", "gamma"); | 
| 300 | + | if (gamstr == NULL)             /* get it from the environment */ | 
| 301 | + | gamstr = getenv("DISPLAY_GAMMA"); | 
| 302 | + | if (gamstr != NULL) | 
| 303 | + | gamcor = atof(gamstr); | 
| 304 |  | /* get best visual for default screen */ | 
| 305 |  | getbestvis(); | 
| 306 |  | /* store image */ | 
| 327 |  | } | 
| 328 |  | } | 
| 329 |  | /* open window */ | 
| 330 | + | i = CWEventMask|CWCursor|CWBackPixel|CWBorderPixel; | 
| 331 |  | ourwinattr.border_pixel = ourwhite; | 
| 332 |  | ourwinattr.background_pixel = ourblack; | 
| 333 | < | ourwinattr.colormap = XCreateColormap(thedisplay, ourroot, | 
| 334 | < | ourvis.visual, AllocNone); | 
| 333 | > | if (ourvis.visual != DefaultVisual(thedisplay,ourscreen)) { | 
| 334 | > | ourwinattr.colormap = newcmap(thedisplay, ourscreen, ourvis.visual); | 
| 335 | > | i |= CWColormap; | 
| 336 | > | } | 
| 337 |  | ourwinattr.event_mask = ExposureMask|KeyPressMask|ButtonPressMask| | 
| 338 |  | ButtonReleaseMask|ButtonMotionMask|StructureNotifyMask; | 
| 339 |  | ourwinattr.cursor = XCreateFontCursor(thedisplay, XC_diamond_cross); | 
| 340 |  | wind = XCreateWindow(thedisplay, ourroot, xszhints.x, xszhints.y, | 
| 341 |  | xszhints.width, xszhints.height, BORWIDTH, | 
| 342 | < | ourvis.depth, InputOutput, ourvis.visual, CWEventMask| | 
| 343 | < | CWCursor|CWBackPixel|CWBorderPixel|CWColormap, &ourwinattr); | 
| 342 | > | ourvis.depth, InputOutput, ourvis.visual, | 
| 343 | > | i, &ourwinattr); | 
| 344 |  | if (wind == 0) | 
| 345 |  | quiterr("cannot create window"); | 
| 346 |  | width = xmax; | 
| 347 |  | height = ymax; | 
| 348 | < | xgcv.foreground = ourblack; | 
| 315 | < | xgcv.background = ourwhite; | 
| 348 | > | /* prepare graphics drawing context */ | 
| 349 |  | if ((xgcv.font = XLoadFont(thedisplay, FONTNAME)) == 0) | 
| 350 |  | quiterr("cannot get font"); | 
| 351 | + | xgcv.foreground = ourblack; | 
| 352 | + | xgcv.background = ourwhite; | 
| 353 |  | ourgc = XCreateGC(thedisplay, wind, GCForeground|GCBackground| | 
| 354 |  | GCFont, &xgcv); | 
| 355 |  | xgcv.function = GXinvert; | 
| 378 |  | XSetWMProtocols(thedisplay, wind, &closedownAtom, 1); | 
| 379 |  |  | 
| 380 |  | XMapWindow(thedisplay, wind); | 
| 346 | – | return; | 
| 381 |  | } /* end of init */ | 
| 382 |  |  | 
| 383 |  |  | 
| 384 |  | quiterr(err)            /* print message and exit */ | 
| 385 |  | char  *err; | 
| 386 |  | { | 
| 387 | < | if (err != NULL) | 
| 387 | > | register int  es; | 
| 388 | > | int  cs; | 
| 389 | > |  | 
| 390 | > | if (es = err != NULL) | 
| 391 |  | fprintf(stderr, "%s: %s: %s\n", progname, | 
| 392 |  | fname==NULL?"<stdin>":fname, err); | 
| 393 | < | if (wind) { | 
| 394 | < | XDestroyWindow(thedisplay, wind); | 
| 395 | < | XFlush(thedisplay); | 
| 396 | < | } | 
| 397 | < | while (parent > 0 && wait(0) != -1)     /* wait for any children */ | 
| 393 | > | if (thedisplay != NULL) | 
| 394 | > | XCloseDisplay(thedisplay); | 
| 395 | > | if (parent < 0 & sigrecv == 0) | 
| 396 | > | kill(getppid(), SIGCONT); | 
| 397 | > | while (parent > 0 && wait(&cs) != -1) { /* wait for any children */ | 
| 398 | > | if (es == 0) | 
| 399 | > | es = cs>>8 & 0xff; | 
| 400 |  | parent--; | 
| 401 | < | exit(err != NULL); | 
| 401 | > | } | 
| 402 | > | exit(es); | 
| 403 |  | } | 
| 404 |  |  | 
| 405 |  |  | 
| 556 |  | xmax, ymax, 8); | 
| 557 |  | if (ourras == NULL) | 
| 558 |  | goto fail; | 
| 559 | < | if (greyscale | ourvis.class == StaticGray) | 
| 559 | > | if (greyscale) | 
| 560 |  | getgrey(); | 
| 561 |  | else | 
| 562 |  | getmapped(); | 
| 586 |  | map_rcolors(ourras, wind); | 
| 587 |  | if (fast) | 
| 588 |  | make_rpixmap(ourras, wind); | 
| 589 | + | if (!sequential & parent < 0 & sigrecv == 0) { | 
| 590 | + | kill(getppid(), SIGCONT); | 
| 591 | + | sigrecv--; | 
| 592 | + | } | 
| 593 |  | break; | 
| 594 |  | case UnmapNotify: | 
| 595 |  | if (!fast) | 
| 602 |  | case ButtonPress: | 
| 603 |  | if (xev.xbutton.state & (ShiftMask|ControlMask)) | 
| 604 |  | moveimage(&xev.xbutton); | 
| 561 | – | else if (xev.xbutton.button == Button2) | 
| 562 | – | traceray(xev.xbutton.x, xev.xbutton.y); | 
| 605 |  | else | 
| 606 | < | getbox(&xev.xbutton); | 
| 606 | > | switch (xev.xbutton.button) { | 
| 607 | > | case Button1: | 
| 608 | > | getbox(&xev.xbutton); | 
| 609 | > | break; | 
| 610 | > | case Button2: | 
| 611 | > | traceray(xev.xbutton.x, xev.xbutton.y); | 
| 612 | > | break; | 
| 613 | > | case Button3: | 
| 614 | > | trackrays(&xev.xbutton); | 
| 615 | > | break; | 
| 616 | > | } | 
| 617 |  | break; | 
| 618 |  | case ClientMessage: | 
| 619 |  | if ((xev.xclient.message_type == wmProtocolsAtom) && | 
| 624 |  | } | 
| 625 |  |  | 
| 626 |  |  | 
| 627 | < | traceray(xpos, ypos)                    /* print ray corresponding to pixel */ | 
| 627 | > | traceray(xpos, ypos)                    /* print requested pixel data */ | 
| 628 |  | int  xpos, ypos; | 
| 629 |  | { | 
| 630 | + | extern char  *index(); | 
| 631 |  | FLOAT  hv[2]; | 
| 632 |  | FVECT  rorg, rdir; | 
| 633 | + | COLOR  cval; | 
| 634 | + | register char  *cp; | 
| 635 |  |  | 
| 636 | < | if (!gotview) {         /* no view, no can do */ | 
| 637 | < | XBell(thedisplay, 0); | 
| 638 | < | return(-1); | 
| 639 | < | } | 
| 636 | > | box.xmin = xpos; box.xsiz = 1; | 
| 637 | > | box.ymin = ypos; box.ysiz = 1; | 
| 638 | > | avgbox(cval); | 
| 639 | > | scalecolor(cval, 1./exposure); | 
| 640 |  | pix2loc(hv, &inpres, xpos-xoff, ypos-yoff); | 
| 641 | < | if (viewray(rorg, rdir, &ourview, hv[0], hv[1]) < 0) | 
| 642 | < | return(-1); | 
| 643 | < | printf("%e %e %e ", rorg[0], rorg[1], rorg[2]); | 
| 644 | < | printf("%e %e %e\n", rdir[0], rdir[1], rdir[2]); | 
| 641 | > | if (!gotview || viewray(rorg, rdir, &ourview, hv[0], hv[1]) < 0) | 
| 642 | > | rorg[0] = rorg[1] = rorg[2] = | 
| 643 | > | rdir[0] = rdir[1] = rdir[2] = 0.; | 
| 644 | > |  | 
| 645 | > | for (cp = tout; *cp; cp++)      /* print what they asked for */ | 
| 646 | > | switch (*cp) { | 
| 647 | > | case 'o':                       /* origin */ | 
| 648 | > | printf("%e %e %e ", rorg[0], rorg[1], rorg[2]); | 
| 649 | > | break; | 
| 650 | > | case 'd':                       /* direction */ | 
| 651 | > | printf("%e %e %e ", rdir[0], rdir[1], rdir[2]); | 
| 652 | > | break; | 
| 653 | > | case 'v':                       /* radiance value */ | 
| 654 | > | printf("%e %e %e ", colval(cval,RED), | 
| 655 | > | colval(cval,GRN), colval(cval,BLU)); | 
| 656 | > | break; | 
| 657 | > | case 'l':                       /* luminance */ | 
| 658 | > | printf("%e ", luminance(cval)); | 
| 659 | > | break; | 
| 660 | > | case 'p':                       /* pixel position */ | 
| 661 | > | printf("%d %d ", (int)(hv[0]*inpres.xr), | 
| 662 | > | (int)(hv[1]*inpres.yr)); | 
| 663 | > | break; | 
| 664 | > | } | 
| 665 | > | putchar('\n'); | 
| 666 |  | fflush(stdout); | 
| 667 |  | return(0); | 
| 668 |  | } | 
| 733 |  | return(0); | 
| 734 |  | case 't':                               /* trace */ | 
| 735 |  | return(traceray(ekey->x, ekey->y)); | 
| 736 | + | case 'a':                               /* auto exposure */ | 
| 737 | + | if (tmflags == TM_F_CAMERA) | 
| 738 | + | return(0); | 
| 739 | + | tmflags = TM_F_CAMERA; | 
| 740 | + | strcpy(buf, "auto exposure..."); | 
| 741 | + | goto remap; | 
| 742 | + | case 'h':                               /* human response */ | 
| 743 | + | if (tmflags == TM_F_HUMAN) | 
| 744 | + | return(0); | 
| 745 | + | tmflags = TM_F_HUMAN; | 
| 746 | + | strcpy(buf, "human exposure..."); | 
| 747 | + | goto remap; | 
| 748 |  | case '=':                               /* adjust exposure */ | 
| 749 |  | case '@':                               /* adaptation level */ | 
| 750 |  | if (avgbox(cval) == -1) | 
| 751 |  | return(-1); | 
| 752 | < | comp = com=='@' | 
| 753 | < | ? 106./pow(1.219+pow(luminance(cval)/exposure,.4),2.5)/exposure | 
| 754 | < | : .5/bright(cval) ; | 
| 752 | > | comp = bright(cval); | 
| 753 | > | if (comp < 1e-20) { | 
| 754 | > | XBell(thedisplay, 0); | 
| 755 | > | return(-1); | 
| 756 | > | } | 
| 757 | > | if (com == '@') | 
| 758 | > | comp = 106./exposure/ | 
| 759 | > | pow(1.219+pow(comp*WHTEFFICACY/exposure,.4),2.5); | 
| 760 | > | else | 
| 761 | > | comp = .5/comp; | 
| 762 |  | comp = log(comp)/.69315 - scale; | 
| 763 |  | n = comp < 0 ? comp-.5 : comp+.5 ;      /* round */ | 
| 764 | < | if (n == 0) | 
| 765 | < | return(0); | 
| 766 | < | scale_rcolors(ourras, pow(2.0, (double)n)); | 
| 764 | > | if (tmflags != -1) | 
| 765 | > | tmflags = -1;           /* turn off tone mapping */ | 
| 766 | > | else { | 
| 767 | > | if (n == 0)             /* else check if any change */ | 
| 768 | > | return(0); | 
| 769 | > | scale_rcolors(ourras, pow(2.0, (double)n)); | 
| 770 | > | } | 
| 771 |  | scale += n; | 
| 772 |  | sprintf(buf, "%+d", scale); | 
| 773 | + | remap: | 
| 774 |  | XDrawImageString(thedisplay, wind, ourgc, | 
| 775 |  | box.xmin, box.ymin+box.ysiz, buf, strlen(buf)); | 
| 776 |  | XFlush(thedisplay); | 
| 865 |  | } | 
| 866 |  |  | 
| 867 |  |  | 
| 868 | + | trackrays(ebut)                         /* trace rays as mouse moves */ | 
| 869 | + | XButtonPressedEvent  *ebut; | 
| 870 | + | { | 
| 871 | + | XEvent  e; | 
| 872 | + | unsigned long   lastrept; | 
| 873 | + |  | 
| 874 | + | traceray(ebut->x, ebut->y); | 
| 875 | + | lastrept = ebut->time; | 
| 876 | + | XMaskEvent(thedisplay, ButtonReleaseMask|ButtonMotionMask, &e); | 
| 877 | + | while (e.type == MotionNotify) { | 
| 878 | + | if (e.xmotion.time >= lastrept + tinterv) { | 
| 879 | + | traceray(e.xmotion.x, e.xmotion.y); | 
| 880 | + | lastrept = e.xmotion.time; | 
| 881 | + | } | 
| 882 | + | XMaskEvent(thedisplay,ButtonReleaseMask|ButtonMotionMask,&e); | 
| 883 | + | } | 
| 884 | + | } | 
| 885 | + |  | 
| 886 | + |  | 
| 887 |  | revbox(x0, y0, x1, y1)                  /* draw box with reversed lines */ | 
| 888 |  | int  x0, y0, x1, y1; | 
| 889 |  | { | 
| 942 |  | } | 
| 943 |  |  | 
| 944 |  |  | 
| 945 | + | make_tonemap()                  /* initialize tone mapping */ | 
| 946 | + | { | 
| 947 | + | int  flags, y; | 
| 948 | + |  | 
| 949 | + | if (tmflags != -1 && fname == NULL) { | 
| 950 | + | fprintf(stderr, "%s: cannot adjust tone of standard input\n", | 
| 951 | + | progname); | 
| 952 | + | tmflags = -1; | 
| 953 | + | } | 
| 954 | + | if (tmflags == -1) { | 
| 955 | + | setcolrcor(pow, 1.0/gamcor); | 
| 956 | + | return; | 
| 957 | + | } | 
| 958 | + | if (tmTop != NULL)              /* already initialized? */ | 
| 959 | + | return; | 
| 960 | + | flags = tmflags; | 
| 961 | + | if (greyscale) flags |= TM_F_BW; | 
| 962 | + | /* initialize tm library */ | 
| 963 | + | if (tmInit(flags, stdprims, gamcor) == NULL) | 
| 964 | + | goto memerr; | 
| 965 | + | if (tmSetSpace(stdprims, WHTEFFICACY/exposure)) | 
| 966 | + | goto tmerr; | 
| 967 | + | /* allocate encoding buffers */ | 
| 968 | + | if ((lscan = (TMbright *)malloc(xmax*sizeof(TMbright))) == NULL) | 
| 969 | + | goto memerr; | 
| 970 | + | if (greyscale) { | 
| 971 | + | cscan = TM_NOCHROM; | 
| 972 | + | if ((pscan = (BYTE *)malloc(sizeof(BYTE)*xmax)) == NULL) | 
| 973 | + | goto memerr; | 
| 974 | + | } else if ((pscan=cscan = (BYTE *)malloc(3*sizeof(BYTE)*xmax)) == NULL) | 
| 975 | + | goto memerr; | 
| 976 | + | /* compute picture histogram */ | 
| 977 | + | for (y = 0; y < ymax; y++) { | 
| 978 | + | getscan(y); | 
| 979 | + | if (tmCvColrs(lscan, TM_NOCHROM, scanline, xmax)) | 
| 980 | + | goto tmerr; | 
| 981 | + | if (tmAddHisto(lscan, xmax, 1)) | 
| 982 | + | goto tmerr; | 
| 983 | + | } | 
| 984 | + | /* compute tone mapping */ | 
| 985 | + | if (tmComputeMapping(gamcor, 0., 0.)) | 
| 986 | + | goto tmerr; | 
| 987 | + | return; | 
| 988 | + | memerr: | 
| 989 | + | quiterr("out of memory in make_tonemap"); | 
| 990 | + | tmerr: | 
| 991 | + | quiterr("tone mapping error"); | 
| 992 | + | } | 
| 993 | + |  | 
| 994 | + |  | 
| 995 | + | tmap_colrs(scn, len)            /* apply tone mapping to scanline */ | 
| 996 | + | register COLR  *scn; | 
| 997 | + | int  len; | 
| 998 | + | { | 
| 999 | + | register BYTE  *ps; | 
| 1000 | + |  | 
| 1001 | + | if (tmflags == -1) { | 
| 1002 | + | if (scale) | 
| 1003 | + | shiftcolrs(scn, len, scale); | 
| 1004 | + | colrs_gambs(scn, len); | 
| 1005 | + | return; | 
| 1006 | + | } | 
| 1007 | + | if (len > xmax) | 
| 1008 | + | quiterr("code error 1 in tmap_colrs"); | 
| 1009 | + | if (tmCvColrs(lscan, cscan, scn, len)) | 
| 1010 | + | goto tmerr; | 
| 1011 | + | if (tmMapPixels(pscan, lscan, cscan, len)) | 
| 1012 | + | goto tmerr; | 
| 1013 | + | ps = pscan; | 
| 1014 | + | if (greyscale) | 
| 1015 | + | while (len--) { | 
| 1016 | + | scn[0][RED] = scn[0][GRN] = scn[0][BLU] = *ps++; | 
| 1017 | + | scn[0][EXP] = COLXS; | 
| 1018 | + | scn++; | 
| 1019 | + | } | 
| 1020 | + | else | 
| 1021 | + | while (len--) { | 
| 1022 | + | scn[0][RED] = *ps++; | 
| 1023 | + | scn[0][GRN] = *ps++; | 
| 1024 | + | scn[0][BLU] = *ps++; | 
| 1025 | + | scn[0][EXP] = COLXS; | 
| 1026 | + | scn++; | 
| 1027 | + | } | 
| 1028 | + | return; | 
| 1029 | + | tmerr: | 
| 1030 | + | quiterr("tone mapping error"); | 
| 1031 | + | } | 
| 1032 | + |  | 
| 1033 | + |  | 
| 1034 |  | getmono()                       /* get monochrome data */ | 
| 1035 |  | { | 
| 1036 |  | register unsigned char  *dp; | 
| 1120 |  | int     y; | 
| 1121 |  | register unsigned int4  *dp; | 
| 1122 |  | register int    x; | 
| 1123 | < | /* set gamma correction */ | 
| 1124 | < | setcolrgam(gamcor); | 
| 1123 | > | /* initialize tone mapping */ | 
| 1124 | > | make_tonemap(); | 
| 1125 |  | /* read and convert file */ | 
| 1126 |  | dp = (unsigned int4 *)ourdata; | 
| 1127 |  | for (y = 0; y < ymax; y++) { | 
| 1128 |  | getscan(y); | 
| 1129 |  | add2icon(y, scanline); | 
| 1130 | < | if (scale) | 
| 923 | < | shiftcolrs(scanline, xmax, scale); | 
| 924 | < | colrs_gambs(scanline, xmax); | 
| 1130 | > | tmap_colrs(scanline, xmax); | 
| 1131 |  | if (ourras->image->blue_mask & 1) | 
| 1132 |  | for (x = 0; x < xmax; x++) | 
| 1133 | < | *dp++ = scanline[x][RED] << 16 | | 
| 1134 | < | scanline[x][GRN] << 8 | | 
| 1135 | < | scanline[x][BLU] ; | 
| 1133 | > | *dp++ = (unsigned int4)scanline[x][RED] << 16 | | 
| 1134 | > | (unsigned int4)scanline[x][GRN] << 8 | | 
| 1135 | > | (unsigned int4)scanline[x][BLU] ; | 
| 1136 |  | else | 
| 1137 |  | for (x = 0; x < xmax; x++) | 
| 1138 | < | *dp++ = scanline[x][RED] | | 
| 1139 | < | scanline[x][GRN] << 8 | | 
| 1140 | < | scanline[x][BLU] << 16 ; | 
| 1138 | > | *dp++ = (unsigned int4)scanline[x][RED] | | 
| 1139 | > | (unsigned int4)scanline[x][GRN] << 8 | | 
| 1140 | > | (unsigned int4)scanline[x][BLU] << 16 ; | 
| 1141 |  | } | 
| 1142 |  | } | 
| 1143 |  |  | 
| 1147 |  | int     y; | 
| 1148 |  | register unsigned char  *dp; | 
| 1149 |  | register int    x; | 
| 1150 | < | /* set gamma correction */ | 
| 1151 | < | setcolrgam(gamcor); | 
| 1150 | > | /* initialize tone mapping */ | 
| 1151 | > | make_tonemap(); | 
| 1152 |  | /* read and convert file */ | 
| 1153 |  | dp = ourdata; | 
| 1154 |  | for (y = 0; y < ymax; y++) { | 
| 1155 |  | getscan(y); | 
| 1156 |  | add2icon(y, scanline); | 
| 1157 | < | if (scale) | 
| 952 | < | shiftcolrs(scanline, xmax, scale); | 
| 953 | < | for (x = 0; x < xmax; x++) | 
| 954 | < | scanline[x][GRN] = normbright(scanline[x]); | 
| 955 | < | colrs_gambs(scanline, xmax); | 
| 1157 | > | tmap_colrs(scanline, xmax); | 
| 1158 |  | if (maxcolors < 256) | 
| 1159 |  | for (x = 0; x < xmax; x++) | 
| 1160 | < | *dp++ = ((long)scanline[x][GRN] * | 
| 1160 | > | *dp++ = ((int4)scanline[x][GRN] * | 
| 1161 |  | maxcolors + maxcolors/2) >> 8; | 
| 1162 |  | else | 
| 1163 |  | for (x = 0; x < xmax; x++) | 
| 1165 |  | } | 
| 1166 |  | for (x = 0; x < maxcolors; x++) | 
| 1167 |  | clrtab[x][RED] = clrtab[x][GRN] = | 
| 1168 | < | clrtab[x][BLU] = ((long)x*256 + 128)/maxcolors; | 
| 1168 | > | clrtab[x][BLU] = ((int4)x*256 + 128)/maxcolors; | 
| 1169 |  | } | 
| 1170 |  |  | 
| 1171 |  |  | 
| 1172 |  | getmapped()                     /* get color-mapped data */ | 
| 1173 |  | { | 
| 1174 |  | int     y; | 
| 1175 | < | /* set gamma correction */ | 
| 1176 | < | setcolrgam(gamcor); | 
| 1175 | > | /* make sure we can do it first */ | 
| 1176 | > | if (fname == NULL) | 
| 1177 | > | quiterr("cannot map colors from standard input"); | 
| 1178 | > | /* initialize tone mapping */ | 
| 1179 | > | make_tonemap(); | 
| 1180 |  | /* make histogram */ | 
| 1181 | < | new_histo(); | 
| 1181 | > | if (new_histo((int4)xmax*ymax) == -1) | 
| 1182 | > | quiterr("cannot initialize histogram"); | 
| 1183 |  | for (y = 0; y < ymax; y++) { | 
| 1184 |  | if (getscan(y) < 0) | 
| 1185 | < | quiterr("seek error in getmapped"); | 
| 1185 | > | break; | 
| 1186 |  | add2icon(y, scanline); | 
| 1187 | < | if (scale) | 
| 982 | < | shiftcolrs(scanline, xmax, scale); | 
| 983 | < | colrs_gambs(scanline, xmax); | 
| 1187 | > | tmap_colrs(scanline, xmax); | 
| 1188 |  | cnt_colrs(scanline, xmax); | 
| 1189 |  | } | 
| 1190 |  | /* map pixels */ | 
| 1191 |  | if (!new_clrtab(maxcolors)) | 
| 1192 |  | quiterr("cannot create color map"); | 
| 1193 |  | for (y = 0; y < ymax; y++) { | 
| 1194 | < | if (getscan(y) < 0) | 
| 1195 | < | quiterr("seek error in getmapped"); | 
| 992 | < | if (scale) | 
| 993 | < | shiftcolrs(scanline, xmax, scale); | 
| 994 | < | colrs_gambs(scanline, xmax); | 
| 1194 | > | getscan(y); | 
| 1195 | > | tmap_colrs(scanline, xmax); | 
| 1196 |  | if (dither) | 
| 1197 |  | dith_colrs(ourdata+y*xmax, scanline, xmax); | 
| 1198 |  | else | 
| 1232 |  | getscan(y) | 
| 1233 |  | int  y; | 
| 1234 |  | { | 
| 1235 | + | static int  trunced = -1;               /* truncated file? */ | 
| 1236 | + | skipit: | 
| 1237 | + | if (trunced >= 0 && y >= trunced) { | 
| 1238 | + | bzero(scanline, xmax*sizeof(COLR)); | 
| 1239 | + | return(-1); | 
| 1240 | + | } | 
| 1241 |  | if (y != cury) { | 
| 1242 |  | if (scanpos == NULL || scanpos[y] == -1) | 
| 1243 |  | return(-1); | 
| 1247 |  | } else if (scanpos != NULL && scanpos[y] == -1) | 
| 1248 |  | scanpos[y] = ftell(fin); | 
| 1249 |  |  | 
| 1250 | < | if (freadcolrs(scanline, xmax, fin) < 0) | 
| 1251 | < | quiterr("read error"); | 
| 1252 | < |  | 
| 1250 | > | if (freadcolrs(scanline, xmax, fin) < 0) { | 
| 1251 | > | fprintf(stderr, "%s: %s: unfinished picture\n", | 
| 1252 | > | progname, fname==NULL?"<stdin>":fname); | 
| 1253 | > | trunced = y; | 
| 1254 | > | goto skipit; | 
| 1255 | > | } | 
| 1256 |  | cury++; | 
| 1257 |  | return(0); | 
| 1258 |  | } |