src/px/x11image.c

/* Copyright (c) 1990 Regents of the University of California */

#ifndef lint
static char SCCSid[] = "$SunId$ LBL";
#endif

/*
 *  x11image.c - driver for X-windows
 *
 *     3/1/90
 */

/*
 *  Modified for X11
 *
 *  January 1990
 *
 *  Anat Grynberg  and Greg Ward
 */


#include  "standard.h"

#include  <X11/Xlib.h>
#include  <X11/cursorfont.h>
#include  <X11/Xutil.h>

#include  "color.h"
#include  "view.h"
#include  "pic.h"
#include  "x11raster.h"
#include  "random.h"

#define  FONTNAME       "8x13"          /* text font we'll use */

#define  CTRL(c)        ('c'-'@')

#define  BORWIDTH       5               /* border width */

#define  ICONSIZ        (8*10)          /* maximum icon dimension (even 8) */

#define  ourscreen      DefaultScreen(thedisplay)
#define  ourblack       BlackPixel(thedisplay,ourscreen)
#define  ourwhite       WhitePixel(thedisplay,ourscreen)
#define  ourroot        RootWindow(thedisplay,ourscreen)
#define  ourgc          DefaultGC(thedisplay,ourscreen)

#define  revline(x0,y0,x1,y1)   XDrawLine(thedisplay,wind,revgc,x0,y0,x1,y1)

#define  redraw(x,y,w,h) patch_raster(wind,(x)-xoff,(y)-yoff,x,y,w,h,ourras)

double  gamcor = 2.2;                   /* gamma correction */

int  dither = 1;                        /* dither colors? */
int  fast = 0;                          /* keep picture in Pixmap? */

Window  wind = 0;                       /* our output window */
Font  fontid;                           /* our font */

int  maxcolors = 0;                     /* maximum colors */
int  greyscale = 0;                     /* in grey */

int  scale = 0;                         /* scalefactor; power of two */

int  xoff = 0;                          /* x image offset */
int  yoff = 0;                          /* y image offset */

VIEW  ourview = STDVIEW;                /* image view parameters */
int  gotview = 0;                       /* got parameters from file */

COLR  *scanline;                        /* scan line buffer */

int  xmax, ymax;                        /* picture resolution */
int  width, height;                     /* window size */
char  *fname = NULL;                    /* input file name */
FILE  *fin = stdin;                     /* input file */
long  *scanpos = NULL;                  /* scan line positions in file */
int  cury = 0;                          /* current scan location */

double  exposure = 1.0;                 /* exposure compensation used */

int  wrongformat = 0;                   /* input in another format? */

GC      revgc;                          /* graphics context with GXinvert */

XRASTER *ourras;                        /* our stored image */
unsigned char   *ourdata;               /* our image data */

struct {
        int  xmin, ymin, xsiz, ysiz;
}  box = {0, 0, 0, 0};                  /* current box */

char  *geometry = NULL;                 /* geometry specification */

char  icondata[ICONSIZ*ICONSIZ/8];      /* icon bitmap data */
int  iconwidth = 0, iconheight = 0;

char  *progname;

char  errmsg[128];

extern long  ftell();

extern char  *malloc(), *calloc();

extern double  atof(), pow(), log();

Display  *thedisplay;

main(argc, argv)
int  argc;
char  *argv[];
{
        int  headline();
        int  i;
        
        progname = argv[0];

        for (i = 1; i < argc; i++)
                if (argv[i][0] == '-')
                        switch (argv[i][1]) {
                        case 'c':
                                maxcolors = atoi(argv[++i]);
                                break;
                        case 'b':
                                greyscale = !greyscale;
                                break;
                        case 'm':
                                maxcolors = 2;
                                break;
                        case 'd':
                                dither = !dither;
                                break;
                        case 'f':
                                fast = !fast;
                                break;
                        case 'e':
                                if (argv[i+1][0] != '+' && argv[i+1][0] != '-')
                                        goto userr;
                                scale = atoi(argv[++i]);
                                break;
                        case 'g':
                                if (!strcmp(argv[i], "-geometry"))
                                        geometry = argv[++i];
                                else
                                        gamcor = atof(argv[++i]);
                                break;
                        default:
                                goto userr;
                        }
                else if (argv[i][0] == '=')
                        geometry = argv[i];
                else
                        break;

        if (i == argc-1) {
                fname = argv[i];
                fin = fopen(fname, "r");
                if (fin == NULL) {
                        sprintf(errmsg, "can't open file \"%s\"", fname);
                        quiterr(errmsg);
                }
        } else if (i != argc)
                goto userr;
                                /* get header */
        getheader(fin, headline, NULL);
                                /* get picture dimensions */
        if (wrongformat || fgetresolu(&xmax, &ymax, fin) != (YMAJOR|YDECR))
                quiterr("bad picture format");
                                /* set view parameters */
        if (gotview && setview(&ourview) != NULL)
                gotview = 0;
        if ((scanline = (COLR *)malloc(xmax*sizeof(COLR))) == NULL)
                quiterr("out of memory");

        init();                 /* get file and open window */

        for ( ; ; )
                getevent();             /* main loop */
userr:
        fprintf(stderr,
        "Usage: %s [-geometry spec][-b][-m][-d][-f][-c ncolors][-e +/-stops] file\n",
                        progname);
        quit(1);
}


headline(s)             /* get relevant info from header */
char  *s;
{
        static char  *altname[] = {"rview","rpict","pinterp",VIEWSTR,NULL};
        register char  **an;
        char  fmt[32];

        if (isexpos(s))
                exposure *= exposval(s);
        else if (isformat(s)) {
                formatval(fmt, s);
                wrongformat = strcmp(fmt, COLRFMT);
        } else
                for (an = altname; *an != NULL; an++)
                        if (!strncmp(*an, s, strlen(*an))) {
                                if (sscanview(&ourview, s+strlen(*an)) > 0)
                                        gotview++;
                                return;
                        }
}


init()                  /* get data and open window */
{
        XWMHints        ourxwmhints;
        XSetWindowAttributes    ourwinattr;
        XSizeHints  oursizhints;
        register int  i;
        
        if (fname != NULL) {
                scanpos = (long *)malloc(ymax*sizeof(long));
                if (scanpos == NULL)
                        goto memerr;
                for (i = 0; i < ymax; i++)
                        scanpos[i] = -1;
        }
        if ((thedisplay = XOpenDisplay(NULL)) == NULL)
                quiterr("can't open display; DISPLAY variable set?");
        if (maxcolors == 0) {           /* get number of available colors */
                i = DisplayPlanes(thedisplay,ourscreen);
                maxcolors = i > 8 ? 256 : 1<<i;
                if (maxcolors > 4) maxcolors -= 2;
        }
                                /* store image */
        getras();
                                /* open window */
        ourwinattr.border_pixel = ourblack;
        ourwinattr.background_pixel = ourwhite;
        wind = XCreateWindow(thedisplay, ourroot, 0, 0, xmax, ymax, BORWIDTH,
                        0, InputOutput, ourras->visual, 
                        CWBackPixel|CWBorderPixel, &ourwinattr);
        if (wind == 0)
                quiterr("can't create window");
        width = xmax;
        height = ymax;
        fontid = XLoadFont(thedisplay, FONTNAME);
        if (fontid == 0)
                quiterr("can't get font");
        XSetFont(thedisplay, ourgc, fontid);
        revgc = XCreateGC(thedisplay, wind, 0, 0);
        XSetFunction(thedisplay, revgc, GXinvert);
        XDefineCursor(thedisplay, wind, XCreateFontCursor(thedisplay, 
                        XC_diamond_cross));
        if (geometry != NULL) {
                bzero((char *)&oursizhints, sizeof(oursizhints));
                i = XParseGeometry(geometry, &oursizhints.x, &oursizhints.y,
                                (unsigned *)&oursizhints.width,
                                (unsigned *)&oursizhints.height);
                if ((i&(WidthValue|HeightValue)) == (WidthValue|HeightValue))
                        oursizhints.flags |= USSize;
                else {
                        oursizhints.width = xmax;
                        oursizhints.height = ymax;
                        oursizhints.flags |= PSize;
                }
                if ((i&(XValue|YValue)) == (XValue|YValue)) {
                        oursizhints.flags |= USPosition;
                        if (i & XNegative)
                                oursizhints.x += DisplayWidth(thedisplay,
                                ourscreen)-1-oursizhints.width-2*BORWIDTH;
                        if (i & YNegative)
                                oursizhints.y += DisplayHeight(thedisplay,
                                ourscreen)-1-oursizhints.height-2*BORWIDTH;
                }
                XSetNormalHints(thedisplay, wind, &oursizhints);
        }
        ourxwmhints.flags = InputHint|IconPixmapHint;
        ourxwmhints.input = True;
        ourxwmhints.icon_pixmap = XCreateBitmapFromData(thedisplay,
                        wind, icondata, iconwidth, iconheight);
        XSetWMHints(thedisplay, wind, &ourxwmhints);
        XSelectInput(thedisplay, wind, ButtonPressMask|ButtonReleaseMask
                        |ButtonMotionMask|StructureNotifyMask
                        |KeyPressMask|ExposureMask);
        XMapWindow(thedisplay, wind);
                                /* store name last as ready signal */
        XStoreName(thedisplay, wind, fname == NULL ? progname : fname);
        return;
memerr:
        quiterr("out of memory");
} /* end of init */


quiterr(err)            /* print message and exit */
char  *err;
{
        if (err != NULL) {
                fprintf(stderr, "%s: %s\n", progname, err);
                exit(1);
        }
        exit(0);
}


eputs(s)
char    *s;
{
        fputs(s, stderr);
}


quit(code)
int  code;
{
        exit(code);
}


getras()                                /* get raster file */
{
        colormap        ourmap;
        XVisualInfo     vinfo;

        if (maxcolors <= 2) {           /* monochrome */
                ourdata = (unsigned char *)malloc(ymax*((xmax+7)/8));
                if (ourdata == NULL)
                        goto fail;
                ourras = make_raster(thedisplay, ourscreen, 1, ourdata,
                                xmax, ymax, 8);
                if (ourras == NULL)
                        goto fail;
                getmono();
        } else if (XMatchVisualInfo(thedisplay,ourscreen,24,TrueColor,&vinfo)
                                                /* kludge for DirectColor */
        || XMatchVisualInfo(thedisplay,ourscreen,24,DirectColor,&vinfo)) {
                ourdata = (unsigned char *)malloc(xmax*ymax*3);
                if (ourdata == NULL)
                        goto fail;
                ourras = make_raster(thedisplay, ourscreen, 24, ourdata,
                                xmax, ymax, 8);
                if (ourras == NULL)
                        goto fail;
                getfull();
        } else {
                ourdata = (unsigned char *)malloc(xmax*ymax);
                if (ourdata == NULL)
                        goto fail;
                ourras = make_raster(thedisplay, ourscreen, 8, ourdata,
                                xmax, ymax, 8);
                if (ourras == NULL)
                        goto fail;
                if (greyscale)
                        biq(dither,maxcolors,1,ourmap);
                else
                        ciq(dither,maxcolors,1,ourmap);
                if (init_rcolors(ourras, ourmap[0], ourmap[1], ourmap[2]) == 0)
                        goto fail;
        }
        return;
fail:
        quiterr("could not create raster image");
}


getevent()                              /* process the next event */
{
        union {
                XEvent  u;
                XConfigureEvent  c;
                XExposeEvent  e;
                XButtonPressedEvent  b;
                XKeyPressedEvent  k;
        } e;

        XNextEvent(thedisplay, &e.u);
        switch (e.u.type) {
        case KeyPress:
                docom(&e.k);
                break;
        case ConfigureNotify:
                width = e.c.width;
                height = e.c.height;
                break;
        case MapNotify:
                map_rcolors(ourras, wind);
                if (fast)
                        make_rpixmap(ourras);
                break;
        case UnmapNotify:
                unmap_rcolors(ourras);
                break;
        case Expose:
                redraw(e.e.x, e.e.y, e.e.width, e.e.height);
                break;
        case ButtonPress:
                if (e.b.state & (ShiftMask|ControlMask))
                        moveimage(&e.b);
                else
                        getbox(&e.b);
                break;
        }
}


docom(ekey)                                     /* execute command */
XKeyPressedEvent  *ekey;
{
        char  buf[80];
        COLOR  cval;
        XColor  cvx;
        int  com, n;
        double  comp;
        FVECT  rorg, rdir;

        n = XLookupString(ekey, buf, sizeof(buf), NULL, NULL); 
        if (n == 0)
                return(0);
        com = buf[0];
        switch (com) {                  /* interpret command */
        case 'q':
        case CTRL(D):                           /* quit */
                quit(0);
        case '\n':
        case '\r':
        case 'l':
        case 'c':                               /* value */
                if (avgbox(cval) == -1)
                        return(-1);
                switch (com) {
                case '\n':
                case '\r':                              /* radiance */
                        sprintf(buf, "%.3f", intens(cval)/exposure);
                        break;
                case 'l':                               /* luminance */
                        sprintf(buf, "%.0fL", luminance(cval)/exposure);
                        break;
                case 'c':                               /* color */
                        comp = pow(2.0, (double)scale);
                        sprintf(buf, "(%.2f,%.2f,%.2f)",
                                        colval(cval,RED)*comp,
                                        colval(cval,GRN)*comp,
                                        colval(cval,BLU)*comp);
                        break;
                }
                XDrawImageString(thedisplay, wind, ourgc,
                                box.xmin, box.ymin+box.ysiz, buf, strlen(buf)); 
                return(0);
        case 'i':                               /* identify (contour) */
                if (ourras->pixels == NULL)
                        return(-1);
                n = ourdata[ekey->x-xoff+xmax*(ekey->y-yoff)];
                n = ourras->pmap[n];
                cvx.pixel = ourras->cdefs[n].pixel;
                cvx.red = random() & 65535;
                cvx.green = random() & 65535;
                cvx.blue = random() & 65535;
                cvx.flags = DoRed|DoGreen|DoBlue;
                XStoreColor(thedisplay, ourras->cmap, &cvx);
                return(0);
        case 'p':                               /* position */
                sprintf(buf, "(%d,%d)", ekey->x-xoff, ymax-1-ekey->y+yoff);
                XDrawImageString(thedisplay, wind, ourgc, ekey->x, ekey->y,
                                        buf, strlen(buf));
                return(0);
        case 't':                               /* trace */
                if (!gotview) {
                        XBell(thedisplay, 0);
                        return(-1);
                }
                if (viewray(rorg, rdir, &ourview,
                                (ekey->x-xoff+.5)/xmax,
                                (ymax-1-ekey->y+yoff+.5)/ymax) < 0)
                        return(-1);
                printf("%e %e %e ", rorg[0], rorg[1], rorg[2]);
                printf("%e %e %e\n", rdir[0], rdir[1], rdir[2]);
                fflush(stdout);
                return(0);
        case '=':                               /* adjust exposure */
                if (avgbox(cval) == -1)
                        return(-1);
                n = log(.5/bright(cval))/.69315 - scale;        /* truncate */
                if (n == 0)
                        return(0);
                scale_rcolors(ourras, pow(2.0, (double)n));
                scale += n;
                sprintf(buf, "%+d", scale);
                XDrawImageString(thedisplay, wind, ourgc,
                                box.xmin, box.ymin+box.ysiz, buf, strlen(buf));
                XFlush(thedisplay);
                free(ourdata);
                free_raster(ourras);
                getras();
        /* fall through */
        case CTRL(R):                           /* redraw */
        case CTRL(L):
                unmap_rcolors(ourras);
                XClearWindow(thedisplay, wind);
                map_rcolors(ourras, wind);
                if (fast)
                        make_rpixmap(ourras);
                redraw(0, 0, width, height);
                return(0);
        case ' ':                               /* clear */
                redraw(box.xmin, box.ymin, box.xsiz, box.ysiz);
                return(0);
        default:
                XBell(thedisplay, 0);
                return(-1);
        }
}


moveimage(ebut)                         /* shift the image */
XButtonPressedEvent  *ebut;
{
        union {
                XEvent  u;
                XButtonReleasedEvent  b;
                XPointerMovedEvent  m;
        }  e;
        int     mxo, myo;

        XMaskEvent(thedisplay, ButtonReleaseMask|ButtonMotionMask, &e.u);
        while (e.u.type == MotionNotify) {
                mxo = e.m.x;
                myo = e.m.y;
                revline(ebut->x, ebut->y, mxo, myo);
                revbox(xoff+mxo-ebut->x, yoff+myo-ebut->y,
                                xoff+mxo-ebut->x+xmax, yoff+myo-ebut->y+ymax);
                XMaskEvent(thedisplay,ButtonReleaseMask|ButtonMotionMask,&e.u);
                revline(ebut->x, ebut->y, mxo, myo);
                revbox(xoff+mxo-ebut->x, yoff+myo-ebut->y,
                                xoff+mxo-ebut->x+xmax, yoff+myo-ebut->y+ymax);
        }
        xoff += e.b.x - ebut->x;
        yoff += e.b.y - ebut->y;
        XClearWindow(thedisplay, wind);
        redraw(0, 0, width, height);
}


getbox(ebut)                            /* get new box */
XButtonPressedEvent  *ebut;
{
        union {
                XEvent  u;
                XButtonReleasedEvent  b;
                XPointerMovedEvent  m;
        }  e;

        XMaskEvent(thedisplay, ButtonReleaseMask|ButtonMotionMask, &e.u);
        while (e.u.type == MotionNotify) {
                revbox(ebut->x, ebut->y, box.xmin = e.m.x, box.ymin = e.m.y);
                XMaskEvent(thedisplay,ButtonReleaseMask|ButtonMotionMask,&e.u);
                revbox(ebut->x, ebut->y, box.xmin, box.ymin);
        }
        box.xmin = e.b.x<0 ? 0 : (e.b.x>=width ? width-1 : e.b.x);
        box.ymin = e.b.y<0 ? 0 : (e.b.y>=height ? height-1 : e.b.y);
        if (box.xmin > ebut->x) {
                box.xsiz = box.xmin - ebut->x + 1;
                box.xmin = ebut->x;
        } else {
                box.xsiz = ebut->x - box.xmin + 1;
        }
        if (box.ymin > ebut->y) {
                box.ysiz = box.ymin - ebut->y + 1;
                box.ymin = ebut->y;
        } else {
                box.ysiz = ebut->y - box.ymin + 1;
        }
}


revbox(x0, y0, x1, y1)                  /* draw box with reversed lines */
int  x0, y0, x1, y1;
{
        revline(x0, y0, x1, y0);
        revline(x0, y1, x1, y1);
        revline(x0, y0, x0, y1);
        revline(x1, y0, x1, y1);
}


avgbox(clr)                             /* average color over current box */
COLOR  clr;
{
        int  left, right, top, bottom;
        int  y;
        double  d;
        COLOR  ctmp;
        register int  x;

        setcolor(clr, 0.0, 0.0, 0.0);
        left = box.xmin - xoff;
        right = left + box.xsiz;
        if (left < 0)
                left = 0;
        if (right > xmax)
                right = xmax;
        if (left >= right)
                return(-1);
        top = box.ymin - yoff;
        bottom = top + box.ysiz;
        if (top < 0)
                top = 0;
        if (bottom > ymax)
                bottom = ymax;
        if (top >= bottom)
                return(-1);
        for (y = top; y < bottom; y++) {
                if (getscan(y) == -1)
                        return(-1);
                for (x = left; x < right; x++) {
                        colr_color(ctmp, scanline[x]);
                        addcolor(clr, ctmp);
                }
        }
        d = 1.0/((right-left)*(bottom-top));
        scalecolor(clr, d);
        return(0);
}


getmono()                       /* get monochrome data */
{
        register unsigned char  *dp;
        register int    x, err;
        int     y;
        short   *cerr;

        if ((cerr = (short *)calloc(xmax,sizeof(short))) == NULL)
                quiterr("out of memory in getmono");
        dp = ourdata - 1;
        for (y = 0; y < ymax; y++) {
                if (getscan(y) < 0)
                        quiterr("seek error in getmono");
                normcolrs(scanline, xmax, scale);
                add2icon(y, scanline);
                err = 0;
                for (x = 0; x < xmax; x++) {
                        if (!(x&7))
                                *++dp = 0;
                        err += normbright(scanline[x]) + cerr[x];
                        if (err > 127)
                                err -= 255;
                        else
                                *dp |= 1<<(7-(x&07));
                        cerr[x] = err >>= 1;
                }
        }
        free((char *)cerr);
}


add2icon(y, scan)               /* add a scanline to our icon data */
int  y;
COLR  *scan;
{
        static short  cerr[ICONSIZ];
        static int  ynext;
        static char  *dp;
        register int  err;
        register int    x, ti;

        if (iconheight == 0) {          /* initialize */
                if (xmax <= ICONSIZ && ymax <= ICONSIZ) {
                        iconwidth = xmax;
                        iconheight = ymax;
                } else if (xmax > ymax) {
                        iconwidth = ICONSIZ;
                        iconheight = ICONSIZ*ymax/xmax;
                } else {
                        iconwidth = ICONSIZ*xmax/ymax;
                        iconheight = ICONSIZ;
                }
                ynext = 0;
                dp = icondata - 1;
        }
        if (y < ynext*ymax/iconheight)  /* skip this one */
                return;
        err = 0;
        for (x = 0; x < iconwidth; x++) {
                if (!(x&7))
                        *++dp = 0;
                ti = x*xmax/iconwidth;
                err += normbright(scan[ti]) + cerr[x];
                if (err > 127)
                        err -= 255;
                else
                        *dp |= 1<<(x&07);
                cerr[x] = err >>= 1;
        }
        ynext++;
}


getfull()                       /* get full (24-bit) data */
{
        int     y;
        register unsigned char  *dp;
        register int    x;
                                        /* set gamma correction */
        setcolrgam(gamcor);
                                        /* read and convert file */
        dp = ourdata;
        for (y = 0; y < ymax; y++) {
                if (getscan(y) < 0)
                        quiterr("seek error in getfull");
                if (scale)
                        shiftcolrs(scanline, xmax, scale);
                colrs_gambs(scanline, xmax);
                add2icon(y, scanline);
                for (x = 0; x < xmax; x++) {
                        *dp++ = scanline[x][RED];
                        *dp++ = scanline[x][GRN];
                        *dp++ = scanline[x][BLU];
                }
        }
}


scale_rcolors(xr, sf)                   /* scale color map */
register XRASTER        *xr;
double  sf;
{
        register int    i;
        long    maxv;

        if (xr->pixels == NULL)
                return;

        sf = pow(sf, 1.0/gamcor);
        maxv = 65535/sf;

        for (i = xr->ncolors; i--; ) {
                xr->cdefs[i].red = xr->cdefs[i].red > maxv ?
                                65535 :
                                xr->cdefs[i].red * sf;
                xr->cdefs[i].green = xr->cdefs[i].green > maxv ?
                                65535 :
                                xr->cdefs[i].green * sf;
                xr->cdefs[i].blue = xr->cdefs[i].blue > maxv ?
                                65535 :
                                xr->cdefs[i].blue * sf;
        }
        XStoreColors(thedisplay, xr->cmap, xr->cdefs, xr->ncolors);
}


getscan(y)
int  y;
{
        if (y != cury) {
                if (scanpos == NULL || scanpos[y] == -1)
                        return(-1);
                if (fseek(fin, scanpos[y], 0) == -1)
                        quiterr("fseek error");
                cury = y;
        } else if (scanpos != NULL)
                scanpos[y] = ftell(fin);

        if (freadcolrs(scanline, xmax, fin) < 0)
                quiterr("read error");

        cury++;
        return(0);
}


picreadline3(y, l3)                     /* read in 3-byte scanline */
int  y;
register rgbpixel  *l3;
{
        register int    i;
                                                        /* read scanline */
        if (getscan(y) < 0)
                quiterr("cannot seek for picreadline");
                                                        /* convert scanline */
        normcolrs(scanline, xmax, scale);
        add2icon(y, scanline);
        for (i = 0; i < xmax; i++) {
                l3[i].r = scanline[i][RED];
                l3[i].g = scanline[i][GRN];
                l3[i].b = scanline[i][BLU];
        }
}


picwriteline(y, l)              /* add 8-bit scanline to image */
int  y;
pixel  *l;
{
        bcopy((char *)l, (char *)ourdata+y*xmax, xmax);
}


picreadcm(map)                  /* do gamma correction */
colormap  map;
{
        extern double  pow();
        register int  i, val;

        for (i = 0; i < 256; i++) {
                val = pow((i+0.5)/256.0, 1.0/gamcor) * 256.0;
                map[0][i] = map[1][i] = map[2][i] = val;
        }
}


picwritecm(map)                 /* handled elsewhere */
colormap  map;
{
#ifdef DEBUG
        register int i;

        for (i = 0; i < 256; i++)
                printf("%d %d %d\n", map[0][i],map[1][i],map[2][i]);
#endif
}
Revision:	1.20
Committed:	Fri May 3 10:54:47 1991 UTC (34 years, 6 months ago) by greg
Content type:	text/plain
Branch:	MAIN
Changes since 1.19:	+8 -9 lines
Log Message:	fixed major bug in add2icon()
#	User	Rev	Content
1	greg	1.1	/* Copyright (c) 1990 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
28			#include "color.h"
29			#include "view.h"
30			#include "pic.h"
31			#include "x11raster.h"
32			#include "random.h"
33
34			#define FONTNAME "8x13" /* text font we'll use */
35
36			#define CTRL(c) ('c'-'@')
37
38			#define BORWIDTH 5 /* border width */
39
40	greg	1.18	#define ICONSIZ (810) / maximum icon dimension (even 8) */
41	greg	1.17
42	greg	1.1	#define ourscreen DefaultScreen(thedisplay)
43			#define ourblack BlackPixel(thedisplay,ourscreen)
44			#define ourwhite WhitePixel(thedisplay,ourscreen)
45			#define ourroot RootWindow(thedisplay,ourscreen)
46			#define ourgc DefaultGC(thedisplay,ourscreen)
47
48	greg	1.5	#define revline(x0,y0,x1,y1) XDrawLine(thedisplay,wind,revgc,x0,y0,x1,y1)
49
50	greg	1.2	#define redraw(x,y,w,h) patch_raster(wind,(x)-xoff,(y)-yoff,x,y,w,h,ourras)
51
52	greg	1.1	double gamcor = 2.2; /* gamma correction */
53
54			int dither = 1; /* dither colors? */
55			int fast = 0; /* keep picture in Pixmap? */
56
57			Window wind = 0; /* our output window */
58			Font fontid; /* our font */
59
60	greg	1.2	int maxcolors = 0; /* maximum colors */
61	greg	1.1	int greyscale = 0; /* in grey */
62
63			int scale = 0; /* scalefactor; power of two */
64
65			int xoff = 0; /* x image offset */
66			int yoff = 0; /* y image offset */
67
68			VIEW ourview = STDVIEW; /* image view parameters */
69			int gotview = 0; /* got parameters from file */
70
71			COLR scanline; / scan line buffer */
72
73			int xmax, ymax; /* picture resolution */
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	greg	1.14	int wrongformat = 0; /* input in another format? */
83
84	greg	1.5	GC revgc; /* graphics context with GXinvert */
85
86	greg	1.1	XRASTER ourras; / our stored image */
87			unsigned char ourdata; / our image data */
88
89			struct {
90			int xmin, ymin, xsiz, ysiz;
91			} box = {0, 0, 0, 0}; /* current box */
92
93			char geometry = NULL; / geometry specification */
94
95	greg	1.18	char icondata[ICONSIZICONSIZ/8]; / icon bitmap data */
96	greg	1.17	int iconwidth = 0, iconheight = 0;
97
98	greg	1.1	char *progname;
99
100			char errmsg[128];
101
102			extern long ftell();
103
104			extern char malloc(), calloc();
105
106			extern double atof(), pow(), log();
107
108			Display *thedisplay;
109
110			main(argc, argv)
111			int argc;
112			char *argv[];
113			{
114			int headline();
115			int i;
116
117			progname = argv[0];
118
119			for (i = 1; i < argc; i++)
120			if (argv[i][0] == '-')
121			switch (argv[i][1]) {
122			case 'c':
123			maxcolors = atoi(argv[++i]);
124			break;
125			case 'b':
126			greyscale = !greyscale;
127			break;
128			case 'm':
129			maxcolors = 2;
130			break;
131			case 'd':
132			dither = !dither;
133			break;
134			case 'f':
135			fast = !fast;
136			break;
137			case 'e':
138			if (argv[i+1][0] != '+' && argv[i+1][0] != '-')
139			goto userr;
140			scale = atoi(argv[++i]);
141			break;
142			case 'g':
143			if (!strcmp(argv[i], "-geometry"))
144			geometry = argv[++i];
145			else
146			gamcor = atof(argv[++i]);
147			break;
148			default:
149			goto userr;
150			}
151	greg	1.3	else if (argv[i][0] == '=')
152			geometry = argv[i];
153	greg	1.1	else
154			break;
155
156	greg	1.3	if (i == argc-1) {
157	greg	1.1	fname = argv[i];
158			fin = fopen(fname, "r");
159			if (fin == NULL) {
160			sprintf(errmsg, "can't open file \"%s\"", fname);
161			quiterr(errmsg);
162			}
163	greg	1.3	} else if (i != argc)
164			goto userr;
165	greg	1.1	/* get header */
166	greg	1.14	getheader(fin, headline, NULL);
167	greg	1.1	/* get picture dimensions */
168	greg	1.14	if (wrongformat \|\| fgetresolu(&xmax, &ymax, fin) != (YMAJOR\|YDECR))
169			quiterr("bad picture format");
170	greg	1.1	/* set view parameters */
171			if (gotview && setview(&ourview) != NULL)
172			gotview = 0;
173			if ((scanline = (COLR )malloc(xmaxsizeof(COLR))) == NULL)
174			quiterr("out of memory");
175
176			init(); /* get file and open window */
177
178			for ( ; ; )
179			getevent(); /* main loop */
180			userr:
181			fprintf(stderr,
182			"Usage: %s [-geometry spec][-b][-m][-d][-f][-c ncolors][-e +/-stops] file\n",
183			progname);
184			quit(1);
185			}
186
187
188			headline(s) /* get relevant info from header */
189			char *s;
190			{
191			static char *altname[] = {"rview","rpict","pinterp",VIEWSTR,NULL};
192			register char **an;
193	greg	1.14	char fmt[32];
194	greg	1.1
195	greg	1.2	if (isexpos(s))
196			exposure *= exposval(s);
197	greg	1.14	else if (isformat(s)) {
198			formatval(fmt, s);
199			wrongformat = strcmp(fmt, COLRFMT);
200			} else
201	greg	1.1	for (an = altname; *an != NULL; an++)
202			if (!strncmp(an, s, strlen(an))) {
203	greg	1.2	if (sscanview(&ourview, s+strlen(*an)) > 0)
204	greg	1.1	gotview++;
205			return;
206			}
207			}
208
209
210			init() /* get data and open window */
211			{
212	greg	1.15	XWMHints ourxwmhints;
213	greg	1.4	XSetWindowAttributes ourwinattr;
214	greg	1.1	XSizeHints oursizhints;
215			register int i;
216
217			if (fname != NULL) {
218			scanpos = (long )malloc(ymaxsizeof(long));
219			if (scanpos == NULL)
220			goto memerr;
221			for (i = 0; i < ymax; i++)
222			scanpos[i] = -1;
223			}
224			if ((thedisplay = XOpenDisplay(NULL)) == NULL)
225			quiterr("can't open display; DISPLAY variable set?");
226			if (maxcolors == 0) { /* get number of available colors */
227	greg	1.2	i = DisplayPlanes(thedisplay,ourscreen);
228			maxcolors = i > 8 ? 256 : 1<<i;
229	greg	1.1	if (maxcolors > 4) maxcolors -= 2;
230			}
231	greg	1.4	/* store image */
232			getras();
233			/* open window */
234			ourwinattr.border_pixel = ourblack;
235			ourwinattr.background_pixel = ourwhite;
236			wind = XCreateWindow(thedisplay, ourroot, 0, 0, xmax, ymax, BORWIDTH,
237			0, InputOutput, ourras->visual,
238			CWBackPixel\|CWBorderPixel, &ourwinattr);
239			if (wind == 0)
240			quiterr("can't create window");
241	greg	1.5	width = xmax;
242			height = ymax;
243	greg	1.1	fontid = XLoadFont(thedisplay, FONTNAME);
244			if (fontid == 0)
245			quiterr("can't get font");
246			XSetFont(thedisplay, ourgc, fontid);
247	greg	1.5	revgc = XCreateGC(thedisplay, wind, 0, 0);
248			XSetFunction(thedisplay, revgc, GXinvert);
249	greg	1.1	XDefineCursor(thedisplay, wind, XCreateFontCursor(thedisplay,
250			XC_diamond_cross));
251			if (geometry != NULL) {
252			bzero((char *)&oursizhints, sizeof(oursizhints));
253			i = XParseGeometry(geometry, &oursizhints.x, &oursizhints.y,
254	greg	1.12	(unsigned *)&oursizhints.width,
255			(unsigned *)&oursizhints.height);
256	greg	1.3	if ((i&(WidthValue\|HeightValue)) == (WidthValue\|HeightValue))
257			oursizhints.flags \|= USSize;
258			else {
259			oursizhints.width = xmax;
260			oursizhints.height = ymax;
261			oursizhints.flags \|= PSize;
262			}
263			if ((i&(XValue\|YValue)) == (XValue\|YValue)) {
264	greg	1.1	oursizhints.flags \|= USPosition;
265			if (i & XNegative)
266	greg	1.3	oursizhints.x += DisplayWidth(thedisplay,
267			ourscreen)-1-oursizhints.width-2*BORWIDTH;
268	greg	1.1	if (i & YNegative)
269	greg	1.3	oursizhints.y += DisplayHeight(thedisplay,
270			ourscreen)-1-oursizhints.height-2*BORWIDTH;
271	greg	1.1	}
272			XSetNormalHints(thedisplay, wind, &oursizhints);
273			}
274	greg	1.15	ourxwmhints.flags = InputHint\|IconPixmapHint;
275			ourxwmhints.input = True;
276			ourxwmhints.icon_pixmap = XCreateBitmapFromData(thedisplay,
277	greg	1.17	wind, icondata, iconwidth, iconheight);
278	greg	1.15	XSetWMHints(thedisplay, wind, &ourxwmhints);
279	greg	1.2	XSelectInput(thedisplay, wind, ButtonPressMask\|ButtonReleaseMask
280			\|ButtonMotionMask\|StructureNotifyMask
281			\|KeyPressMask\|ExposureMask);
282	greg	1.1	XMapWindow(thedisplay, wind);
283	greg	1.19	/* store name last as ready signal */
284			XStoreName(thedisplay, wind, fname == NULL ? progname : fname);
285	greg	1.1	return;
286			memerr:
287			quiterr("out of memory");
288			} /* end of init */
289
290
291			quiterr(err) /* print message and exit */
292			char *err;
293			{
294			if (err != NULL) {
295			fprintf(stderr, "%s: %s\n", progname, err);
296			exit(1);
297			}
298			exit(0);
299			}
300
301
302			eputs(s)
303			char *s;
304			{
305			fputs(s, stderr);
306			}
307
308
309			quit(code)
310			int code;
311			{
312			exit(code);
313			}
314
315
316			getras() /* get raster file */
317			{
318			colormap ourmap;
319			XVisualInfo vinfo;
320
321			if (maxcolors <= 2) { /* monochrome */
322			ourdata = (unsigned char )malloc(ymax((xmax+7)/8));
323			if (ourdata == NULL)
324			goto fail;
325			ourras = make_raster(thedisplay, ourscreen, 1, ourdata,
326			xmax, ymax, 8);
327			if (ourras == NULL)
328			goto fail;
329			getmono();
330	greg	1.7	} else if (XMatchVisualInfo(thedisplay,ourscreen,24,TrueColor,&vinfo)
331			/* kludge for DirectColor */
332			\|\| XMatchVisualInfo(thedisplay,ourscreen,24,DirectColor,&vinfo)) {
333	greg	1.1	ourdata = (unsigned char )malloc(xmaxymax*3);
334			if (ourdata == NULL)
335			goto fail;
336			ourras = make_raster(thedisplay, ourscreen, 24, ourdata,
337			xmax, ymax, 8);
338			if (ourras == NULL)
339			goto fail;
340			getfull();
341			} else {
342			ourdata = (unsigned char )malloc(xmaxymax);
343			if (ourdata == NULL)
344			goto fail;
345			ourras = make_raster(thedisplay, ourscreen, 8, ourdata,
346			xmax, ymax, 8);
347			if (ourras == NULL)
348			goto fail;
349			if (greyscale)
350			biq(dither,maxcolors,1,ourmap);
351			else
352			ciq(dither,maxcolors,1,ourmap);
353	greg	1.3	if (init_rcolors(ourras, ourmap[0], ourmap[1], ourmap[2]) == 0)
354	greg	1.1	goto fail;
355			}
356			return;
357			fail:
358	greg	1.9	quiterr("could not create raster image");
359	greg	1.1	}
360
361
362			getevent() /* process the next event */
363			{
364			union {
365			XEvent u;
366			XConfigureEvent c;
367			XExposeEvent e;
368			XButtonPressedEvent b;
369			XKeyPressedEvent k;
370			} e;
371
372			XNextEvent(thedisplay, &e.u);
373			switch (e.u.type) {
374			case KeyPress:
375			docom(&e.k);
376			break;
377			case ConfigureNotify:
378			width = e.c.width;
379			height = e.c.height;
380			break;
381			case MapNotify:
382			map_rcolors(ourras, wind);
383			if (fast)
384			make_rpixmap(ourras);
385			break;
386			case UnmapNotify:
387			unmap_rcolors(ourras);
388			break;
389			case Expose:
390			redraw(e.e.x, e.e.y, e.e.width, e.e.height);
391			break;
392			case ButtonPress:
393			if (e.b.state & (ShiftMask\|ControlMask))
394			moveimage(&e.b);
395			else
396			getbox(&e.b);
397			break;
398			}
399			}
400
401
402			docom(ekey) /* execute command */
403			XKeyPressedEvent *ekey;
404			{
405			char buf[80];
406			COLOR cval;
407			XColor cvx;
408			int com, n;
409			double comp;
410			FVECT rorg, rdir;
411
412			n = XLookupString(ekey, buf, sizeof(buf), NULL, NULL);
413			if (n == 0)
414			return(0);
415			com = buf[0];
416			switch (com) { /* interpret command */
417			case 'q':
418			case CTRL(D): /* quit */
419			quit(0);
420			case '\n':
421			case '\r':
422			case 'l':
423			case 'c': /* value */
424			if (avgbox(cval) == -1)
425			return(-1);
426			switch (com) {
427			case '\n':
428			case '\r': /* radiance */
429			sprintf(buf, "%.3f", intens(cval)/exposure);
430			break;
431			case 'l': /* luminance */
432	greg	1.16	sprintf(buf, "%.0fL", luminance(cval)/exposure);
433	greg	1.1	break;
434			case 'c': /* color */
435			comp = pow(2.0, (double)scale);
436			sprintf(buf, "(%.2f,%.2f,%.2f)",
437			colval(cval,RED)*comp,
438			colval(cval,GRN)*comp,
439			colval(cval,BLU)*comp);
440			break;
441			}
442	greg	1.3	XDrawImageString(thedisplay, wind, ourgc,
443			box.xmin, box.ymin+box.ysiz, buf, strlen(buf));
444	greg	1.1	return(0);
445			case 'i': /* identify (contour) */
446			if (ourras->pixels == NULL)
447			return(-1);
448			n = ourdata[ekey->x-xoff+xmax*(ekey->y-yoff)];
449			n = ourras->pmap[n];
450			cvx.pixel = ourras->cdefs[n].pixel;
451			cvx.red = random() & 65535;
452			cvx.green = random() & 65535;
453			cvx.blue = random() & 65535;
454	greg	1.2	cvx.flags = DoRed\|DoGreen\|DoBlue;
455			XStoreColor(thedisplay, ourras->cmap, &cvx);
456	greg	1.1	return(0);
457			case 'p': /* position */
458			sprintf(buf, "(%d,%d)", ekey->x-xoff, ymax-1-ekey->y+yoff);
459			XDrawImageString(thedisplay, wind, ourgc, ekey->x, ekey->y,
460			buf, strlen(buf));
461			return(0);
462			case 't': /* trace */
463			if (!gotview) {
464			XBell(thedisplay, 0);
465			return(-1);
466			}
467	greg	1.8	if (viewray(rorg, rdir, &ourview,
468	greg	1.1	(ekey->x-xoff+.5)/xmax,
469	greg	1.8	(ymax-1-ekey->y+yoff+.5)/ymax) < 0)
470			return(-1);
471	greg	1.1	printf("%e %e %e ", rorg[0], rorg[1], rorg[2]);
472			printf("%e %e %e\n", rdir[0], rdir[1], rdir[2]);
473			fflush(stdout);
474			return(0);
475			case '=': /* adjust exposure */
476			if (avgbox(cval) == -1)
477			return(-1);
478			n = log(.5/bright(cval))/.69315 - scale; /* truncate */
479			if (n == 0)
480			return(0);
481			scale_rcolors(ourras, pow(2.0, (double)n));
482			scale += n;
483			sprintf(buf, "%+d", scale);
484	greg	1.3	XDrawImageString(thedisplay, wind, ourgc,
485			box.xmin, box.ymin+box.ysiz, buf, strlen(buf));
486	greg	1.1	XFlush(thedisplay);
487			free(ourdata);
488			free_raster(ourras);
489			getras();
490			/* fall through */
491			case CTRL(R): /* redraw */
492			case CTRL(L):
493			unmap_rcolors(ourras);
494			XClearWindow(thedisplay, wind);
495			map_rcolors(ourras, wind);
496			if (fast)
497			make_rpixmap(ourras);
498			redraw(0, 0, width, height);
499			return(0);
500			case ' ': /* clear */
501			redraw(box.xmin, box.ymin, box.xsiz, box.ysiz);
502			return(0);
503			default:
504	greg	1.2	XBell(thedisplay, 0);
505	greg	1.1	return(-1);
506			}
507			}
508
509
510	greg	1.2	moveimage(ebut) /* shift the image */
511			XButtonPressedEvent *ebut;
512	greg	1.1	{
513	greg	1.2	union {
514	greg	1.3	XEvent u;
515			XButtonReleasedEvent b;
516			XPointerMovedEvent m;
517			} e;
518	greg	1.5	int mxo, myo;
519	greg	1.1
520	greg	1.3	XMaskEvent(thedisplay, ButtonReleaseMask\|ButtonMotionMask, &e.u);
521			while (e.u.type == MotionNotify) {
522	greg	1.5	mxo = e.m.x;
523			myo = e.m.y;
524			revline(ebut->x, ebut->y, mxo, myo);
525			revbox(xoff+mxo-ebut->x, yoff+myo-ebut->y,
526			xoff+mxo-ebut->x+xmax, yoff+myo-ebut->y+ymax);
527	greg	1.3	XMaskEvent(thedisplay,ButtonReleaseMask\|ButtonMotionMask,&e.u);
528	greg	1.5	revline(ebut->x, ebut->y, mxo, myo);
529			revbox(xoff+mxo-ebut->x, yoff+myo-ebut->y,
530			xoff+mxo-ebut->x+xmax, yoff+myo-ebut->y+ymax);
531	greg	1.3	}
532	greg	1.2	xoff += e.b.x - ebut->x;
533			yoff += e.b.y - ebut->y;
534	greg	1.1	XClearWindow(thedisplay, wind);
535			redraw(0, 0, width, height);
536			}
537
538
539			getbox(ebut) /* get new box */
540			XButtonPressedEvent *ebut;
541			{
542			union {
543	greg	1.2	XEvent u;
544	greg	1.1	XButtonReleasedEvent b;
545			XPointerMovedEvent m;
546			} e;
547
548	greg	1.2	XMaskEvent(thedisplay, ButtonReleaseMask\|ButtonMotionMask, &e.u);
549			while (e.u.type == MotionNotify) {
550	greg	1.1	revbox(ebut->x, ebut->y, box.xmin = e.m.x, box.ymin = e.m.y);
551	greg	1.2	XMaskEvent(thedisplay,ButtonReleaseMask\|ButtonMotionMask,&e.u);
552	greg	1.1	revbox(ebut->x, ebut->y, box.xmin, box.ymin);
553			}
554			box.xmin = e.b.x<0 ? 0 : (e.b.x>=width ? width-1 : e.b.x);
555			box.ymin = e.b.y<0 ? 0 : (e.b.y>=height ? height-1 : e.b.y);
556			if (box.xmin > ebut->x) {
557			box.xsiz = box.xmin - ebut->x + 1;
558			box.xmin = ebut->x;
559			} else {
560			box.xsiz = ebut->x - box.xmin + 1;
561			}
562			if (box.ymin > ebut->y) {
563			box.ysiz = box.ymin - ebut->y + 1;
564			box.ymin = ebut->y;
565			} else {
566			box.ysiz = ebut->y - box.ymin + 1;
567			}
568			}
569
570
571			revbox(x0, y0, x1, y1) /* draw box with reversed lines */
572			int x0, y0, x1, y1;
573			{
574	greg	1.5	revline(x0, y0, x1, y0);
575			revline(x0, y1, x1, y1);
576			revline(x0, y0, x0, y1);
577			revline(x1, y0, x1, y1);
578			}
579	greg	1.1
580
581			avgbox(clr) /* average color over current box */
582			COLOR clr;
583			{
584			int left, right, top, bottom;
585			int y;
586			double d;
587			COLOR ctmp;
588			register int x;
589
590			setcolor(clr, 0.0, 0.0, 0.0);
591			left = box.xmin - xoff;
592			right = left + box.xsiz;
593			if (left < 0)
594			left = 0;
595			if (right > xmax)
596			right = xmax;
597			if (left >= right)
598			return(-1);
599			top = box.ymin - yoff;
600			bottom = top + box.ysiz;
601			if (top < 0)
602			top = 0;
603			if (bottom > ymax)
604			bottom = ymax;
605			if (top >= bottom)
606			return(-1);
607			for (y = top; y < bottom; y++) {
608			if (getscan(y) == -1)
609			return(-1);
610			for (x = left; x < right; x++) {
611			colr_color(ctmp, scanline[x]);
612			addcolor(clr, ctmp);
613			}
614			}
615			d = 1.0/((right-left)*(bottom-top));
616			scalecolor(clr, d);
617			return(0);
618			}
619
620
621			getmono() /* get monochrome data */
622			{
623			register unsigned char *dp;
624			register int x, err;
625			int y;
626			short *cerr;
627
628	greg	1.10	if ((cerr = (short *)calloc(xmax,sizeof(short))) == NULL)
629	greg	1.9	quiterr("out of memory in getmono");
630	greg	1.1	dp = ourdata - 1;
631			for (y = 0; y < ymax; y++) {
632	greg	1.10	if (getscan(y) < 0)
633			quiterr("seek error in getmono");
634			normcolrs(scanline, xmax, scale);
635	greg	1.17	add2icon(y, scanline);
636	greg	1.1	err = 0;
637			for (x = 0; x < xmax; x++) {
638			if (!(x&7))
639			*++dp = 0;
640	greg	1.10	err += normbright(scanline[x]) + cerr[x];
641	greg	1.1	if (err > 127)
642			err -= 255;
643			else
644	greg	1.3	*dp \|= 1<<(7-(x&07));
645	greg	1.1	cerr[x] = err >>= 1;
646			}
647			}
648			free((char *)cerr);
649			}
650
651
652	greg	1.17	add2icon(y, scan) /* add a scanline to our icon data */
653			int y;
654			COLR *scan;
655			{
656			static short cerr[ICONSIZ];
657	greg	1.20	static int ynext;
658			static char *dp;
659	greg	1.17	register int err;
660	greg	1.20	register int x, ti;
661	greg	1.17
662			if (iconheight == 0) { /* initialize */
663	greg	1.18	if (xmax <= ICONSIZ && ymax <= ICONSIZ) {
664	greg	1.17	iconwidth = xmax;
665			iconheight = ymax;
666			} else if (xmax > ymax) {
667			iconwidth = ICONSIZ;
668			iconheight = ICONSIZ*ymax/xmax;
669			} else {
670			iconwidth = ICONSIZ*xmax/ymax;
671			iconheight = ICONSIZ;
672			}
673	greg	1.20	ynext = 0;
674	greg	1.17	dp = icondata - 1;
675			}
676	greg	1.20	if (y < ynextymax/iconheight) / skip this one */
677	greg	1.17	return;
678			err = 0;
679			for (x = 0; x < iconwidth; x++) {
680			if (!(x&7))
681			*++dp = 0;
682	greg	1.20	ti = x*xmax/iconwidth;
683			err += normbright(scan[ti]) + cerr[x];
684	greg	1.17	if (err > 127)
685			err -= 255;
686			else
687			*dp \|= 1<<(x&07);
688			cerr[x] = err >>= 1;
689			}
690	greg	1.20	ynext++;
691	greg	1.17	}
692
693
694	greg	1.1	getfull() /* get full (24-bit) data */
695			{
696			int y;
697	greg	1.10	register unsigned char *dp;
698			register int x;
699			/* set gamma correction */
700			setcolrgam(gamcor);
701			/* read and convert file */
702			dp = ourdata;
703			for (y = 0; y < ymax; y++) {
704			if (getscan(y) < 0)
705			quiterr("seek error in getfull");
706			if (scale)
707			shiftcolrs(scanline, xmax, scale);
708			colrs_gambs(scanline, xmax);
709	greg	1.17	add2icon(y, scanline);
710	greg	1.10	for (x = 0; x < xmax; x++) {
711			*dp++ = scanline[x][RED];
712			*dp++ = scanline[x][GRN];
713			*dp++ = scanline[x][BLU];
714			}
715			}
716	greg	1.1	}
717
718
719			scale_rcolors(xr, sf) /* scale color map */
720			register XRASTER *xr;
721			double sf;
722			{
723			register int i;
724			long maxv;
725
726			if (xr->pixels == NULL)
727			return;
728
729			sf = pow(sf, 1.0/gamcor);
730			maxv = 65535/sf;
731
732			for (i = xr->ncolors; i--; ) {
733			xr->cdefs[i].red = xr->cdefs[i].red > maxv ?
734			65535 :
735			xr->cdefs[i].red * sf;
736			xr->cdefs[i].green = xr->cdefs[i].green > maxv ?
737			65535 :
738			xr->cdefs[i].green * sf;
739			xr->cdefs[i].blue = xr->cdefs[i].blue > maxv ?
740			65535 :
741			xr->cdefs[i].blue * sf;
742			}
743			XStoreColors(thedisplay, xr->cmap, xr->cdefs, xr->ncolors);
744			}
745
746
747			getscan(y)
748			int y;
749			{
750			if (y != cury) {
751			if (scanpos == NULL \|\| scanpos[y] == -1)
752			return(-1);
753			if (fseek(fin, scanpos[y], 0) == -1)
754	greg	1.9	quiterr("fseek error");
755	greg	1.1	cury = y;
756			} else if (scanpos != NULL)
757			scanpos[y] = ftell(fin);
758
759			if (freadcolrs(scanline, xmax, fin) < 0)
760			quiterr("read error");
761
762			cury++;
763			return(0);
764			}
765
766
767			picreadline3(y, l3) /* read in 3-byte scanline */
768			int y;
769			register rgbpixel *l3;
770			{
771			register int i;
772			/* read scanline */
773			if (getscan(y) < 0)
774			quiterr("cannot seek for picreadline");
775			/* convert scanline */
776			normcolrs(scanline, xmax, scale);
777	greg	1.17	add2icon(y, scanline);
778	greg	1.1	for (i = 0; i < xmax; i++) {
779			l3[i].r = scanline[i][RED];
780			l3[i].g = scanline[i][GRN];
781			l3[i].b = scanline[i][BLU];
782			}
783			}
784
785
786			picwriteline(y, l) /* add 8-bit scanline to image */
787			int y;
788			pixel *l;
789			{
790			bcopy((char )l, (char )ourdata+y*xmax, xmax);
791			}
792
793
794			picreadcm(map) /* do gamma correction */
795			colormap map;
796			{
797			extern double pow();
798			register int i, val;
799
800			for (i = 0; i < 256; i++) {
801	greg	1.11	val = pow((i+0.5)/256.0, 1.0/gamcor) * 256.0;
802	greg	1.1	map[0][i] = map[1][i] = map[2][i] = val;
803			}
804			}
805
806
807			picwritecm(map) /* handled elsewhere */
808			colormap map;
809			{
810			#ifdef DEBUG
811			register int i;
812
813			for (i = 0; i < 256; i++)
814			printf("%d %d %d\n", map[0][i],map[1][i],map[2][i]);
815			#endif
816			}