home *** CD-ROM | disk | FTP | other *** search
/ Club Amiga de Montreal - CAM / CAM_CD_1.iso / files / 309.lha / PBM_PLUS / ppm / ppmarith.c < prev    next >
C/C++ Source or Header  |  1980-12-04  |  3KB  |  136 lines
  1. /* ppmarith.c - perform arithmetic on two portable pixmaps
  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. #else    SYSV
  17. #include <strings.h>
  18. #endif    SYSV
  19. #include "ppm.h"
  20.  
  21. #define max(a,b) ((a) > (b) ? (a) : (b))
  22.  
  23. main( argc, argv )
  24. int argc;
  25. char *argv[];
  26.     {
  27.     FILE *ifd1, *ifd2;
  28.     register pixel *pixrow1, *pixrow2, *p1P, *p2P;
  29.     pixval maxval1, maxval2, maxval3;
  30.     int argn, rows1, cols1, format1, rows2, cols2, format2, row, col;
  31.     char function;
  32.     char *usage = "-add|-subtract|-multiply ppmfile1 ppmfile2";
  33.  
  34.     pm_progname = argv[0];
  35.  
  36.     argn = 1;
  37.     function = ' ';
  38.  
  39.     /* Check for flags. */
  40.     if ( argn < argc && argv[argn][0] == '-' )
  41.     {
  42.     if ( strncmp(argv[argn],"-add",max(strlen(argv[argn]),2)) == 0 )
  43.         function = '+';
  44.     else if ( strncmp(argv[argn],"-subtract",max(strlen(argv[argn]),2)) == 0 )
  45.         function = '-';
  46.     else if ( strncmp(argv[argn],"-multiply",max(strlen(argv[argn]),2)) == 0 )
  47.         function = '*';
  48.     else
  49.         pm_usage( usage );
  50.     argn++;
  51.     }
  52.  
  53.     if ( function == ' ' )
  54.     pm_usage( usage );
  55.  
  56.     if ( argn == argc )
  57.     pm_usage( usage );
  58.     ifd1 = pm_openr( argv[argn] );
  59.     argn++;
  60.  
  61.     if ( argn == argc )
  62.     pm_usage( usage );
  63.     ifd2 = pm_openr( argv[argn] );
  64.     argn++;
  65.  
  66.     if ( argn != argc )
  67.     pm_usage( usage );
  68.  
  69.     ppm_readppminit( ifd1, &cols1, &rows1, &maxval1, &format1 );
  70.     pixrow1 = ppm_allocrow( cols1 );
  71.     ppm_readppminit( ifd2, &cols2, &rows2, &maxval2, &format2 );
  72.     if ( cols2 != cols1 || rows2 != rows1 )
  73.     pm_error(
  74.         "the two pixmaps must be the same width and height", 0,0,0,0,0 );
  75.     pixrow2 = ppm_allocrow( cols1 );
  76.  
  77.     maxval3 = max( maxval1, maxval2 );
  78.     ppm_writeppminit( stdout, cols1, rows1, maxval3 );
  79.     for ( row = 0; row < rows1; row++ )
  80.     {
  81.     ppm_readppmrow( ifd1, pixrow1, cols1, maxval1, format1 );
  82.     ppm_readppmrow( ifd2, pixrow2, cols1, maxval2, format2 );
  83.         for ( col = 0, p1P = pixrow1, p2P = pixrow2; col < cols1; col++, p1P++, p2P++ )
  84.         {
  85.         int r1, g1, b1, r2, g2, b2;
  86.  
  87.         if ( maxval1 != maxval3 )
  88.         PPM_CSCALE( *p1P, *p1P, maxval1, maxval3 );
  89.         r1 = PPM_GETR( *p1P );
  90.         g1 = PPM_GETG( *p1P );
  91.         b1 = PPM_GETB( *p1P );
  92.         if ( maxval2 != maxval3 )
  93.         PPM_CSCALE( *p2P, *p2P, maxval2, maxval3 );
  94.         r2 = PPM_GETR( *p2P );
  95.         g2 = PPM_GETG( *p2P );
  96.         b2 = PPM_GETB( *p2P );
  97.         switch ( function )
  98.         {
  99.         case '+':
  100.         r1 += r2;
  101.         g1 += g2;
  102.         b1 += b2;
  103.         break;
  104.  
  105.         case '-':
  106.         r1 -= r2;
  107.         g1 -= g2;
  108.         b1 -= b2;
  109.         break;
  110.  
  111.         case '*':
  112.         r1 = r1 * r2 / maxval3;
  113.         g1 = g1 * g2 / maxval3;
  114.         b1 = b1 * b2 / maxval3;
  115.         break;
  116.  
  117.         default:
  118.         pm_error( "can't happen", 0,0,0,0,0 );
  119.         }
  120.         if ( r1 < 0 ) r1 = 0;
  121.         else if ( r1 > maxval3 ) r1 = maxval3;
  122.         if ( g1 < 0 ) g1 = 0;
  123.         else if ( g1 > maxval3 ) g1 = maxval3;
  124.         if ( b1 < 0 ) b1 = 0;
  125.         else if ( b1 > maxval3 ) b1 = maxval3;
  126.         PPM_ASSIGN( *p1P, r1, g1, b1 );
  127.         }
  128.     ppm_writeppmrow( stdout, pixrow1, cols1, maxval3 );
  129.     }
  130.  
  131.     pm_close( ifd1 );
  132.     pm_close( ifd2 );
  133.  
  134.     exit( 0 );
  135.     }
  136.