ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/px/x11image.c
(Generate patch)

Comparing ray/src/px/x11image.c (file contents):
Revision 2.41 by greg, Wed Jan 25 10:53:02 1995 UTC vs.
Revision 2.63 by schorsch, Mon Jul 21 22:30:19 2003 UTC

# Line 1 | Line 1
1 /* Copyright (c) 1995 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   *
# Line 21 | Line 18 | static char SCCSid[] = "$SunId$ LBL";
18  
19   #include  "standard.h"
20  
21 + #include  <string.h>
22   #include  <signal.h>
23   #include  <X11/Xlib.h>
24   #include  <X11/cursorfont.h>
# Line 28 | Line 26 | static char SCCSid[] = "$SunId$ LBL";
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)-'@')
# Line 48 | Line 39 | static char SCCSid[] = "$SunId$ LBL";
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  
# Line 76 | Line 69 | int  yoff = 0;                         /* y image offset */
69   int  parent = 0;                        /* number of children, -1 if child */
70   int  sequential = 0;                    /* display images in sequence */
71  
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  
# Line 103 | Line 104 | unsigned char  *ourdata;               /* our image data */
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  
# Line 125 | Line 126 | Atom  closedownAtom, wmProtocolsAtom;
126  
127   int  sigrecv;
128  
129 < int  onsig() { sigrecv++; }
129 > void  onsig(int i) { sigrecv++; }
130  
131  
132   main(argc, argv)
# Line 137 | Line 138 | char  *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] == '-')
145                          switch (argv[i][1]) {
146 <                        case 'c':
146 >                        case 'c':                       /* number of colors */
147                                  maxcolors = atoi(argv[++i]);
148                                  break;
149 <                        case 'b':
149 >                        case 'b':                       /* greyscale only */
150                                  greyscale = !greyscale;
151                                  break;
152 <                        case 'm':
152 >                        case 'm':                       /* monochrome */
153                                  greyscale = 1;
154                                  maxcolors = 2;
155                                  break;
156 <                        case 'd':
156 >                        case 'd':                       /* display or dither */
157                                  if (argv[i][2] == 'i')
158                                          dispname = argv[++i];
159                                  else
160                                          dither = !dither;
161                                  break;
162 <                        case 'f':
162 >                        case 'f':                       /* save pixmap */
163                                  fast = !fast;
164                                  break;
165 <                        case 's':
165 >                        case 's':                       /* one at a time */
166                                  sequential = !sequential;
167                                  break;
168 <                        case 'e':
169 <                                if (argv[i+1][0] != '+' && argv[i+1][0] != '-')
168 >                        case 'o':                       /* 't' output */
169 >                                tout = argv[i]+2;
170 >                                break;
171 >                        case 't':                       /* msec interval */
172 >                                tinterv = atoi(argv[++i]);
173 >                                break;
174 >                        case 'e':                       /* exposure comp. */
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':
188 >                        case 'g':                       /* gamma comp. */
189                                  if (argv[i][2] == 'e')
190                                          geometry = argv[++i];
191                                  else
# Line 223 | Line 240 | char  *argv[];
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] 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   {
# Line 236 | Line 254 | char  *s;
254  
255          if (isexpos(s))
256                  exposure *= exposval(s);
257 <        else if (isformat(s)) {
240 <                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  
# Line 396 | Line 414 | register XVisualInfo   *v1, *v2;
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);
# Line 475 | Line 494 | static char  vistype[][12] = {
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);
# Line 519 | Line 538 | getras()                               /* get raster file */
538                          goto fail;
539                  getmono();
540          } else if (ourvis.class == TrueColor | ourvis.class == DirectColor) {
541 <                ourdata = (unsigned char *)malloc(sizeof(int4)*xmax*ymax);
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();
# Line 535 | Line 555 | getras()                               /* get raster file */
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();
# Line 581 | Line 601 | getevent()                             /* process the next event */
601          case ButtonPress:
602                  if (xev.xbutton.state & (ShiftMask|ControlMask))
603                          moveimage(&xev.xbutton);
584                else if (xev.xbutton.button == Button2)
585                        traceray(xev.xbutton.x, xev.xbutton.y);
604                  else
605 <                        getbox(&xev.xbutton);
605 >                        switch (xev.xbutton.button) {
606 >                        case Button1:
607 >                                getbox(&xev.xbutton);
608 >                                break;
609 >                        case Button2:
610 >                                traceray(xev.xbutton.x, xev.xbutton.y);
611 >                                break;
612 >                        case Button3:
613 >                                trackrays(&xev.xbutton);
614 >                                break;
615 >                        }
616                  break;
617          case ClientMessage:
618                  if ((xev.xclient.message_type == wmProtocolsAtom) &&
# Line 595 | Line 623 | getevent()                             /* process the next event */
623   }
624  
625  
626 < traceray(xpos, ypos)                    /* print ray corresponding to pixel */
626 > traceray(xpos, ypos)                    /* print requested pixel data */
627   int  xpos, ypos;
628   {
629 <        FLOAT  hv[2];
629 >        RREAL  hv[2];
630          FVECT  rorg, rdir;
631 +        COLOR  cval;
632 +        register char  *cp;
633  
634 <        if (!gotview) {         /* no view, no can do */
635 <                XBell(thedisplay, 0);
636 <                return(-1);
637 <        }
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);
639 <        if (viewray(rorg, rdir, &ourview, hv[0], hv[1]) < 0)
640 <                return(-1);
641 <        printf("%e %e %e ", rorg[0], rorg[1], rorg[2]);
642 <        printf("%e %e %e\n", rdir[0], rdir[1], rdir[2]);
639 >        if (!gotview || viewray(rorg, rdir, &ourview, hv[0], hv[1]) < 0)
640 >                rorg[0] = rorg[1] = rorg[2] =
641 >                rdir[0] = rdir[1] = rdir[2] = 0.;
642 >
643 >        for (cp = tout; *cp; cp++)      /* print what they asked for */
644 >                switch (*cp) {
645 >                case 'o':                       /* origin */
646 >                        printf("%e %e %e ", rorg[0], rorg[1], rorg[2]);
647 >                        break;
648 >                case 'd':                       /* direction */
649 >                        printf("%e %e %e ", rdir[0], rdir[1], rdir[2]);
650 >                        break;
651 >                case 'v':                       /* radiance value */
652 >                        printf("%e %e %e ", colval(cval,RED),
653 >                                        colval(cval,GRN), colval(cval,BLU));
654 >                        break;
655 >                case 'l':                       /* luminance */
656 >                        printf("%e ", luminance(cval));
657 >                        break;
658 >                case 'p':                       /* pixel position */
659 >                        printf("%d %d ", (int)(hv[0]*inpres.xr),
660 >                                        (int)(hv[1]*inpres.yr));
661 >                        break;
662 >                }
663 >        putchar('\n');
664          fflush(stdout);
665          return(0);
666   }
# Line 623 | Line 674 | XKeyPressedEvent  *ekey;
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)
# Line 638 | Line 689 | XKeyPressedEvent  *ekey;
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':
# Line 646 | Line 697 | XKeyPressedEvent  *ekey;
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);
# Line 657 | Line 708 | XKeyPressedEvent  *ekey;
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)
# Line 680 | Line 731 | XKeyPressedEvent  *ekey;
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) {
# Line 696 | Line 759 | XKeyPressedEvent  *ekey;
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 */
# Line 733 | Line 801 | XKeyPressedEvent  *ekey;
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);
# Line 767 | Line 835 | XButtonPressedEvent  *ebut;
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  
865  
866 < revbox(x0, y0, x1, y1)                  /* draw box with reversed lines */
866 > trackrays(ebut)                         /* trace rays as mouse moves */
867 > XButtonPressedEvent  *ebut;
868 > {
869 >        XEvent  e;
870 >        unsigned long   lastrept;
871 >
872 >        traceray(ebut->x, ebut->y);
873 >        lastrept = ebut->time;
874 >        XMaskEvent(thedisplay, ButtonReleaseMask|ButtonMotionMask, &e);
875 >        while (e.type == MotionNotify) {
876 >                if (e.xmotion.time >= lastrept + tinterv) {
877 >                        traceray(e.xmotion.x, e.xmotion.y);
878 >                        lastrept = e.xmotion.time;
879 >                }
880 >                XMaskEvent(thedisplay,ButtonReleaseMask|ButtonMotionMask,&e);
881 >        }
882 > }
883 >
884 >
885 > revbox(x0, y0, x1, y1)                  /* draw bbox with reversed lines */
886   int  x0, y0, x1, y1;
887   {
888          revline(x0, y0, x1, y0);
# Line 805 | Line 892 | int  x0, y0, x1, y1;
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;
815        double  d;
816        COLOR  ctmp;
817        register int  x;
934  
935 <        setcolor(clr, 0.0, 0.0, 0.0);
936 <        left = box.xmin - xoff;
821 <        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)
835                return(-1);
836        if (left == ll && right == lr && top == lt && bottom == lb) {
837                copycolor(clr, lc);
950                  return(0);
839        }
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;
# Line 881 | Line 1094 | getmono()                      /* get monochrome data */
1094                          cerr[x] = err + errp;
1095                  }
1096          }
1097 <        free((char *)cerr);
1097 >        free((void *)cerr);
1098   }
1099  
1100  
# Line 924 | Line 1137 | COLR  *scan;
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;
# Line 940 | Line 1153 | COLR  *scan;
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);
955 <                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  
# Line 971 | Line 1222 | getgrey()                      /* get greyscale data */
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)
982 <                        shiftcolrs(scanline, xmax, scale);
983 <                for (x = 0; x < xmax; x++)
984 <                        scanline[x][GRN] = normbright(scanline[x]);
985 <                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++)
# Line 993 | Line 1240 | getgrey()                      /* get greyscale data */
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  
# Line 1003 | Line 1250 | getmapped()                    /* get color-mapped data */
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)
1016 <                        shiftcolrs(scanline, xmax, scale);
1017 <                colrs_gambs(scanline, xmax);
1262 >                tmap_colrs(scanline, xmax);
1263                  cnt_colrs(scanline, xmax);
1264          }
1265                                          /* map pixels */
# Line 1022 | Line 1267 | getmapped()                    /* get color-mapped data */
1267                  quiterr("cannot create color map");
1268          for (y = 0; y < ymax; y++) {
1269                  getscan(y);
1270 <                if (scale)
1026 <                        shiftcolrs(scanline, xmax, scale);
1027 <                colrs_gambs(scanline, xmax);
1270 >                tmap_colrs(scanline, xmax);
1271                  if (dither)
1272                          dith_colrs(ourdata+y*xmax, scanline, xmax);
1273                  else
# Line 1067 | Line 1310 | int  y;
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) {

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines