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.35 by greg, Mon Nov 22 11:33:18 1993 UTC vs.
Revision 2.57 by greg, Mon Mar 10 17:26:26 2003 UTC

# Line 1 | Line 1
1 – /* Copyright (c) 1993 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 28 | Line 25 | static char SCCSid[] = "$SunId$ LBL";
25   #include  <X11/Xatom.h>
26  
27   #include  "color.h"
28 + #include  "tonemap.h"
29   #include  "view.h"
30   #include  "x11raster.h"
31   #include  "random.h"
34 – #include  "resolu.h"
32  
36 – #ifdef  __alpha
37 – #define  int4           int
38 – #endif
39 – #ifndef  int4
40 – #define  int4           long
41 – #endif
42 –
33   #define  FONTNAME       "8x13"          /* text font we'll use */
34  
35   #define  CTRL(c)        ((c)-'@')
# Line 48 | Line 38 | static char SCCSid[] = "$SunId$ LBL";
38  
39   #define  ICONSIZ        (8*10)          /* maximum icon dimension (even 8) */
40  
41 + #define  FIXWEIGHT      20              /* weight to add for fixation points */
42 +
43   #define  ourscreen      DefaultScreen(thedisplay)
44   #define  ourroot        RootWindow(thedisplay,ourscreen)
45  
# Line 74 | Line 66 | int  xoff = 0;                         /* x image offset */
66   int  yoff = 0;                          /* y image offset */
67  
68   int  parent = 0;                        /* number of children, -1 if child */
69 + int  sequential = 0;                    /* display images in sequence */
70  
71 + char  *tout = "od";                     /* output of 't' command */
72 + int  tinterv = 0;                       /* interval between mouse reports */
73 +
74 + int  tmflags = TM_F_LINEAR;             /* tone mapping flags */
75 +
76   VIEW  ourview = STDVIEW;                /* image view parameters */
77   int  gotview = 0;                       /* got parameters from file */
78  
79   COLR  *scanline;                        /* scan line buffer */
80 + TMbright  *lscan;                       /* encoded luminance scanline */
81 + BYTE  *cscan;                           /* encoded chroma scanline */
82 + BYTE  *pscan;                           /* compute pixel scanline */
83  
84   RESOLU  inpres;                         /* input resolution and ordering */
85   int  xmax, ymax;                        /* picture dimensions */
86   int  width, height;                     /* window size */
87   char  *fname = NULL;                    /* input file name */
88 < FILE  *fin = stdin;                     /* input file */
88 > FILE  *fin = NULL;                      /* input file */
89   long  *scanpos = NULL;                  /* scan line positions in file */
90   int  cury = 0;                          /* current scan location */
91  
# Line 102 | Line 103 | unsigned char  *ourdata;               /* our image data */
103  
104   struct {
105          int  xmin, ymin, xsiz, ysiz;
106 < }  box = {0, 0, 0, 0};                  /* current box */
106 > }  bbox = {0, 0, 0, 0};                 /* current bbox */
107  
108   char  *geometry = NULL;                 /* geometry specification */
109  
# Line 113 | Line 114 | char  *progname;
114  
115   char  errmsg[128];
116  
117 < extern BYTE  clrtab[256][3];            /* global color map */
117 > BYTE  clrtab[256][3];                   /* global color map */
118  
119   extern long  ftell();
120  
# Line 136 | Line 137 | char  *argv[];
137          int  pid;
138          
139          progname = argv[0];
140 +        fin = stdin;
141  
142          for (i = 1; i < argc; i++)
143                  if (argv[i][0] == '-')
144                          switch (argv[i][1]) {
145 <                        case 'c':
145 >                        case 'c':                       /* number of colors */
146                                  maxcolors = atoi(argv[++i]);
147                                  break;
148 <                        case 'b':
148 >                        case 'b':                       /* greyscale only */
149                                  greyscale = !greyscale;
150                                  break;
151 <                        case 'm':
151 >                        case 'm':                       /* monochrome */
152 >                                greyscale = 1;
153                                  maxcolors = 2;
154                                  break;
155 <                        case 'd':
155 >                        case 'd':                       /* display or dither */
156                                  if (argv[i][2] == 'i')
157                                          dispname = argv[++i];
158                                  else
159                                          dither = !dither;
160                                  break;
161 <                        case 'f':
161 >                        case 'f':                       /* save pixmap */
162                                  fast = !fast;
163                                  break;
164 <                        case 'e':
165 <                                if (argv[i+1][0] != '+' && argv[i+1][0] != '-')
164 >                        case 's':                       /* one at a time */
165 >                                sequential = !sequential;
166 >                                break;
167 >                        case 'o':                       /* 't' output */
168 >                                tout = argv[i]+2;
169 >                                break;
170 >                        case 't':                       /* msec interval */
171 >                                tinterv = atoi(argv[++i]);
172 >                                break;
173 >                        case 'e':                       /* exposure comp. */
174 >                                i++;
175 >                                if (argv[i][0] == 'a') {
176 >                                        tmflags = TM_F_CAMERA;
177 >                                        break;
178 >                                }
179 >                                if (argv[i][0] == 'h') {
180 >                                        tmflags = TM_F_HUMAN;
181 >                                        break;
182 >                                }
183 >                                if (argv[i][0] != '+' && argv[i][0] != '-')
184                                          goto userr;
185 <                                scale = atoi(argv[++i]);
185 >                                scale = atoi(argv[i]);
186                                  break;
187 <                        case 'g':
187 >                        case 'g':                       /* gamma comp. */
188                                  if (argv[i][2] == 'e')
189                                          geometry = argv[++i];
190                                  else
# Line 214 | Line 235 | char  *argv[];
235  
236          init(argc, argv);                       /* get file and open window */
237  
217 –        if (parent < 0) {
218 –                kill(getppid(), SIGCONT);       /* signal parent if child */
219 –                sigrecv--;
220 –        }
238          for ( ; ; )
239                  getevent();             /* main loop */
240   userr:
241          fprintf(stderr,
242 < "Usage: %s [-di disp][[-ge] spec][-b][-m][-d][-f][-c nclrs][-e +/-stops][-g gamcor] pic ..\n",
242 > "Usage: %s [-di disp][[-ge] spec][-b][-m][-d][-f][-c nclrs][-e spec][-g gamcor][-s][-ospec][-t intvl] pic ..\n",
243                          progname);
244          exit(1);
245   }
246  
247  
248 + int
249   headline(s)             /* get relevant info from header */
250   char  *s;
251   {
# Line 235 | Line 253 | char  *s;
253  
254          if (isexpos(s))
255                  exposure *= exposval(s);
256 <        else if (isformat(s)) {
239 <                formatval(fmt, s);
256 >        else if (formatval(fmt, s))
257                  wrongformat = strcmp(fmt, COLRFMT);
258 <        } else if (isview(s) && sscanview(&ourview, s) > 0)
258 >        else if (isview(s) && sscanview(&ourview, s) > 0)
259                  gotview++;
260 +        return(0);
261   }
262  
263  
# Line 276 | Line 294 | char **argv;
294          if (gamstr == NULL)             /* get it from the X server */
295                  gamstr = XGetDefault(thedisplay, "radiance", "gamma");
296          if (gamstr == NULL)             /* get it from the environment */
297 <                gamstr = getenv("GAMMA");
297 >                gamstr = getenv("DISPLAY_GAMMA");
298          if (gamstr != NULL)
299                  gamcor = atof(gamstr);
300                                  /* get best visual for default screen */
# Line 305 | Line 323 | char **argv;
323                  }
324          }
325          /* open window */
326 +        i = CWEventMask|CWCursor|CWBackPixel|CWBorderPixel;
327          ourwinattr.border_pixel = ourwhite;
328          ourwinattr.background_pixel = ourblack;
329 <        ourwinattr.colormap = XCreateColormap(thedisplay, ourroot,
330 <                        ourvis.visual, AllocNone);
329 >        if (ourvis.visual != DefaultVisual(thedisplay,ourscreen)) {
330 >                ourwinattr.colormap = newcmap(thedisplay, ourscreen, ourvis.visual);
331 >                i |= CWColormap;
332 >        }
333          ourwinattr.event_mask = ExposureMask|KeyPressMask|ButtonPressMask|
334                          ButtonReleaseMask|ButtonMotionMask|StructureNotifyMask;
335          ourwinattr.cursor = XCreateFontCursor(thedisplay, XC_diamond_cross);
336          wind = XCreateWindow(thedisplay, ourroot, xszhints.x, xszhints.y,
337                          xszhints.width, xszhints.height, BORWIDTH,
338 <                        ourvis.depth, InputOutput, ourvis.visual, CWEventMask|
339 <                        CWCursor|CWBackPixel|CWBorderPixel|CWColormap, &ourwinattr);
338 >                        ourvis.depth, InputOutput, ourvis.visual,
339 >                        i, &ourwinattr);
340          if (wind == 0)
341                  quiterr("cannot create window");
342          width = xmax;
# Line 353 | Line 374 | char **argv;
374          XSetWMProtocols(thedisplay, wind, &closedownAtom, 1);
375  
376          XMapWindow(thedisplay, wind);
356 –        return;
377   } /* end of init */
378  
379  
# Line 366 | Line 386 | char  *err;
386          if (es = err != NULL)
387                  fprintf(stderr, "%s: %s: %s\n", progname,
388                                  fname==NULL?"<stdin>":fname, err);
389 <        if (parent > 0 & wind != 0) {
390 <                XDestroyWindow(thedisplay, wind);
391 <                XFlush(thedisplay);
372 <        } else if (parent < 0 & sigrecv == 0)
389 >        if (thedisplay != NULL)
390 >                XCloseDisplay(thedisplay);
391 >        if (parent < 0 & sigrecv == 0)
392                  kill(getppid(), SIGCONT);
393          while (parent > 0 && wait(&cs) != -1) { /* wait for any children */
394                  if (es == 0)
# Line 394 | Line 413 | register XVisualInfo   *v1, *v2;
413                                  return(-1);
414                          if (v1->depth == 32 && v2->depth == 24)
415                                  return(1);
416 <                        return(0);
416 >                                        /* go for maximum depth otherwise */
417 >                        return(v2->depth - v1->depth);
418                  }
419                                          /* don't be too greedy */
420                  if (maxcolors <= 1<<v1->depth && maxcolors <= 1<<v2->depth)
421                          return(v1->depth - v2->depth);
422                  return(v2->depth - v1->depth);
423          }
424 <                                        /* prefer Pseudo when < 24-bit */
424 >                                        /* prefer Pseudo when < 15-bit */
425          if ((v1->class == TrueColor || v1->class == DirectColor) &&
426 <                        v1->depth < 24)
426 >                        v1->depth < 15)
427                  bad1 = 1;
428          if ((v2->class == TrueColor || v2->class == DirectColor) &&
429 <                        v2->depth < 24)
429 >                        v2->depth < 15)
430                  bad2 = -1;
431          if (bad1 | bad2)
432                  return(bad1+bad2);
# Line 517 | Line 537 | getras()                               /* get raster file */
537                          goto fail;
538                  getmono();
539          } else if (ourvis.class == TrueColor | ourvis.class == DirectColor) {
540 <                ourdata = (unsigned char *)malloc(sizeof(int4)*xmax*ymax);
540 >                int  datsiz = ourvis.depth>16 ? sizeof(int4) : sizeof(int2);
541 >                ourdata = (unsigned char *)malloc(datsiz*xmax*ymax);
542                  if (ourdata == NULL)
543                          goto fail;
544 <                ourras = make_raster(thedisplay, &ourvis, sizeof(int4)*8,
545 <                                ourdata, xmax, ymax, 32);
544 >                ourras = make_raster(thedisplay, &ourvis, datsiz*8,
545 >                                ourdata, xmax, ymax, datsiz*8);
546                  if (ourras == NULL)
547                          goto fail;
548                  getfull();
# Line 533 | Line 554 | getras()                               /* get raster file */
554                                  xmax, ymax, 8);
555                  if (ourras == NULL)
556                          goto fail;
557 <                if (greyscale | ourvis.class == StaticGray)
557 >                if (greyscale)
558                          getgrey();
559                  else
560                          getmapped();
# Line 563 | Line 584 | getevent()                             /* process the next event */
584                  map_rcolors(ourras, wind);
585                  if (fast)
586                          make_rpixmap(ourras, wind);
587 +                if (!sequential & parent < 0 & sigrecv == 0) {
588 +                        kill(getppid(), SIGCONT);
589 +                        sigrecv--;
590 +                }
591                  break;
592          case UnmapNotify:
593                  if (!fast)
# Line 575 | Line 600 | getevent()                             /* process the next event */
600          case ButtonPress:
601                  if (xev.xbutton.state & (ShiftMask|ControlMask))
602                          moveimage(&xev.xbutton);
578 –                else if (xev.xbutton.button == Button2)
579 –                        traceray(xev.xbutton.x, xev.xbutton.y);
603                  else
604 <                        getbox(&xev.xbutton);
604 >                        switch (xev.xbutton.button) {
605 >                        case Button1:
606 >                                getbox(&xev.xbutton);
607 >                                break;
608 >                        case Button2:
609 >                                traceray(xev.xbutton.x, xev.xbutton.y);
610 >                                break;
611 >                        case Button3:
612 >                                trackrays(&xev.xbutton);
613 >                                break;
614 >                        }
615                  break;
616          case ClientMessage:
617                  if ((xev.xclient.message_type == wmProtocolsAtom) &&
# Line 589 | Line 622 | getevent()                             /* process the next event */
622   }
623  
624  
625 < traceray(xpos, ypos)                    /* print ray corresponding to pixel */
625 > traceray(xpos, ypos)                    /* print requested pixel data */
626   int  xpos, ypos;
627   {
628          FLOAT  hv[2];
629          FVECT  rorg, rdir;
630 +        COLOR  cval;
631 +        register char  *cp;
632  
633 <        if (!gotview) {         /* no view, no can do */
634 <                XBell(thedisplay, 0);
635 <                return(-1);
636 <        }
633 >        bbox.xmin = xpos; bbox.xsiz = 1;
634 >        bbox.ymin = ypos; bbox.ysiz = 1;
635 >        avgbox(cval);
636 >        scalecolor(cval, 1./exposure);
637          pix2loc(hv, &inpres, xpos-xoff, ypos-yoff);
638 <        if (viewray(rorg, rdir, &ourview, hv[0], hv[1]) < 0)
639 <                return(-1);
640 <        printf("%e %e %e ", rorg[0], rorg[1], rorg[2]);
641 <        printf("%e %e %e\n", rdir[0], rdir[1], rdir[2]);
638 >        if (!gotview || viewray(rorg, rdir, &ourview, hv[0], hv[1]) < 0)
639 >                rorg[0] = rorg[1] = rorg[2] =
640 >                rdir[0] = rdir[1] = rdir[2] = 0.;
641 >
642 >        for (cp = tout; *cp; cp++)      /* print what they asked for */
643 >                switch (*cp) {
644 >                case 'o':                       /* origin */
645 >                        printf("%e %e %e ", rorg[0], rorg[1], rorg[2]);
646 >                        break;
647 >                case 'd':                       /* direction */
648 >                        printf("%e %e %e ", rdir[0], rdir[1], rdir[2]);
649 >                        break;
650 >                case 'v':                       /* radiance value */
651 >                        printf("%e %e %e ", colval(cval,RED),
652 >                                        colval(cval,GRN), colval(cval,BLU));
653 >                        break;
654 >                case 'l':                       /* luminance */
655 >                        printf("%e ", luminance(cval));
656 >                        break;
657 >                case 'p':                       /* pixel position */
658 >                        printf("%d %d ", (int)(hv[0]*inpres.xr),
659 >                                        (int)(hv[1]*inpres.yr));
660 >                        break;
661 >                }
662 >        putchar('\n');
663          fflush(stdout);
664          return(0);
665   }
# Line 632 | Line 688 | XKeyPressedEvent  *ekey;
688          case '\r':
689          case 'l':
690          case 'c':                               /* value */
691 <                if (avgbox(cval) == -1)
691 >                if (!avgbox(cval))
692                          return(-1);
693                  switch (com) {
694                  case '\n':
# Line 640 | Line 696 | XKeyPressedEvent  *ekey;
696                          sprintf(buf, "%.3f", intens(cval)/exposure);
697                          break;
698                  case 'l':                               /* luminance */
699 <                        sprintf(buf, "%.0fL", luminance(cval)/exposure);
699 >                        sprintf(buf, "%.1fL", luminance(cval)/exposure);
700                          break;
701                  case 'c':                               /* color */
702                          comp = pow(2.0, (double)scale);
# Line 651 | Line 707 | XKeyPressedEvent  *ekey;
707                          break;
708                  }
709                  XDrawImageString(thedisplay, wind, ourgc,
710 <                                box.xmin, box.ymin+box.ysiz, buf, strlen(buf));
710 >                                bbox.xmin, bbox.ymin+bbox.ysiz, buf, strlen(buf));
711                  return(0);
712          case 'i':                               /* identify (contour) */
713                  if (ourras->pixels == NULL)
# Line 674 | Line 730 | XKeyPressedEvent  *ekey;
730                  return(0);
731          case 't':                               /* trace */
732                  return(traceray(ekey->x, ekey->y));
733 +        case 'a':                               /* auto exposure */
734 +                if (fname == NULL)
735 +                        return(-1);
736 +                tmflags = TM_F_CAMERA;
737 +                strcpy(buf, "auto exposure...");
738 +                goto remap;
739 +        case 'h':                               /* human response */
740 +                if (fname == NULL)
741 +                        return(-1);
742 +                tmflags = TM_F_HUMAN;
743 +                strcpy(buf, "human exposure...");
744 +                goto remap;
745          case '=':                               /* adjust exposure */
746          case '@':                               /* adaptation level */
747 <                if (avgbox(cval) == -1)
747 >                if (!avgbox(cval))
748                          return(-1);
749                  comp = bright(cval);
750                  if (comp < 1e-20) {
# Line 690 | Line 758 | XKeyPressedEvent  *ekey;
758                          comp = .5/comp;
759                  comp = log(comp)/.69315 - scale;
760                  n = comp < 0 ? comp-.5 : comp+.5 ;      /* round */
761 <                if (n == 0)
762 <                        return(0);
763 <                scale_rcolors(ourras, pow(2.0, (double)n));
761 >                if (tmflags != TM_F_LINEAR)
762 >                        tmflags = TM_F_LINEAR;  /* turn off tone mapping */
763 >                else {
764 >                        if (n == 0)             /* else check if any change */
765 >                                return(0);
766 >                        scale_rcolors(ourras, pow(2.0, (double)n));
767 >                }
768                  scale += n;
769                  sprintf(buf, "%+d", scale);
770 +        remap:
771                  XDrawImageString(thedisplay, wind, ourgc,
772 <                                box.xmin, box.ymin+box.ysiz, buf, strlen(buf));
772 >                                bbox.xmin, bbox.ymin+bbox.ysiz, buf, strlen(buf));
773                  XFlush(thedisplay);
774 <                free(ourdata);
774 >                /* free(ourdata);       This is done in XDestroyImage()! */
775                  free_raster(ourras);
776                  getras();
777          /* fall through */
# Line 727 | Line 800 | XKeyPressedEvent  *ekey;
800                  redraw(0, 0, width, height);
801                  return(0);
802          case ' ':                               /* clear */
803 <                redraw(box.xmin, box.ymin, box.xsiz, box.ysiz);
803 >                redraw(bbox.xmin, bbox.ymin, bbox.xsiz, bbox.ysiz);
804                  return(0);
805          default:
806                  XBell(thedisplay, 0);
# Line 761 | Line 834 | XButtonPressedEvent  *ebut;
834   }
835  
836  
837 < getbox(ebut)                            /* get new box */
837 > getbox(ebut)                            /* get new bbox */
838   XButtonPressedEvent  *ebut;
839   {
840          XEvent  e;
841  
842          XMaskEvent(thedisplay, ButtonReleaseMask|ButtonMotionMask, &e);
843          while (e.type == MotionNotify) {
844 <                revbox(ebut->x, ebut->y, box.xmin = e.xmotion.x, box.ymin = e.xmotion.y);
844 >                revbox(ebut->x, ebut->y, bbox.xmin = e.xmotion.x, bbox.ymin = e.xmotion.y);
845                  XMaskEvent(thedisplay,ButtonReleaseMask|ButtonMotionMask,&e);
846 <                revbox(ebut->x, ebut->y, box.xmin, box.ymin);
846 >                revbox(ebut->x, ebut->y, bbox.xmin, bbox.ymin);
847          }
848 <        box.xmin = e.xbutton.x<0 ? 0 : (e.xbutton.x>=width ? width-1 : e.xbutton.x);
849 <        box.ymin = e.xbutton.y<0 ? 0 : (e.xbutton.y>=height ? height-1 : e.xbutton.y);
850 <        if (box.xmin > ebut->x) {
851 <                box.xsiz = box.xmin - ebut->x + 1;
852 <                box.xmin = ebut->x;
848 >        bbox.xmin = e.xbutton.x<0 ? 0 : (e.xbutton.x>=width ? width-1 : e.xbutton.x);
849 >        bbox.ymin = e.xbutton.y<0 ? 0 : (e.xbutton.y>=height ? height-1 : e.xbutton.y);
850 >        if (bbox.xmin > ebut->x) {
851 >                bbox.xsiz = bbox.xmin - ebut->x + 1;
852 >                bbox.xmin = ebut->x;
853          } else {
854 <                box.xsiz = ebut->x - box.xmin + 1;
854 >                bbox.xsiz = ebut->x - bbox.xmin + 1;
855          }
856 <        if (box.ymin > ebut->y) {
857 <                box.ysiz = box.ymin - ebut->y + 1;
858 <                box.ymin = ebut->y;
856 >        if (bbox.ymin > ebut->y) {
857 >                bbox.ysiz = bbox.ymin - ebut->y + 1;
858 >                bbox.ymin = ebut->y;
859          } else {
860 <                box.ysiz = ebut->y - box.ymin + 1;
860 >                bbox.ysiz = ebut->y - bbox.ymin + 1;
861          }
862   }
863  
864  
865 < revbox(x0, y0, x1, y1)                  /* draw box with reversed lines */
865 > trackrays(ebut)                         /* trace rays as mouse moves */
866 > XButtonPressedEvent  *ebut;
867 > {
868 >        XEvent  e;
869 >        unsigned long   lastrept;
870 >
871 >        traceray(ebut->x, ebut->y);
872 >        lastrept = ebut->time;
873 >        XMaskEvent(thedisplay, ButtonReleaseMask|ButtonMotionMask, &e);
874 >        while (e.type == MotionNotify) {
875 >                if (e.xmotion.time >= lastrept + tinterv) {
876 >                        traceray(e.xmotion.x, e.xmotion.y);
877 >                        lastrept = e.xmotion.time;
878 >                }
879 >                XMaskEvent(thedisplay,ButtonReleaseMask|ButtonMotionMask,&e);
880 >        }
881 > }
882 >
883 >
884 > revbox(x0, y0, x1, y1)                  /* draw bbox with reversed lines */
885   int  x0, y0, x1, y1;
886   {
887          revline(x0, y0, x1, y0);
# Line 799 | Line 891 | int  x0, y0, x1, y1;
891   }
892  
893  
894 < avgbox(clr)                             /* average color over current box */
895 < COLOR  clr;
894 > int
895 > colavg(scn, n, cavg)
896 > register COLR   *scn;
897 > register int    n;
898 > COLOR   cavg;
899   {
900 <        static COLOR  lc;
901 <        static int  ll, lr, lt, lb;
900 >        COLOR   col;
901 >
902 >        while (n--) {
903 >                colr_color(col, *scn++);
904 >                addcolor(cavg, col);
905 >        }
906 > }
907 >
908 >
909 > int
910 > avgbox(cavg)                            /* average color over current bbox */
911 > COLOR   cavg;
912 > {
913 >        double  d;
914 >        register int    rval;
915 >
916 >        setcolor(cavg, 0., 0., 0.);
917 >        rval = dobox(colavg, (char *)cavg);
918 >        if (rval > 0) {
919 >                d = 1./rval;
920 >                scalecolor(cavg, d);
921 >        }
922 >        return(rval);
923 > }
924 >
925 >
926 > int
927 > dobox(f, p)                             /* run function over bbox */
928 > int     (*f)();                 /* function to call for each subscan */
929 > char    *p;                     /* pointer to private data */
930 > {
931          int  left, right, top, bottom;
932          int  y;
809 –        double  d;
810 –        COLOR  ctmp;
811 –        register int  x;
933  
934 <        setcolor(clr, 0.0, 0.0, 0.0);
935 <        left = box.xmin - xoff;
815 <        right = left + box.xsiz;
934 >        left = bbox.xmin - xoff;
935 >        right = left + bbox.xsiz;
936          if (left < 0)
937                  left = 0;
938          if (right > xmax)
939                  right = xmax;
940          if (left >= right)
941 <                return(-1);
942 <        top = box.ymin - yoff;
943 <        bottom = top + box.ysiz;
941 >                return(0);
942 >        top = bbox.ymin - yoff;
943 >        bottom = top + bbox.ysiz;
944          if (top < 0)
945                  top = 0;
946          if (bottom > ymax)
947                  bottom = ymax;
948          if (top >= bottom)
829 –                return(-1);
830 –        if (left == ll && right == lr && top == lt && bottom == lb) {
831 –                copycolor(clr, lc);
949                  return(0);
833 –        }
950          for (y = top; y < bottom; y++) {
951                  if (getscan(y) == -1)
952                          return(-1);
953 <                for (x = left; x < right; x++) {
954 <                        colr_color(ctmp, scanline[x]);
955 <                        addcolor(clr, ctmp);
953 >                (*f)(scanline+left, right-left, p);
954 >        }
955 >        return((right-left)*(bottom-top));
956 > }
957 >
958 >
959 > int
960 > addfix(scn, n)                  /* add fixation points to histogram */
961 > COLR    *scn;
962 > int     n;
963 > {
964 >        if (tmCvColrs(lscan, TM_NOCHROM, scn, n))
965 >                goto tmerr;
966 >        if (tmAddHisto(lscan, n, FIXWEIGHT))
967 >                goto tmerr;
968 >        return;
969 > tmerr:
970 >        quiterr("tone mapping error");
971 > }
972 >
973 >
974 > make_tonemap()                  /* initialize tone mapping */
975 > {
976 >        int  flags, y;
977 >
978 >        if (tmflags != TM_F_LINEAR && fname == NULL) {
979 >                fprintf(stderr, "%s: cannot adjust tone of standard input\n",
980 >                                progname);
981 >                tmflags = TM_F_LINEAR;
982 >        }
983 >        if (tmflags == TM_F_LINEAR) {   /* linear with clamping */
984 >                setcolrcor(pow, 1.0/gamcor);
985 >                return;
986 >        }
987 >        flags = tmflags;                /* histogram adjustment */
988 >        if (greyscale) flags |= TM_F_BW;
989 >        if (tmTop != NULL) {            /* reuse old histogram if one */
990 >                tmDone(tmTop);
991 >                tmTop->flags = flags;
992 >        } else {                        /* else initialize */
993 >                if ((lscan = (TMbright *)malloc(xmax*sizeof(TMbright))) == NULL)
994 >                        goto memerr;
995 >                if (greyscale) {
996 >                        cscan = TM_NOCHROM;
997 >                        if ((pscan = (BYTE *)malloc(sizeof(BYTE)*xmax)) == NULL)
998 >                                goto memerr;
999 >                } else if ((pscan=cscan = (BYTE *)malloc(3*sizeof(BYTE)*xmax))
1000 >                                == NULL)
1001 >                        goto memerr;
1002 >                                                /* initialize tm library */
1003 >                if (tmInit(flags, stdprims, gamcor) == NULL)
1004 >                        goto memerr;
1005 >                if (tmSetSpace(stdprims, WHTEFFICACY/exposure))
1006 >                        goto tmerr;
1007 >                                                /* compute picture histogram */
1008 >                for (y = 0; y < ymax; y++) {
1009 >                        getscan(y);
1010 >                        if (tmCvColrs(lscan, TM_NOCHROM, scanline, xmax))
1011 >                                goto tmerr;
1012 >                        if (tmAddHisto(lscan, xmax, 1))
1013 >                                goto tmerr;
1014                  }
1015          }
1016 <        d = 1.0/((right-left)*(bottom-top));
1017 <        scalecolor(clr, d);
1018 <        ll = left; lr = right; lt = top; lb = bottom;
1019 <        copycolor(lc, clr);
1020 <        return(0);
1016 >        tmDup();                        /* add fixations to duplicate map */
1017 >        dobox(addfix, NULL);
1018 >                                        /* (re)compute tone mapping */
1019 >        if (tmComputeMapping(gamcor, 0., 0.))
1020 >                goto tmerr;
1021 >        return;
1022 > memerr:
1023 >        quiterr("out of memory in make_tonemap");
1024 > tmerr:
1025 >        quiterr("tone mapping error");
1026   }
1027  
1028  
1029 + tmap_colrs(scn, len)            /* apply tone mapping to scanline */
1030 + register COLR  *scn;
1031 + int  len;
1032 + {
1033 +        register BYTE  *ps;
1034 +
1035 +        if (tmflags == TM_F_LINEAR) {
1036 +                if (scale)
1037 +                        shiftcolrs(scn, len, scale);
1038 +                colrs_gambs(scn, len);
1039 +                return;
1040 +        }
1041 +        if (len > xmax)
1042 +                quiterr("code error 1 in tmap_colrs");
1043 +        if (tmCvColrs(lscan, cscan, scn, len))
1044 +                goto tmerr;
1045 +        if (tmMapPixels(pscan, lscan, cscan, len))
1046 +                goto tmerr;
1047 +        ps = pscan;
1048 +        if (greyscale)
1049 +                while (len--) {
1050 +                        scn[0][RED] = scn[0][GRN] = scn[0][BLU] = *ps++;
1051 +                        scn[0][EXP] = COLXS;
1052 +                        scn++;
1053 +                }
1054 +        else
1055 +                while (len--) {
1056 +                        scn[0][RED] = *ps++;
1057 +                        scn[0][GRN] = *ps++;
1058 +                        scn[0][BLU] = *ps++;
1059 +                        scn[0][EXP] = COLXS;
1060 +                        scn++;
1061 +                }
1062 +        return;
1063 + tmerr:
1064 +        quiterr("tone mapping error");
1065 + }
1066 +
1067 +
1068   getmono()                       /* get monochrome data */
1069   {
1070          register unsigned char  *dp;
# Line 875 | Line 1093 | getmono()                      /* get monochrome data */
1093                          cerr[x] = err + errp;
1094                  }
1095          }
1096 <        free((char *)cerr);
1096 >        free((void *)cerr);
1097   }
1098  
1099  
# Line 918 | Line 1136 | COLR  *scan;
1136                  errp = err;
1137                  ti = x*xmax/iconwidth;
1138                  copycolr(clr, scan[ti]);
1139 <                normcolrs(clr, 1, scale);
1139 >                normcolrs(&clr, 1, scale);
1140                  err += normbright(clr) + cerr[x];
1141                  if (err > 127)
1142                          err -= 255;
# Line 935 | Line 1153 | getfull()                      /* get full (24-bit) data */
1153   {
1154          int     y;
1155          register unsigned int4  *dp;
1156 +        register unsigned int2  *dph;
1157          register int    x;
1158 <                                        /* set gamma correction */
1159 <        setcolrgam(gamcor);
1158 >                                        /* initialize tone mapping */
1159 >        make_tonemap();
1160                                          /* read and convert file */
1161          dp = (unsigned int4 *)ourdata;
1162 +        dph = (unsigned int2 *)ourdata;
1163          for (y = 0; y < ymax; y++) {
1164                  getscan(y);
1165                  add2icon(y, scanline);
1166 <                if (scale)
1167 <                        shiftcolrs(scanline, xmax, scale);
1168 <                colrs_gambs(scanline, xmax);
949 <                if (ourras->image->blue_mask & 1)
1166 >                tmap_colrs(scanline, xmax);
1167 >                switch (ourras->image->blue_mask) {
1168 >                case 0xff:              /* 24-bit RGB */
1169                          for (x = 0; x < xmax; x++)
1170                                  *dp++ = (unsigned int4)scanline[x][RED] << 16 |
1171                                          (unsigned int4)scanline[x][GRN] << 8 |
1172                                          (unsigned int4)scanline[x][BLU] ;
1173 <                else
1173 >                        break;
1174 >                case 0xff0000:          /* 24-bit BGR */
1175                          for (x = 0; x < xmax; x++)
1176                                  *dp++ = (unsigned int4)scanline[x][RED] |
1177                                          (unsigned int4)scanline[x][GRN] << 8 |
1178                                          (unsigned int4)scanline[x][BLU] << 16 ;
1179 +                        break;
1180 + #if 0
1181 +                case 0x1f:              /* 15-bit RGB */
1182 +                        for (x = 0; x < xmax; x++)
1183 +                                *dph++ =        (scanline[x][RED] << 7 & 0x7c00) |
1184 +                                        (scanline[x][GRN] << 2 & 0x3e0) |
1185 +                                        (unsigned)scanline[x][BLU] >> 3 ;
1186 +                        break;
1187 +                case 0x7c00:            /* 15-bit BGR */
1188 +                        for (x = 0; x < xmax; x++)
1189 +                                *dph++ =        (unsigned)scanline[x][RED] >> 3 |
1190 +                                        (scanline[x][GRN] << 2 & 0x3e0) |
1191 +                                        (scanline[x][BLU] << 7 & 0x7c00) ;
1192 +                        break;
1193 + #endif
1194 +                default:                /* unknown */
1195 +                        if (ourvis.depth > 16)
1196 +                                for (x = 0; x < xmax; x++) {
1197 +                                        *dp = ourras->image->red_mask &
1198 +                                                ourras->image->red_mask*scanline[x][RED]/255;
1199 +                                        *dp |= ourras->image->green_mask &
1200 +                                                ourras->image->green_mask*scanline[x][GRN]/255;
1201 +                                        *dp++ |= ourras->image->blue_mask &
1202 +                                                ourras->image->blue_mask*scanline[x][BLU]/255;
1203 +                                }
1204 +                        else
1205 +                                for (x = 0; x < xmax; x++) {
1206 +                                        *dph = ourras->image->red_mask &
1207 +                                                ourras->image->red_mask*scanline[x][RED]/255;
1208 +                                        *dph |= ourras->image->green_mask &
1209 +                                                ourras->image->green_mask*scanline[x][GRN]/255;
1210 +                                        *dph++ |= ourras->image->blue_mask &
1211 +                                                ourras->image->blue_mask*scanline[x][BLU]/255;
1212 +                                }
1213 +                        break;
1214 +                }
1215          }
1216   }
1217  
# Line 965 | Line 1221 | getgrey()                      /* get greyscale data */
1221          int     y;
1222          register unsigned char  *dp;
1223          register int    x;
1224 <                                        /* set gamma correction */
1225 <        setcolrgam(gamcor);
1224 >                                        /* initialize tone mapping */
1225 >        make_tonemap();
1226                                          /* read and convert file */
1227          dp = ourdata;
1228          for (y = 0; y < ymax; y++) {
1229                  getscan(y);
1230                  add2icon(y, scanline);
1231 <                if (scale)
976 <                        shiftcolrs(scanline, xmax, scale);
977 <                for (x = 0; x < xmax; x++)
978 <                        scanline[x][GRN] = normbright(scanline[x]);
979 <                colrs_gambs(scanline, xmax);
1231 >                tmap_colrs(scanline, xmax);
1232                  if (maxcolors < 256)
1233                          for (x = 0; x < xmax; x++)
1234 <                                *dp++ = ((long)scanline[x][GRN] *
1234 >                                *dp++ = ((int4)scanline[x][GRN] *
1235                                          maxcolors + maxcolors/2) >> 8;
1236                  else
1237                          for (x = 0; x < xmax; x++)
# Line 987 | Line 1239 | getgrey()                      /* get greyscale data */
1239          }
1240          for (x = 0; x < maxcolors; x++)
1241                  clrtab[x][RED] = clrtab[x][GRN] =
1242 <                        clrtab[x][BLU] = ((long)x*256 + 128)/maxcolors;
1242 >                        clrtab[x][BLU] = ((int4)x*256 + 128)/maxcolors;
1243   }
1244  
1245  
# Line 997 | Line 1249 | getmapped()                    /* get color-mapped data */
1249                                          /* make sure we can do it first */
1250          if (fname == NULL)
1251                  quiterr("cannot map colors from standard input");
1252 <                                        /* set gamma correction */
1253 <        setcolrgam(gamcor);
1252 >                                        /* initialize tone mapping */
1253 >        make_tonemap();
1254                                          /* make histogram */
1255 <        new_histo();
1255 >        if (new_histo((int4)xmax*ymax) == -1)
1256 >                quiterr("cannot initialize histogram");
1257          for (y = 0; y < ymax; y++) {
1258                  if (getscan(y) < 0)
1259                          break;
1260                  add2icon(y, scanline);
1261 <                if (scale)
1009 <                        shiftcolrs(scanline, xmax, scale);
1010 <                colrs_gambs(scanline, xmax);
1261 >                tmap_colrs(scanline, xmax);
1262                  cnt_colrs(scanline, xmax);
1263          }
1264                                          /* map pixels */
# Line 1015 | Line 1266 | getmapped()                    /* get color-mapped data */
1266                  quiterr("cannot create color map");
1267          for (y = 0; y < ymax; y++) {
1268                  getscan(y);
1269 <                if (scale)
1019 <                        shiftcolrs(scanline, xmax, scale);
1020 <                colrs_gambs(scanline, xmax);
1269 >                tmap_colrs(scanline, xmax);
1270                  if (dither)
1271                          dith_colrs(ourdata+y*xmax, scanline, xmax);
1272                  else

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines