/* pnmtorast.c - read a portable anymap and produce a Sun rasterfile
**
** Copyright (C) 1989, 1991 by Jef Poskanzer.
**
** 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. This software is provided "as is" without express or
** implied warranty.
*/
#include "pnm.h"
#include "rast.h"
#include "mallocvar.h"
#define MAXCOLORS 256
static colormap_t *
alloc_pr_colormap(void) {
colormap_t* pr_colormapP;
MALLOCVAR(pr_colormapP);
if ( pr_colormapP == NULL )
pm_error( "out of memory" );
pr_colormapP->type = RMT_EQUAL_RGB;
pr_colormapP->length = MAXCOLORS;
MALLOCARRAY(pr_colormapP->map[0], MAXCOLORS);
MALLOCARRAY(pr_colormapP->map[1], MAXCOLORS);
MALLOCARRAY(pr_colormapP->map[2], MAXCOLORS);
if ( pr_colormapP->map[0] == NULL ||
pr_colormapP->map[1] == NULL ||
pr_colormapP->map[2] == NULL )
pm_error( "out of memory" );
return pr_colormapP;
}
static colormap_t*
make_pr_colormap(colorhist_vector const chv,
int const colors) {
colormap_t* pr_colormapP;
int i;
pr_colormapP = alloc_pr_colormap( );
for ( i = 0; i < colors; ++i )
{
pr_colormapP->map[0][i] = PPM_GETR( chv[i].color );
pr_colormapP->map[1][i] = PPM_GETG( chv[i].color );
pr_colormapP->map[2][i] = PPM_GETB( chv[i].color );
}
for ( ; i < MAXCOLORS; ++i )
pr_colormapP->map[0][i] = pr_colormapP->map[1][i] =
pr_colormapP->map[2][i] = 0;
return pr_colormapP;
}
static colormap_t *
make_gray_pr_colormap(void) {
colormap_t* pr_colormapP;
int i;
pr_colormapP = alloc_pr_colormap( );
for ( i = 0; i < MAXCOLORS; ++i )
{
pr_colormapP->map[0][i] = i;
pr_colormapP->map[1][i] = i;
pr_colormapP->map[2][i] = i;
}
return pr_colormapP;
}
int
main(int argc, char ** argv) {
FILE* ifp;
xel** xels;
xel* xelrow;
xel p;
register xel* xP;
colorhist_vector chv;
colorhash_table cht;
colormap_t* pr_colormapP;
int argn, pr_type, rows, cols, format, i;
int depth, colors, linesize, row;
register int col, bitcount;
xelval maxval;
struct pixrect* pr;
unsigned char* data;
register unsigned char* byteP;
const char* const usage = "[-standard|-rle] [pnmfile]";
pnm_init( &argc, argv );
argn = 1;
pr_type = RT_BYTE_ENCODED;
while ( argn < argc && argv[argn][0] == '-' && argv[argn][1] != '\0' )
{
if ( pm_keymatch( argv[argn], "-standard", 2 ) )
pr_type = RT_STANDARD;
else if ( pm_keymatch( argv[argn], "-rle", 2 ) )
pr_type = RT_BYTE_ENCODED;
else
pm_usage( usage );
++argn;
}
if ( argn != argc )
{
ifp = pm_openr( argv[argn] );
++argn;
}
else
ifp = stdin;
if ( argn != argc )
pm_usage( usage );
xels = pnm_readpnm( ifp, &cols, &rows, &maxval, &format );
pm_close( ifp );
/* Figure out the proper depth and colormap. */
switch ( PNM_FORMAT_TYPE(format) )
{
case PPM_TYPE:
pm_message( "computing colormap..." );
chv = ppm_computecolorhist( xels, cols, rows, MAXCOLORS, &colors );
if ( chv == (colorhist_vector) 0 )
{
pm_message(
"Too many colors - proceeding to write a 24-bit non-mapped" );
pm_message(
"rasterfile. If you want 8 bits, try doing a 'pnmquant %d'.",
MAXCOLORS );
depth = 24;
pr_type = RT_STANDARD;
pr_colormapP = (colormap_t*) 0;
}
else
{
pm_message( "%d colors found", colors );
if ( maxval != 255 )
for ( i = 0; i < colors; ++i )
PPM_DEPTH( chv[i].color, chv[i].color, maxval, 255 );
/* Force white to slot 0 and black to slot 1, if possible. */
PPM_ASSIGN( p, 255, 255, 255 );
ppm_addtocolorhist( chv, &colors, MAXCOLORS, &p, 0, 0 );
PPM_ASSIGN( p, 0, 0, 0 );
ppm_addtocolorhist( chv, &colors, MAXCOLORS, &p, 0, 1 );
if ( colors == 2 )
{
/* Monochrome. */
depth = 1;
pr_colormapP = (colormap_t*) 0;
}
else
{
/* Turn the ppm colormap into the appropriate Sun colormap. */
depth = 8;
pr_colormapP = make_pr_colormap( chv, colors );
}
cht = ppm_colorhisttocolorhash( chv, colors );
ppm_freecolorhist( chv );
}
break;
case PGM_TYPE:
depth = 8;
pr_colormapP = make_gray_pr_colormap( );
break;
default:
depth = 1;
pr_colormapP = (colormap_t*) 0;
break;
}
if ( maxval > 255 && depth != 1 )
pm_message(
"maxval is not 255 - automatically rescaling colors" );
/* Allocate space for the Sun-format image. */
if ( (pr = mem_create(cols, rows, depth)) == (struct pixrect*) 0 )
pm_error( "unable to create new pixrect" );
data = ( (struct mpr_data*) pr->pr_data )->md_image;
linesize = ( (struct mpr_data*) pr->pr_data )->md_linebytes;
/* And compute the Sun image. The variables at this point are:
** cht is null or not
** depth is 1, 8, or 24
*/
for ( row = 0; row < rows; ++row )
{
xelrow = xels[row];
byteP = data;
switch ( depth )
{
case 1:
*byteP = 0;
bitcount = 7;
for ( col = 0, xP = xelrow; col < cols; ++col, ++xP )
{
register int color;
switch ( PNM_FORMAT_TYPE(format) )
{
case PPM_TYPE:
if ( maxval != 255 )
PPM_DEPTH( *xP, *xP, maxval, 255 );
color = ppm_lookupcolor( cht, xP );
if ( color == -1 )
pm_error(
"color not found?!? row=%d col=%d r=%d g=%d b=%d",
row, col, PPM_GETR(*xP), PPM_GETG(*xP),
PPM_GETB(*xP) );
if ( color )
*byteP |= 1 << bitcount;
break;
default:
color = PNM_GET1( *xP );
if ( ! color )
*byteP |= 1 << bitcount;
break;
}
--bitcount;
if ( bitcount < 0 )
{
++byteP;
*byteP = 0;
bitcount = 7;
}
}
break;
case 8:
for ( col = 0, xP = xelrow; col < cols; ++col, ++xP )
{
register int color;
switch ( PNM_FORMAT_TYPE(format) )
{
case PPM_TYPE:
if ( maxval != 255 )
PPM_DEPTH( *xP, *xP, maxval, 255 );
color = ppm_lookupcolor( cht, xP );
if ( color == -1 )
pm_error(
"color not found?!? row=%d col=%d r=%d g=%d b=%d",
row, col, PPM_GETR(*xP), PPM_GETG(*xP),
PPM_GETB(*xP) );
break;
case PGM_TYPE:
color = PNM_GET1( *xP );
if ( maxval != 255 )
color = color * 255 / maxval;
break;
default:
color = PNM_GET1( *xP );
}
*byteP++ = color;
}
break;
case 24:
/* If depth is 24, we do NOT have a valid cht. */
for ( col = 0, xP = xelrow; col < cols; ++col, ++xP )
{
if ( maxval != 255 )
PPM_DEPTH( *xP, *xP, maxval, 255 );
*byteP++ = PPM_GETB( *xP );
*byteP++ = PPM_GETG( *xP );
*byteP++ = PPM_GETR( *xP );
}
break;
default:
pm_error( "can't happen" );
}
data += linesize;
}
pnm_freearray( xels, rows );
/* Finally, write the sucker out. */
if ( pr_dump( pr, stdout, pr_colormapP, pr_type, 0 ) == PIX_ERR )
pm_error( "error writing rasterfile" );
exit( 0 );
}