tiff2ps.c

Go to the documentation of this file.

/* -*- Mode: C; tab-width: 8; indent-tabs-mode: t; c-basic-offset: 8; c-indent-level: 8 -*- */
/* $Id$ */

/*
 * Copyright (c) 1988-1997 Sam Leffler
 * Copyright (c) 1991-1997 Silicon Graphics, Inc.
 *
 * Permission to use, copy, modify, distribute, and sell this software and 
 * its documentation for any purpose is hereby granted without fee, provided
 * that (i) the above copyright notices and this permission notice appear in
 * all copies of the software and related documentation, and (ii) the names of
 * Sam Leffler and Silicon Graphics may not be used in any advertising or
 * publicity relating to the software without the specific, prior written
 * permission of Sam Leffler and Silicon Graphics.
 * 
 * THE SOFTWARE IS PROVIDED "AS-IS" AND WITHOUT WARRANTY OF ANY KIND, 
 * EXPRESS, IMPLIED OR OTHERWISE, INCLUDING WITHOUT LIMITATION, ANY 
 * WARRANTY OF MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE.  
 * 
 * IN NO EVENT SHALL SAM LEFFLER OR SILICON GRAPHICS BE LIABLE FOR
 * ANY SPECIAL, INCIDENTAL, INDIRECT OR CONSEQUENTIAL DAMAGES OF ANY KIND,
 * OR ANY DAMAGES WHATSOEVER RESULTING FROM LOSS OF USE, DATA OR PROFITS,
 * WHETHER OR NOT ADVISED OF THE POSSIBILITY OF DAMAGE, AND ON ANY THEORY OF 
 * LIABILITY, ARISING OUT OF OR IN CONNECTION WITH THE USE OR PERFORMANCE 
 * OF THIS SOFTWARE.
 */

/*
 * Modified for use as Evince TIFF ps exporter by
 * Matthew S. Wilson <msw@rpath.com>
 * Modifications Copyright (C) 2005 rpath, Inc.
 *
 */

#include <config.h>
#include <stdio.h>
#include <stdlib.h>         /* for atof */
#include <math.h>
#include <time.h>
#include <string.h>
#include <unistd.h>

#include <glib.h>
#include <glib/gstdio.h>

#include "tiff2ps.h"

/*
 * Revision history
 *
 * 2001-Mar-21
 *    I (Bruce A. Mallett) added this revision history comment ;)
 *
 *    Fixed PS_Lvl2page() code which outputs non-ASCII85 raw
 *    data.  Moved test for when to output a line break to
 *    *after* the output of a character.  This just serves
 *    to fix an eye-nuisance where the first line of raw
 *    data was one character shorter than subsequent lines.
 *
 *    Added an experimental ASCII85 encoder which can be used
 *    only when there is a single buffer of bytes to be encoded.
 *    This version is much faster at encoding a straight-line
 *    buffer of data because it can avoid alot of the loop
 *    overhead of the byte-by-bye version.  To use this version
 *    you need to define EXP_ASCII85ENCODER (experimental ...).
 *
 *    Added bug fix given by Michael Schmidt to PS_Lvl2page()
 *    in which an end-of-data marker ('>') was not being output
 *    when producing non-ASCII85 encoded PostScript Level 2
 *    data.
 *
 *    Fixed PS_Lvl2colorspace() so that it no longer assumes that
 *    a TIFF having more than 2 planes is a CMYK.  This routine
 *    no longer looks at the samples per pixel but instead looks
 *    at the "photometric" value.  This change allows support of
 *    CMYK TIFFs.
 *
 *    Modified the PostScript L2 imaging loop so as to test if
 *    the input stream is still open before attempting to do a
 *    flushfile on it.  This was done because some RIPs close
 *    the stream after doing the image operation.
 *
 *    Got rid of the realloc() being done inside a loop in the
 *    PSRawDataBW() routine.  The code now walks through the
 *    byte-size array outside the loop to determine the largest
 *    size memory block that will be needed.
 *
 *    Added "-m" switch to ask tiff2ps to, where possible, use the
 *    "imagemask" operator instead of the "image" operator.
 *
 *    Added the "-i #" switch to allow interpolation to be disabled.
 *
 *    Unrolled a loop or two to improve performance.
 */

/*
 * Define EXP_ASCII85ENCODER if you want to use an experimental
 * version of the ASCII85 encoding routine.  The advantage of
 * using this routine is that tiff2ps will convert to ASCII85
 * encoding at between 3 and 4 times the speed as compared to
 * using the old (non-experimental) encoder.  The disadvantage
 * is that you will be using a new (and unproven) encoding
 * routine.  So user beware, you have been warned!
 */

#define EXP_ASCII85ENCODER

/*
 * NB: this code assumes uint32 works with printf's %l[ud].
 */

struct _TIFF2PSContext
{
    char *filename;     /* input filename */
    FILE *fd;       /* output file stream */
    int ascii85;        /* use ASCII85 encoding */
    int interpolate;    /* interpolate level2 image */
    int level2;     /* generate PostScript level 2 */
    int level3;     /* generate PostScript level 3 */
    int generateEPSF;   /* generate Encapsulated PostScript */
    int PSduplex;       /* enable duplex printing */
    int PStumble;       /* enable top edge binding */
    int PSavoiddeadzone;    /* enable avoiding printer deadzone */
    double maxPageHeight;   /* maximum size to fit on page */
    double splitOverlap;    /* amount for split pages to overlag */
    int rotate;     /* rotate image by 180 degrees */
    int useImagemask;   /* Use imagemask instead of image operator */
    uint16 res_unit;    /* Resolution units: 2 - inches, 3 - cm */
    int npages;     /* number of pages processed */

    tsize_t tf_bytesperrow;
    tsize_t ps_bytesperrow;
    tsize_t tf_rowsperstrip;
    tsize_t tf_numberstrips;

    /*
     * ASCII85 Encoding Support.
     */
    unsigned char ascii85buf[10];
    int ascii85count;
    int ascii85breaklen;
    uint16 samplesperpixel;
    uint16 bitspersample;
    uint16 planarconfiguration;
    uint16 photometric;
    uint16 compression;
    uint16 extrasamples;
    int alpha;
};

static void PSpage(TIFF2PSContext*, TIFF*, uint32, uint32);
static void PSColorContigPreamble(TIFF2PSContext*, uint32, uint32, int);
static void PSColorSeparatePreamble(TIFF2PSContext*, uint32, uint32, int);
static void PSDataColorContig(TIFF2PSContext*, TIFF*, uint32, uint32, int);
static void PSDataColorSeparate(TIFF2PSContext*, TIFF*, uint32, uint32, int);
static void PSDataPalette(TIFF2PSContext*, TIFF*, uint32, uint32);
static void PSDataBW(TIFF2PSContext*, TIFF*, uint32, uint32);
static void Ascii85Init(TIFF2PSContext*);
static void Ascii85Put(TIFF2PSContext*, unsigned char);
static void Ascii85Flush(TIFF2PSContext*);
static void PSHead(TIFF2PSContext*, TIFF*, uint32, uint32,
           double, double, double, double);
static void PSTail(TIFF2PSContext*);

#if defined( EXP_ASCII85ENCODER )
static int Ascii85EncodeBlock(TIFF2PSContext*, uint8 * ascii85_p,
                  unsigned f_eod, const uint8 * raw_p, int raw_l);
#endif

#define IMAGEOP(ctx) ((ctx)->useImagemask && ((ctx)->bitspersample == 1)) ? "imagemask" : "image"

TIFF2PSContext* tiff2ps_context_new(const gchar *filename) {
    TIFF2PSContext* ctx;

    ctx = g_new0(TIFF2PSContext, 1);
    ctx->filename = g_strdup(filename);
    ctx->fd = g_fopen(ctx->filename, "w");
    if (ctx->fd == NULL) {
        g_free (ctx->filename);
        g_free (ctx);
        return NULL;
    }
    ctx->interpolate = TRUE;     /* interpolate level2 image */
    ctx->PSavoiddeadzone = TRUE; /* enable avoiding printer deadzone */
    return ctx;
}

void tiff2ps_context_finalize(TIFF2PSContext *ctx) {
    PSTail(ctx);
    fclose(ctx->fd);
    g_free(ctx->filename);
    g_free(ctx);
}

static int
checkImage(TIFF2PSContext *ctx, TIFF* tif)
{
    switch (ctx->photometric) {
    case PHOTOMETRIC_YCBCR:
        if ((ctx->compression == COMPRESSION_JPEG
             || ctx->compression == COMPRESSION_OJPEG)
            && ctx->planarconfiguration == PLANARCONFIG_CONTIG) {
            /* can rely on libjpeg to convert to RGB */
            TIFFSetField(tif, TIFFTAG_JPEGCOLORMODE,
                     JPEGCOLORMODE_RGB);
            ctx->photometric = PHOTOMETRIC_RGB;
        } else {
            if (ctx->level2 || ctx->level3)
                break;
            TIFFError(ctx->filename, "Can not handle image with %s",
                "Ctx->PhotometricInterpretation=YCbCr");
            return (0);
        }
        /* fall thru... */
    case PHOTOMETRIC_RGB:
        if (ctx->alpha && ctx->bitspersample != 8) {
            TIFFError(ctx->filename,
                "Can not handle %d-bit/sample RGB image with ctx->alpha",
                ctx->bitspersample);
            return (0);
        }
        /* fall thru... */
    case PHOTOMETRIC_SEPARATED:
    case PHOTOMETRIC_PALETTE:
    case PHOTOMETRIC_MINISBLACK:
    case PHOTOMETRIC_MINISWHITE:
        break;
    case PHOTOMETRIC_LOGL:
    case PHOTOMETRIC_LOGLUV:
        if (ctx->compression != COMPRESSION_SGILOG &&
            ctx->compression != COMPRESSION_SGILOG24) {
            TIFFError(ctx->filename,
            "Can not handle %s data with ctx->compression other than SGILog",
                (ctx->photometric == PHOTOMETRIC_LOGL) ?
                "LogL" : "LogLuv"
            );
            return (0);
        }
        /* rely on library to convert to RGB/greyscale */
        TIFFSetField(tif, TIFFTAG_SGILOGDATAFMT, SGILOGDATAFMT_8BIT);
        ctx->photometric = (ctx->photometric == PHOTOMETRIC_LOGL) ?
            PHOTOMETRIC_MINISBLACK : PHOTOMETRIC_RGB;
        ctx->bitspersample = 8;
        break;
    case PHOTOMETRIC_CIELAB:
        /* fall thru... */
    default:
        TIFFError(ctx->filename,
            "Can not handle image with Ctx->PhotometricInterpretation=%d",
            ctx->photometric);
        return (0);
    }
    switch (ctx->bitspersample) {
    case 1: case 2:
    case 4: case 8:
        break;
    default:
        TIFFError(ctx->filename, "Can not handle %d-bit/sample image",
            ctx->bitspersample);
        return (0);
    }
    if (ctx->planarconfiguration == PLANARCONFIG_SEPARATE &&
        ctx->extrasamples > 0)
        TIFFWarning(ctx->filename, "Ignoring extra samples");
    return (1);
}

#define PS_UNIT_SIZE    72.0F
#define PSUNITS(npix,res)   ((npix) * (PS_UNIT_SIZE / (res)))

static  char RGBcolorimage[] = "\
/bwproc {\n\
    rgbproc\n\
    dup length 3 idiv string 0 3 0\n\
    5 -1 roll {\n\
    add 2 1 roll 1 sub dup 0 eq {\n\
        pop 3 idiv\n\
        3 -1 roll\n\
        dup 4 -1 roll\n\
        dup 3 1 roll\n\
        5 -1 roll put\n\
        1 add 3 0\n\
    } { 2 1 roll } ifelse\n\
    } forall\n\
    pop pop pop\n\
} def\n\
/colorimage where {pop} {\n\
    /colorimage {pop pop /rgbproc exch def {bwproc} image} bind def\n\
} ifelse\n\
";

/*
 * Adobe Photoshop requires a comment line of the form:
 *
 * %ImageData: <cols> <rows> <depth>  <main channels> <pad channels>
 *  <block size> <1 for binary|2 for hex> "data start"
 *
 * It is claimed to be part of some future revision of the EPS spec.
 */
G_GNUC_PRINTF (6, 7) static void
PhotoshopBanner(TIFF2PSContext* ctx, uint32 w, uint32 h, int bs, int nc,
        const char* startline, ...)
{
    va_list args;
    fprintf(ctx->fd, "%%ImageData: %ld %ld %d %d 0 %d 2 \"",
        (long) w, (long) h, ctx->bitspersample, nc, bs);

    va_start(args, startline);
    vfprintf(ctx->fd, startline, args);
    va_end(args);

    fprintf(ctx->fd, "\"\n");
}

/*
 *   pw : image width in pixels
 *   ph : image height in pixels
 * pprw : image width in PS units (72 dpi)
 * pprh : image height in PS units (72 dpi)
 */
static void
setupPageState(TIFF2PSContext *ctx, TIFF* tif, uint32* pw, uint32* ph,
           double* pprw, double* pprh)
{
    float xres = 0.0F, yres = 0.0F;

    TIFFGetField(tif, TIFFTAG_IMAGEWIDTH, pw);
    TIFFGetField(tif, TIFFTAG_IMAGELENGTH, ph);
    if (ctx->res_unit == 0)
        TIFFGetFieldDefaulted(tif, TIFFTAG_RESOLUTIONUNIT, &ctx->res_unit);
    /*
     * Calculate printable area.
     */
    if (!TIFFGetField(tif, TIFFTAG_XRESOLUTION, &xres)
            || fabs(xres) < 0.0000001)
        xres = PS_UNIT_SIZE;
    if (!TIFFGetField(tif, TIFFTAG_YRESOLUTION, &yres)
            || fabs(yres) < 0.0000001)
        yres = PS_UNIT_SIZE;
    switch (ctx->res_unit) {
    case RESUNIT_CENTIMETER:
        xres *= 2.54F, yres *= 2.54F;
        break;
    case RESUNIT_INCH:
        break;
    case RESUNIT_NONE:
    default:
        xres *= PS_UNIT_SIZE, yres *= PS_UNIT_SIZE;
        break;
    }
    *pprh = PSUNITS(*ph, yres);
    *pprw = PSUNITS(*pw, xres);
}

static int
isCCITTCompression(TIFF* tif)
{
    uint16 compress;
    TIFFGetField(tif, TIFFTAG_COMPRESSION, &compress);
    return (compress == COMPRESSION_CCITTFAX3 ||
        compress == COMPRESSION_CCITTFAX4 ||
        compress == COMPRESSION_CCITTRLE ||
        compress == COMPRESSION_CCITTRLEW);
}

static  char *hex = "0123456789abcdef";

/*
 * imagewidth & imageheight are 1/72 inches
 * pagewidth & pageheight are inches
 */
static int
PlaceImage(TIFF2PSContext *ctx, double pagewidth, double pageheight,
       double imagewidth, double imageheight, int splitpage,
       double lm, double bm, int cnt)
{
    double xtran = 0;
    double ytran = 0;
    double xscale = 1;
    double yscale = 1;
    double left_offset = lm * PS_UNIT_SIZE;
    double bottom_offset = bm * PS_UNIT_SIZE;
    double subimageheight;
    double splitheight;
    double overlap;
    /* buffers for locale-insitive number formatting */
    gchar buf[2][G_ASCII_DTOSTR_BUF_SIZE];

    pagewidth *= PS_UNIT_SIZE;
    pageheight *= PS_UNIT_SIZE;

    if (ctx->maxPageHeight==0)
        splitheight = 0;
    else
        splitheight = ctx->maxPageHeight * PS_UNIT_SIZE;
    overlap = ctx->splitOverlap * PS_UNIT_SIZE;

    /*
     * WIDTH:
     *      if too wide, scrunch to fit
     *      else leave it alone
     */
    if (imagewidth <= pagewidth) {
        xscale = imagewidth;
    } else {
        xscale = pagewidth;
    }

    /* HEIGHT:
     *      if too long, scrunch to fit
     *      if too short, move to top of page
     */
    if (imageheight <= pageheight) {
        yscale = imageheight;
        ytran = pageheight - imageheight;
    } else if (imageheight > pageheight &&
        (splitheight == 0 || imageheight <= splitheight)) {
        yscale = pageheight;
    } else /* imageheight > splitheight */ {
        subimageheight = imageheight - (pageheight-overlap)*splitpage;
        if (subimageheight <= pageheight) {
            yscale = imageheight;
            ytran = pageheight - subimageheight;
            splitpage = 0;
        } else if ( subimageheight > pageheight && subimageheight <= splitheight) {
            yscale = imageheight * pageheight / subimageheight;
            ytran = 0;
            splitpage = 0;
        } else /* sumimageheight > splitheight */ {
            yscale = imageheight;
            ytran = pageheight - subimageheight;
            splitpage++;
        }
    }

    bottom_offset += ytran / (cnt?2:1);
    if (cnt)
        left_offset += xtran / 2;

    fprintf(ctx->fd, "%s %s translate\n",
        g_ascii_dtostr(buf[0], sizeof(buf[0]), left_offset),
        g_ascii_dtostr(buf[1], sizeof(buf[1]), bottom_offset));
    fprintf(ctx->fd, "%s %s scale\n",
        g_ascii_dtostr(buf[0], sizeof(buf[0]), xscale),
        g_ascii_dtostr(buf[1], sizeof(buf[1]), yscale));
    if (ctx->rotate)
        fputs ("1 1 translate 180 ctx->rotate\n", ctx->fd);

    return splitpage;
}


void
tiff2ps_process_page(TIFF2PSContext* ctx, TIFF* tif, double pw, double ph,
             double lm, double bm, gboolean cnt)
{
    uint32 w, h;
    float ox, oy;
        double prw, prh;
    double scale = 1.0;
    double left_offset = lm * PS_UNIT_SIZE;
    double bottom_offset = bm * PS_UNIT_SIZE;
    uint16* sampleinfo;
    int split;
    /* buffers for locale-insitive number formatting */
    gchar buf[2][G_ASCII_DTOSTR_BUF_SIZE];

    if (!TIFFGetField(tif, TIFFTAG_XPOSITION, &ox))
        ox = 0;
    if (!TIFFGetField(tif, TIFFTAG_YPOSITION, &oy))
        oy = 0;
    setupPageState(ctx, tif, &w, &h, &prw, &prh);

    ctx->tf_numberstrips = TIFFNumberOfStrips(tif);
    TIFFGetFieldDefaulted(tif, TIFFTAG_ROWSPERSTRIP,
                  &ctx->tf_rowsperstrip);
    setupPageState(ctx, tif, &w, &h, &prw, &prh);
    if (!ctx->npages)
        PSHead(ctx, tif, w, h, prw, prh, ox, oy);
    TIFFGetFieldDefaulted(tif, TIFFTAG_BITSPERSAMPLE,
                  &ctx->bitspersample);
    TIFFGetFieldDefaulted(tif, TIFFTAG_SAMPLESPERPIXEL,
                  &ctx->samplesperpixel);
    TIFFGetFieldDefaulted(tif, TIFFTAG_PLANARCONFIG,
                  &ctx->planarconfiguration);
    TIFFGetField(tif, TIFFTAG_COMPRESSION, &ctx->compression);
    TIFFGetFieldDefaulted(tif, TIFFTAG_EXTRASAMPLES,
                  &ctx->extrasamples, &sampleinfo);
    ctx->alpha = (ctx->extrasamples == 1 &&
              sampleinfo[0] == EXTRASAMPLE_ASSOCALPHA);
    if (!TIFFGetField(tif, TIFFTAG_PHOTOMETRIC, &ctx->photometric)) {
        switch (ctx->samplesperpixel - ctx->extrasamples) {
        case 1:
            if (isCCITTCompression(tif))
                ctx->photometric = PHOTOMETRIC_MINISWHITE;
            else
                ctx->photometric = PHOTOMETRIC_MINISBLACK;
            break;
        case 3:
            ctx->photometric = PHOTOMETRIC_RGB;
            break;
        case 4:
            ctx->photometric = PHOTOMETRIC_SEPARATED;
            break;
        }
    }
    if (checkImage(ctx, tif)) {
        ctx->tf_bytesperrow = TIFFScanlineSize(tif);
        ctx->npages++;
        fprintf(ctx->fd, "%%%%Page: %d %d\n", ctx->npages,
            ctx->npages);
        if (!ctx->generateEPSF && ( ctx->level2 || ctx->level3 )) {
            double psw = 0.0, psh = 0.0;
            if (psw != 0.0) {
                psw = pw * PS_UNIT_SIZE;
                if (ctx->res_unit == RESUNIT_CENTIMETER)
                    psw *= 2.54F;
            } else
                psw=ctx->rotate ? prh:prw;
            if (psh != 0.0) {
                psh = ph * PS_UNIT_SIZE;
                if (ctx->res_unit == RESUNIT_CENTIMETER)
                    psh *= 2.54F;
            } else
                psh=ctx->rotate ? prw:prh;
            fprintf(ctx->fd,
                "1 dict begin /PageSize [ %s %s ] def currentdict end setpagedevice\n",
                g_ascii_dtostr(buf[0], sizeof(buf[0]), psw),
                g_ascii_dtostr(buf[1], sizeof(buf[1]), psh));
            fputs(
                  "<<\n  /Policies <<\n    /PageSize 3\n  >>\n>> setpagedevice\n",
                  ctx->fd);
        }
        fprintf(ctx->fd, "gsave\n");
        fprintf(ctx->fd, "100 dict begin\n");
        if (pw != 0 || ph != 0) {
            if (!pw)
                pw = prw;
            if (!ph)
                ph = prh;
            if (ctx->maxPageHeight) { /* used -H option */
                split = PlaceImage(ctx,pw,ph,prw,prh,
                           0,lm,bm,cnt);
                while( split ) {
                    PSpage(ctx, tif, w, h);
                    fprintf(ctx->fd, "end\n");
                    fprintf(ctx->fd, "grestore\n");
                    fprintf(ctx->fd, "showpage\n");
                    ctx->npages++;
                    fprintf(ctx->fd, "%%%%Page: %d %d\n",
                        ctx->npages, ctx->npages);
                    fprintf(ctx->fd, "gsave\n");
                    fprintf(ctx->fd, "100 dict begin\n");
                    split = PlaceImage(ctx,pw,ph,prw,prh,
                               split,lm,bm,cnt);
                }
            } else {
                pw *= PS_UNIT_SIZE;
                ph *= PS_UNIT_SIZE;

                /* NB: maintain image aspect ratio */
                scale = pw/prw < ph/prh ?
                    pw/prw : ph/prh;
                if (scale > 1.0)
                    scale = 1.0;
                if (cnt) {
                    bottom_offset +=
                        (ph - prh * scale) / 2;
                    left_offset +=
                        (pw - prw * scale) / 2;
                }
                fprintf(ctx->fd, "%s %s translate\n",
                    g_ascii_dtostr(buf[0], sizeof(buf[0]), left_offset),
                    g_ascii_dtostr(buf[1], sizeof(buf[1]), bottom_offset));
                fprintf(ctx->fd, "%s %s scale\n",
                    g_ascii_dtostr(buf[0], sizeof(buf[0]), prw * scale),
                    g_ascii_dtostr(buf[1], sizeof(buf[1]), prh * scale));
                if (ctx->rotate)
                    fputs ("1 1 translate 180 ctx->rotate\n", ctx->fd);
            }
        } else {
            fprintf(ctx->fd, "%s %s scale\n",
                g_ascii_dtostr(buf[0], sizeof(buf[0]), prw),
                g_ascii_dtostr(buf[1], sizeof(buf[1]), prh));
            if (ctx->rotate)
                fputs ("1 1 translate 180 ctx->rotate\n", ctx->fd);
        }
        PSpage(ctx, tif, w, h);
        fprintf(ctx->fd, "end\n");
        fprintf(ctx->fd, "grestore\n");
        fprintf(ctx->fd, "showpage\n");
    }
}


static char DuplexPreamble[] = "\
%%BeginFeature: *Duplex True\n\
systemdict begin\n\
  /languagelevel where { pop languagelevel } { 1 } ifelse\n\
  2 ge { 1 dict dup /Duplex true put setpagedevice }\n\
  { statusdict /setduplex known { statusdict begin setduplex true end } if\n\
  } ifelse\n\
end\n\
%%EndFeature\n\
";

static char TumblePreamble[] = "\
%%BeginFeature: *Tumble True\n\
systemdict begin\n\
  /languagelevel where { pop languagelevel } { 1 } ifelse\n\
  2 ge { 1 dict dup /Tumble true put setpagedevice }\n\
  { statusdict /settumble known { statusdict begin true settumble end } if\n\
  } ifelse\n\
end\n\
%%EndFeature\n\
";

static char AvoidDeadZonePreamble[] = "\
gsave newpath clippath pathbbox grestore\n\
  4 2 roll 2 copy translate\n\
  exch 3 1 roll sub 3 1 roll sub exch\n\
  currentpagedevice /PageSize get aload pop\n\
  exch 3 1 roll div 3 1 roll div abs exch abs\n\
  2 copy gt { exch } if pop\n\
  dup 1 lt { dup scale } { pop } ifelse\n\
";

void
PSHead(TIFF2PSContext *ctx, TIFF *tif, uint32 w, uint32 h,
       double pw, double ph, double ox, double oy)
{
    time_t t;

    (void) tif; (void) w; (void) h;
    t = time(0);
    fprintf(ctx->fd, "%%!PS-Adobe-3.0%s\n",
        ctx->generateEPSF ? " EPSF-3.0" : "");
    fprintf(ctx->fd, "%%%%Creator: Evince\n");
    fprintf(ctx->fd, "%%%%CreationDate: %s", ctime(&t));
    fprintf(ctx->fd, "%%%%DocumentData: Clean7Bit\n");
    fprintf(ctx->fd, "%%%%Origin: %ld %ld\n", (long) ox, (long) oy);
    /* NB: should use PageBoundingBox */
    fprintf(ctx->fd, "%%%%BoundingBox: 0 0 %ld %ld\n",
        (long) ceil(pw), (long) ceil(ph));
    fprintf(ctx->fd, "%%%%LanguageLevel: %d\n",
        (ctx->level3 ? 3 : (ctx->level2 ? 2 : 1)));
    fprintf(ctx->fd, "%%%%Pages: (atend)\n");
    fprintf(ctx->fd, "%%%%EndComments\n");
    fprintf(ctx->fd, "%%%%BeginSetup\n");
    if (ctx->PSduplex)
        fprintf(ctx->fd, "%s", DuplexPreamble);
    if (ctx->PStumble)
        fprintf(ctx->fd, "%s", TumblePreamble);
    if (ctx->PSavoiddeadzone && (ctx->level2 || ctx->level3))
        fprintf(ctx->fd, "%s", AvoidDeadZonePreamble);
    fprintf(ctx->fd, "%%%%EndSetup\n");
}

static void
PSTail(TIFF2PSContext *ctx)
{
    if (!ctx->npages)
        return;
    fprintf(ctx->fd, "%%%%Trailer\n");
    fprintf(ctx->fd, "%%%%Pages: %d\n", ctx->npages);
    fprintf(ctx->fd, "%%%%EOF\n");
}

static int
checkcmap(TIFF2PSContext* ctx, TIFF* tif, int n,
      uint16* r, uint16* g, uint16* b)
{
    (void) tif;
    while (n-- > 0)
        if (*r++ >= 256 || *g++ >= 256 || *b++ >= 256)
            return (16);
    TIFFWarning(ctx->filename, "Assuming 8-bit colormap");
    return (8);
}

static void
PS_Lvl2colorspace(TIFF2PSContext* ctx, TIFF* tif)
{
    uint16 *rmap, *gmap, *bmap;
    int i, num_colors;
    const char * colorspace_p;

    switch ( ctx->photometric )
    {
    case PHOTOMETRIC_SEPARATED:
        colorspace_p = "CMYK";
        break;

    case PHOTOMETRIC_RGB:
        colorspace_p = "RGB";
        break;

    default:
        colorspace_p = "Gray";
    }

    /*
     * Set up PostScript Level 2 colorspace according to
     * section 4.8 in the PostScript refenence manual.
     */
    fputs("% PostScript Level 2 only.\n", ctx->fd);
    if (ctx->photometric != PHOTOMETRIC_PALETTE) {
        if (ctx->photometric == PHOTOMETRIC_YCBCR) {
            /* MORE CODE HERE */
        }
        fprintf(ctx->fd, "/Device%s setcolorspace\n", colorspace_p );
        return;
    }

    /*
     * Set up an indexed/palette colorspace
     */
    num_colors = (1 << ctx->bitspersample);
    if (!TIFFGetField(tif, TIFFTAG_COLORMAP, &rmap, &gmap, &bmap)) {
        TIFFError(ctx->filename,
            "Palette image w/o \"Colormap\" tag");
        return;
    }
    if (checkcmap(ctx, tif, num_colors, rmap, gmap, bmap) == 16) {
        /*
         * Convert colormap to 8-bits values.
         */
#define CVT(x)      (((x) * 255) / ((1L<<16)-1))
        for (i = 0; i < num_colors; i++) {
            rmap[i] = CVT(rmap[i]);
            gmap[i] = CVT(gmap[i]);
            bmap[i] = CVT(bmap[i]);
        }
#undef CVT
    }
    fprintf(ctx->fd, "[ /Indexed /DeviceRGB %d", num_colors - 1);
    if (ctx->ascii85) {
        Ascii85Init(ctx);
        fputs("\n<~", ctx->fd);
        ctx->ascii85breaklen -= 2;
    } else
        fputs(" <", ctx->fd);
    for (i = 0; i < num_colors; i++) {
        if (ctx->ascii85) {
            Ascii85Put(ctx, (unsigned char)rmap[i]);
            Ascii85Put(ctx, (unsigned char)gmap[i]);
            Ascii85Put(ctx, (unsigned char)bmap[i]);
        } else {
            fputs((i % 8) ? " " : "\n  ", ctx->fd);
            fprintf(ctx->fd, "%02x%02x%02x",
                rmap[i], gmap[i], bmap[i]);
        }
    }
    if (ctx->ascii85)
        Ascii85Flush(ctx);
    else
        fputs(">\n", ctx->fd);
    fputs("] setcolorspace\n", ctx->fd);
}

static int
PS_Lvl2ImageDict(TIFF2PSContext* ctx, TIFF* tif, uint32 w, uint32 h)
{
    int use_rawdata;
    uint32 tile_width, tile_height;
    uint16 predictor, minsamplevalue, maxsamplevalue;
    int repeat_count;
    char im_h[64], im_x[64], im_y[64];

    (void)strcpy(im_x, "0");
    (void)sprintf(im_y, "%lu", (long) h);
    (void)sprintf(im_h, "%lu", (long) h);
    tile_width = w;
    tile_height = h;
    if (TIFFIsTiled(tif)) {
        repeat_count = TIFFNumberOfTiles(tif);
        TIFFGetField(tif, TIFFTAG_TILEWIDTH, &tile_width);
        TIFFGetField(tif, TIFFTAG_TILELENGTH, &tile_height);
        if (tile_width > w || tile_height > h ||
            (w % tile_width) != 0 || (h % tile_height != 0)) {
            /*
             * The tiles does not fit image width and height.
             * Set up a clip rectangle for the image unit square.
             */
            fputs("0 0 1 1 rectclip\n", ctx->fd);
        }
        if (tile_width < w) {
            fputs("/im_x 0 def\n", ctx->fd);
            (void)strcpy(im_x, "im_x neg");
        }
        if (tile_height < h) {
            fputs("/im_y 0 def\n", ctx->fd);
            (void)sprintf(im_y, "%lu im_y sub", (unsigned long) h);
        }
    } else {
        repeat_count = ctx->tf_numberstrips;
        tile_height = ctx->tf_rowsperstrip;
        if (tile_height > h)
            tile_height = h;
        if (repeat_count > 1) {
            fputs("/im_y 0 def\n", ctx->fd);
            fprintf(ctx->fd, "/im_h %lu def\n",
                (unsigned long) tile_height);
            (void)strcpy(im_h, "im_h");
            (void)sprintf(im_y, "%lu im_y sub", (unsigned long) h);
        }
    }

    /*
     * Output start of exec block
     */
    fputs("{ % exec\n", ctx->fd);

    if (repeat_count > 1)
        fprintf(ctx->fd, "%d { %% repeat\n", repeat_count);

    /*
     * Output filter options and image dictionary.
     */
    if (ctx->ascii85)
        fputs(" /im_stream currentfile /ASCII85Decode filter def\n",
              ctx->fd);
    fputs(" <<\n", ctx->fd);
    fputs("  /ImageType 1\n", ctx->fd);
    fprintf(ctx->fd, "  /Width %lu\n", (unsigned long) tile_width);
    /*
     * Workaround for some software that may crash when last strip
     * of image contains fewer number of scanlines than specified
     * by the `/Height' variable. So for stripped images with multiple
     * strips we will set `/Height' as `im_h', because one is 
     * recalculated for each strip - including the (smaller) final strip.
     * For tiled images and images with only one strip `/Height' will
     * contain number of scanlines in tile (or image height in case of
     * one-stripped image).
     */
    if (TIFFIsTiled(tif) || ctx->tf_numberstrips == 1)
        fprintf(ctx->fd, "  /Height %lu\n", (unsigned long) tile_height);
    else
        fprintf(ctx->fd, "  /Height im_h\n");

    if (ctx->planarconfiguration == PLANARCONFIG_SEPARATE && ctx->samplesperpixel > 1)
        fputs("  /MultipleDataSources true\n", ctx->fd);
    fprintf(ctx->fd, "  /ImageMatrix [ %lu 0 0 %ld %s %s ]\n",
        (unsigned long) w, - (long)h, im_x, im_y);
    fprintf(ctx->fd, "  /BitsPerComponent %d\n", ctx->bitspersample);
    fprintf(ctx->fd, "  /Ctx->Interpolate %s\n", ctx->interpolate ? "true" : "false");

    switch (ctx->samplesperpixel - ctx->extrasamples) {
    case 1:
        switch (ctx->photometric) {
        case PHOTOMETRIC_MINISBLACK:
            fputs("  /Decode [0 1]\n", ctx->fd);
            break;
        case PHOTOMETRIC_MINISWHITE:
            switch (ctx->compression) {
            case COMPRESSION_CCITTRLE:
            case COMPRESSION_CCITTRLEW:
            case COMPRESSION_CCITTFAX3:
            case COMPRESSION_CCITTFAX4:
                /*
                 * Manage inverting with /Blackis1 flag
                 * since there migth be uncompressed parts
                 */
                fputs("  /Decode [0 1]\n", ctx->fd);
                break;
            default:
                /*
                 * ERROR...
                 */
                fputs("  /Decode [1 0]\n", ctx->fd);
                break;
            }
            break;
        case PHOTOMETRIC_PALETTE:
            TIFFGetFieldDefaulted(tif, TIFFTAG_MINSAMPLEVALUE,
                &minsamplevalue);
            TIFFGetFieldDefaulted(tif, TIFFTAG_MAXSAMPLEVALUE,
                &maxsamplevalue);
            fprintf(ctx->fd, "  /Decode [%u %u]\n",
                    minsamplevalue, maxsamplevalue);
            break;
        default:
            /*
             * ERROR ?
             */
            fputs("  /Decode [0 1]\n", ctx->fd);
            break;
        }
        break;
    case 3:
        switch (ctx->photometric) {
        case PHOTOMETRIC_RGB:
            fputs("  /Decode [0 1 0 1 0 1]\n", ctx->fd);
            break;
        case PHOTOMETRIC_MINISWHITE:
        case PHOTOMETRIC_MINISBLACK:
        default:
            /*
             * ERROR??
             */
            fputs("  /Decode [0 1 0 1 0 1]\n", ctx->fd);
            break;
        }
        break;
    case 4:
        /*
         * ERROR??
         */
        fputs("  /Decode [0 1 0 1 0 1 0 1]\n", ctx->fd);
        break;
    }
    fputs("  /DataSource", ctx->fd);
    if (ctx->planarconfiguration == PLANARCONFIG_SEPARATE &&
        ctx->samplesperpixel > 1)
        fputs(" [", ctx->fd);
    if (ctx->ascii85)
        fputs(" im_stream", ctx->fd);
    else
        fputs(" currentfile /ASCIIHexDecode filter", ctx->fd);

    use_rawdata = TRUE;
    switch (ctx->compression) {
    case COMPRESSION_NONE:      /* 1: uncompressed */
        break;
    case COMPRESSION_CCITTRLE:  /* 2: CCITT modified Huffman RLE */
    case COMPRESSION_CCITTRLEW: /* 32771: #1 w/ word alignment */
    case COMPRESSION_CCITTFAX3: /* 3: CCITT Group 3 fax encoding */
    case COMPRESSION_CCITTFAX4: /* 4: CCITT Group 4 fax encoding */
        fputs("\n\t<<\n", ctx->fd);
        if (ctx->compression == COMPRESSION_CCITTFAX3) {
            uint32 g3_options;

            fputs("\t /EndOfLine true\n", ctx->fd);
            fputs("\t /EndOfBlock false\n", ctx->fd);
            if (!TIFFGetField(tif, TIFFTAG_GROUP3OPTIONS,
                        &g3_options))
                g3_options = 0;
            if (g3_options & GROUP3OPT_2DENCODING)
                fprintf(ctx->fd, "\t /K %s\n", im_h);
            if (g3_options & GROUP3OPT_UNCOMPRESSED)
                fputs("\t /Uncompressed true\n", ctx->fd);
            if (g3_options & GROUP3OPT_FILLBITS)
                fputs("\t /EncodedByteAlign true\n", ctx->fd);
        }
        if (ctx->compression == COMPRESSION_CCITTFAX4) {
            uint32 g4_options;

            fputs("\t /K -1\n", ctx->fd);
            TIFFGetFieldDefaulted(tif, TIFFTAG_GROUP4OPTIONS,
                           &g4_options);
            if (g4_options & GROUP4OPT_UNCOMPRESSED)
                fputs("\t /Uncompressed true\n", ctx->fd);
        }
        if (!(tile_width == w && w == 1728U))
            fprintf(ctx->fd, "\t /Columns %lu\n",
                (unsigned long) tile_width);
        fprintf(ctx->fd, "\t /Rows %s\n", im_h);
        if (ctx->compression == COMPRESSION_CCITTRLE ||
            ctx->compression == COMPRESSION_CCITTRLEW) {
            fputs("\t /EncodedByteAlign true\n", ctx->fd);
            fputs("\t /EndOfBlock false\n", ctx->fd);
        }
        if (ctx->photometric == PHOTOMETRIC_MINISBLACK)
            fputs("\t /BlackIs1 true\n", ctx->fd);
        fprintf(ctx->fd, "\t>> /CCITTFaxDecode filter");
        break;
    case COMPRESSION_LZW:   /* 5: Lempel-Ziv & Welch */
        TIFFGetFieldDefaulted(tif, TIFFTAG_PREDICTOR, &predictor);
        if (predictor == 2) {
            fputs("\n\t<<\n", ctx->fd);
            fprintf(ctx->fd, "\t /Predictor %u\n", predictor);
            fprintf(ctx->fd, "\t /Columns %lu\n",
                (unsigned long) tile_width);
            fprintf(ctx->fd, "\t /Colors %u\n", ctx->samplesperpixel);
            fprintf(ctx->fd, "\t /BitsPerComponent %u\n",
                ctx->bitspersample);
            fputs("\t>>", ctx->fd);
        }
        fputs(" /LZWDecode filter", ctx->fd);
        break;
    case COMPRESSION_DEFLATE:   /* 5: ZIP */
    case COMPRESSION_ADOBE_DEFLATE:
        if ( ctx->level3 ) {
             TIFFGetFieldDefaulted(tif, TIFFTAG_PREDICTOR, &predictor);
             if (predictor > 1) {
                fprintf(ctx->fd, "\t %% PostScript Level 3 only.");
                fputs("\n\t<<\n", ctx->fd);
                fprintf(ctx->fd, "\t /Predictor %u\n", predictor);
                fprintf(ctx->fd, "\t /Columns %lu\n",
                    (unsigned long) tile_width);
                fprintf(ctx->fd, "\t /Colors %u\n", ctx->samplesperpixel);
                    fprintf(ctx->fd, "\t /BitsPerComponent %u\n",
                    ctx->bitspersample);
                fputs("\t>>", ctx->fd);
             }
             fputs(" /FlateDecode filter", ctx->fd);
        } else {
            use_rawdata = FALSE ;
        }
        break;
    case COMPRESSION_PACKBITS:  /* 32773: Macintosh RLE */
        fputs(" /RunLengthDecode filter", ctx->fd);
        use_rawdata = TRUE;
        break;
    case COMPRESSION_OJPEG:     /* 6: !6.0 JPEG */
    case COMPRESSION_JPEG:      /* 7: %JPEG DCT ctx->compression */
#ifdef notdef
        /*
         * Code not tested yet
         */
        fputs(" /DCTDecode filter", ctx->fd);
        use_rawdata = TRUE;
#else
        use_rawdata = FALSE;
#endif
        break;
    case COMPRESSION_NEXT:      /* 32766: NeXT 2-bit RLE */
    case COMPRESSION_THUNDERSCAN:   /* 32809: ThunderScan RLE */
    case COMPRESSION_PIXARFILM: /* 32908: Pixar companded 10bit LZW */
    case COMPRESSION_JBIG:      /* 34661: ISO JBIG */
        use_rawdata = FALSE;
        break;
    case COMPRESSION_SGILOG:    /* 34676: SGI LogL or LogLuv */
    case COMPRESSION_SGILOG24:  /* 34677: SGI 24-bit LogLuv */
        use_rawdata = FALSE;
        break;
    default:
        /*
         * ERROR...
         */
        use_rawdata = FALSE;
        break;
    }
    if (ctx->planarconfiguration == PLANARCONFIG_SEPARATE &&
        ctx->samplesperpixel > 1) {
        uint16 i;

        /*
         * NOTE: This code does not work yet...
         */
        for (i = 1; i < ctx->samplesperpixel; i++)
            fputs(" dup", ctx->fd);
        fputs(" ]", ctx->fd);
    }

    fprintf( ctx->fd, "\n >> %s\n", IMAGEOP(ctx) );
    if (ctx->ascii85)
        fputs(" im_stream status { im_stream flushfile } if\n", ctx->fd);
    if (repeat_count > 1) {
        if (tile_width < w) {
            fprintf(ctx->fd, " /im_x im_x %lu add def\n",
                (unsigned long) tile_width);
            if (tile_height < h) {
                fprintf(ctx->fd, " im_x %lu ge {\n",
                    (unsigned long) w);
                fputs("  /im_x 0 def\n", ctx->fd);
                fprintf(ctx->fd, " /im_y im_y %lu add def\n",
                    (unsigned long) tile_height);
                fputs(" } if\n", ctx->fd);
            }
        }
        if (tile_height < h) {
            if (tile_width >= w) {
                fprintf(ctx->fd, " /im_y im_y %lu add def\n",
                    (unsigned long) tile_height);
                if (!TIFFIsTiled(tif)) {
                    fprintf(ctx->fd, " /im_h %lu im_y sub",
                        (unsigned long) h);
                    fprintf(ctx->fd, " dup %lu gt { pop",
                        (unsigned long) tile_height);
                    fprintf(ctx->fd, " %lu } if def\n",
                        (unsigned long) tile_height);
                }
            }
        }
        fputs("} repeat\n", ctx->fd);
    }
    /*
     * End of exec function
     */
    fputs("}\n", ctx->fd);

    return(use_rawdata);
}

#define MAXLINE     36

static int
PS_Lvl2page(TIFF2PSContext* ctx, TIFF* tif, uint32 w, uint32 h)
{
    uint16 fillorder;
    int use_rawdata, tiled_image, breaklen = MAXLINE;
    uint32 chunk_no, num_chunks, *bc;
    unsigned char *buf_data, *cp;
    tsize_t chunk_size, byte_count;

#if defined( EXP_ASCII85ENCODER )
    int         ascii85_l;  /* Length, in bytes, of ascii85_p[] data */
    uint8       *   ascii85_p = 0;  /* Holds ASCII85 encoded data */
#endif

    PS_Lvl2colorspace(ctx, tif);
    use_rawdata = PS_Lvl2ImageDict(ctx, tif, w, h);

/* See http://bugzilla.remotesensing.org/show_bug.cgi?id=80 */
#ifdef ENABLE_BROKEN_BEGINENDDATA
    fputs("%%BeginData:\n", ctx->fd);
#endif
    fputs("exec\n", ctx->fd);

    tiled_image = TIFFIsTiled(tif);
    if (tiled_image) {
        num_chunks = TIFFNumberOfTiles(tif);
        TIFFGetField(tif, TIFFTAG_TILEBYTECOUNTS, &bc);
    } else {
        num_chunks = TIFFNumberOfStrips(tif);
        TIFFGetField(tif, TIFFTAG_STRIPBYTECOUNTS, &bc);
    }

    if (use_rawdata) {
        chunk_size = (tsize_t) bc[0];
        for (chunk_no = 1; chunk_no < num_chunks; chunk_no++)
            if ((tsize_t) bc[chunk_no] > chunk_size)
                chunk_size = (tsize_t) bc[chunk_no];
    } else {
        if (tiled_image)
            chunk_size = TIFFTileSize(tif);
        else
            chunk_size = TIFFStripSize(tif);
    }
    buf_data = (unsigned char *)_TIFFmalloc(chunk_size);
    if (!buf_data) {
        TIFFError(ctx->filename, "Can't alloc %u bytes for %s.",
            (unsigned int) chunk_size, tiled_image ? "tiles" : "strips");
        return(FALSE);
    }

#if defined( EXP_ASCII85ENCODER )
    if ( ctx->ascii85 ) {
        /*
         * Allocate a buffer to hold the ASCII85 encoded data.  Note
         * that it is allocated with sufficient room to hold the
         * encoded data (5*chunk_size/4) plus the EOD marker (+8)
         * and formatting line breaks.  The line breaks are more
         * than taken care of by using 6*chunk_size/4 rather than
         * 5*chunk_size/4.
         */

        ascii85_p = _TIFFmalloc( (chunk_size+(chunk_size/2)) + 8 );

        if ( !ascii85_p ) {
        _TIFFfree( buf_data );

        TIFFError( ctx->filename,
               "Cannot allocate ASCII85 encoding buffer." );
        return ( FALSE );
        }
    }
#endif

    TIFFGetFieldDefaulted(tif, TIFFTAG_FILLORDER, &fillorder);
    for (chunk_no = 0; chunk_no < num_chunks; chunk_no++) {
        if (ctx->ascii85)
            Ascii85Init(ctx);
        else
            breaklen = MAXLINE;
        if (use_rawdata) {
            if (tiled_image)
                byte_count = TIFFReadRawTile(tif, chunk_no,
                          buf_data, chunk_size);
            else
                byte_count = TIFFReadRawStrip(tif, chunk_no,
                          buf_data, chunk_size);
            if (fillorder == FILLORDER_LSB2MSB)
                TIFFReverseBits(buf_data, byte_count);
        } else {
            if (tiled_image)
                byte_count = TIFFReadEncodedTile(tif,
                        chunk_no, buf_data,
                        chunk_size);
            else
                byte_count = TIFFReadEncodedStrip(tif,
                        chunk_no, buf_data,
                        chunk_size);
        }
        if (byte_count < 0) {
            TIFFError(ctx->filename, "Can't read %s %d.",
                tiled_image ? "tile" : "strip", chunk_no);
            if (ctx->ascii85)
                Ascii85Put(ctx, '\0');
        }
        /*
         * For images with ctx->alpha, matte against a white background;
         * i.e. Cback * (1 - Aimage) where Cback = 1. We will fill the
         * lower part of the buffer with the modified values.
         *
         * XXX: needs better solution
         */
        if (ctx->alpha) {
            int adjust, i, j = 0;
            int ncomps = ctx->samplesperpixel - ctx->extrasamples;
            for (i = 0; i < byte_count; i+=ctx->samplesperpixel) {
                adjust = 255 - buf_data[i + ncomps];
                switch (ncomps) {
                    case 1:
                        buf_data[j++] = buf_data[i] + adjust;
                        break;
                    case 2:
                        buf_data[j++] = buf_data[i] + adjust;
                        buf_data[j++] = buf_data[i+1] + adjust;
                        break;
                    case 3:
                        buf_data[j++] = buf_data[i] + adjust;
                        buf_data[j++] = buf_data[i+1] + adjust;
                        buf_data[j++] = buf_data[i+2] + adjust;
                        break;
                }
            }
            byte_count -= j;
        }

        if (ctx->ascii85) {
#if defined( EXP_ASCII85ENCODER )
            ascii85_l = Ascii85EncodeBlock(ctx, ascii85_p, 1,
                               buf_data, byte_count);

            if ( ascii85_l > 0 )
                fwrite( ascii85_p, ascii85_l, 1, ctx->fd );
#else
            for (cp = buf_data; byte_count > 0; byte_count--)
                Ascii85Put(ctx, *cp++);
#endif
        }
        else
        {
            for (cp = buf_data; byte_count > 0; byte_count--) {
                putc(hex[((*cp)>>4)&0xf], ctx->fd);
                putc(hex[(*cp)&0xf], ctx->fd);
                cp++;

                if (--breaklen <= 0) {
                    putc('\n', ctx->fd);
                    breaklen = MAXLINE;
                }
            }
        }

        if ( !ctx->ascii85 ) {
            if ( ctx->level2 || ctx->level3 )
                putc( '>', ctx->fd );
            putc('\n', ctx->fd);
        }
#if !defined( EXP_ASCII85ENCODER )
        else
            Ascii85Flush(ctx);
#endif
    }

#if defined( EXP_ASCII85ENCODER )
    if ( ascii85_p )
        _TIFFfree( ascii85_p );
#endif
    _TIFFfree(buf_data);
#ifdef ENABLE_BROKEN_BEGINENDDATA
    fputs("%%EndData\n", ctx->fd);
#endif
    return(TRUE);
}

void
PSpage(TIFF2PSContext* ctx, TIFF* tif, uint32 w, uint32 h)
{
    if ((ctx->level2 || ctx->level3) && PS_Lvl2page(ctx, tif, w, h))
        return;
    ctx->ps_bytesperrow = ctx->tf_bytesperrow - (ctx->extrasamples * ctx->bitspersample / 8)*w;
    switch (ctx->photometric) {
    case PHOTOMETRIC_RGB:
        if (ctx->planarconfiguration == PLANARCONFIG_CONTIG) {
            fprintf(ctx->fd, "%s", RGBcolorimage);
            PSColorContigPreamble(ctx, w, h, 3);
            PSDataColorContig(ctx, tif, w, h, 3);
        } else {
            PSColorSeparatePreamble(ctx, w, h, 3);
            PSDataColorSeparate(ctx, tif, w, h, 3);
        }
        break;
    case PHOTOMETRIC_SEPARATED:
        /* XXX should emit CMYKcolorimage */
        if (ctx->planarconfiguration == PLANARCONFIG_CONTIG) {
            PSColorContigPreamble(ctx, w, h, 4);
            PSDataColorContig(ctx, tif, w, h, 4);
        } else {
            PSColorSeparatePreamble(ctx, w, h, 4);
            PSDataColorSeparate(ctx, tif, w, h, 4);
        }
        break;
    case PHOTOMETRIC_PALETTE:
        fprintf(ctx->fd, "%s", RGBcolorimage);
        PhotoshopBanner(ctx, w, h, 1, 3, "false 3 colorimage");
        fprintf(ctx->fd, "/scanLine %ld string def\n",
            (long) ctx->ps_bytesperrow * 3L);
        fprintf(ctx->fd, "%lu %lu 8\n",
            (unsigned long) w, (unsigned long) h);
        fprintf(ctx->fd, "[%lu 0 0 -%lu 0 %lu]\n",
            (unsigned long) w, (unsigned long) h,
            (unsigned long) h);
        fprintf(ctx->fd,
            "{currentfile scanLine readhexstring pop} bind\n");
        fprintf(ctx->fd, "false 3 colorimage\n");
        PSDataPalette(ctx, tif, w, h);
        break;
    case PHOTOMETRIC_MINISBLACK:
    case PHOTOMETRIC_MINISWHITE:
        PhotoshopBanner(ctx, w, h, 1, 1, IMAGEOP(ctx));
        fprintf(ctx->fd, "/scanLine %ld string def\n",
            (long) ctx->ps_bytesperrow);
        fprintf(ctx->fd, "%lu %lu %d\n",
            (unsigned long) w, (unsigned long) h, ctx->bitspersample);
        fprintf(ctx->fd, "[%lu 0 0 -%lu 0 %lu]\n",
            (unsigned long) w, (unsigned long) h, (unsigned long) h);
        fprintf(ctx->fd,
            "{currentfile scanLine readhexstring pop} bind\n");
        fprintf(ctx->fd, "%s\n", IMAGEOP(ctx));
        PSDataBW(ctx, tif, w, h);
        break;
    }
    putc('\n', ctx->fd);
}

void
PSColorContigPreamble(TIFF2PSContext* ctx, uint32 w, uint32 h, int nc)
{
    ctx->ps_bytesperrow = nc * (ctx->tf_bytesperrow / ctx->samplesperpixel);
    PhotoshopBanner(ctx, w, h, 1, nc, "false %d colorimage", nc);
    fprintf(ctx->fd, "/line %ld string def\n", (long) ctx->ps_bytesperrow);
    fprintf(ctx->fd, "%lu %lu %d\n",
        (unsigned long) w, (unsigned long) h, ctx->bitspersample);
    fprintf(ctx->fd, "[%lu 0 0 -%lu 0 %lu]\n",
        (unsigned long) w, (unsigned long) h, (unsigned long) h);
    fprintf(ctx->fd, "{currentfile line readhexstring pop} bind\n");
    fprintf(ctx->fd, "false %d colorimage\n", nc);
}

void
PSColorSeparatePreamble(TIFF2PSContext* ctx, uint32 w, uint32 h, int nc)
{
    int i;

    PhotoshopBanner(ctx, w, h, ctx->ps_bytesperrow, nc, "true %d colorimage", nc);
    for (i = 0; i < nc; i++)
        fprintf(ctx->fd, "/line%d %ld string def\n",
            i, (long) ctx->ps_bytesperrow);
    fprintf(ctx->fd, "%lu %lu %d\n",
        (unsigned long) w, (unsigned long) h, ctx->bitspersample);
    fprintf(ctx->fd, "[%lu 0 0 -%lu 0 %lu] \n",
        (unsigned long) w, (unsigned long) h, (unsigned long) h);
    for (i = 0; i < nc; i++)
        fprintf(ctx->fd, "{currentfile line%d readhexstring pop}bind\n", i);
    fprintf(ctx->fd, "true %d colorimage\n", nc);
}

#define DOBREAK(len, howmany, fd) \
    if (((len) -= (howmany)) <= 0) {    \
        putc('\n', fd);         \
        (len) = MAXLINE-(howmany);  \
    }
#define PUTHEX(c,fd)    putc(hex[((c)>>4)&0xf],fd); putc(hex[(c)&0xf],fd)

void
PSDataColorContig(TIFF2PSContext* ctx, TIFF* tif, uint32 w, uint32 h, int nc)
{
    uint32 row;
    int breaklen = MAXLINE, cc, es = ctx->samplesperpixel - nc;
    unsigned char *tf_buf;
    unsigned char *cp, c;

    (void) w;
    tf_buf = (unsigned char *) _TIFFmalloc(ctx->tf_bytesperrow);
    if (tf_buf == NULL) {
        TIFFError(ctx->filename, "No space for scanline buffer");
        return;
    }
    for (row = 0; row < h; row++) {
        if (TIFFReadScanline(tif, tf_buf, row, 0) < 0)
            break;
        cp = tf_buf;
        if (ctx->alpha) {
            int adjust;
            cc = 0;
            for (; cc < ctx->tf_bytesperrow; cc += ctx->samplesperpixel) {
                DOBREAK(breaklen, nc, ctx->fd);
                /*
                 * For images with ctx->alpha, matte against
                 * a white background; i.e.
                 *    Cback * (1 - Aimage)
                 * where Cback = 1.
                 */
                adjust = 255 - cp[nc];
                switch (nc) {
                case 4: c = *cp++ + adjust; PUTHEX(c,ctx->fd);
                case 3: c = *cp++ + adjust; PUTHEX(c,ctx->fd);
                case 2: c = *cp++ + adjust; PUTHEX(c,ctx->fd);
                case 1: c = *cp++ + adjust; PUTHEX(c,ctx->fd);
                }
                cp += es;
            }
        } else {
            cc = 0;
            for (; cc < ctx->tf_bytesperrow; cc += ctx->samplesperpixel) {
                DOBREAK(breaklen, nc, ctx->fd);
                switch (nc) {
                case 4: c = *cp++; PUTHEX(c,ctx->fd);
                case 3: c = *cp++; PUTHEX(c,ctx->fd);
                case 2: c = *cp++; PUTHEX(c,ctx->fd);
                case 1: c = *cp++; PUTHEX(c,ctx->fd);
                }
                cp += es;
            }
        }
    }
    _TIFFfree((char *) tf_buf);
}

void
PSDataColorSeparate(TIFF2PSContext* ctx, TIFF* tif, uint32 w, uint32 h, int nc)
{
    uint32 row;
    int breaklen = MAXLINE, cc;
    tsample_t s, maxs;
    unsigned char *tf_buf;
    unsigned char *cp, c;

    (void) w;
    tf_buf = (unsigned char *) _TIFFmalloc(ctx->tf_bytesperrow);
    if (tf_buf == NULL) {
        TIFFError(ctx->filename, "No space for scanline buffer");
        return;
    }
    maxs = (ctx->samplesperpixel > nc ? nc : ctx->samplesperpixel);
    for (row = 0; row < h; row++) {
        for (s = 0; s < maxs; s++) {
            if (TIFFReadScanline(tif, tf_buf, row, s) < 0)
                break;
            for (cp = tf_buf, cc = 0; cc < ctx->tf_bytesperrow; cc++) {
                DOBREAK(breaklen, 1, ctx->fd);
                c = *cp++;
                PUTHEX(c,ctx->fd);
            }
        }
    }
    _TIFFfree((char *) tf_buf);
}

#define PUTRGBHEX(c,fd) \
    PUTHEX(rmap[c],fd); PUTHEX(gmap[c],fd); PUTHEX(bmap[c],fd)

void
PSDataPalette(TIFF2PSContext* ctx, TIFF* tif, uint32 w, uint32 h)
{
    uint16 *rmap, *gmap, *bmap;
    uint32 row;
    int breaklen = MAXLINE, cc, nc;
    unsigned char *tf_buf;
    unsigned char *cp, c;

    (void) w;
    if (!TIFFGetField(tif, TIFFTAG_COLORMAP, &rmap, &gmap, &bmap)) {
        TIFFError(ctx->filename, "Palette image w/o \"Colormap\" tag");
        return;
    }
    switch (ctx->bitspersample) {
    case 8: case 4: case 2: case 1:
        break;
    default:
        TIFFError(ctx->filename, "Depth %d not supported", ctx->bitspersample);
        return;
    }
    nc = 3 * (8 / ctx->bitspersample);
    tf_buf = (unsigned char *) _TIFFmalloc(ctx->tf_bytesperrow);
    if (tf_buf == NULL) {
        TIFFError(ctx->filename, "No space for scanline buffer");
        return;
    }
    if (checkcmap(ctx, tif, 1<<ctx->bitspersample, rmap, gmap, bmap) == 16) {
        int i;
#define CVT(x)      ((unsigned short) (((x) * 255) / ((1U<<16)-1)))
        for (i = (1<<ctx->bitspersample)-1; i >= 0; i--) {
            rmap[i] = CVT(rmap[i]);
            gmap[i] = CVT(gmap[i]);
            bmap[i] = CVT(bmap[i]);
        }
#undef CVT
    }
    for (row = 0; row < h; row++) {
        if (TIFFReadScanline(tif, tf_buf, row, 0) < 0)
            break;
        for (cp = tf_buf, cc = 0; cc < ctx->tf_bytesperrow; cc++) {
            DOBREAK(breaklen, nc, ctx->fd);
            switch (ctx->bitspersample) {
            case 8:
                c = *cp++; PUTRGBHEX(c, ctx->fd);
                break;
            case 4:
                c = *cp++; PUTRGBHEX(c&0xf, ctx->fd);
                c >>= 4;   PUTRGBHEX(c, ctx->fd);
                break;
            case 2:
                c = *cp++; PUTRGBHEX(c&0x3, ctx->fd);
                c >>= 2;   PUTRGBHEX(c&0x3, ctx->fd);
                c >>= 2;   PUTRGBHEX(c&0x3, ctx->fd);
                c >>= 2;   PUTRGBHEX(c, ctx->fd);
                break;
            case 1:
                c = *cp++; PUTRGBHEX(c&0x1, ctx->fd);
                c >>= 1;   PUTRGBHEX(c&0x1, ctx->fd);
                c >>= 1;   PUTRGBHEX(c&0x1, ctx->fd);
                c >>= 1;   PUTRGBHEX(c&0x1, ctx->fd);
                c >>= 1;   PUTRGBHEX(c&0x1, ctx->fd);
                c >>= 1;   PUTRGBHEX(c&0x1, ctx->fd);
                c >>= 1;   PUTRGBHEX(c&0x1, ctx->fd);
                c >>= 1;   PUTRGBHEX(c, ctx->fd);
                break;
            }
        }
    }
    _TIFFfree((char *) tf_buf);
}

void
PSDataBW(TIFF2PSContext* ctx, TIFF* tif, uint32 w, uint32 h)
{
    int breaklen = MAXLINE;
    unsigned char* tf_buf;
    unsigned char* cp;
    tsize_t stripsize = TIFFStripSize(tif);
    tstrip_t s;

#if defined( EXP_ASCII85ENCODER )
    int ascii85_l;      /* Length, in bytes, of ascii85_p[] data */
    uint8   *ascii85_p = 0;     /* Holds ASCII85 encoded data */
#endif

    (void) w; (void) h;
    tf_buf = (unsigned char *) _TIFFmalloc(stripsize);
        memset(tf_buf, 0, stripsize);
    if (tf_buf == NULL) {
        TIFFError(ctx->filename, "No space for scanline buffer");
        return;
    }

#if defined( EXP_ASCII85ENCODER )
    if ( ctx->ascii85 ) {
        /*
         * Allocate a buffer to hold the ASCII85 encoded data.  Note
         * that it is allocated with sufficient room to hold the
         * encoded data (5*stripsize/4) plus the EOD marker (+8)
         * and formatting line breaks.  The line breaks are more
         * than taken care of by using 6*stripsize/4 rather than
         * 5*stripsize/4.
         */

        ascii85_p = _TIFFmalloc( (stripsize+(stripsize/2)) + 8 );

        if ( !ascii85_p ) {
        _TIFFfree( tf_buf );

        TIFFError( ctx->filename,
               "Cannot allocate ASCII85 encoding buffer." );
        return;
        }
    }
#endif

    if (ctx->ascii85)
        Ascii85Init(ctx);

    for (s = 0; s < TIFFNumberOfStrips(tif); s++) {
        int cc = TIFFReadEncodedStrip(tif, s, tf_buf, stripsize);
        if (cc < 0) {
            TIFFError(ctx->filename, "Can't read strip");
            break;
        }
        cp = tf_buf;
        if (ctx->photometric == PHOTOMETRIC_MINISWHITE) {
            for (cp += cc; --cp >= tf_buf;)
                *cp = ~*cp;
            cp++;
        }
        if (ctx->ascii85) {
#if defined( EXP_ASCII85ENCODER )
            if (ctx->alpha) {
                int adjust, i;
                for (i = 0; i < cc; i+=2) {
                    adjust = 255 - cp[i + 1];
                    cp[i / 2] = cp[i] + adjust;
                }
                cc /= 2;
            }

            ascii85_l = Ascii85EncodeBlock(ctx, ascii85_p, 1, cp,
                               cc);

            if ( ascii85_l > 0 )
                fwrite( ascii85_p, ascii85_l, 1, ctx->fd );
#else
            while (cc-- > 0)
                Ascii85Put(ctx, *cp++);
#endif /* EXP_ASCII85_ENCODER */
        } else {
            unsigned char c;

            if (ctx->alpha) {
                int adjust;
                while (cc-- > 0) {
                    DOBREAK(breaklen, 1, ctx->fd);
                    /*
                     * For images with ctx->alpha, matte against
                     * a white background; i.e.
                     *    Cback * (1 - Aimage)
                     * where Cback = 1.
                     */
                    adjust = 255 - cp[1];
                    c = *cp++ + adjust; PUTHEX(c,ctx->fd);
                    cp++, cc--;
                }
            } else {
                while (cc-- > 0) {
                    c = *cp++;
                    DOBREAK(breaklen, 1, ctx->fd);
                    PUTHEX(c, ctx->fd);
                }
            }
        }
    }

    if ( !ctx->ascii85 )
    {
        if ( ctx->level2 || ctx->level3)
        fputs(">\n", ctx->fd);
    }
#if !defined( EXP_ASCII85ENCODER )
    else
        Ascii85Flush(ctx);
#else
    if ( ascii85_p )
        _TIFFfree( ascii85_p );
#endif

    _TIFFfree(tf_buf);
}

static void
Ascii85Init(TIFF2PSContext *ctx)
{
    ctx->ascii85breaklen = 2*MAXLINE;
    ctx->ascii85count = 0;
}

static void
Ascii85Encode(unsigned char* raw, char *buf)
{
    uint32 word;

    word = (((raw[0]<<8)+raw[1])<<16) + (raw[2]<<8) + raw[3];
    if (word != 0L) {
        uint32 q;
        uint16 w1;

        q = word / (85L*85*85*85);  /* actually only a byte */
        buf[0] = (char) (q + '!');

        word -= q * (85L*85*85*85); q = word / (85L*85*85);
        buf[1] = (char) (q + '!');

        word -= q * (85L*85*85); q = word / (85*85);
        buf[2] = (char) (q + '!');

        w1 = (uint16) (word - q*(85L*85));
        buf[3] = (char) ((w1 / 85) + '!');
        buf[4] = (char) ((w1 % 85) + '!');
        buf[5] = '\0';
    } else
        buf[0] = 'z', buf[1] = '\0';
}

void
Ascii85Put(TIFF2PSContext *ctx, unsigned char code)
{
    ctx->ascii85buf[ctx->ascii85count++] = code;
    if (ctx->ascii85count >= 4) {
        unsigned char* p;
        int n;
        char buf[6];

        for (n = ctx->ascii85count, p = ctx->ascii85buf;
             n >= 4; n -= 4, p += 4) {
            char* cp;
            Ascii85Encode(p, buf);
            for (cp = buf; *cp; cp++) {
                putc(*cp, ctx->fd);
                if (--ctx->ascii85breaklen == 0) {
                    putc('\n', ctx->fd);
                    ctx->ascii85breaklen = 2*MAXLINE;
                }
            }
        }
        _TIFFmemcpy(ctx->ascii85buf, p, n);
        ctx->ascii85count = n;
    }
}

void
Ascii85Flush(TIFF2PSContext* ctx)
{
    if (ctx->ascii85count > 0) {
        char res[6];
        _TIFFmemset(&ctx->ascii85buf[ctx->ascii85count], 0, 3);
        Ascii85Encode(ctx->ascii85buf, res);
        fwrite(res[0] == 'z' ? "!!!!" : res, ctx->ascii85count + 1, 1, ctx->fd);
    }
    fputs("~>\n", ctx->fd);
}
#if defined( EXP_ASCII85ENCODER)

#define A85BREAKCNTR    ctx->ascii85breaklen
#define A85BREAKLEN     (2*MAXLINE)

/*****************************************************************************
*
* Name:         Ascii85EncodeBlock( ascii85_p, f_eod, raw_p, raw_l )
*
* Description:  This routine will encode the raw data in the buffer described
*               by raw_p and raw_l into ASCII85 format and store the encoding
*               in the buffer given by ascii85_p.
*
* Parameters:   ctx         -   TIFF2PS context
*               ascii85_p   -   A buffer supplied by the caller which will
*                               contain the encoded ASCII85 data.
*               f_eod       -   Flag: Nz means to end the encoded buffer with
*                               an End-Of-Data marker.
*               raw_p       -   Pointer to the buffer of data to be encoded
*               raw_l       -   Number of bytes in raw_p[] to be encoded
*
* Returns:      (int)   <   0   Error, see errno
*                       >=  0   Number of bytes written to ascii85_p[].
*
* Notes:        An external variable given by A85BREAKCNTR is used to
*               determine when to insert newline characters into the
*               encoded data.  As each byte is placed into ascii85_p this
*               external is decremented.  If the variable is decrement to
*               or past zero then a newline is inserted into ascii85_p
*               and the A85BREAKCNTR is then reset to A85BREAKLEN.
*                   Note:  for efficiency reasons the A85BREAKCNTR variable
*                          is not actually checked on *every* character
*                          placed into ascii85_p but often only for every
*                          5 characters.
*
*               THE CALLER IS RESPONSIBLE FOR ENSURING THAT ASCII85_P[] IS
*               SUFFICIENTLY LARGE TO THE ENCODED DATA!
*                   You will need at least 5 * (raw_l/4) bytes plus space for
*                   newline characters and space for an EOD marker (if
*                   requested).  A safe calculation is to use 6*(raw_l/4) + 8
*                   to size ascii85_p.
*
*****************************************************************************/

int Ascii85EncodeBlock( TIFF2PSContext *ctx, uint8 * ascii85_p,
            unsigned f_eod, const uint8 * raw_p, int raw_l )

{
    char                        ascii85[5];     /* Encoded 5 tuple */
    int                         ascii85_l;      /* Number of bytes written to ascii85_p[] */
    int                         rc;             /* Return code */
    uint32                      val32;          /* Unencoded 4 tuple */

    ascii85_l = 0;                              /* Nothing written yet */

    if ( raw_p )
    {
        --raw_p;                                /* Prepare for pre-increment fetches */

        for ( ; raw_l > 3; raw_l -= 4 )
        {
            val32  = *(++raw_p) << 24;
            val32 += *(++raw_p) << 16;
            val32 += *(++raw_p) <<  8;
            val32 += *(++raw_p);

            if ( val32 == 0 )                   /* Special case */
            {
                ascii85_p[ascii85_l] = 'z';
                rc = 1;
            }

            else
            {
                ascii85[4] = (char) ((val32 % 85) + 33);
                val32 /= 85;

                ascii85[3] = (char) ((val32 % 85) + 33);
                val32 /= 85;

                ascii85[2] = (char) ((val32 % 85) + 33);
                val32 /= 85;

                ascii85[1] = (char) ((val32 % 85) + 33);
                ascii85[0] = (char) ((val32 / 85) + 33);

                _TIFFmemcpy( &ascii85_p[ascii85_l], ascii85, sizeof(ascii85) );
                rc = sizeof(ascii85);
            }

            ascii85_l += rc;

            if ( (A85BREAKCNTR -= rc) <= 0 )
            {
                ascii85_p[ascii85_l] = '\n';
                ++ascii85_l;
                A85BREAKCNTR = A85BREAKLEN;
            }
        }

        /*
         * Output any straggler bytes:
         */

        if ( raw_l > 0 )
        {
            int             len;                /* Output this many bytes */

            len = raw_l + 1;
            val32 = *++raw_p << 24;             /* Prime the pump */

            if ( --raw_l > 0 )  val32 += *(++raw_p) << 16;
            if ( --raw_l > 0 )  val32 += *(++raw_p) <<  8;

            val32 /= 85;

            ascii85[3] = (char) ((val32 % 85) + 33);
            val32 /= 85;

            ascii85[2] = (char) ((val32 % 85) + 33);
            val32 /= 85;

            ascii85[1] = (char) ((val32 % 85) + 33);
            ascii85[0] = (char) ((val32 / 85) + 33);

            _TIFFmemcpy( &ascii85_p[ascii85_l], ascii85, len );
            ascii85_l += len;
        }
    }

    /*
     * If requested add an ASCII85 End Of Data marker:
     */

    if ( f_eod )
    {
        ascii85_p[ascii85_l++] = '~';
        ascii85_p[ascii85_l++] = '>';
        ascii85_p[ascii85_l++] = '\n';
    }

    return ( ascii85_l );

}   /* Ascii85EncodeBlock() */

#endif  /* EXP_ASCII85ENCODER */

/* vim: set ts=8 sts=8 sw=8 noet: */