ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/px/x11image.c
Revision: 2.26
Committed: Wed May 5 10:27:58 1993 UTC (30 years, 11 months ago) by greg
Content type: text/plain
Branch: MAIN
Changes since 2.25: +1 -1 lines
Log Message:
corrected factor in '@' calculation

File Contents

# Content
1 /* Copyright (c) 1993 Regents of the University of California */
2
3 #ifndef lint
4 static char SCCSid[] = "$SunId$ LBL";
5 #endif
6
7 /*
8 * x11image.c - driver for X-windows
9 *
10 * 3/1/90
11 */
12
13 /*
14 * Modified for X11
15 *
16 * January 1990
17 *
18 * Anat Grynberg and Greg Ward
19 */
20
21
22 #include "standard.h"
23
24 #include <X11/Xlib.h>
25 #include <X11/cursorfont.h>
26 #include <X11/Xutil.h>
27 #include <X11/Xatom.h>
28
29 #include "color.h"
30 #include "view.h"
31 #include "x11raster.h"
32 #include "random.h"
33 #include "resolu.h"
34
35 #define FONTNAME "8x13" /* text font we'll use */
36
37 #define CTRL(c) ((c)-'@')
38
39 #define BORWIDTH 5 /* border width */
40
41 #define ICONSIZ (8*10) /* maximum icon dimension (even 8) */
42
43 #define ourscreen DefaultScreen(thedisplay)
44 #define ourroot RootWindow(thedisplay,ourscreen)
45
46 #define revline(x0,y0,x1,y1) XDrawLine(thedisplay,wind,revgc,x0,y0,x1,y1)
47
48 #define redraw(x,y,w,h) patch_raster(wind,(x)-xoff,(y)-yoff,x,y,w,h,ourras)
49
50 double gamcor = 2.2; /* gamma correction */
51
52 int dither = 1; /* dither colors? */
53 int fast = 0; /* keep picture in Pixmap? */
54
55 char *dispname = NULL; /* our display name */
56
57 Window wind = 0; /* our output window */
58 unsigned long ourblack=0, ourwhite=1; /* black and white for this visual */
59 int maxcolors = 0; /* maximum colors */
60 int greyscale = 0; /* in grey */
61
62 int scale = 0; /* scalefactor; power of two */
63
64 int xoff = 0; /* x image offset */
65 int yoff = 0; /* y image offset */
66
67 VIEW ourview = STDVIEW; /* image view parameters */
68 int gotview = 0; /* got parameters from file */
69
70 COLR *scanline; /* scan line buffer */
71
72 RESOLU inpres; /* input resolution and ordering */
73 int xmax, ymax; /* picture dimensions */
74 int width, height; /* window size */
75 char *fname = NULL; /* input file name */
76 FILE *fin = stdin; /* input file */
77 long *scanpos = NULL; /* scan line positions in file */
78 int cury = 0; /* current scan location */
79
80 double exposure = 1.0; /* exposure compensation used */
81
82 int wrongformat = 0; /* input in another format? */
83
84 GC ourgc; /* standard graphics context */
85 GC revgc; /* graphics context with GXinvert */
86
87 int *ourrank; /* our visual class ranking */
88 XVisualInfo ourvis; /* our visual */
89 XRASTER *ourras; /* our stored image */
90 unsigned char *ourdata; /* our image data */
91
92 struct {
93 int xmin, ymin, xsiz, ysiz;
94 } box = {0, 0, 0, 0}; /* current box */
95
96 char *geometry = NULL; /* geometry specification */
97
98 char icondata[ICONSIZ*ICONSIZ/8]; /* icon bitmap data */
99 int iconwidth = 0, iconheight = 0;
100
101 char *progname;
102
103 char errmsg[128];
104
105 extern BYTE clrtab[256][3]; /* global color map */
106
107 extern long ftell();
108
109 Display *thedisplay;
110 Atom closedownAtom, wmProtocolsAtom;
111
112
113 main(argc, argv)
114 int argc;
115 char *argv[];
116 {
117 extern char *getenv();
118 char *gv;
119 int headline();
120 int i;
121
122 progname = argv[0];
123 if ((gv = getenv("GAMMA")) != NULL)
124 gamcor = atof(gv);
125
126 for (i = 1; i < argc; i++)
127 if (argv[i][0] == '-')
128 switch (argv[i][1]) {
129 case 'c':
130 maxcolors = atoi(argv[++i]);
131 break;
132 case 'b':
133 greyscale = !greyscale;
134 break;
135 case 'm':
136 maxcolors = 2;
137 break;
138 case 'd':
139 if (argv[i][2] == 'i')
140 dispname = argv[++i];
141 else
142 dither = !dither;
143 break;
144 case 'f':
145 fast = !fast;
146 break;
147 case 'e':
148 if (argv[i+1][0] != '+' && argv[i+1][0] != '-')
149 goto userr;
150 scale = atoi(argv[++i]);
151 break;
152 case 'g':
153 if (argv[i][2] == 'e')
154 geometry = argv[++i];
155 else
156 gamcor = atof(argv[++i]);
157 break;
158 default:
159 goto userr;
160 }
161 else if (argv[i][0] == '=')
162 geometry = argv[i];
163 else
164 break;
165
166 if (i == argc-1) {
167 fname = argv[i];
168 fin = fopen(fname, "r");
169 if (fin == NULL) {
170 sprintf(errmsg, "cannot open file \"%s\"", fname);
171 quiterr(errmsg);
172 }
173 } else if (i != argc)
174 goto userr;
175 /* get header */
176 getheader(fin, headline, NULL);
177 /* get picture dimensions */
178 if (wrongformat || !fgetsresolu(&inpres, fin))
179 quiterr("bad picture format");
180 xmax = scanlen(&inpres);
181 ymax = numscans(&inpres);
182 /* set view parameters */
183 if (gotview && setview(&ourview) != NULL)
184 gotview = 0;
185 if ((scanline = (COLR *)malloc(xmax*sizeof(COLR))) == NULL)
186 quiterr("out of memory");
187
188 init(argc, argv); /* get file and open window */
189
190 for ( ; ; )
191 getevent(); /* main loop */
192 userr:
193 fprintf(stderr,
194 "Usage: %s [-di disp][[-ge] spec][-b][-m][-d][-f][-c nclrs][-e +/-stops] pic\n",
195 progname);
196 exit(1);
197 }
198
199
200 headline(s) /* get relevant info from header */
201 char *s;
202 {
203 char fmt[32];
204
205 if (isexpos(s))
206 exposure *= exposval(s);
207 else if (isformat(s)) {
208 formatval(fmt, s);
209 wrongformat = strcmp(fmt, COLRFMT);
210 } else if (isview(s) && sscanview(&ourview, s) > 0)
211 gotview++;
212 }
213
214
215 init(argc, argv) /* get data and open window */
216 int argc;
217 char **argv;
218 {
219 XSetWindowAttributes ourwinattr;
220 XClassHint xclshints;
221 XWMHints xwmhints;
222 XSizeHints xszhints;
223 XTextProperty windowName, iconName;
224 XGCValues xgcv;
225 char *name;
226 register int i;
227
228 if (fname != NULL) {
229 scanpos = (long *)malloc(ymax*sizeof(long));
230 if (scanpos == NULL)
231 quiterr("out of memory");
232 for (i = 0; i < ymax; i++)
233 scanpos[i] = -1;
234 name = fname;
235 } else
236 name = progname;
237 /* remove directory prefix from name */
238 for (i = strlen(name); i-- > 0; )
239 if (name[i] == '/')
240 break;
241 name += i+1;
242 if ((thedisplay = XOpenDisplay(dispname)) == NULL)
243 quiterr("cannot open display");
244 /* get best visual for default screen */
245 getbestvis();
246 /* store image */
247 getras();
248 /* get size and position */
249 xszhints.flags = 0;
250 xszhints.width = xmax; xszhints.height = ymax;
251 if (geometry != NULL) {
252 i = XParseGeometry(geometry, &xszhints.x, &xszhints.y,
253 (unsigned *)&xszhints.width,
254 (unsigned *)&xszhints.height);
255 if ((i&(WidthValue|HeightValue)) == (WidthValue|HeightValue))
256 xszhints.flags |= USSize;
257 else
258 xszhints.flags |= PSize;
259 if ((i&(XValue|YValue)) == (XValue|YValue)) {
260 xszhints.flags |= USPosition;
261 if (i & XNegative)
262 xszhints.x += DisplayWidth(thedisplay,
263 ourscreen)-1-xszhints.width-2*BORWIDTH;
264 if (i & YNegative)
265 xszhints.y += DisplayHeight(thedisplay,
266 ourscreen)-1-xszhints.height-2*BORWIDTH;
267 }
268 }
269 /* open window */
270 ourwinattr.border_pixel = ourwhite;
271 ourwinattr.background_pixel = ourblack;
272 ourwinattr.colormap = XCreateColormap(thedisplay, ourroot,
273 ourvis.visual, AllocNone);
274 ourwinattr.event_mask = ExposureMask|KeyPressMask|ButtonPressMask|
275 ButtonReleaseMask|ButtonMotionMask|StructureNotifyMask;
276 ourwinattr.cursor = XCreateFontCursor(thedisplay, XC_diamond_cross);
277 wind = XCreateWindow(thedisplay, ourroot, xszhints.x, xszhints.y,
278 xszhints.width, xszhints.height, BORWIDTH,
279 ourvis.depth, InputOutput, ourvis.visual, CWEventMask|
280 CWCursor|CWBackPixel|CWBorderPixel|CWColormap, &ourwinattr);
281 if (wind == 0)
282 quiterr("cannot create window");
283 width = xmax;
284 height = ymax;
285 xgcv.foreground = ourblack;
286 xgcv.background = ourwhite;
287 if ((xgcv.font = XLoadFont(thedisplay, FONTNAME)) == 0)
288 quiterr("cannot get font");
289 ourgc = XCreateGC(thedisplay, wind, GCForeground|GCBackground|
290 GCFont, &xgcv);
291 xgcv.function = GXinvert;
292 revgc = XCreateGC(thedisplay, wind, GCForeground|GCBackground|
293 GCFunction, &xgcv);
294
295 /* set up the window manager */
296 xwmhints.flags = InputHint|IconPixmapHint;
297 xwmhints.input = True;
298 xwmhints.icon_pixmap = XCreateBitmapFromData(thedisplay,
299 wind, icondata, iconwidth, iconheight);
300
301 windowName.encoding = iconName.encoding = XA_STRING;
302 windowName.format = iconName.format = 8;
303 windowName.value = (u_char *)name;
304 windowName.nitems = strlen(windowName.value);
305 iconName.value = (u_char *)name;
306 iconName.nitems = strlen(windowName.value);
307
308 xclshints.res_name = NULL;
309 xclshints.res_class = "Ximage";
310 XSetWMProperties(thedisplay, wind, &windowName, &iconName,
311 argv, argc, &xszhints, &xwmhints, &xclshints);
312 closedownAtom = XInternAtom(thedisplay, "WM_DELETE_WINDOW", False);
313 wmProtocolsAtom = XInternAtom(thedisplay, "WM_PROTOCOLS", False);
314 XSetWMProtocols(thedisplay, wind, &closedownAtom, 1);
315
316 XMapWindow(thedisplay, wind);
317 return;
318 } /* end of init */
319
320
321 quiterr(err) /* print message and exit */
322 char *err;
323 {
324 if (err != NULL) {
325 fprintf(stderr, "%s: %s\n", progname, err);
326 exit(1);
327 }
328 exit(0);
329 }
330
331
332 static int
333 viscmp(v1,v2) /* compare visual to see which is better, descending */
334 register XVisualInfo *v1, *v2;
335 {
336 int bad1 = 0, bad2 = 0;
337 register int *rp;
338
339 if (v1->class == v2->class) {
340 if (v1->class == TrueColor || v1->class == DirectColor) {
341 /* prefer 24-bit to 32-bit */
342 if (v1->depth == 24 && v2->depth == 32)
343 return(-1);
344 if (v1->depth == 32 && v2->depth == 24)
345 return(1);
346 return(0);
347 }
348 /* don't be too greedy */
349 if (maxcolors <= 1<<v1->depth && maxcolors <= 1<<v2->depth)
350 return(v1->depth - v2->depth);
351 return(v2->depth - v1->depth);
352 }
353 /* prefer Pseudo when < 24-bit */
354 if ((v1->class == TrueColor || v1->class == DirectColor) &&
355 v1->depth < 24)
356 bad1 = 1;
357 if ((v2->class == TrueColor || v2->class == DirectColor) &&
358 v2->depth < 24)
359 bad2 = -1;
360 if (bad1 | bad2)
361 return(bad1+bad2);
362 /* otherwise, use class ranking */
363 for (rp = ourrank; *rp != -1; rp++) {
364 if (v1->class == *rp)
365 return(-1);
366 if (v2->class == *rp)
367 return(1);
368 }
369 return(0);
370 }
371
372
373 getbestvis() /* get the best visual for this screen */
374 {
375 #ifdef DEBUG
376 static char vistype[][12] = {
377 "StaticGray",
378 "GrayScale",
379 "StaticColor",
380 "PseudoColor",
381 "TrueColor",
382 "DirectColor"
383 };
384 #endif
385 static int rankings[3][6] = {
386 {TrueColor,DirectColor,PseudoColor,GrayScale,StaticGray,-1},
387 {PseudoColor,GrayScale,StaticGray,-1},
388 {PseudoColor,GrayScale,StaticGray,-1}
389 };
390 XVisualInfo *xvi;
391 int vismatched;
392 register int i, j;
393
394 if (greyscale) {
395 ourrank = rankings[2];
396 if (maxcolors < 2) maxcolors = 256;
397 } else if (maxcolors >= 2 && maxcolors <= 256)
398 ourrank = rankings[1];
399 else {
400 ourrank = rankings[0];
401 maxcolors = 256;
402 }
403 /* find best visual */
404 ourvis.screen = ourscreen;
405 xvi = XGetVisualInfo(thedisplay,VisualScreenMask,&ourvis,&vismatched);
406 if (xvi == NULL)
407 quiterr("no visuals for this screen!");
408 #ifdef DEBUG
409 fprintf(stderr, "Supported visuals:\n");
410 for (i = 0; i < vismatched; i++)
411 fprintf(stderr, "\ttype %s, depth %d\n",
412 vistype[xvi[i].class], xvi[i].depth);
413 #endif
414 for (i = 0, j = 1; j < vismatched; j++)
415 if (viscmp(&xvi[i],&xvi[j]) > 0)
416 i = j;
417 /* compare to least acceptable */
418 for (j = 0; ourrank[j++] != -1; )
419 ;
420 ourvis.class = ourrank[--j];
421 ourvis.depth = 1;
422 if (viscmp(&xvi[i],&ourvis) > 0)
423 quiterr("inadequate visuals on this screen");
424 /* OK, we'll use it */
425 copystruct(&ourvis, &xvi[i]);
426 #ifdef DEBUG
427 fprintf(stderr, "Selected visual type %s, depth %d\n",
428 vistype[ourvis.class], ourvis.depth);
429 #endif
430 /* make appropriate adjustments */
431 if (ourvis.class == GrayScale || ourvis.class == StaticGray)
432 greyscale = 1;
433 if (ourvis.depth <= 8 && ourvis.colormap_size < maxcolors)
434 maxcolors = ourvis.colormap_size;
435 if (ourvis.class == StaticGray) {
436 ourblack = 0;
437 ourwhite = 255;
438 } else if (ourvis.class == PseudoColor) {
439 ourblack = BlackPixel(thedisplay,ourscreen);
440 ourwhite = WhitePixel(thedisplay,ourscreen);
441 if ((ourblack|ourwhite) & ~255L) {
442 ourblack = 0;
443 ourwhite = 1;
444 }
445 if (maxcolors > 4)
446 maxcolors -= 2;
447 } else {
448 ourblack = 0;
449 ourwhite = ourvis.red_mask|ourvis.green_mask|ourvis.blue_mask;
450 }
451 XFree((char *)xvi);
452 }
453
454
455 getras() /* get raster file */
456 {
457 XVisualInfo vinfo;
458
459 if (maxcolors <= 2) { /* monochrome */
460 ourdata = (unsigned char *)malloc(ymax*((xmax+7)/8));
461 if (ourdata == NULL)
462 goto fail;
463 ourras = make_raster(thedisplay, &ourvis, 1, ourdata,
464 xmax, ymax, 8);
465 if (ourras == NULL)
466 goto fail;
467 getmono();
468 } else if (ourvis.class == TrueColor | ourvis.class == DirectColor) {
469 ourdata = (unsigned char *)malloc(4*xmax*ymax);
470 if (ourdata == NULL)
471 goto fail;
472 ourras = make_raster(thedisplay, &ourvis, 32,
473 ourdata, xmax, ymax, 32);
474 if (ourras == NULL)
475 goto fail;
476 getfull();
477 } else {
478 ourdata = (unsigned char *)malloc(xmax*ymax);
479 if (ourdata == NULL)
480 goto fail;
481 ourras = make_raster(thedisplay, &ourvis, 8, ourdata,
482 xmax, ymax, 8);
483 if (ourras == NULL)
484 goto fail;
485 if (greyscale | ourvis.class == StaticGray)
486 getgrey();
487 else
488 getmapped();
489 if (ourvis.class != StaticGray && !init_rcolors(ourras,clrtab))
490 goto fail;
491 }
492 return;
493 fail:
494 quiterr("could not create raster image");
495 }
496
497
498 getevent() /* process the next event */
499 {
500 XEvent xev;
501
502 XNextEvent(thedisplay, &xev);
503 switch ((int)xev.type) {
504 case KeyPress:
505 docom(&xev.xkey);
506 break;
507 case ConfigureNotify:
508 width = xev.xconfigure.width;
509 height = xev.xconfigure.height;
510 break;
511 case MapNotify:
512 map_rcolors(ourras, wind);
513 if (fast)
514 make_rpixmap(ourras, wind);
515 break;
516 case UnmapNotify:
517 if (!fast)
518 unmap_rcolors(ourras);
519 break;
520 case Expose:
521 redraw(xev.xexpose.x, xev.xexpose.y,
522 xev.xexpose.width, xev.xexpose.height);
523 break;
524 case ButtonPress:
525 if (xev.xbutton.state & (ShiftMask|ControlMask))
526 moveimage(&xev.xbutton);
527 else if (xev.xbutton.button == Button2)
528 traceray(xev.xbutton.x, xev.xbutton.y);
529 else
530 getbox(&xev.xbutton);
531 break;
532 case ClientMessage:
533 if ((xev.xclient.message_type == wmProtocolsAtom) &&
534 (xev.xclient.data.l[0] == closedownAtom))
535 quiterr(NULL);
536 break;
537 }
538 }
539
540
541 traceray(xpos, ypos) /* print ray corresponding to pixel */
542 int xpos, ypos;
543 {
544 FLOAT hv[2];
545 FVECT rorg, rdir;
546
547 if (!gotview) { /* no view, no can do */
548 XBell(thedisplay, 0);
549 return(-1);
550 }
551 pix2loc(hv, &inpres, xpos-xoff, ypos-yoff);
552 if (viewray(rorg, rdir, &ourview, hv[0], hv[1]) < 0)
553 return(-1);
554 printf("%e %e %e ", rorg[0], rorg[1], rorg[2]);
555 printf("%e %e %e\n", rdir[0], rdir[1], rdir[2]);
556 fflush(stdout);
557 return(0);
558 }
559
560
561 docom(ekey) /* execute command */
562 XKeyPressedEvent *ekey;
563 {
564 char buf[80];
565 COLOR cval;
566 XColor cvx;
567 int com, n;
568 double comp;
569 FLOAT hv[2];
570
571 n = XLookupString(ekey, buf, sizeof(buf), NULL, NULL);
572 if (n == 0)
573 return(0);
574 com = buf[0];
575 switch (com) { /* interpret command */
576 case 'q':
577 case 'Q':
578 case CTRL('D'): /* quit */
579 quiterr(NULL);
580 case '\n':
581 case '\r':
582 case 'l':
583 case 'c': /* value */
584 if (avgbox(cval) == -1)
585 return(-1);
586 switch (com) {
587 case '\n':
588 case '\r': /* radiance */
589 sprintf(buf, "%.3f", intens(cval)/exposure);
590 break;
591 case 'l': /* luminance */
592 sprintf(buf, "%.0fL", luminance(cval)/exposure);
593 break;
594 case 'c': /* color */
595 comp = pow(2.0, (double)scale);
596 sprintf(buf, "(%.2f,%.2f,%.2f)",
597 colval(cval,RED)*comp,
598 colval(cval,GRN)*comp,
599 colval(cval,BLU)*comp);
600 break;
601 }
602 XDrawImageString(thedisplay, wind, ourgc,
603 box.xmin, box.ymin+box.ysiz, buf, strlen(buf));
604 return(0);
605 case 'i': /* identify (contour) */
606 if (ourras->pixels == NULL)
607 return(-1);
608 n = ourdata[ekey->x-xoff+xmax*(ekey->y-yoff)];
609 n = ourras->pmap[n];
610 cvx.pixel = ourras->cdefs[n].pixel;
611 cvx.red = random() & 65535;
612 cvx.green = random() & 65535;
613 cvx.blue = random() & 65535;
614 cvx.flags = DoRed|DoGreen|DoBlue;
615 XStoreColor(thedisplay, ourras->cmap, &cvx);
616 return(0);
617 case 'p': /* position */
618 pix2loc(hv, &inpres, ekey->x-xoff, ekey->y-yoff);
619 sprintf(buf, "(%d,%d)", (int)(hv[0]*inpres.xr),
620 (int)(hv[1]*inpres.yr));
621 XDrawImageString(thedisplay, wind, ourgc, ekey->x, ekey->y,
622 buf, strlen(buf));
623 return(0);
624 case 't': /* trace */
625 return(traceray(ekey->x, ekey->y));
626 case '=': /* adjust exposure */
627 case '@': /* adaptation level */
628 if (avgbox(cval) == -1)
629 return(-1);
630 comp = com=='@'
631 ? 106./pow(1.219+pow(luminance(cval)/exposure,.4),2.5)/exposure
632 : .5/bright(cval) ;
633 comp = log(comp)/.69315 - scale;
634 n = comp < 0 ? comp-.5 : comp+.5 ; /* round */
635 if (n == 0)
636 return(0);
637 scale_rcolors(ourras, pow(2.0, (double)n));
638 scale += n;
639 sprintf(buf, "%+d", scale);
640 XDrawImageString(thedisplay, wind, ourgc,
641 box.xmin, box.ymin+box.ysiz, buf, strlen(buf));
642 XFlush(thedisplay);
643 free(ourdata);
644 free_raster(ourras);
645 getras();
646 /* fall through */
647 case CTRL('R'): /* redraw */
648 case CTRL('L'):
649 unmap_rcolors(ourras);
650 XClearWindow(thedisplay, wind);
651 map_rcolors(ourras, wind);
652 if (fast)
653 make_rpixmap(ourras, wind);
654 redraw(0, 0, width, height);
655 return(0);
656 case '0': /* recenter origin */
657 if (xoff == 0 & yoff == 0)
658 return(0);
659 xoff = yoff = 0;
660 XClearWindow(thedisplay, wind);
661 redraw(0, 0, width, height);
662 return(0);
663 case ' ': /* clear */
664 redraw(box.xmin, box.ymin, box.xsiz, box.ysiz);
665 return(0);
666 default:
667 XBell(thedisplay, 0);
668 return(-1);
669 }
670 }
671
672
673 moveimage(ebut) /* shift the image */
674 XButtonPressedEvent *ebut;
675 {
676 XEvent e;
677 int mxo, myo;
678
679 XMaskEvent(thedisplay, ButtonReleaseMask|ButtonMotionMask, &e);
680 while (e.type == MotionNotify) {
681 mxo = e.xmotion.x;
682 myo = e.xmotion.y;
683 revline(ebut->x, ebut->y, mxo, myo);
684 revbox(xoff+mxo-ebut->x, yoff+myo-ebut->y,
685 xoff+mxo-ebut->x+xmax, yoff+myo-ebut->y+ymax);
686 XMaskEvent(thedisplay,ButtonReleaseMask|ButtonMotionMask,&e);
687 revline(ebut->x, ebut->y, mxo, myo);
688 revbox(xoff+mxo-ebut->x, yoff+myo-ebut->y,
689 xoff+mxo-ebut->x+xmax, yoff+myo-ebut->y+ymax);
690 }
691 xoff += e.xbutton.x - ebut->x;
692 yoff += e.xbutton.y - ebut->y;
693 XClearWindow(thedisplay, wind);
694 redraw(0, 0, width, height);
695 }
696
697
698 getbox(ebut) /* get new box */
699 XButtonPressedEvent *ebut;
700 {
701 XEvent e;
702
703 XMaskEvent(thedisplay, ButtonReleaseMask|ButtonMotionMask, &e);
704 while (e.type == MotionNotify) {
705 revbox(ebut->x, ebut->y, box.xmin = e.xmotion.x, box.ymin = e.xmotion.y);
706 XMaskEvent(thedisplay,ButtonReleaseMask|ButtonMotionMask,&e);
707 revbox(ebut->x, ebut->y, box.xmin, box.ymin);
708 }
709 box.xmin = e.xbutton.x<0 ? 0 : (e.xbutton.x>=width ? width-1 : e.xbutton.x);
710 box.ymin = e.xbutton.y<0 ? 0 : (e.xbutton.y>=height ? height-1 : e.xbutton.y);
711 if (box.xmin > ebut->x) {
712 box.xsiz = box.xmin - ebut->x + 1;
713 box.xmin = ebut->x;
714 } else {
715 box.xsiz = ebut->x - box.xmin + 1;
716 }
717 if (box.ymin > ebut->y) {
718 box.ysiz = box.ymin - ebut->y + 1;
719 box.ymin = ebut->y;
720 } else {
721 box.ysiz = ebut->y - box.ymin + 1;
722 }
723 }
724
725
726 revbox(x0, y0, x1, y1) /* draw box with reversed lines */
727 int x0, y0, x1, y1;
728 {
729 revline(x0, y0, x1, y0);
730 revline(x0, y1, x1, y1);
731 revline(x0, y0, x0, y1);
732 revline(x1, y0, x1, y1);
733 }
734
735
736 avgbox(clr) /* average color over current box */
737 COLOR clr;
738 {
739 static COLOR lc;
740 static int ll, lr, lt, lb;
741 int left, right, top, bottom;
742 int y;
743 double d;
744 COLOR ctmp;
745 register int x;
746
747 setcolor(clr, 0.0, 0.0, 0.0);
748 left = box.xmin - xoff;
749 right = left + box.xsiz;
750 if (left < 0)
751 left = 0;
752 if (right > xmax)
753 right = xmax;
754 if (left >= right)
755 return(-1);
756 top = box.ymin - yoff;
757 bottom = top + box.ysiz;
758 if (top < 0)
759 top = 0;
760 if (bottom > ymax)
761 bottom = ymax;
762 if (top >= bottom)
763 return(-1);
764 if (left == ll && right == lr && top == lt && bottom == lb) {
765 copycolor(clr, lc);
766 return(0);
767 }
768 for (y = top; y < bottom; y++) {
769 if (getscan(y) == -1)
770 return(-1);
771 for (x = left; x < right; x++) {
772 colr_color(ctmp, scanline[x]);
773 addcolor(clr, ctmp);
774 }
775 }
776 d = 1.0/((right-left)*(bottom-top));
777 scalecolor(clr, d);
778 ll = left; lr = right; lt = top; lb = bottom;
779 copycolor(lc, clr);
780 return(0);
781 }
782
783
784 getmono() /* get monochrome data */
785 {
786 register unsigned char *dp;
787 register int x, err;
788 int y, errp;
789 short *cerr;
790
791 if ((cerr = (short *)calloc(xmax,sizeof(short))) == NULL)
792 quiterr("out of memory in getmono");
793 dp = ourdata - 1;
794 for (y = 0; y < ymax; y++) {
795 getscan(y);
796 add2icon(y, scanline);
797 normcolrs(scanline, xmax, scale);
798 err = 0;
799 for (x = 0; x < xmax; x++) {
800 if (!(x&7))
801 *++dp = 0;
802 errp = err;
803 err += normbright(scanline[x]) + cerr[x];
804 if (err > 127)
805 err -= 255;
806 else
807 *dp |= 1<<(7-(x&07));
808 err /= 3;
809 cerr[x] = err + errp;
810 }
811 }
812 free((char *)cerr);
813 }
814
815
816 add2icon(y, scan) /* add a scanline to our icon data */
817 int y;
818 COLR *scan;
819 {
820 static short cerr[ICONSIZ];
821 static int ynext;
822 static char *dp;
823 COLR clr;
824 register int err;
825 register int x, ti;
826 int errp;
827
828 if (iconheight == 0) { /* initialize */
829 if (xmax <= ICONSIZ && ymax <= ICONSIZ) {
830 iconwidth = xmax;
831 iconheight = ymax;
832 } else if (xmax > ymax) {
833 iconwidth = ICONSIZ;
834 iconheight = ICONSIZ*ymax/xmax;
835 if (iconheight < 1)
836 iconheight = 1;
837 } else {
838 iconwidth = ICONSIZ*xmax/ymax;
839 if (iconwidth < 1)
840 iconwidth = 1;
841 iconheight = ICONSIZ;
842 }
843 ynext = 0;
844 dp = icondata - 1;
845 }
846 if (y < ynext*ymax/iconheight) /* skip this one */
847 return;
848 err = 0;
849 for (x = 0; x < iconwidth; x++) {
850 if (!(x&7))
851 *++dp = 0;
852 errp = err;
853 ti = x*xmax/iconwidth;
854 copycolr(clr, scan[ti]);
855 normcolrs(clr, 1, scale);
856 err += normbright(clr) + cerr[x];
857 if (err > 127)
858 err -= 255;
859 else
860 *dp |= 1<<(x&07);
861 err /= 3;
862 cerr[x] = err + errp;
863 }
864 ynext++;
865 }
866
867
868 getfull() /* get full (24-bit) data */
869 {
870 int y;
871 register unsigned long *dp;
872 register int x;
873 /* set gamma correction */
874 setcolrgam(gamcor);
875 /* read and convert file */
876 dp = (unsigned long *)ourdata;
877 for (y = 0; y < ymax; y++) {
878 getscan(y);
879 add2icon(y, scanline);
880 if (scale)
881 shiftcolrs(scanline, xmax, scale);
882 colrs_gambs(scanline, xmax);
883 if (ourras->image->blue_mask & 1)
884 for (x = 0; x < xmax; x++)
885 *dp++ = scanline[x][RED] << 16 |
886 scanline[x][GRN] << 8 |
887 scanline[x][BLU] ;
888 else
889 for (x = 0; x < xmax; x++)
890 *dp++ = scanline[x][RED] |
891 scanline[x][GRN] << 8 |
892 scanline[x][BLU] << 16 ;
893 }
894 }
895
896
897 getgrey() /* get greyscale data */
898 {
899 int y;
900 register unsigned char *dp;
901 register int x;
902 /* set gamma correction */
903 setcolrgam(gamcor);
904 /* read and convert file */
905 dp = ourdata;
906 for (y = 0; y < ymax; y++) {
907 getscan(y);
908 add2icon(y, scanline);
909 if (scale)
910 shiftcolrs(scanline, xmax, scale);
911 for (x = 0; x < xmax; x++)
912 scanline[x][GRN] = normbright(scanline[x]);
913 colrs_gambs(scanline, xmax);
914 if (maxcolors < 256)
915 for (x = 0; x < xmax; x++)
916 *dp++ = ((long)scanline[x][GRN] *
917 maxcolors + maxcolors/2) >> 8;
918 else
919 for (x = 0; x < xmax; x++)
920 *dp++ = scanline[x][GRN];
921 }
922 for (x = 0; x < maxcolors; x++)
923 clrtab[x][RED] = clrtab[x][GRN] =
924 clrtab[x][BLU] = ((long)x*256 + 128)/maxcolors;
925 }
926
927
928 getmapped() /* get color-mapped data */
929 {
930 int y;
931 /* set gamma correction */
932 setcolrgam(gamcor);
933 /* make histogram */
934 new_histo();
935 for (y = 0; y < ymax; y++) {
936 if (getscan(y) < 0)
937 quiterr("seek error in getmapped");
938 add2icon(y, scanline);
939 if (scale)
940 shiftcolrs(scanline, xmax, scale);
941 colrs_gambs(scanline, xmax);
942 cnt_colrs(scanline, xmax);
943 }
944 /* map pixels */
945 if (!new_clrtab(maxcolors))
946 quiterr("cannot create color map");
947 for (y = 0; y < ymax; y++) {
948 if (getscan(y) < 0)
949 quiterr("seek error in getmapped");
950 if (scale)
951 shiftcolrs(scanline, xmax, scale);
952 colrs_gambs(scanline, xmax);
953 if (dither)
954 dith_colrs(ourdata+y*xmax, scanline, xmax);
955 else
956 map_colrs(ourdata+y*xmax, scanline, xmax);
957 }
958 }
959
960
961 scale_rcolors(xr, sf) /* scale color map */
962 register XRASTER *xr;
963 double sf;
964 {
965 register int i;
966 long maxv;
967
968 if (xr->pixels == NULL)
969 return;
970
971 sf = pow(sf, 1.0/gamcor);
972 maxv = 65535/sf;
973
974 for (i = xr->ncolors; i--; ) {
975 xr->cdefs[i].red = xr->cdefs[i].red > maxv ?
976 65535 :
977 xr->cdefs[i].red * sf;
978 xr->cdefs[i].green = xr->cdefs[i].green > maxv ?
979 65535 :
980 xr->cdefs[i].green * sf;
981 xr->cdefs[i].blue = xr->cdefs[i].blue > maxv ?
982 65535 :
983 xr->cdefs[i].blue * sf;
984 }
985 XStoreColors(thedisplay, xr->cmap, xr->cdefs, xr->ncolors);
986 }
987
988
989 getscan(y)
990 int y;
991 {
992 if (y != cury) {
993 if (scanpos == NULL || scanpos[y] == -1)
994 return(-1);
995 if (fseek(fin, scanpos[y], 0) == -1)
996 quiterr("fseek error");
997 cury = y;
998 } else if (scanpos != NULL && scanpos[y] == -1)
999 scanpos[y] = ftell(fin);
1000
1001 if (freadcolrs(scanline, xmax, fin) < 0)
1002 quiterr("read error");
1003
1004 cury++;
1005 return(0);
1006 }