conv.c

   1 #include "image.h"
   2
   3 /*
   4   General convolution for 2d decoupled filters
   5   end effects are acounted for by increasing
   6   scaling the result with the sum of used coefficients
   7
   8   coeff: (float array) coefficients for filter
   9     len: length of filter.. number of coefficients
  10            note that this has to be an odd number
  11            (since the filter is even);
  12 */
  13
  14 void
  15 i_conv(i_img *im,float *coeff,int len) {
  16   int i,l,c,ch,center;
  17   float pc;
  18   i_color rcolor;
  19   float res[11];
  20   i_img timg;
  21
  22   mm_log((1,"i_conv(im %p, coeff %p, len %d)\n",im,coeff,len));
  23
  24   i_img_empty_ch(&timg,im->xsize,im->ysize,im->channels);
  25
  26   center=(len-1)/2;
  27
  28
  29   for(l=0;l<im->ysize;l++) {
  30     for(i=0;i<im->xsize;i++) {
  31       pc=0.0;
  32       for(ch=0;ch<im->channels;ch++) res[ch]=0;
  33       for(c=0;c<len;c++)
  34         if (i_gpix(im,i+c-center,l,&rcolor)!=-1) {
  35           for(ch=0;ch<im->channels;ch++)
  36             res[ch]+=(float)(rcolor.channel[ch])*coeff[c];
  37           pc+=coeff[c];
  38         }
  39       for(ch=0;ch<im->channels;ch++) {
  40         double temp = res[ch]/pc;
  41         rcolor.channel[ch] =
  42           temp < 0 ? 0 : temp > 255 ? 255 : (unsigned char)temp;
  43       }
  44       i_ppix(&timg,i,l,&rcolor);
  45     }
  46   }
  47
  48   for(l=0;l<im->xsize;l++)
  49     {
  50       for(i=0;i<im->ysize;i++)
  51         {
  52           pc=0.0;
  53           for(ch=0;ch<im->channels;ch++) res[ch]=0;
  54           for(c=0;c<len;c++)
  55             if (i_gpix(&timg,l,i+c-center,&rcolor)!=-1)
  56               {
  57                 for(ch=0;ch<im->channels;ch++)
  58                   res[ch]+=(float)(rcolor.channel[ch])*coeff[c];
  59                 pc+=coeff[c];
  60               }
  61           for(ch=0;ch<im->channels;ch++) {
  62             double temp = res[ch]/pc;
  63             rcolor.channel[ch]=
  64               temp < 0 ? 0 : temp > 255 ? 255 : (unsigned char)temp;
  65           }
  66           i_ppix(im,l,i,&rcolor);
  67         }
  68     }
  69   i_img_exorcise(&timg);
  70 }
  71
  72
  73
  74
  75
  76
  77
  78
  79