home *** CD-ROM | disk | FTP | other *** search
/ Club Amiga de Montreal - CAM / CAM_CD_1.iso / files / 309.lha / PBM_PLUS / pbm / pbmreduce.c < prev    next >
C/C++ Source or Header  |  1980-12-04  |  5KB  |  203 lines
  1. /* pbmreduce.c - read a portable bitmap and reduce it N times
  2. **
  3. ** Copyright (C) 1989 by Jef Poskanzer.
  4. **
  5. ** Permission to use, copy, modify, and distribute this software and its
  6. ** documentation for any purpose and without fee is hereby granted, provided
  7. ** that the above copyright notice appear in all copies and that both that
  8. ** copyright notice and this permission notice appear in supporting
  9. ** documentation.  This software is provided "as is" without express or
  10. ** implied warranty.
  11. */
  12.  
  13. #include <stdio.h>
  14. #ifdef    SYSV
  15. #include <string.h>
  16. #define srandom srand
  17. #define random rand
  18. #else    SYSV
  19. #include <strings.h>
  20. #endif    SYSV
  21. #include "pbm.h"
  22.  
  23. #define max(a,b) ((a) > (b) ? (a) : (b))
  24.  
  25. main( argc, argv )
  26. int argc;
  27. char *argv[];
  28.     {
  29.     FILE *ifd;
  30.     register bit **bitslice, *newbitrow, *nbP;
  31.     int argn, n, rows, cols, format, newrows, newcols;
  32.     int row, col, limitcol, subrow, subcol, count, direction;
  33.     char *usage = "[-floyd|-fs | -threshold ] [-value <val>] N [pbmfile]";
  34.     int halftone;
  35. #define QT_FS 1
  36. #define QT_THRESH 2
  37. #define SCALE 1024
  38. #define HALFSCALE 512
  39.     long threshval, sum, *thiserr, *nexterr, *temperr;
  40.  
  41.     pm_progname = argv[0];
  42.  
  43.     argn = 1;
  44.     halftone = QT_FS;
  45.     threshval = HALFSCALE;
  46.  
  47.     while ( argn < argc && argv[argn][0] == '-' )
  48.     {
  49.     if ( strncmp(argv[argn],"-fs",max(strlen(argv[argn]),2)) == 0 ||
  50.          strncmp(argv[argn],"-floyd",max(strlen(argv[argn]),2)) == 0 )
  51.         halftone = QT_FS;
  52.     else if ( strncmp(argv[argn],"-threshold",max(strlen(argv[argn]),2)) == 0 )
  53.         halftone = QT_THRESH;
  54.     else if ( strncmp(argv[argn],"-value",max(strlen(argv[argn]),2)) == 0 )
  55.         {
  56.         float f;
  57.  
  58.         argn++;
  59.         if ( argn == argc || sscanf( argv[argn], "%g", &f ) != 1 ||
  60.          f < 0.0 || f > 1.0 )
  61.         pm_usage( usage );
  62.         threshval = f * SCALE;
  63.         }
  64.     else
  65.         pm_usage( usage );
  66.     argn++;
  67.     }
  68.  
  69.     if ( argn == argc )
  70.     pm_usage( usage );
  71.     if ( sscanf( argv[argn], "%d", &n ) != 1 )
  72.     pm_usage( usage );
  73.     if ( n < 2 )
  74.     pm_error( "N must be greater than 1", 0,0,0,0,0 );
  75.     argn++;
  76.  
  77.     if ( argn == argc )
  78.     ifd = stdin;
  79.     else
  80.     {
  81.     ifd = pm_openr( argv[argn] );
  82.     argn++;
  83.     }
  84.  
  85.     if ( argn != argc )
  86.     pm_usage( usage );
  87.  
  88.     pbm_readpbminit( ifd, &cols, &rows, &format );
  89.     bitslice = pbm_allocarray( cols, n );
  90.  
  91.     newrows = rows / n;
  92.     newcols = cols / n;
  93.     pbm_writepbminit( stdout, newcols, newrows );
  94.     newbitrow = pbm_allocrow( newcols );
  95.  
  96.     if ( halftone == QT_FS )
  97.     {
  98.     /* Initialize Floyd-Steinberg. */
  99.     thiserr = (long *) malloc( ( newcols + 2 ) * sizeof(long) );
  100.     nexterr = (long *) malloc( ( newcols + 2 ) * sizeof(long) );
  101.     if ( thiserr == 0 || nexterr == 0 )
  102.         pm_error( "out of memory", 0,0,0,0,0 );
  103.  
  104.     srandom( (int) time( 0 ) );
  105.     for ( col = 0; col < newcols + 2; col++ )
  106.         thiserr[col] = ( random( ) % SCALE - HALFSCALE ) / 4;
  107.         /* (random errors in [-SCALE/8 .. SCALE/8]) */
  108.     }
  109.     direction = 1;
  110.  
  111.     for ( row = 0; row < newrows; row++ )
  112.     {
  113.     for ( subrow = 0; subrow < n; subrow++ )
  114.         pbm_readpbmrow( ifd, bitslice[subrow], cols, format );
  115.  
  116.     if ( halftone == QT_FS )
  117.         for ( col = 0; col < newcols + 2; col++ )
  118.         nexterr[col] = 0;
  119.     if ( direction )
  120.         {
  121.         col = 0;
  122.         limitcol = newcols;
  123.         nbP = newbitrow;
  124.         }
  125.     else
  126.         {
  127.         col = newcols - 1;
  128.         limitcol = -1;
  129.         nbP = &(newbitrow[col]);
  130.         }
  131.  
  132.     do
  133.         {
  134.         sum = 0;
  135.         count = 0;
  136.         for ( subrow = 0; subrow < n; subrow++ )
  137.         for ( subcol = 0; subcol < n; subcol++ )
  138.             if ( row * n + subrow < rows && col * n + subcol < cols )
  139.             {
  140.             count += 1;
  141.             if ( bitslice[subrow][col * n + subcol] == PBM_WHITE )
  142.                 sum += 1;
  143.             }
  144.         sum = ( sum * SCALE ) / count;
  145.  
  146.         if ( halftone == QT_FS )
  147.         sum += thiserr[col + 1];
  148.  
  149.         if ( sum >= threshval )
  150.         {
  151.         *nbP = PBM_WHITE;
  152.         if ( halftone == QT_FS )
  153.             sum = sum - threshval - HALFSCALE;
  154.         }
  155.         else
  156.         *nbP = PBM_BLACK;
  157.  
  158.         if ( halftone == QT_FS )
  159.         {
  160.         if ( direction )
  161.             {
  162.             thiserr[col + 2] += ( sum * 7 ) / 16;
  163.             nexterr[col    ] += ( sum * 3 ) / 16;
  164.             nexterr[col + 1] += ( sum * 5 ) / 16;
  165.             nexterr[col + 2] += ( sum     ) / 16;
  166.             }
  167.         else
  168.             {
  169.             thiserr[col    ] += ( sum * 7 ) / 16;
  170.             nexterr[col + 2] += ( sum * 3 ) / 16;
  171.             nexterr[col + 1] += ( sum * 5 ) / 16;
  172.             nexterr[col    ] += ( sum     ) / 16;
  173.             }
  174.         }
  175.         if ( direction )
  176.         {
  177.         col++;
  178.         nbP++;
  179.         }
  180.         else
  181.         {
  182.         col--;
  183.         nbP--;
  184.         }
  185.         }
  186.     while ( col != limitcol );
  187.  
  188.     pbm_writepbmrow( stdout, newbitrow, newcols );
  189.  
  190.     if ( halftone == QT_FS )
  191.         {
  192.         temperr = thiserr;
  193.         thiserr = nexterr;
  194.         nexterr = temperr;
  195.         direction = ! direction;
  196.         }
  197.     }
  198.  
  199.     pm_close( ifd );
  200.  
  201.     exit( 0 );
  202.     }
  203.