| 1 | – | /* Copyright (c) 1997 Regents of the University of California */ | 
| 2 | – |  | 
| 1 |  | #ifndef lint | 
| 2 | < | static char SCCSid[] = "$SunId$ LBL"; | 
| 2 | > | static const char RCSid[] = "$Id"; | 
| 3 |  | #endif | 
| 6 | – |  | 
| 4 |  | /* | 
| 5 |  | *  x11image.c - driver for X-windows | 
| 6 |  | * | 
| 18 |  |  | 
| 19 |  | #include  "standard.h" | 
| 20 |  |  | 
| 21 | + | #include  <string.h> | 
| 22 |  | #include  <signal.h> | 
| 23 |  | #include  <X11/Xlib.h> | 
| 24 |  | #include  <X11/cursorfont.h> | 
| 26 |  | #include  <X11/Xatom.h> | 
| 27 |  |  | 
| 28 |  | #include  "color.h" | 
| 29 | + | #include  "tonemap.h" | 
| 30 |  | #include  "view.h" | 
| 31 |  | #include  "x11raster.h" | 
| 32 |  | #include  "random.h" | 
| 34 | – | #include  "resolu.h" | 
| 33 |  |  | 
| 36 | – | #ifdef  __alpha | 
| 37 | – | #define  int4           int | 
| 38 | – | #endif | 
| 39 | – | #ifndef  int4 | 
| 40 | – | #define  int4           long | 
| 41 | – | #endif | 
| 42 | – |  | 
| 34 |  | #define  FONTNAME       "8x13"          /* text font we'll use */ | 
| 35 |  |  | 
| 36 |  | #define  CTRL(c)        ((c)-'@') | 
| 39 |  |  | 
| 40 |  | #define  ICONSIZ        (8*10)          /* maximum icon dimension (even 8) */ | 
| 41 |  |  | 
| 42 | + | #define  FIXWEIGHT      20              /* weight to add for fixation points */ | 
| 43 | + |  | 
| 44 |  | #define  ourscreen      DefaultScreen(thedisplay) | 
| 45 |  | #define  ourroot        RootWindow(thedisplay,ourscreen) | 
| 46 |  |  | 
| 72 |  | char  *tout = "od";                     /* output of 't' command */ | 
| 73 |  | int  tinterv = 0;                       /* interval between mouse reports */ | 
| 74 |  |  | 
| 75 | + | int  tmflags = TM_F_LINEAR;             /* tone mapping flags */ | 
| 76 | + |  | 
| 77 |  | VIEW  ourview = STDVIEW;                /* image view parameters */ | 
| 78 |  | int  gotview = 0;                       /* got parameters from file */ | 
| 79 |  |  | 
| 80 |  | COLR  *scanline;                        /* scan line buffer */ | 
| 81 | + | TMbright  *lscan;                       /* encoded luminance scanline */ | 
| 82 | + | BYTE  *cscan;                           /* encoded chroma scanline */ | 
| 83 | + | BYTE  *pscan;                           /* compute pixel scanline */ | 
| 84 |  |  | 
| 85 |  | RESOLU  inpres;                         /* input resolution and ordering */ | 
| 86 |  | int  xmax, ymax;                        /* picture dimensions */ | 
| 87 |  | int  width, height;                     /* window size */ | 
| 88 |  | char  *fname = NULL;                    /* input file name */ | 
| 89 | < | FILE  *fin = stdin;                     /* input file */ | 
| 89 | > | FILE  *fin = NULL;                      /* input file */ | 
| 90 |  | long  *scanpos = NULL;                  /* scan line positions in file */ | 
| 91 |  | int  cury = 0;                          /* current scan location */ | 
| 92 |  |  | 
| 104 |  |  | 
| 105 |  | struct { | 
| 106 |  | int  xmin, ymin, xsiz, ysiz; | 
| 107 | < | }  box = {0, 0, 0, 0};                  /* current box */ | 
| 107 | > | }  bbox = {0, 0, 0, 0};                 /* current bbox */ | 
| 108 |  |  | 
| 109 |  | char  *geometry = NULL;                 /* geometry specification */ | 
| 110 |  |  | 
| 126 |  |  | 
| 127 |  | int  sigrecv; | 
| 128 |  |  | 
| 129 | < | int  onsig() { sigrecv++; } | 
| 129 | > | void  onsig(int i) { sigrecv++; } | 
| 130 |  |  | 
| 131 |  |  | 
| 132 |  | main(argc, argv) | 
| 138 |  | int  pid; | 
| 139 |  |  | 
| 140 |  | progname = argv[0]; | 
| 141 | + | fin = stdin; | 
| 142 |  |  | 
| 143 |  | for (i = 1; i < argc; i++) | 
| 144 |  | if (argv[i][0] == '-') | 
| 172 |  | tinterv = atoi(argv[++i]); | 
| 173 |  | break; | 
| 174 |  | case 'e':                       /* exposure comp. */ | 
| 175 | < | if (argv[i+1][0] != '+' && argv[i+1][0] != '-') | 
| 175 | > | i++; | 
| 176 | > | if (argv[i][0] == 'a') { | 
| 177 | > | tmflags = TM_F_CAMERA; | 
| 178 | > | break; | 
| 179 | > | } | 
| 180 | > | if (argv[i][0] == 'h') { | 
| 181 | > | tmflags = TM_F_HUMAN; | 
| 182 | > | break; | 
| 183 | > | } | 
| 184 | > | if (argv[i][0] != '+' && argv[i][0] != '-') | 
| 185 |  | goto userr; | 
| 186 | < | scale = atoi(argv[++i]); | 
| 186 | > | scale = atoi(argv[i]); | 
| 187 |  | break; | 
| 188 |  | case 'g':                       /* gamma comp. */ | 
| 189 |  | if (argv[i][2] == 'e') | 
| 240 |  | getevent();             /* main loop */ | 
| 241 |  | userr: | 
| 242 |  | fprintf(stderr, | 
| 243 | < | "Usage: %s [-di disp][[-ge] spec][-b][-m][-d][-f][-c nclrs][-e +/-stops][-g gamcor][-s][-ospec][-t intvl] pic ..\n", | 
| 243 | > | "Usage: %s [-di disp][[-ge] spec][-b][-m][-d][-f][-c nclrs][-e spec][-g gamcor][-s][-ospec][-t intvl] pic ..\n", | 
| 244 |  | progname); | 
| 245 |  | exit(1); | 
| 246 |  | } | 
| 247 |  |  | 
| 248 |  |  | 
| 249 | + | int | 
| 250 |  | headline(s)             /* get relevant info from header */ | 
| 251 |  | char  *s; | 
| 252 |  | { | 
| 254 |  |  | 
| 255 |  | if (isexpos(s)) | 
| 256 |  | exposure *= exposval(s); | 
| 257 | < | else if (isformat(s)) { | 
| 249 | < | formatval(fmt, s); | 
| 257 | > | else if (formatval(fmt, s)) | 
| 258 |  | wrongformat = strcmp(fmt, COLRFMT); | 
| 259 | < | } else if (isview(s) && sscanview(&ourview, s) > 0) | 
| 259 | > | else if (isview(s) && sscanview(&ourview, s) > 0) | 
| 260 |  | gotview++; | 
| 261 | + | return(0); | 
| 262 |  | } | 
| 263 |  |  | 
| 264 |  |  | 
| 384 |  | register int  es; | 
| 385 |  | int  cs; | 
| 386 |  |  | 
| 387 | < | if (es = err != NULL) | 
| 387 | > | if ( (es = err != NULL) ) | 
| 388 |  | fprintf(stderr, "%s: %s: %s\n", progname, | 
| 389 |  | fname==NULL?"<stdin>":fname, err); | 
| 390 |  | if (thedisplay != NULL) | 
| 391 |  | XCloseDisplay(thedisplay); | 
| 392 | < | if (parent < 0 & sigrecv == 0) | 
| 392 | > | if ((parent < 0) & (sigrecv == 0)) | 
| 393 |  | kill(getppid(), SIGCONT); | 
| 394 |  | while (parent > 0 && wait(&cs) != -1) { /* wait for any children */ | 
| 395 |  | if (es == 0) | 
| 414 |  | return(-1); | 
| 415 |  | if (v1->depth == 32 && v2->depth == 24) | 
| 416 |  | return(1); | 
| 417 | < | return(0); | 
| 417 | > | /* go for maximum depth otherwise */ | 
| 418 | > | return(v2->depth - v1->depth); | 
| 419 |  | } | 
| 420 |  | /* don't be too greedy */ | 
| 421 |  | if (maxcolors <= 1<<v1->depth && maxcolors <= 1<<v2->depth) | 
| 422 |  | return(v1->depth - v2->depth); | 
| 423 |  | return(v2->depth - v1->depth); | 
| 424 |  | } | 
| 425 | < | /* prefer Pseudo when < 24-bit */ | 
| 425 | > | /* prefer Pseudo when < 15-bit */ | 
| 426 |  | if ((v1->class == TrueColor || v1->class == DirectColor) && | 
| 427 | < | v1->depth < 24) | 
| 427 | > | v1->depth < 15) | 
| 428 |  | bad1 = 1; | 
| 429 |  | if ((v2->class == TrueColor || v2->class == DirectColor) && | 
| 430 | < | v2->depth < 24) | 
| 430 | > | v2->depth < 15) | 
| 431 |  | bad2 = -1; | 
| 432 |  | if (bad1 | bad2) | 
| 433 |  | return(bad1+bad2); | 
| 494 |  | if (viscmp(&xvi[i],&ourvis) > 0) | 
| 495 |  | quiterr("inadequate visuals on this screen"); | 
| 496 |  | /* OK, we'll use it */ | 
| 497 | < | copystruct(&ourvis, &xvi[i]); | 
| 497 | > | ourvis = xvi[i]; | 
| 498 |  | #ifdef DEBUG | 
| 499 |  | fprintf(stderr, "Selected visual type %s, depth %d\n", | 
| 500 |  | vistype[ourvis.class], ourvis.depth); | 
| 537 |  | if (ourras == NULL) | 
| 538 |  | goto fail; | 
| 539 |  | getmono(); | 
| 540 | < | } else if (ourvis.class == TrueColor | ourvis.class == DirectColor) { | 
| 541 | < | ourdata = (unsigned char *)malloc(sizeof(int4)*xmax*ymax); | 
| 540 | > | } else if ((ourvis.class == TrueColor) | (ourvis.class == DirectColor)) { | 
| 541 | > | int  datsiz = ourvis.depth>16 ? sizeof(int32) : sizeof(int16); | 
| 542 | > | ourdata = (unsigned char *)malloc(datsiz*xmax*ymax); | 
| 543 |  | if (ourdata == NULL) | 
| 544 |  | goto fail; | 
| 545 | < | ourras = make_raster(thedisplay, &ourvis, sizeof(int4)*8, | 
| 546 | < | ourdata, xmax, ymax, 32); | 
| 545 | > | ourras = make_raster(thedisplay, &ourvis, datsiz*8, | 
| 546 | > | ourdata, xmax, ymax, datsiz*8); | 
| 547 |  | if (ourras == NULL) | 
| 548 |  | goto fail; | 
| 549 |  | getfull(); | 
| 555 |  | xmax, ymax, 8); | 
| 556 |  | if (ourras == NULL) | 
| 557 |  | goto fail; | 
| 558 | < | if (greyscale | ourvis.class == StaticGray) | 
| 558 | > | if (greyscale) | 
| 559 |  | getgrey(); | 
| 560 |  | else | 
| 561 |  | getmapped(); | 
| 585 |  | map_rcolors(ourras, wind); | 
| 586 |  | if (fast) | 
| 587 |  | make_rpixmap(ourras, wind); | 
| 588 | < | if (!sequential & parent < 0 & sigrecv == 0) { | 
| 588 | > | if ((!sequential) & (parent < 0) & (sigrecv == 0)) { | 
| 589 |  | kill(getppid(), SIGCONT); | 
| 590 |  | sigrecv--; | 
| 591 |  | } | 
| 626 |  | traceray(xpos, ypos)                    /* print requested pixel data */ | 
| 627 |  | int  xpos, ypos; | 
| 628 |  | { | 
| 629 | < | extern char  *index(); | 
| 619 | < | FLOAT  hv[2]; | 
| 629 | > | RREAL  hv[2]; | 
| 630 |  | FVECT  rorg, rdir; | 
| 631 |  | COLOR  cval; | 
| 632 |  | register char  *cp; | 
| 633 |  |  | 
| 634 | < | box.xmin = xpos; box.xsiz = 1; | 
| 635 | < | box.ymin = ypos; box.ysiz = 1; | 
| 634 | > | bbox.xmin = xpos; bbox.xsiz = 1; | 
| 635 | > | bbox.ymin = ypos; bbox.ysiz = 1; | 
| 636 |  | avgbox(cval); | 
| 637 |  | scalecolor(cval, 1./exposure); | 
| 638 |  | pix2loc(hv, &inpres, xpos-xoff, ypos-yoff); | 
| 674 |  | XColor  cvx; | 
| 675 |  | int  com, n; | 
| 676 |  | double  comp; | 
| 677 | < | FLOAT  hv[2]; | 
| 677 | > | RREAL  hv[2]; | 
| 678 |  |  | 
| 679 |  | n = XLookupString(ekey, buf, sizeof(buf), NULL, NULL); | 
| 680 |  | if (n == 0) | 
| 689 |  | case '\r': | 
| 690 |  | case 'l': | 
| 691 |  | case 'c':                               /* value */ | 
| 692 | < | if (avgbox(cval) == -1) | 
| 692 | > | if (!avgbox(cval)) | 
| 693 |  | return(-1); | 
| 694 |  | switch (com) { | 
| 695 |  | case '\n': | 
| 697 |  | sprintf(buf, "%.3f", intens(cval)/exposure); | 
| 698 |  | break; | 
| 699 |  | case 'l':                               /* luminance */ | 
| 700 | < | sprintf(buf, "%.0fL", luminance(cval)/exposure); | 
| 700 | > | sprintf(buf, "%.1fL", luminance(cval)/exposure); | 
| 701 |  | break; | 
| 702 |  | case 'c':                               /* color */ | 
| 703 |  | comp = pow(2.0, (double)scale); | 
| 708 |  | break; | 
| 709 |  | } | 
| 710 |  | XDrawImageString(thedisplay, wind, ourgc, | 
| 711 | < | box.xmin, box.ymin+box.ysiz, buf, strlen(buf)); | 
| 711 | > | bbox.xmin, bbox.ymin+bbox.ysiz, buf, strlen(buf)); | 
| 712 |  | return(0); | 
| 713 |  | case 'i':                               /* identify (contour) */ | 
| 714 |  | if (ourras->pixels == NULL) | 
| 731 |  | return(0); | 
| 732 |  | case 't':                               /* trace */ | 
| 733 |  | return(traceray(ekey->x, ekey->y)); | 
| 734 | + | case 'a':                               /* auto exposure */ | 
| 735 | + | if (fname == NULL) | 
| 736 | + | return(-1); | 
| 737 | + | tmflags = TM_F_CAMERA; | 
| 738 | + | strcpy(buf, "auto exposure..."); | 
| 739 | + | goto remap; | 
| 740 | + | case 'h':                               /* human response */ | 
| 741 | + | if (fname == NULL) | 
| 742 | + | return(-1); | 
| 743 | + | tmflags = TM_F_HUMAN; | 
| 744 | + | strcpy(buf, "human exposure..."); | 
| 745 | + | goto remap; | 
| 746 |  | case '=':                               /* adjust exposure */ | 
| 747 |  | case '@':                               /* adaptation level */ | 
| 748 | < | if (avgbox(cval) == -1) | 
| 748 | > | if (!avgbox(cval)) | 
| 749 |  | return(-1); | 
| 750 |  | comp = bright(cval); | 
| 751 |  | if (comp < 1e-20) { | 
| 759 |  | comp = .5/comp; | 
| 760 |  | comp = log(comp)/.69315 - scale; | 
| 761 |  | n = comp < 0 ? comp-.5 : comp+.5 ;      /* round */ | 
| 762 | < | if (n == 0) | 
| 763 | < | return(0); | 
| 764 | < | scale_rcolors(ourras, pow(2.0, (double)n)); | 
| 762 | > | if (tmflags != TM_F_LINEAR) | 
| 763 | > | tmflags = TM_F_LINEAR;  /* turn off tone mapping */ | 
| 764 | > | else { | 
| 765 | > | if (n == 0)             /* else check if any change */ | 
| 766 | > | return(0); | 
| 767 | > | scale_rcolors(ourras, pow(2.0, (double)n)); | 
| 768 | > | } | 
| 769 |  | scale += n; | 
| 770 |  | sprintf(buf, "%+d", scale); | 
| 771 | + | remap: | 
| 772 |  | XDrawImageString(thedisplay, wind, ourgc, | 
| 773 | < | box.xmin, box.ymin+box.ysiz, buf, strlen(buf)); | 
| 773 | > | bbox.xmin, bbox.ymin+bbox.ysiz, buf, strlen(buf)); | 
| 774 |  | XFlush(thedisplay); | 
| 775 | < | free(ourdata); | 
| 775 | > | /* free(ourdata);       This is done in XDestroyImage()! */ | 
| 776 |  | free_raster(ourras); | 
| 777 |  | getras(); | 
| 778 |  | /* fall through */ | 
| 794 |  | free_rpixmap(ourras); | 
| 795 |  | return(0); | 
| 796 |  | case '0':                               /* recenter origin */ | 
| 797 | < | if (xoff == 0 & yoff == 0) | 
| 797 | > | if ((xoff == 0) & (yoff == 0)) | 
| 798 |  | return(0); | 
| 799 |  | xoff = yoff = 0; | 
| 800 |  | XClearWindow(thedisplay, wind); | 
| 801 |  | redraw(0, 0, width, height); | 
| 802 |  | return(0); | 
| 803 |  | case ' ':                               /* clear */ | 
| 804 | < | redraw(box.xmin, box.ymin, box.xsiz, box.ysiz); | 
| 804 | > | redraw(bbox.xmin, bbox.ymin, bbox.xsiz, bbox.ysiz); | 
| 805 |  | return(0); | 
| 806 |  | default: | 
| 807 |  | XBell(thedisplay, 0); | 
| 835 |  | } | 
| 836 |  |  | 
| 837 |  |  | 
| 838 | < | getbox(ebut)                            /* get new box */ | 
| 838 | > | getbox(ebut)                            /* get new bbox */ | 
| 839 |  | XButtonPressedEvent  *ebut; | 
| 840 |  | { | 
| 841 |  | XEvent  e; | 
| 842 |  |  | 
| 843 |  | XMaskEvent(thedisplay, ButtonReleaseMask|ButtonMotionMask, &e); | 
| 844 |  | while (e.type == MotionNotify) { | 
| 845 | < | revbox(ebut->x, ebut->y, box.xmin = e.xmotion.x, box.ymin = e.xmotion.y); | 
| 845 | > | revbox(ebut->x, ebut->y, bbox.xmin = e.xmotion.x, bbox.ymin = e.xmotion.y); | 
| 846 |  | XMaskEvent(thedisplay,ButtonReleaseMask|ButtonMotionMask,&e); | 
| 847 | < | revbox(ebut->x, ebut->y, box.xmin, box.ymin); | 
| 847 | > | revbox(ebut->x, ebut->y, bbox.xmin, bbox.ymin); | 
| 848 |  | } | 
| 849 | < | box.xmin = e.xbutton.x<0 ? 0 : (e.xbutton.x>=width ? width-1 : e.xbutton.x); | 
| 850 | < | box.ymin = e.xbutton.y<0 ? 0 : (e.xbutton.y>=height ? height-1 : e.xbutton.y); | 
| 851 | < | if (box.xmin > ebut->x) { | 
| 852 | < | box.xsiz = box.xmin - ebut->x + 1; | 
| 853 | < | box.xmin = ebut->x; | 
| 849 | > | bbox.xmin = e.xbutton.x<0 ? 0 : (e.xbutton.x>=width ? width-1 : e.xbutton.x); | 
| 850 | > | bbox.ymin = e.xbutton.y<0 ? 0 : (e.xbutton.y>=height ? height-1 : e.xbutton.y); | 
| 851 | > | if (bbox.xmin > ebut->x) { | 
| 852 | > | bbox.xsiz = bbox.xmin - ebut->x + 1; | 
| 853 | > | bbox.xmin = ebut->x; | 
| 854 |  | } else { | 
| 855 | < | box.xsiz = ebut->x - box.xmin + 1; | 
| 855 | > | bbox.xsiz = ebut->x - bbox.xmin + 1; | 
| 856 |  | } | 
| 857 | < | if (box.ymin > ebut->y) { | 
| 858 | < | box.ysiz = box.ymin - ebut->y + 1; | 
| 859 | < | box.ymin = ebut->y; | 
| 857 | > | if (bbox.ymin > ebut->y) { | 
| 858 | > | bbox.ysiz = bbox.ymin - ebut->y + 1; | 
| 859 | > | bbox.ymin = ebut->y; | 
| 860 |  | } else { | 
| 861 | < | box.ysiz = ebut->y - box.ymin + 1; | 
| 861 | > | bbox.ysiz = ebut->y - bbox.ymin + 1; | 
| 862 |  | } | 
| 863 |  | } | 
| 864 |  |  | 
| 882 |  | } | 
| 883 |  |  | 
| 884 |  |  | 
| 885 | < | revbox(x0, y0, x1, y1)                  /* draw box with reversed lines */ | 
| 885 | > | revbox(x0, y0, x1, y1)                  /* draw bbox with reversed lines */ | 
| 886 |  | int  x0, y0, x1, y1; | 
| 887 |  | { | 
| 888 |  | revline(x0, y0, x1, y0); | 
| 892 |  | } | 
| 893 |  |  | 
| 894 |  |  | 
| 895 | < | avgbox(clr)                             /* average color over current box */ | 
| 896 | < | COLOR  clr; | 
| 895 | > | void | 
| 896 | > | colavg(scn, n, cavg) | 
| 897 | > | register COLR   *scn; | 
| 898 | > | register int    n; | 
| 899 | > | COLOR   cavg; | 
| 900 |  | { | 
| 901 | < | static COLOR  lc; | 
| 902 | < | static int  ll, lr, lt, lb; | 
| 901 | > | COLOR   col; | 
| 902 | > |  | 
| 903 | > | while (n--) { | 
| 904 | > | colr_color(col, *scn++); | 
| 905 | > | addcolor(cavg, col); | 
| 906 | > | } | 
| 907 | > | } | 
| 908 | > |  | 
| 909 | > |  | 
| 910 | > | int | 
| 911 | > | avgbox(cavg)                            /* average color over current bbox */ | 
| 912 | > | COLOR   cavg; | 
| 913 | > | { | 
| 914 | > | double  d; | 
| 915 | > | register int    rval; | 
| 916 | > |  | 
| 917 | > | setcolor(cavg, 0., 0., 0.); | 
| 918 | > | rval = dobox(colavg, (char *)cavg); | 
| 919 | > | if (rval > 0) { | 
| 920 | > | d = 1./rval; | 
| 921 | > | scalecolor(cavg, d); | 
| 922 | > | } | 
| 923 | > | return(rval); | 
| 924 | > | } | 
| 925 | > |  | 
| 926 | > |  | 
| 927 | > | int | 
| 928 | > | dobox(f, p)                             /* run function over bbox */ | 
| 929 | > | void    (*f)();                 /* function to call for each subscan */ | 
| 930 | > | char    *p;                     /* pointer to private data */ | 
| 931 | > | { | 
| 932 |  | int  left, right, top, bottom; | 
| 933 |  | int  y; | 
| 875 | – | double  d; | 
| 876 | – | COLOR  ctmp; | 
| 877 | – | register int  x; | 
| 934 |  |  | 
| 935 | < | setcolor(clr, 0.0, 0.0, 0.0); | 
| 936 | < | left = box.xmin - xoff; | 
| 881 | < | right = left + box.xsiz; | 
| 935 | > | left = bbox.xmin - xoff; | 
| 936 | > | right = left + bbox.xsiz; | 
| 937 |  | if (left < 0) | 
| 938 |  | left = 0; | 
| 939 |  | if (right > xmax) | 
| 940 |  | right = xmax; | 
| 941 |  | if (left >= right) | 
| 942 | < | return(-1); | 
| 943 | < | top = box.ymin - yoff; | 
| 944 | < | bottom = top + box.ysiz; | 
| 942 | > | return(0); | 
| 943 | > | top = bbox.ymin - yoff; | 
| 944 | > | bottom = top + bbox.ysiz; | 
| 945 |  | if (top < 0) | 
| 946 |  | top = 0; | 
| 947 |  | if (bottom > ymax) | 
| 948 |  | bottom = ymax; | 
| 949 |  | if (top >= bottom) | 
| 895 | – | return(-1); | 
| 896 | – | if (left == ll && right == lr && top == lt && bottom == lb) { | 
| 897 | – | copycolor(clr, lc); | 
| 950 |  | return(0); | 
| 899 | – | } | 
| 951 |  | for (y = top; y < bottom; y++) { | 
| 952 |  | if (getscan(y) == -1) | 
| 953 |  | return(-1); | 
| 954 | < | for (x = left; x < right; x++) { | 
| 955 | < | colr_color(ctmp, scanline[x]); | 
| 956 | < | addcolor(clr, ctmp); | 
| 954 | > | (*f)(scanline+left, right-left, p); | 
| 955 | > | } | 
| 956 | > | return((right-left)*(bottom-top)); | 
| 957 | > | } | 
| 958 | > |  | 
| 959 | > |  | 
| 960 | > | void | 
| 961 | > | addfix(scn, n)                  /* add fixation points to histogram */ | 
| 962 | > | COLR    *scn; | 
| 963 | > | int     n; | 
| 964 | > | { | 
| 965 | > | if (tmCvColrs(lscan, TM_NOCHROM, scn, n)) | 
| 966 | > | goto tmerr; | 
| 967 | > | if (tmAddHisto(lscan, n, FIXWEIGHT)) | 
| 968 | > | goto tmerr; | 
| 969 | > | return; | 
| 970 | > | tmerr: | 
| 971 | > | quiterr("tone mapping error"); | 
| 972 | > | } | 
| 973 | > |  | 
| 974 | > |  | 
| 975 | > | make_tonemap()                  /* initialize tone mapping */ | 
| 976 | > | { | 
| 977 | > | int  flags, y; | 
| 978 | > |  | 
| 979 | > | if (tmflags != TM_F_LINEAR && fname == NULL) { | 
| 980 | > | fprintf(stderr, "%s: cannot adjust tone of standard input\n", | 
| 981 | > | progname); | 
| 982 | > | tmflags = TM_F_LINEAR; | 
| 983 | > | } | 
| 984 | > | if (tmflags == TM_F_LINEAR) {   /* linear with clamping */ | 
| 985 | > | setcolrcor(pow, 1.0/gamcor); | 
| 986 | > | return; | 
| 987 | > | } | 
| 988 | > | flags = tmflags;                /* histogram adjustment */ | 
| 989 | > | if (greyscale) flags |= TM_F_BW; | 
| 990 | > | if (tmTop != NULL) {            /* reuse old histogram if one */ | 
| 991 | > | tmDone(tmTop); | 
| 992 | > | tmTop->flags = flags; | 
| 993 | > | } else {                        /* else initialize */ | 
| 994 | > | if ((lscan = (TMbright *)malloc(xmax*sizeof(TMbright))) == NULL) | 
| 995 | > | goto memerr; | 
| 996 | > | if (greyscale) { | 
| 997 | > | cscan = TM_NOCHROM; | 
| 998 | > | if ((pscan = (BYTE *)malloc(sizeof(BYTE)*xmax)) == NULL) | 
| 999 | > | goto memerr; | 
| 1000 | > | } else if ((pscan=cscan = (BYTE *)malloc(3*sizeof(BYTE)*xmax)) | 
| 1001 | > | == NULL) | 
| 1002 | > | goto memerr; | 
| 1003 | > | /* initialize tm library */ | 
| 1004 | > | if (tmInit(flags, stdprims, gamcor) == NULL) | 
| 1005 | > | goto memerr; | 
| 1006 | > | if (tmSetSpace(stdprims, WHTEFFICACY/exposure)) | 
| 1007 | > | goto tmerr; | 
| 1008 | > | /* compute picture histogram */ | 
| 1009 | > | for (y = 0; y < ymax; y++) { | 
| 1010 | > | getscan(y); | 
| 1011 | > | if (tmCvColrs(lscan, TM_NOCHROM, scanline, xmax)) | 
| 1012 | > | goto tmerr; | 
| 1013 | > | if (tmAddHisto(lscan, xmax, 1)) | 
| 1014 | > | goto tmerr; | 
| 1015 |  | } | 
| 1016 |  | } | 
| 1017 | < | d = 1.0/((right-left)*(bottom-top)); | 
| 1018 | < | scalecolor(clr, d); | 
| 1019 | < | ll = left; lr = right; lt = top; lb = bottom; | 
| 1020 | < | copycolor(lc, clr); | 
| 1021 | < | return(0); | 
| 1017 | > | tmDup();                        /* add fixations to duplicate map */ | 
| 1018 | > | dobox(addfix, NULL); | 
| 1019 | > | /* (re)compute tone mapping */ | 
| 1020 | > | if (tmComputeMapping(gamcor, 0., 0.)) | 
| 1021 | > | goto tmerr; | 
| 1022 | > | return; | 
| 1023 | > | memerr: | 
| 1024 | > | quiterr("out of memory in make_tonemap"); | 
| 1025 | > | tmerr: | 
| 1026 | > | quiterr("tone mapping error"); | 
| 1027 |  | } | 
| 1028 |  |  | 
| 1029 |  |  | 
| 1030 | + | tmap_colrs(scn, len)            /* apply tone mapping to scanline */ | 
| 1031 | + | register COLR  *scn; | 
| 1032 | + | int  len; | 
| 1033 | + | { | 
| 1034 | + | register BYTE  *ps; | 
| 1035 | + |  | 
| 1036 | + | if (tmflags == TM_F_LINEAR) { | 
| 1037 | + | if (scale) | 
| 1038 | + | shiftcolrs(scn, len, scale); | 
| 1039 | + | colrs_gambs(scn, len); | 
| 1040 | + | return; | 
| 1041 | + | } | 
| 1042 | + | if (len > xmax) | 
| 1043 | + | quiterr("code error 1 in tmap_colrs"); | 
| 1044 | + | if (tmCvColrs(lscan, cscan, scn, len)) | 
| 1045 | + | goto tmerr; | 
| 1046 | + | if (tmMapPixels(pscan, lscan, cscan, len)) | 
| 1047 | + | goto tmerr; | 
| 1048 | + | ps = pscan; | 
| 1049 | + | if (greyscale) | 
| 1050 | + | while (len--) { | 
| 1051 | + | scn[0][RED] = scn[0][GRN] = scn[0][BLU] = *ps++; | 
| 1052 | + | scn[0][EXP] = COLXS; | 
| 1053 | + | scn++; | 
| 1054 | + | } | 
| 1055 | + | else | 
| 1056 | + | while (len--) { | 
| 1057 | + | scn[0][RED] = *ps++; | 
| 1058 | + | scn[0][GRN] = *ps++; | 
| 1059 | + | scn[0][BLU] = *ps++; | 
| 1060 | + | scn[0][EXP] = COLXS; | 
| 1061 | + | scn++; | 
| 1062 | + | } | 
| 1063 | + | return; | 
| 1064 | + | tmerr: | 
| 1065 | + | quiterr("tone mapping error"); | 
| 1066 | + | } | 
| 1067 | + |  | 
| 1068 | + |  | 
| 1069 |  | getmono()                       /* get monochrome data */ | 
| 1070 |  | { | 
| 1071 |  | register unsigned char  *dp; | 
| 1094 |  | cerr[x] = err + errp; | 
| 1095 |  | } | 
| 1096 |  | } | 
| 1097 | < | free((char *)cerr); | 
| 1097 | > | free((void *)cerr); | 
| 1098 |  | } | 
| 1099 |  |  | 
| 1100 |  |  | 
| 1137 |  | errp = err; | 
| 1138 |  | ti = x*xmax/iconwidth; | 
| 1139 |  | copycolr(clr, scan[ti]); | 
| 1140 | < | normcolrs(clr, 1, scale); | 
| 1140 | > | normcolrs(&clr, 1, scale); | 
| 1141 |  | err += normbright(clr) + cerr[x]; | 
| 1142 |  | if (err > 127) | 
| 1143 |  | err -= 255; | 
| 1153 |  | getfull()                       /* get full (24-bit) data */ | 
| 1154 |  | { | 
| 1155 |  | int     y; | 
| 1156 | < | register unsigned int4  *dp; | 
| 1156 | > | register uint32 *dp; | 
| 1157 | > | register uint16 *dph; | 
| 1158 |  | register int    x; | 
| 1159 | < | /* set gamma correction */ | 
| 1160 | < | setcolrgam(gamcor); | 
| 1159 | > | /* initialize tone mapping */ | 
| 1160 | > | make_tonemap(); | 
| 1161 |  | /* read and convert file */ | 
| 1162 | < | dp = (unsigned int4 *)ourdata; | 
| 1162 | > | dp = (uint32 *)ourdata; | 
| 1163 | > | dph = (uint16 *)ourdata; | 
| 1164 |  | for (y = 0; y < ymax; y++) { | 
| 1165 |  | getscan(y); | 
| 1166 |  | add2icon(y, scanline); | 
| 1167 | < | if (scale) | 
| 1168 | < | shiftcolrs(scanline, xmax, scale); | 
| 1169 | < | colrs_gambs(scanline, xmax); | 
| 1015 | < | if (ourras->image->blue_mask & 1) | 
| 1167 | > | tmap_colrs(scanline, xmax); | 
| 1168 | > | switch (ourras->image->blue_mask) { | 
| 1169 | > | case 0xff:              /* 24-bit RGB */ | 
| 1170 |  | for (x = 0; x < xmax; x++) | 
| 1171 | < | *dp++ = (unsigned int4)scanline[x][RED] << 16 | | 
| 1172 | < | (unsigned int4)scanline[x][GRN] << 8 | | 
| 1173 | < | (unsigned int4)scanline[x][BLU] ; | 
| 1174 | < | else | 
| 1171 | > | *dp++ = (uint32)scanline[x][RED] << 16 | | 
| 1172 | > | (uint32)scanline[x][GRN] << 8 | | 
| 1173 | > | (uint32)scanline[x][BLU] ; | 
| 1174 | > | break; | 
| 1175 | > | case 0xff0000:          /* 24-bit BGR */ | 
| 1176 |  | for (x = 0; x < xmax; x++) | 
| 1177 | < | *dp++ = (unsigned int4)scanline[x][RED] | | 
| 1178 | < | (unsigned int4)scanline[x][GRN] << 8 | | 
| 1179 | < | (unsigned int4)scanline[x][BLU] << 16 ; | 
| 1177 | > | *dp++ = (uint32)scanline[x][RED] | | 
| 1178 | > | (uint32)scanline[x][GRN] << 8 | | 
| 1179 | > | (uint32)scanline[x][BLU] << 16 ; | 
| 1180 | > | break; | 
| 1181 | > | #if 0 | 
| 1182 | > | case 0x1f:              /* 15-bit RGB */ | 
| 1183 | > | for (x = 0; x < xmax; x++) | 
| 1184 | > | *dph++ =        (scanline[x][RED] << 7 & 0x7c00) | | 
| 1185 | > | (scanline[x][GRN] << 2 & 0x3e0) | | 
| 1186 | > | (unsigned)scanline[x][BLU] >> 3 ; | 
| 1187 | > | break; | 
| 1188 | > | case 0x7c00:            /* 15-bit BGR */ | 
| 1189 | > | for (x = 0; x < xmax; x++) | 
| 1190 | > | *dph++ =        (unsigned)scanline[x][RED] >> 3 | | 
| 1191 | > | (scanline[x][GRN] << 2 & 0x3e0) | | 
| 1192 | > | (scanline[x][BLU] << 7 & 0x7c00) ; | 
| 1193 | > | break; | 
| 1194 | > | #endif | 
| 1195 | > | default:                /* unknown */ | 
| 1196 | > | if (ourvis.depth > 16) | 
| 1197 | > | for (x = 0; x < xmax; x++) { | 
| 1198 | > | *dp = ourras->image->red_mask & | 
| 1199 | > | ourras->image->red_mask*scanline[x][RED]/255; | 
| 1200 | > | *dp |= ourras->image->green_mask & | 
| 1201 | > | ourras->image->green_mask*scanline[x][GRN]/255; | 
| 1202 | > | *dp++ |= ourras->image->blue_mask & | 
| 1203 | > | ourras->image->blue_mask*scanline[x][BLU]/255; | 
| 1204 | > | } | 
| 1205 | > | else | 
| 1206 | > | for (x = 0; x < xmax; x++) { | 
| 1207 | > | *dph = ourras->image->red_mask & | 
| 1208 | > | ourras->image->red_mask*scanline[x][RED]/255; | 
| 1209 | > | *dph |= ourras->image->green_mask & | 
| 1210 | > | ourras->image->green_mask*scanline[x][GRN]/255; | 
| 1211 | > | *dph++ |= ourras->image->blue_mask & | 
| 1212 | > | ourras->image->blue_mask*scanline[x][BLU]/255; | 
| 1213 | > | } | 
| 1214 | > | break; | 
| 1215 | > | } | 
| 1216 |  | } | 
| 1217 |  | } | 
| 1218 |  |  | 
| 1222 |  | int     y; | 
| 1223 |  | register unsigned char  *dp; | 
| 1224 |  | register int    x; | 
| 1225 | < | /* set gamma correction */ | 
| 1226 | < | setcolrgam(gamcor); | 
| 1225 | > | /* initialize tone mapping */ | 
| 1226 | > | make_tonemap(); | 
| 1227 |  | /* read and convert file */ | 
| 1228 |  | dp = ourdata; | 
| 1229 |  | for (y = 0; y < ymax; y++) { | 
| 1230 |  | getscan(y); | 
| 1231 |  | add2icon(y, scanline); | 
| 1232 | < | if (scale) | 
| 1042 | < | shiftcolrs(scanline, xmax, scale); | 
| 1043 | < | for (x = 0; x < xmax; x++) | 
| 1044 | < | scanline[x][GRN] = normbright(scanline[x]); | 
| 1045 | < | colrs_gambs(scanline, xmax); | 
| 1232 | > | tmap_colrs(scanline, xmax); | 
| 1233 |  | if (maxcolors < 256) | 
| 1234 |  | for (x = 0; x < xmax; x++) | 
| 1235 | < | *dp++ = ((long)scanline[x][GRN] * | 
| 1235 | > | *dp++ = ((int32)scanline[x][GRN] * | 
| 1236 |  | maxcolors + maxcolors/2) >> 8; | 
| 1237 |  | else | 
| 1238 |  | for (x = 0; x < xmax; x++) | 
| 1240 |  | } | 
| 1241 |  | for (x = 0; x < maxcolors; x++) | 
| 1242 |  | clrtab[x][RED] = clrtab[x][GRN] = | 
| 1243 | < | clrtab[x][BLU] = ((long)x*256 + 128)/maxcolors; | 
| 1243 | > | clrtab[x][BLU] = ((int32)x*256 + 128)/maxcolors; | 
| 1244 |  | } | 
| 1245 |  |  | 
| 1246 |  |  | 
| 1250 |  | /* make sure we can do it first */ | 
| 1251 |  | if (fname == NULL) | 
| 1252 |  | quiterr("cannot map colors from standard input"); | 
| 1253 | < | /* set gamma correction */ | 
| 1254 | < | setcolrgam(gamcor); | 
| 1253 | > | /* initialize tone mapping */ | 
| 1254 | > | make_tonemap(); | 
| 1255 |  | /* make histogram */ | 
| 1256 | < | if (new_histo((long)xmax*ymax) == -1) | 
| 1256 | > | if (new_histo((int32)xmax*ymax) == -1) | 
| 1257 |  | quiterr("cannot initialize histogram"); | 
| 1258 |  | for (y = 0; y < ymax; y++) { | 
| 1259 |  | if (getscan(y) < 0) | 
| 1260 |  | break; | 
| 1261 |  | add2icon(y, scanline); | 
| 1262 | < | if (scale) | 
| 1076 | < | shiftcolrs(scanline, xmax, scale); | 
| 1077 | < | colrs_gambs(scanline, xmax); | 
| 1262 | > | tmap_colrs(scanline, xmax); | 
| 1263 |  | cnt_colrs(scanline, xmax); | 
| 1264 |  | } | 
| 1265 |  | /* map pixels */ | 
| 1267 |  | quiterr("cannot create color map"); | 
| 1268 |  | for (y = 0; y < ymax; y++) { | 
| 1269 |  | getscan(y); | 
| 1270 | < | if (scale) | 
| 1086 | < | shiftcolrs(scanline, xmax, scale); | 
| 1087 | < | colrs_gambs(scanline, xmax); | 
| 1270 | > | tmap_colrs(scanline, xmax); | 
| 1271 |  | if (dither) | 
| 1272 |  | dith_colrs(ourdata+y*xmax, scanline, xmax); | 
| 1273 |  | else | 
| 1310 |  | static int  trunced = -1;               /* truncated file? */ | 
| 1311 |  | skipit: | 
| 1312 |  | if (trunced >= 0 && y >= trunced) { | 
| 1313 | < | bzero(scanline, xmax*sizeof(COLR)); | 
| 1313 | > | memset(scanline, '\0', xmax*sizeof(COLR)); | 
| 1314 |  | return(-1); | 
| 1315 |  | } | 
| 1316 |  | if (y != cury) { |