ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/px/x11image.c
Revision: 2.27
Committed: Fri Jun 18 10:22:26 1993 UTC (30 years, 10 months ago) by greg
Content type: text/plain
Branch: MAIN
Changes since 2.26: +8 -0 lines
Log Message:
added 'f' and 'F' commands to turn on and off fast redraw

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 'f': /* turn on fast redraw */
657 fast = 1;
658 make_rpixmap(ourras, wind);
659 return(0);
660 case 'F': /* turn off fast redraw */
661 fast = 0;
662 free_rpixmap(ourras);
663 return(0);
664 case '0': /* recenter origin */
665 if (xoff == 0 & yoff == 0)
666 return(0);
667 xoff = yoff = 0;
668 XClearWindow(thedisplay, wind);
669 redraw(0, 0, width, height);
670 return(0);
671 case ' ': /* clear */
672 redraw(box.xmin, box.ymin, box.xsiz, box.ysiz);
673 return(0);
674 default:
675 XBell(thedisplay, 0);
676 return(-1);
677 }
678 }
679
680
681 moveimage(ebut) /* shift the image */
682 XButtonPressedEvent *ebut;
683 {
684 XEvent e;
685 int mxo, myo;
686
687 XMaskEvent(thedisplay, ButtonReleaseMask|ButtonMotionMask, &e);
688 while (e.type == MotionNotify) {
689 mxo = e.xmotion.x;
690 myo = e.xmotion.y;
691 revline(ebut->x, ebut->y, mxo, myo);
692 revbox(xoff+mxo-ebut->x, yoff+myo-ebut->y,
693 xoff+mxo-ebut->x+xmax, yoff+myo-ebut->y+ymax);
694 XMaskEvent(thedisplay,ButtonReleaseMask|ButtonMotionMask,&e);
695 revline(ebut->x, ebut->y, mxo, myo);
696 revbox(xoff+mxo-ebut->x, yoff+myo-ebut->y,
697 xoff+mxo-ebut->x+xmax, yoff+myo-ebut->y+ymax);
698 }
699 xoff += e.xbutton.x - ebut->x;
700 yoff += e.xbutton.y - ebut->y;
701 XClearWindow(thedisplay, wind);
702 redraw(0, 0, width, height);
703 }
704
705
706 getbox(ebut) /* get new box */
707 XButtonPressedEvent *ebut;
708 {
709 XEvent e;
710
711 XMaskEvent(thedisplay, ButtonReleaseMask|ButtonMotionMask, &e);
712 while (e.type == MotionNotify) {
713 revbox(ebut->x, ebut->y, box.xmin = e.xmotion.x, box.ymin = e.xmotion.y);
714 XMaskEvent(thedisplay,ButtonReleaseMask|ButtonMotionMask,&e);
715 revbox(ebut->x, ebut->y, box.xmin, box.ymin);
716 }
717 box.xmin = e.xbutton.x<0 ? 0 : (e.xbutton.x>=width ? width-1 : e.xbutton.x);
718 box.ymin = e.xbutton.y<0 ? 0 : (e.xbutton.y>=height ? height-1 : e.xbutton.y);
719 if (box.xmin > ebut->x) {
720 box.xsiz = box.xmin - ebut->x + 1;
721 box.xmin = ebut->x;
722 } else {
723 box.xsiz = ebut->x - box.xmin + 1;
724 }
725 if (box.ymin > ebut->y) {
726 box.ysiz = box.ymin - ebut->y + 1;
727 box.ymin = ebut->y;
728 } else {
729 box.ysiz = ebut->y - box.ymin + 1;
730 }
731 }
732
733
734 revbox(x0, y0, x1, y1) /* draw box with reversed lines */
735 int x0, y0, x1, y1;
736 {
737 revline(x0, y0, x1, y0);
738 revline(x0, y1, x1, y1);
739 revline(x0, y0, x0, y1);
740 revline(x1, y0, x1, y1);
741 }
742
743
744 avgbox(clr) /* average color over current box */
745 COLOR clr;
746 {
747 static COLOR lc;
748 static int ll, lr, lt, lb;
749 int left, right, top, bottom;
750 int y;
751 double d;
752 COLOR ctmp;
753 register int x;
754
755 setcolor(clr, 0.0, 0.0, 0.0);
756 left = box.xmin - xoff;
757 right = left + box.xsiz;
758 if (left < 0)
759 left = 0;
760 if (right > xmax)
761 right = xmax;
762 if (left >= right)
763 return(-1);
764 top = box.ymin - yoff;
765 bottom = top + box.ysiz;
766 if (top < 0)
767 top = 0;
768 if (bottom > ymax)
769 bottom = ymax;
770 if (top >= bottom)
771 return(-1);
772 if (left == ll && right == lr && top == lt && bottom == lb) {
773 copycolor(clr, lc);
774 return(0);
775 }
776 for (y = top; y < bottom; y++) {
777 if (getscan(y) == -1)
778 return(-1);
779 for (x = left; x < right; x++) {
780 colr_color(ctmp, scanline[x]);
781 addcolor(clr, ctmp);
782 }
783 }
784 d = 1.0/((right-left)*(bottom-top));
785 scalecolor(clr, d);
786 ll = left; lr = right; lt = top; lb = bottom;
787 copycolor(lc, clr);
788 return(0);
789 }
790
791
792 getmono() /* get monochrome data */
793 {
794 register unsigned char *dp;
795 register int x, err;
796 int y, errp;
797 short *cerr;
798
799 if ((cerr = (short *)calloc(xmax,sizeof(short))) == NULL)
800 quiterr("out of memory in getmono");
801 dp = ourdata - 1;
802 for (y = 0; y < ymax; y++) {
803 getscan(y);
804 add2icon(y, scanline);
805 normcolrs(scanline, xmax, scale);
806 err = 0;
807 for (x = 0; x < xmax; x++) {
808 if (!(x&7))
809 *++dp = 0;
810 errp = err;
811 err += normbright(scanline[x]) + cerr[x];
812 if (err > 127)
813 err -= 255;
814 else
815 *dp |= 1<<(7-(x&07));
816 err /= 3;
817 cerr[x] = err + errp;
818 }
819 }
820 free((char *)cerr);
821 }
822
823
824 add2icon(y, scan) /* add a scanline to our icon data */
825 int y;
826 COLR *scan;
827 {
828 static short cerr[ICONSIZ];
829 static int ynext;
830 static char *dp;
831 COLR clr;
832 register int err;
833 register int x, ti;
834 int errp;
835
836 if (iconheight == 0) { /* initialize */
837 if (xmax <= ICONSIZ && ymax <= ICONSIZ) {
838 iconwidth = xmax;
839 iconheight = ymax;
840 } else if (xmax > ymax) {
841 iconwidth = ICONSIZ;
842 iconheight = ICONSIZ*ymax/xmax;
843 if (iconheight < 1)
844 iconheight = 1;
845 } else {
846 iconwidth = ICONSIZ*xmax/ymax;
847 if (iconwidth < 1)
848 iconwidth = 1;
849 iconheight = ICONSIZ;
850 }
851 ynext = 0;
852 dp = icondata - 1;
853 }
854 if (y < ynext*ymax/iconheight) /* skip this one */
855 return;
856 err = 0;
857 for (x = 0; x < iconwidth; x++) {
858 if (!(x&7))
859 *++dp = 0;
860 errp = err;
861 ti = x*xmax/iconwidth;
862 copycolr(clr, scan[ti]);
863 normcolrs(clr, 1, scale);
864 err += normbright(clr) + cerr[x];
865 if (err > 127)
866 err -= 255;
867 else
868 *dp |= 1<<(x&07);
869 err /= 3;
870 cerr[x] = err + errp;
871 }
872 ynext++;
873 }
874
875
876 getfull() /* get full (24-bit) data */
877 {
878 int y;
879 register unsigned long *dp;
880 register int x;
881 /* set gamma correction */
882 setcolrgam(gamcor);
883 /* read and convert file */
884 dp = (unsigned long *)ourdata;
885 for (y = 0; y < ymax; y++) {
886 getscan(y);
887 add2icon(y, scanline);
888 if (scale)
889 shiftcolrs(scanline, xmax, scale);
890 colrs_gambs(scanline, xmax);
891 if (ourras->image->blue_mask & 1)
892 for (x = 0; x < xmax; x++)
893 *dp++ = scanline[x][RED] << 16 |
894 scanline[x][GRN] << 8 |
895 scanline[x][BLU] ;
896 else
897 for (x = 0; x < xmax; x++)
898 *dp++ = scanline[x][RED] |
899 scanline[x][GRN] << 8 |
900 scanline[x][BLU] << 16 ;
901 }
902 }
903
904
905 getgrey() /* get greyscale data */
906 {
907 int y;
908 register unsigned char *dp;
909 register int x;
910 /* set gamma correction */
911 setcolrgam(gamcor);
912 /* read and convert file */
913 dp = ourdata;
914 for (y = 0; y < ymax; y++) {
915 getscan(y);
916 add2icon(y, scanline);
917 if (scale)
918 shiftcolrs(scanline, xmax, scale);
919 for (x = 0; x < xmax; x++)
920 scanline[x][GRN] = normbright(scanline[x]);
921 colrs_gambs(scanline, xmax);
922 if (maxcolors < 256)
923 for (x = 0; x < xmax; x++)
924 *dp++ = ((long)scanline[x][GRN] *
925 maxcolors + maxcolors/2) >> 8;
926 else
927 for (x = 0; x < xmax; x++)
928 *dp++ = scanline[x][GRN];
929 }
930 for (x = 0; x < maxcolors; x++)
931 clrtab[x][RED] = clrtab[x][GRN] =
932 clrtab[x][BLU] = ((long)x*256 + 128)/maxcolors;
933 }
934
935
936 getmapped() /* get color-mapped data */
937 {
938 int y;
939 /* set gamma correction */
940 setcolrgam(gamcor);
941 /* make histogram */
942 new_histo();
943 for (y = 0; y < ymax; y++) {
944 if (getscan(y) < 0)
945 quiterr("seek error in getmapped");
946 add2icon(y, scanline);
947 if (scale)
948 shiftcolrs(scanline, xmax, scale);
949 colrs_gambs(scanline, xmax);
950 cnt_colrs(scanline, xmax);
951 }
952 /* map pixels */
953 if (!new_clrtab(maxcolors))
954 quiterr("cannot create color map");
955 for (y = 0; y < ymax; y++) {
956 if (getscan(y) < 0)
957 quiterr("seek error in getmapped");
958 if (scale)
959 shiftcolrs(scanline, xmax, scale);
960 colrs_gambs(scanline, xmax);
961 if (dither)
962 dith_colrs(ourdata+y*xmax, scanline, xmax);
963 else
964 map_colrs(ourdata+y*xmax, scanline, xmax);
965 }
966 }
967
968
969 scale_rcolors(xr, sf) /* scale color map */
970 register XRASTER *xr;
971 double sf;
972 {
973 register int i;
974 long maxv;
975
976 if (xr->pixels == NULL)
977 return;
978
979 sf = pow(sf, 1.0/gamcor);
980 maxv = 65535/sf;
981
982 for (i = xr->ncolors; i--; ) {
983 xr->cdefs[i].red = xr->cdefs[i].red > maxv ?
984 65535 :
985 xr->cdefs[i].red * sf;
986 xr->cdefs[i].green = xr->cdefs[i].green > maxv ?
987 65535 :
988 xr->cdefs[i].green * sf;
989 xr->cdefs[i].blue = xr->cdefs[i].blue > maxv ?
990 65535 :
991 xr->cdefs[i].blue * sf;
992 }
993 XStoreColors(thedisplay, xr->cmap, xr->cdefs, xr->ncolors);
994 }
995
996
997 getscan(y)
998 int y;
999 {
1000 if (y != cury) {
1001 if (scanpos == NULL || scanpos[y] == -1)
1002 return(-1);
1003 if (fseek(fin, scanpos[y], 0) == -1)
1004 quiterr("fseek error");
1005 cury = y;
1006 } else if (scanpos != NULL && scanpos[y] == -1)
1007 scanpos[y] = ftell(fin);
1008
1009 if (freadcolrs(scanline, xmax, fin) < 0)
1010 quiterr("read error");
1011
1012 cury++;
1013 return(0);
1014 }