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

Comparing ray/src/px/x11image.c (file contents):
Revision 2.43 by greg, Mon Jan 27 15:07:03 1997 UTC vs.
Revision 2.51 by gwlarson, Mon Feb 9 14:37:09 1998 UTC

# Line 1 | Line 1
1 < /* Copyright (c) 1997 Regents of the University of California */
1 > /* Copyright (c) 1998 Silicon Graphics, Inc. */
2  
3   #ifndef lint
4 < static char SCCSid[] = "$SunId$ LBL";
4 > static char SCCSid[] = "$SunId$ SGI";
5   #endif
6  
7   /*
# Line 27 | Line 27 | static char SCCSid[] = "$SunId$ LBL";
27   #include  <X11/Xutil.h>
28   #include  <X11/Xatom.h>
29  
30 < #include  "color.h"
30 > #undef  NOPROTO
31 > #define NOPROTO 1
32 > #include  "tonemap.h"
33   #include  "view.h"
34   #include  "x11raster.h"
35   #include  "random.h"
36   #include  "resolu.h"
37  
36 – #ifdef  __alpha
37 – #define  int4           int
38 – #endif
38   #ifndef  int4
39 < #define  int4           long
39 > #define  int4           int             /* most int's are 32-bit */
40   #endif
41  
42   #define  FONTNAME       "8x13"          /* text font we'll use */
# Line 48 | Line 47 | static char SCCSid[] = "$SunId$ LBL";
47  
48   #define  ICONSIZ        (8*10)          /* maximum icon dimension (even 8) */
49  
50 + #define  FIXWEIGHT      20              /* weight to add for fixation points */
51 +
52   #define  ourscreen      DefaultScreen(thedisplay)
53   #define  ourroot        RootWindow(thedisplay,ourscreen)
54  
# Line 79 | Line 80 | int  sequential = 0;                   /* display images in sequence *
80   char  *tout = "od";                     /* output of 't' command */
81   int  tinterv = 0;                       /* interval between mouse reports */
82  
83 + int  tmflags = TM_F_LINEAR;             /* tone mapping flags */
84 +
85   VIEW  ourview = STDVIEW;                /* image view parameters */
86   int  gotview = 0;                       /* got parameters from file */
87  
88   COLR  *scanline;                        /* scan line buffer */
89 + TMbright  *lscan;                       /* encoded luminance scanline */
90 + BYTE  *cscan;                           /* encoded chroma scanline */
91 + BYTE  *pscan;                           /* compute pixel scanline */
92  
93   RESOLU  inpres;                         /* input resolution and ordering */
94   int  xmax, ymax;                        /* picture dimensions */
# Line 173 | Line 179 | char  *argv[];
179                                  tinterv = atoi(argv[++i]);
180                                  break;
181                          case 'e':                       /* exposure comp. */
182 <                                if (argv[i+1][0] != '+' && argv[i+1][0] != '-')
182 >                                i++;
183 >                                if (argv[i][0] == 'a') {
184 >                                        tmflags = TM_F_CAMERA;
185 >                                        break;
186 >                                }
187 >                                if (argv[i][0] == 'h') {
188 >                                        tmflags = TM_F_HUMAN;
189 >                                        break;
190 >                                }
191 >                                if (argv[i][0] != '+' && argv[i][0] != '-')
192                                          goto userr;
193 <                                scale = atoi(argv[++i]);
193 >                                scale = atoi(argv[i]);
194                                  break;
195                          case 'g':                       /* gamma comp. */
196                                  if (argv[i][2] == 'e')
# Line 232 | Line 247 | char  *argv[];
247                  getevent();             /* main loop */
248   userr:
249          fprintf(stderr,
250 < "Usage: %s [-di disp][[-ge] spec][-b][-m][-d][-f][-c nclrs][-e +/-stops][-g gamcor][-s][-ospec][-t intvl] pic ..\n",
250 > "Usage: %s [-di disp][[-ge] spec][-b][-m][-d][-f][-c nclrs][-e spec][-g gamcor][-s][-ospec][-t intvl] pic ..\n",
251                          progname);
252          exit(1);
253   }
# Line 245 | Line 260 | char  *s;
260  
261          if (isexpos(s))
262                  exposure *= exposval(s);
263 <        else if (isformat(s)) {
249 <                formatval(fmt, s);
263 >        else if (formatval(fmt, s))
264                  wrongformat = strcmp(fmt, COLRFMT);
265 <        } else if (isview(s) && sscanview(&ourview, s) > 0)
265 >        else if (isview(s) && sscanview(&ourview, s) > 0)
266                  gotview++;
267   }
268  
# Line 544 | Line 558 | getras()                               /* get raster file */
558                                  xmax, ymax, 8);
559                  if (ourras == NULL)
560                          goto fail;
561 <                if (greyscale | ourvis.class == StaticGray)
561 >                if (greyscale)
562                          getgrey();
563                  else
564                          getmapped();
# Line 679 | Line 693 | XKeyPressedEvent  *ekey;
693          case '\r':
694          case 'l':
695          case 'c':                               /* value */
696 <                if (avgbox(cval) == -1)
696 >                if (!avgbox(cval))
697                          return(-1);
698                  switch (com) {
699                  case '\n':
# Line 721 | Line 735 | XKeyPressedEvent  *ekey;
735                  return(0);
736          case 't':                               /* trace */
737                  return(traceray(ekey->x, ekey->y));
738 +        case 'a':                               /* auto exposure */
739 +                if (fname == NULL)
740 +                        return(-1);
741 +                tmflags = TM_F_CAMERA;
742 +                strcpy(buf, "auto exposure...");
743 +                goto remap;
744 +        case 'h':                               /* human response */
745 +                if (fname == NULL)
746 +                        return(-1);
747 +                tmflags = TM_F_HUMAN;
748 +                strcpy(buf, "human exposure...");
749 +                goto remap;
750          case '=':                               /* adjust exposure */
751          case '@':                               /* adaptation level */
752 <                if (avgbox(cval) == -1)
752 >                if (!avgbox(cval))
753                          return(-1);
754                  comp = bright(cval);
755                  if (comp < 1e-20) {
# Line 737 | Line 763 | XKeyPressedEvent  *ekey;
763                          comp = .5/comp;
764                  comp = log(comp)/.69315 - scale;
765                  n = comp < 0 ? comp-.5 : comp+.5 ;      /* round */
766 <                if (n == 0)
767 <                        return(0);
768 <                scale_rcolors(ourras, pow(2.0, (double)n));
766 >                if (tmflags != TM_F_LINEAR)
767 >                        tmflags = TM_F_LINEAR;  /* turn off tone mapping */
768 >                else {
769 >                        if (n == 0)             /* else check if any change */
770 >                                return(0);
771 >                        scale_rcolors(ourras, pow(2.0, (double)n));
772 >                }
773                  scale += n;
774                  sprintf(buf, "%+d", scale);
775 +        remap:
776                  XDrawImageString(thedisplay, wind, ourgc,
777                                  box.xmin, box.ymin+box.ysiz, buf, strlen(buf));
778                  XFlush(thedisplay);
# Line 865 | Line 896 | int  x0, y0, x1, y1;
896   }
897  
898  
899 < avgbox(clr)                             /* average color over current box */
900 < COLOR  clr;
899 > int
900 > colavg(scn, n, cavg)
901 > register COLR   *scn;
902 > register int    n;
903 > COLOR   cavg;
904   {
905 <        static COLOR  lc;
906 <        static int  ll, lr, lt, lb;
905 >        COLOR   col;
906 >
907 >        while (n--) {
908 >                colr_color(col, scn++);
909 >                addcolor(cavg, col);
910 >        }
911 > }
912 >
913 >
914 > int
915 > avgbox(cavg)                            /* average color over current box */
916 > COLOR   cavg;
917 > {
918 >        double  d;
919 >        register int    rval;
920 >
921 >        setcolor(cavg, 0., 0., 0.);
922 >        rval = dobox(colavg, (char *)cavg);
923 >        if (rval > 0) {
924 >                d = 1./rval;
925 >                scalecolor(cavg, d);
926 >        }
927 >        return(rval);
928 > }
929 >
930 >
931 > int
932 > dobox(f, p)                             /* run function over box */
933 > int     (*f)();                 /* function to call for each subscan */
934 > char    *p;                     /* pointer to private data */
935 > {
936          int  left, right, top, bottom;
937          int  y;
875 –        double  d;
876 –        COLOR  ctmp;
877 –        register int  x;
938  
879 –        setcolor(clr, 0.0, 0.0, 0.0);
939          left = box.xmin - xoff;
940          right = left + box.xsiz;
941          if (left < 0)
# Line 884 | Line 943 | COLOR  clr;
943          if (right > xmax)
944                  right = xmax;
945          if (left >= right)
946 <                return(-1);
946 >                return(0);
947          top = box.ymin - yoff;
948          bottom = top + box.ysiz;
949          if (top < 0)
# Line 892 | Line 951 | COLOR  clr;
951          if (bottom > ymax)
952                  bottom = ymax;
953          if (top >= bottom)
895 –                return(-1);
896 –        if (left == ll && right == lr && top == lt && bottom == lb) {
897 –                copycolor(clr, lc);
954                  return(0);
899 –        }
955          for (y = top; y < bottom; y++) {
956                  if (getscan(y) == -1)
957                          return(-1);
958 <                for (x = left; x < right; x++) {
959 <                        colr_color(ctmp, scanline[x]);
960 <                        addcolor(clr, ctmp);
958 >                (*f)(scanline+left, right-left, p);
959 >        }
960 >        return((right-left)*(bottom-top));
961 > }
962 >
963 >
964 > int
965 > addfix(scn, n)                  /* add fixation points to histogram */
966 > COLR    *scn;
967 > int     n;
968 > {
969 >        if (tmCvColrs(lscan, TM_NOCHROM, scn, n))
970 >                goto tmerr;
971 >        if (tmAddHisto(lscan, n, FIXWEIGHT))
972 >                goto tmerr;
973 >        return;
974 > tmerr:
975 >        quiterr("tone mapping error");
976 > }
977 >
978 >
979 > make_tonemap()                  /* initialize tone mapping */
980 > {
981 >        int  flags, y;
982 >
983 >        if (tmflags != TM_F_LINEAR && fname == NULL) {
984 >                fprintf(stderr, "%s: cannot adjust tone of standard input\n",
985 >                                progname);
986 >                tmflags = TM_F_LINEAR;
987 >        }
988 >        if (tmflags == TM_F_LINEAR) {   /* linear with clamping */
989 >                setcolrcor(pow, 1.0/gamcor);
990 >                return;
991 >        }
992 >        flags = tmflags;                /* histogram adjustment */
993 >        if (greyscale) flags |= TM_F_BW;
994 >        if (tmTop != NULL) {            /* reuse old histogram if one */
995 >                tmDone(tmTop);
996 >                tmTop->flags = flags;
997 >        } else {                        /* else initialize */
998 >                if ((lscan = (TMbright *)malloc(xmax*sizeof(TMbright))) == NULL)
999 >                        goto memerr;
1000 >                if (greyscale) {
1001 >                        cscan = TM_NOCHROM;
1002 >                        if ((pscan = (BYTE *)malloc(sizeof(BYTE)*xmax)) == NULL)
1003 >                                goto memerr;
1004 >                } else if ((pscan=cscan = (BYTE *)malloc(3*sizeof(BYTE)*xmax))
1005 >                                == NULL)
1006 >                        goto memerr;
1007 >                                                /* initialize tm library */
1008 >                if (tmInit(flags, stdprims, gamcor) == NULL)
1009 >                        goto memerr;
1010 >                if (tmSetSpace(stdprims, WHTEFFICACY/exposure))
1011 >                        goto tmerr;
1012 >                                                /* compute picture histogram */
1013 >                for (y = 0; y < ymax; y++) {
1014 >                        getscan(y);
1015 >                        if (tmCvColrs(lscan, TM_NOCHROM, scanline, xmax))
1016 >                                goto tmerr;
1017 >                        if (tmAddHisto(lscan, xmax, 1))
1018 >                                goto tmerr;
1019                  }
1020          }
1021 <        d = 1.0/((right-left)*(bottom-top));
1022 <        scalecolor(clr, d);
1023 <        ll = left; lr = right; lt = top; lb = bottom;
1024 <        copycolor(lc, clr);
1025 <        return(0);
1021 >        tmDup();                        /* add fixations to duplicate map */
1022 >        dobox(addfix, NULL);
1023 >                                        /* (re)compute tone mapping */
1024 >        if (tmComputeMapping(gamcor, 0., 0.))
1025 >                goto tmerr;
1026 >        return;
1027 > memerr:
1028 >        quiterr("out of memory in make_tonemap");
1029 > tmerr:
1030 >        quiterr("tone mapping error");
1031   }
1032  
1033  
1034 + tmap_colrs(scn, len)            /* apply tone mapping to scanline */
1035 + register COLR  *scn;
1036 + int  len;
1037 + {
1038 +        register BYTE  *ps;
1039 +
1040 +        if (tmflags == TM_F_LINEAR) {
1041 +                if (scale)
1042 +                        shiftcolrs(scn, len, scale);
1043 +                colrs_gambs(scn, len);
1044 +                return;
1045 +        }
1046 +        if (len > xmax)
1047 +                quiterr("code error 1 in tmap_colrs");
1048 +        if (tmCvColrs(lscan, cscan, scn, len))
1049 +                goto tmerr;
1050 +        if (tmMapPixels(pscan, lscan, cscan, len))
1051 +                goto tmerr;
1052 +        ps = pscan;
1053 +        if (greyscale)
1054 +                while (len--) {
1055 +                        scn[0][RED] = scn[0][GRN] = scn[0][BLU] = *ps++;
1056 +                        scn[0][EXP] = COLXS;
1057 +                        scn++;
1058 +                }
1059 +        else
1060 +                while (len--) {
1061 +                        scn[0][RED] = *ps++;
1062 +                        scn[0][GRN] = *ps++;
1063 +                        scn[0][BLU] = *ps++;
1064 +                        scn[0][EXP] = COLXS;
1065 +                        scn++;
1066 +                }
1067 +        return;
1068 + tmerr:
1069 +        quiterr("tone mapping error");
1070 + }
1071 +
1072 +
1073   getmono()                       /* get monochrome data */
1074   {
1075          register unsigned char  *dp;
# Line 1002 | Line 1159 | getfull()                      /* get full (24-bit) data */
1159          int     y;
1160          register unsigned int4  *dp;
1161          register int    x;
1162 <                                        /* set gamma correction */
1163 <        setcolrgam(gamcor);
1162 >                                        /* initialize tone mapping */
1163 >        make_tonemap();
1164                                          /* read and convert file */
1165          dp = (unsigned int4 *)ourdata;
1166          for (y = 0; y < ymax; y++) {
1167                  getscan(y);
1168                  add2icon(y, scanline);
1169 <                if (scale)
1013 <                        shiftcolrs(scanline, xmax, scale);
1014 <                colrs_gambs(scanline, xmax);
1169 >                tmap_colrs(scanline, xmax);
1170                  if (ourras->image->blue_mask & 1)
1171                          for (x = 0; x < xmax; x++)
1172                                  *dp++ = (unsigned int4)scanline[x][RED] << 16 |
# Line 1031 | Line 1186 | getgrey()                      /* get greyscale data */
1186          int     y;
1187          register unsigned char  *dp;
1188          register int    x;
1189 <                                        /* set gamma correction */
1190 <        setcolrgam(gamcor);
1189 >                                        /* initialize tone mapping */
1190 >        make_tonemap();
1191                                          /* read and convert file */
1192          dp = ourdata;
1193          for (y = 0; y < ymax; y++) {
1194                  getscan(y);
1195                  add2icon(y, scanline);
1196 <                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);
1196 >                tmap_colrs(scanline, xmax);
1197                  if (maxcolors < 256)
1198                          for (x = 0; x < xmax; x++)
1199 <                                *dp++ = ((long)scanline[x][GRN] *
1199 >                                *dp++ = ((int4)scanline[x][GRN] *
1200                                          maxcolors + maxcolors/2) >> 8;
1201                  else
1202                          for (x = 0; x < xmax; x++)
# Line 1053 | Line 1204 | getgrey()                      /* get greyscale data */
1204          }
1205          for (x = 0; x < maxcolors; x++)
1206                  clrtab[x][RED] = clrtab[x][GRN] =
1207 <                        clrtab[x][BLU] = ((long)x*256 + 128)/maxcolors;
1207 >                        clrtab[x][BLU] = ((int4)x*256 + 128)/maxcolors;
1208   }
1209  
1210  
# Line 1063 | Line 1214 | getmapped()                    /* get color-mapped data */
1214                                          /* make sure we can do it first */
1215          if (fname == NULL)
1216                  quiterr("cannot map colors from standard input");
1217 <                                        /* set gamma correction */
1218 <        setcolrgam(gamcor);
1217 >                                        /* initialize tone mapping */
1218 >        make_tonemap();
1219                                          /* make histogram */
1220 <        if (new_histo((long)xmax*ymax) == -1)
1220 >        if (new_histo((int4)xmax*ymax) == -1)
1221                  quiterr("cannot initialize histogram");
1222          for (y = 0; y < ymax; y++) {
1223                  if (getscan(y) < 0)
1224                          break;
1225                  add2icon(y, scanline);
1226 <                if (scale)
1076 <                        shiftcolrs(scanline, xmax, scale);
1077 <                colrs_gambs(scanline, xmax);
1226 >                tmap_colrs(scanline, xmax);
1227                  cnt_colrs(scanline, xmax);
1228          }
1229                                          /* map pixels */
# Line 1082 | Line 1231 | getmapped()                    /* get color-mapped data */
1231                  quiterr("cannot create color map");
1232          for (y = 0; y < ymax; y++) {
1233                  getscan(y);
1234 <                if (scale)
1086 <                        shiftcolrs(scanline, xmax, scale);
1087 <                colrs_gambs(scanline, xmax);
1234 >                tmap_colrs(scanline, xmax);
1235                  if (dither)
1236                          dith_colrs(ourdata+y*xmax, scanline, xmax);
1237                  else

Diff Legend

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