src/px/x11image.c

/* Copyright (c) 1991 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));
        XStoreName(thedisplay, wind, fname == NULL ? progname : fname);
        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);
        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;
{
        static COLOR  lc;
        static int  ll, lr, lt, lb;
        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);
        if (left == ll && right == lr && top == lt && bottom == lb) {
                copycolor(clr, lc);
                return;
        }
        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);
        ll = left; lr = right; lt = top; lb = bottom;
        copycolor(lc, clr);
        return(0);
}


getmono()                       /* get monochrome data */
{
        register unsigned char  *dp;
        register int    x, err;
        int     y, errp;
        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;
                        errp = err;
                        err += normbright(scanline[x]) + cerr[x];
                        if (err > 127)
                                err -= 255;
                        else
                                *dp |= 1<<(7-(x&07));
                        err /= 3;
                        cerr[x] = err + errp;
                }
        }
        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;
        int  errp;

        if (iconheight == 0) {          /* initialize */
                if (xmax <= ICONSIZ && ymax <= ICONSIZ) {
                        iconwidth = xmax;
                        iconheight = ymax;
                } else if (xmax > ymax) {
                        iconwidth = ICONSIZ;
                        iconheight = ICONSIZ*ymax/xmax;
                        if (iconheight < 1)
                                iconheight = 1;
                } else {
                        iconwidth = ICONSIZ*xmax/ymax;
                        if (iconwidth < 1)
                                iconwidth = 1;
                        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;
                errp = err;
                ti = x*xmax/iconwidth;
                err += normbright(scan[ti]) + cerr[x];
                if (err > 127)
                        err -= 255;
                else
                        *dp |= 1<<(x&07);
                err /= 3;
                cerr[x] = err + errp;
        }
        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.25
Committed:	Mon Aug 26 10:15:01 1991 UTC (32 years, 8 months ago) by greg
Content type:	text/plain
Branch:	MAIN
Changes since 1.24:	+9 -1 lines
Log Message:	added check for duplicate box measurement
#	User	Rev	Content
1	greg	1.25	/* Copyright (c) 1991 Regents of the University of California */
2	greg	1.1
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	greg	1.22	XStoreName(thedisplay, wind, fname == NULL ? progname : fname);
252	greg	1.1	if (geometry != NULL) {
253			bzero((char *)&oursizhints, sizeof(oursizhints));
254			i = XParseGeometry(geometry, &oursizhints.x, &oursizhints.y,
255	greg	1.12	(unsigned *)&oursizhints.width,
256			(unsigned *)&oursizhints.height);
257	greg	1.3	if ((i&(WidthValue\|HeightValue)) == (WidthValue\|HeightValue))
258			oursizhints.flags \|= USSize;
259			else {
260			oursizhints.width = xmax;
261			oursizhints.height = ymax;
262			oursizhints.flags \|= PSize;
263			}
264			if ((i&(XValue\|YValue)) == (XValue\|YValue)) {
265	greg	1.1	oursizhints.flags \|= USPosition;
266			if (i & XNegative)
267	greg	1.3	oursizhints.x += DisplayWidth(thedisplay,
268			ourscreen)-1-oursizhints.width-2*BORWIDTH;
269	greg	1.1	if (i & YNegative)
270	greg	1.3	oursizhints.y += DisplayHeight(thedisplay,
271			ourscreen)-1-oursizhints.height-2*BORWIDTH;
272	greg	1.1	}
273			XSetNormalHints(thedisplay, wind, &oursizhints);
274			}
275	greg	1.15	ourxwmhints.flags = InputHint\|IconPixmapHint;
276			ourxwmhints.input = True;
277			ourxwmhints.icon_pixmap = XCreateBitmapFromData(thedisplay,
278	greg	1.17	wind, icondata, iconwidth, iconheight);
279	greg	1.15	XSetWMHints(thedisplay, wind, &ourxwmhints);
280	greg	1.2	XSelectInput(thedisplay, wind, ButtonPressMask\|ButtonReleaseMask
281			\|ButtonMotionMask\|StructureNotifyMask
282			\|KeyPressMask\|ExposureMask);
283	greg	1.1	XMapWindow(thedisplay, wind);
284			return;
285			memerr:
286			quiterr("out of memory");
287			} /* end of init */
288
289
290			quiterr(err) /* print message and exit */
291			char *err;
292			{
293			if (err != NULL) {
294			fprintf(stderr, "%s: %s\n", progname, err);
295			exit(1);
296			}
297			exit(0);
298			}
299
300
301			eputs(s)
302			char *s;
303			{
304			fputs(s, stderr);
305			}
306
307
308			quit(code)
309			int code;
310			{
311			exit(code);
312			}
313
314
315			getras() /* get raster file */
316			{
317			colormap ourmap;
318			XVisualInfo vinfo;
319
320			if (maxcolors <= 2) { /* monochrome */
321			ourdata = (unsigned char )malloc(ymax((xmax+7)/8));
322			if (ourdata == NULL)
323			goto fail;
324			ourras = make_raster(thedisplay, ourscreen, 1, ourdata,
325			xmax, ymax, 8);
326			if (ourras == NULL)
327			goto fail;
328			getmono();
329	greg	1.7	} else if (XMatchVisualInfo(thedisplay,ourscreen,24,TrueColor,&vinfo)
330			/* kludge for DirectColor */
331			\|\| XMatchVisualInfo(thedisplay,ourscreen,24,DirectColor,&vinfo)) {
332	greg	1.1	ourdata = (unsigned char )malloc(xmaxymax*3);
333			if (ourdata == NULL)
334			goto fail;
335			ourras = make_raster(thedisplay, ourscreen, 24, ourdata,
336			xmax, ymax, 8);
337			if (ourras == NULL)
338			goto fail;
339			getfull();
340			} else {
341			ourdata = (unsigned char )malloc(xmaxymax);
342			if (ourdata == NULL)
343			goto fail;
344			ourras = make_raster(thedisplay, ourscreen, 8, ourdata,
345			xmax, ymax, 8);
346			if (ourras == NULL)
347			goto fail;
348			if (greyscale)
349			biq(dither,maxcolors,1,ourmap);
350			else
351			ciq(dither,maxcolors,1,ourmap);
352	greg	1.3	if (init_rcolors(ourras, ourmap[0], ourmap[1], ourmap[2]) == 0)
353	greg	1.1	goto fail;
354			}
355			return;
356			fail:
357	greg	1.9	quiterr("could not create raster image");
358	greg	1.1	}
359
360
361			getevent() /* process the next event */
362			{
363			union {
364			XEvent u;
365			XConfigureEvent c;
366			XExposeEvent e;
367			XButtonPressedEvent b;
368			XKeyPressedEvent k;
369			} e;
370
371			XNextEvent(thedisplay, &e.u);
372			switch (e.u.type) {
373			case KeyPress:
374			docom(&e.k);
375			break;
376			case ConfigureNotify:
377			width = e.c.width;
378			height = e.c.height;
379			break;
380			case MapNotify:
381			map_rcolors(ourras, wind);
382			if (fast)
383			make_rpixmap(ourras);
384			break;
385			case UnmapNotify:
386			unmap_rcolors(ourras);
387			break;
388			case Expose:
389			redraw(e.e.x, e.e.y, e.e.width, e.e.height);
390			break;
391			case ButtonPress:
392			if (e.b.state & (ShiftMask\|ControlMask))
393			moveimage(&e.b);
394			else
395			getbox(&e.b);
396			break;
397			}
398			}
399
400
401			docom(ekey) /* execute command */
402			XKeyPressedEvent *ekey;
403			{
404			char buf[80];
405			COLOR cval;
406			XColor cvx;
407			int com, n;
408			double comp;
409			FVECT rorg, rdir;
410
411			n = XLookupString(ekey, buf, sizeof(buf), NULL, NULL);
412			if (n == 0)
413			return(0);
414			com = buf[0];
415			switch (com) { /* interpret command */
416			case 'q':
417			case CTRL(D): /* quit */
418			quit(0);
419			case '\n':
420			case '\r':
421			case 'l':
422			case 'c': /* value */
423			if (avgbox(cval) == -1)
424			return(-1);
425			switch (com) {
426			case '\n':
427			case '\r': /* radiance */
428			sprintf(buf, "%.3f", intens(cval)/exposure);
429			break;
430			case 'l': /* luminance */
431	greg	1.16	sprintf(buf, "%.0fL", luminance(cval)/exposure);
432	greg	1.1	break;
433			case 'c': /* color */
434			comp = pow(2.0, (double)scale);
435			sprintf(buf, "(%.2f,%.2f,%.2f)",
436			colval(cval,RED)*comp,
437			colval(cval,GRN)*comp,
438			colval(cval,BLU)*comp);
439			break;
440			}
441	greg	1.3	XDrawImageString(thedisplay, wind, ourgc,
442			box.xmin, box.ymin+box.ysiz, buf, strlen(buf));
443	greg	1.1	return(0);
444			case 'i': /* identify (contour) */
445			if (ourras->pixels == NULL)
446			return(-1);
447			n = ourdata[ekey->x-xoff+xmax*(ekey->y-yoff)];
448			n = ourras->pmap[n];
449			cvx.pixel = ourras->cdefs[n].pixel;
450			cvx.red = random() & 65535;
451			cvx.green = random() & 65535;
452			cvx.blue = random() & 65535;
453	greg	1.2	cvx.flags = DoRed\|DoGreen\|DoBlue;
454			XStoreColor(thedisplay, ourras->cmap, &cvx);
455	greg	1.1	return(0);
456			case 'p': /* position */
457			sprintf(buf, "(%d,%d)", ekey->x-xoff, ymax-1-ekey->y+yoff);
458			XDrawImageString(thedisplay, wind, ourgc, ekey->x, ekey->y,
459			buf, strlen(buf));
460			return(0);
461			case 't': /* trace */
462			if (!gotview) {
463			XBell(thedisplay, 0);
464			return(-1);
465			}
466	greg	1.8	if (viewray(rorg, rdir, &ourview,
467	greg	1.1	(ekey->x-xoff+.5)/xmax,
468	greg	1.8	(ymax-1-ekey->y+yoff+.5)/ymax) < 0)
469			return(-1);
470	greg	1.1	printf("%e %e %e ", rorg[0], rorg[1], rorg[2]);
471			printf("%e %e %e\n", rdir[0], rdir[1], rdir[2]);
472			fflush(stdout);
473			return(0);
474			case '=': /* adjust exposure */
475			if (avgbox(cval) == -1)
476			return(-1);
477			n = log(.5/bright(cval))/.69315 - scale; /* truncate */
478			if (n == 0)
479			return(0);
480			scale_rcolors(ourras, pow(2.0, (double)n));
481			scale += n;
482			sprintf(buf, "%+d", scale);
483	greg	1.3	XDrawImageString(thedisplay, wind, ourgc,
484			box.xmin, box.ymin+box.ysiz, buf, strlen(buf));
485	greg	1.1	XFlush(thedisplay);
486			free(ourdata);
487			free_raster(ourras);
488			getras();
489			/* fall through */
490			case CTRL(R): /* redraw */
491			case CTRL(L):
492			unmap_rcolors(ourras);
493			XClearWindow(thedisplay, wind);
494			map_rcolors(ourras, wind);
495			if (fast)
496			make_rpixmap(ourras);
497			redraw(0, 0, width, height);
498			return(0);
499			case ' ': /* clear */
500			redraw(box.xmin, box.ymin, box.xsiz, box.ysiz);
501			return(0);
502			default:
503	greg	1.2	XBell(thedisplay, 0);
504	greg	1.1	return(-1);
505			}
506			}
507
508
509	greg	1.2	moveimage(ebut) /* shift the image */
510			XButtonPressedEvent *ebut;
511	greg	1.1	{
512	greg	1.2	union {
513	greg	1.3	XEvent u;
514			XButtonReleasedEvent b;
515			XPointerMovedEvent m;
516			} e;
517	greg	1.5	int mxo, myo;
518	greg	1.1
519	greg	1.3	XMaskEvent(thedisplay, ButtonReleaseMask\|ButtonMotionMask, &e.u);
520			while (e.u.type == MotionNotify) {
521	greg	1.5	mxo = e.m.x;
522			myo = e.m.y;
523			revline(ebut->x, ebut->y, mxo, myo);
524			revbox(xoff+mxo-ebut->x, yoff+myo-ebut->y,
525			xoff+mxo-ebut->x+xmax, yoff+myo-ebut->y+ymax);
526	greg	1.3	XMaskEvent(thedisplay,ButtonReleaseMask\|ButtonMotionMask,&e.u);
527	greg	1.5	revline(ebut->x, ebut->y, mxo, myo);
528			revbox(xoff+mxo-ebut->x, yoff+myo-ebut->y,
529			xoff+mxo-ebut->x+xmax, yoff+myo-ebut->y+ymax);
530	greg	1.3	}
531	greg	1.2	xoff += e.b.x - ebut->x;
532			yoff += e.b.y - ebut->y;
533	greg	1.1	XClearWindow(thedisplay, wind);
534			redraw(0, 0, width, height);
535			}
536
537
538			getbox(ebut) /* get new box */
539			XButtonPressedEvent *ebut;
540			{
541			union {
542	greg	1.2	XEvent u;
543	greg	1.1	XButtonReleasedEvent b;
544			XPointerMovedEvent m;
545			} e;
546
547	greg	1.2	XMaskEvent(thedisplay, ButtonReleaseMask\|ButtonMotionMask, &e.u);
548			while (e.u.type == MotionNotify) {
549	greg	1.1	revbox(ebut->x, ebut->y, box.xmin = e.m.x, box.ymin = e.m.y);
550	greg	1.2	XMaskEvent(thedisplay,ButtonReleaseMask\|ButtonMotionMask,&e.u);
551	greg	1.1	revbox(ebut->x, ebut->y, box.xmin, box.ymin);
552			}
553			box.xmin = e.b.x<0 ? 0 : (e.b.x>=width ? width-1 : e.b.x);
554			box.ymin = e.b.y<0 ? 0 : (e.b.y>=height ? height-1 : e.b.y);
555			if (box.xmin > ebut->x) {
556			box.xsiz = box.xmin - ebut->x + 1;
557			box.xmin = ebut->x;
558			} else {
559			box.xsiz = ebut->x - box.xmin + 1;
560			}
561			if (box.ymin > ebut->y) {
562			box.ysiz = box.ymin - ebut->y + 1;
563			box.ymin = ebut->y;
564			} else {
565			box.ysiz = ebut->y - box.ymin + 1;
566			}
567			}
568
569
570			revbox(x0, y0, x1, y1) /* draw box with reversed lines */
571			int x0, y0, x1, y1;
572			{
573	greg	1.5	revline(x0, y0, x1, y0);
574			revline(x0, y1, x1, y1);
575			revline(x0, y0, x0, y1);
576			revline(x1, y0, x1, y1);
577			}
578	greg	1.1
579
580			avgbox(clr) /* average color over current box */
581			COLOR clr;
582			{
583	greg	1.25	static COLOR lc;
584			static int ll, lr, lt, lb;
585	greg	1.1	int left, right, top, bottom;
586			int y;
587			double d;
588			COLOR ctmp;
589			register int x;
590
591			setcolor(clr, 0.0, 0.0, 0.0);
592			left = box.xmin - xoff;
593			right = left + box.xsiz;
594			if (left < 0)
595			left = 0;
596			if (right > xmax)
597			right = xmax;
598			if (left >= right)
599			return(-1);
600			top = box.ymin - yoff;
601			bottom = top + box.ysiz;
602			if (top < 0)
603			top = 0;
604			if (bottom > ymax)
605			bottom = ymax;
606			if (top >= bottom)
607			return(-1);
608	greg	1.25	if (left == ll && right == lr && top == lt && bottom == lb) {
609			copycolor(clr, lc);
610			return;
611			}
612	greg	1.1	for (y = top; y < bottom; y++) {
613			if (getscan(y) == -1)
614			return(-1);
615			for (x = left; x < right; x++) {
616			colr_color(ctmp, scanline[x]);
617			addcolor(clr, ctmp);
618			}
619			}
620			d = 1.0/((right-left)*(bottom-top));
621			scalecolor(clr, d);
622	greg	1.25	ll = left; lr = right; lt = top; lb = bottom;
623			copycolor(lc, clr);
624	greg	1.1	return(0);
625			}
626
627
628			getmono() /* get monochrome data */
629			{
630			register unsigned char *dp;
631			register int x, err;
632	greg	1.23	int y, errp;
633	greg	1.1	short *cerr;
634
635	greg	1.10	if ((cerr = (short *)calloc(xmax,sizeof(short))) == NULL)
636	greg	1.9	quiterr("out of memory in getmono");
637	greg	1.1	dp = ourdata - 1;
638			for (y = 0; y < ymax; y++) {
639	greg	1.10	if (getscan(y) < 0)
640			quiterr("seek error in getmono");
641			normcolrs(scanline, xmax, scale);
642	greg	1.17	add2icon(y, scanline);
643	greg	1.1	err = 0;
644			for (x = 0; x < xmax; x++) {
645			if (!(x&7))
646			*++dp = 0;
647	greg	1.23	errp = err;
648	greg	1.10	err += normbright(scanline[x]) + cerr[x];
649	greg	1.1	if (err > 127)
650			err -= 255;
651			else
652	greg	1.3	*dp \|= 1<<(7-(x&07));
653	greg	1.23	err /= 3;
654			cerr[x] = err + errp;
655	greg	1.1	}
656			}
657			free((char *)cerr);
658			}
659
660
661	greg	1.17	add2icon(y, scan) /* add a scanline to our icon data */
662			int y;
663			COLR *scan;
664			{
665			static short cerr[ICONSIZ];
666	greg	1.20	static int ynext;
667			static char *dp;
668	greg	1.17	register int err;
669	greg	1.20	register int x, ti;
670	greg	1.23	int errp;
671	greg	1.17
672			if (iconheight == 0) { /* initialize */
673	greg	1.18	if (xmax <= ICONSIZ && ymax <= ICONSIZ) {
674	greg	1.17	iconwidth = xmax;
675			iconheight = ymax;
676			} else if (xmax > ymax) {
677			iconwidth = ICONSIZ;
678			iconheight = ICONSIZ*ymax/xmax;
679	greg	1.24	if (iconheight < 1)
680			iconheight = 1;
681	greg	1.17	} else {
682			iconwidth = ICONSIZ*xmax/ymax;
683	greg	1.24	if (iconwidth < 1)
684			iconwidth = 1;
685	greg	1.17	iconheight = ICONSIZ;
686			}
687	greg	1.20	ynext = 0;
688	greg	1.17	dp = icondata - 1;
689			}
690	greg	1.20	if (y < ynextymax/iconheight) / skip this one */
691	greg	1.17	return;
692			err = 0;
693			for (x = 0; x < iconwidth; x++) {
694			if (!(x&7))
695			*++dp = 0;
696	greg	1.23	errp = err;
697	greg	1.20	ti = x*xmax/iconwidth;
698			err += normbright(scan[ti]) + cerr[x];
699	greg	1.17	if (err > 127)
700			err -= 255;
701			else
702			*dp \|= 1<<(x&07);
703	greg	1.23	err /= 3;
704			cerr[x] = err + errp;
705	greg	1.17	}
706	greg	1.20	ynext++;
707	greg	1.17	}
708
709
710	greg	1.1	getfull() /* get full (24-bit) data */
711			{
712			int y;
713	greg	1.10	register unsigned char *dp;
714			register int x;
715			/* set gamma correction */
716			setcolrgam(gamcor);
717			/* read and convert file */
718			dp = ourdata;
719			for (y = 0; y < ymax; y++) {
720			if (getscan(y) < 0)
721			quiterr("seek error in getfull");
722			if (scale)
723			shiftcolrs(scanline, xmax, scale);
724			colrs_gambs(scanline, xmax);
725	greg	1.17	add2icon(y, scanline);
726	greg	1.10	for (x = 0; x < xmax; x++) {
727			*dp++ = scanline[x][RED];
728			*dp++ = scanline[x][GRN];
729			*dp++ = scanline[x][BLU];
730			}
731			}
732	greg	1.1	}
733
734
735			scale_rcolors(xr, sf) /* scale color map */
736			register XRASTER *xr;
737			double sf;
738			{
739			register int i;
740			long maxv;
741
742			if (xr->pixels == NULL)
743			return;
744
745			sf = pow(sf, 1.0/gamcor);
746			maxv = 65535/sf;
747
748			for (i = xr->ncolors; i--; ) {
749			xr->cdefs[i].red = xr->cdefs[i].red > maxv ?
750			65535 :
751			xr->cdefs[i].red * sf;
752			xr->cdefs[i].green = xr->cdefs[i].green > maxv ?
753			65535 :
754			xr->cdefs[i].green * sf;
755			xr->cdefs[i].blue = xr->cdefs[i].blue > maxv ?
756			65535 :
757			xr->cdefs[i].blue * sf;
758			}
759			XStoreColors(thedisplay, xr->cmap, xr->cdefs, xr->ncolors);
760			}
761
762
763			getscan(y)
764			int y;
765			{
766			if (y != cury) {
767			if (scanpos == NULL \|\| scanpos[y] == -1)
768			return(-1);
769			if (fseek(fin, scanpos[y], 0) == -1)
770	greg	1.9	quiterr("fseek error");
771	greg	1.1	cury = y;
772			} else if (scanpos != NULL)
773			scanpos[y] = ftell(fin);
774
775			if (freadcolrs(scanline, xmax, fin) < 0)
776			quiterr("read error");
777
778			cury++;
779			return(0);
780			}
781
782
783			picreadline3(y, l3) /* read in 3-byte scanline */
784			int y;
785			register rgbpixel *l3;
786			{
787			register int i;
788			/* read scanline */
789			if (getscan(y) < 0)
790			quiterr("cannot seek for picreadline");
791			/* convert scanline */
792			normcolrs(scanline, xmax, scale);
793	greg	1.17	add2icon(y, scanline);
794	greg	1.1	for (i = 0; i < xmax; i++) {
795			l3[i].r = scanline[i][RED];
796			l3[i].g = scanline[i][GRN];
797			l3[i].b = scanline[i][BLU];
798			}
799			}
800
801
802			picwriteline(y, l) /* add 8-bit scanline to image */
803			int y;
804			pixel *l;
805			{
806			bcopy((char )l, (char )ourdata+y*xmax, xmax);
807			}
808
809
810			picreadcm(map) /* do gamma correction */
811			colormap map;
812			{
813			extern double pow();
814			register int i, val;
815
816			for (i = 0; i < 256; i++) {
817	greg	1.11	val = pow((i+0.5)/256.0, 1.0/gamcor) * 256.0;
818	greg	1.1	map[0][i] = map[1][i] = map[2][i] = val;
819			}
820			}
821
822
823			picwritecm(map) /* handled elsewhere */
824			colormap map;
825			{
826			#ifdef DEBUG
827			register int i;
828
829			for (i = 0; i < 256; i++)
830			printf("%d %d %d\n", map[0][i],map[1][i],map[2][i]);
831			#endif
832			}