/* $Id: tiff-grayscale.c,v 1.6 2010-06-08 18:55:15 bfriesen Exp $ */

/*

 * tiff-grayscale.c -- create a Class G (grayscale) TIFF file

 *      with a gray response curve in linear optical density

 *

 * Copyright 1990 by Digital Equipment Corporation, Maynard, Massachusetts.

 *

 *                        All Rights Reserved

 *

 * Permission to use, copy, modify, and distribute this software and its

 * documentation for any purpose and without fee is hereby granted,

 * provided that the above copyright notice appear in all copies and that

 * both that copyright notice and this permission notice appear in

 * supporting documentation, and that the name of Digital not be

 * used in advertising or publicity pertaining to distribution of the

 * software without specific, written prior permission.

 *

 * DIGITAL DISCLAIMS ALL WARRANTIES WITH REGARD TO THIS SOFTWARE, INCLUDING

 * ALL IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS, IN NO EVENT SHALL

 * DIGITAL BE LIABLE FOR ANY SPECIAL, INDIRECT OR CONSEQUENTIAL DAMAGES OR

 * ANY DAMAGES WHATSOEVER RESULTING FROM LOSS OF USE, DATA OR PROFITS,

 * WHETHER IN AN ACTION OF CONTRACT, NEGLIGENCE OR OTHER TORTIOUS ACTION,

 * ARISING OUT OF OR IN CONNECTION WITH THE USE OR PERFORMANCE OF THIS

 * SOFTWARE.

 */

#include <math.h></math.h>

#include <stdio.h></stdio.h>

#include <stdlib.h></stdlib.h>

#include <string.h></string.h>

#include "tiffio.h"

#define WIDTH       512

#define HEIGHT      WIDTH

char *              programName;

void                Usage();

int main(int argc, char **argv)

{

    int             bits_per_pixel = 8, cmsize, i, j, k,

                    gray_index, chunk_size = 32, nchunks = 16;

    unsigned char * scan_line;

    uint16 *        gray;

    float           refblackwhite[2*1];

    TIFF *          tif;

    programName = argv[0];

    if (argc != 4)

        Usage();

    if (!strcmp(argv[1], "-depth"))

         bits_per_pixel = atoi(argv[2]);

    else

         Usage();

    switch (bits_per_pixel) {

        case 8:

            nchunks = 16;

            chunk_size = 32;

            break;

        case 4:

            nchunks = 4;

            chunk_size = 128;

            break;

        case 2:

            nchunks = 2;

            chunk_size = 256;

            break;

        default:

            Usage();

    }

    cmsize = nchunks * nchunks;

    gray = (uint16 *) malloc(cmsize * sizeof(uint16));

    gray[0] = 3000;

    for (i = 1; i < cmsize; i++)

        gray[i] = (uint16) (-log10((double) i / (cmsize - 1)) * 1000);

    refblackwhite[0] = 0.0;

    refblackwhite[1] = (float)((1L<

    if ((tif = TIFFOpen(argv[3], "w")) == NULL) {

        fprintf(stderr, "can't open %s as a TIFF file\n", argv[3]);

		free(gray);

        return 0;

    }

    TIFFSetField(tif, TIFFTAG_IMAGEWIDTH, WIDTH);

    TIFFSetField(tif, TIFFTAG_IMAGELENGTH, HEIGHT);

    TIFFSetField(tif, TIFFTAG_BITSPERSAMPLE, bits_per_pixel);

    TIFFSetField(tif, TIFFTAG_COMPRESSION, COMPRESSION_NONE);

    TIFFSetField(tif, TIFFTAG_PHOTOMETRIC, PHOTOMETRIC_MINISBLACK);

    TIFFSetField(tif, TIFFTAG_SAMPLESPERPIXEL, 1);

    TIFFSetField(tif, TIFFTAG_ROWSPERSTRIP, 1);

    TIFFSetField(tif, TIFFTAG_PLANARCONFIG, PLANARCONFIG_CONTIG);

    TIFFSetField(tif, TIFFTAG_REFERENCEBLACKWHITE, refblackwhite);

    TIFFSetField(tif, TIFFTAG_TRANSFERFUNCTION, gray);

    TIFFSetField(tif, TIFFTAG_RESOLUTIONUNIT, RESUNIT_NONE);

    scan_line = (unsigned char *) malloc(WIDTH / (8 / bits_per_pixel));

    for (i = 0; i < HEIGHT; i++) {

        for (j = 0, k = 0; j < WIDTH;) {

            gray_index = (j / chunk_size) + ((i / chunk_size) * nchunks);

            switch (bits_per_pixel) {

            case 8:

                scan_line[k++] = gray_index;

                j++;

                break;

            case 4:

                scan_line[k++] = (gray_index << 4) + gray_index;

                j += 2;

                break;

            case 2:

                scan_line[k++] = (gray_index << 6) + (gray_index << 4)

                    + (gray_index << 2) + gray_index;

                j += 4;

                break;

            }

        }

        TIFFWriteScanline(tif, scan_line, i, 0);

    }

    free(scan_line);

    TIFFClose(tif);

    return 0;

}

void

Usage()

{

    fprintf(stderr, "Usage: %s -depth (8 | 4 | 2) tiff-image\n", programName);

    exit(0);

}

/*

 * Local Variables:

 * mode: c

 * c-basic-offset: 8

 * fill-column: 78

 * End:

 */