ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/px/x11image.c
Revision: 2.40
Committed: Wed Nov 23 14:24:34 1994 UTC (29 years, 5 months ago) by greg
Content type: text/plain
Branch: MAIN
Changes since 2.39: +8 -4 lines
Log Message:
fixed unnecessary color table allocation

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