implement unsharp mask
[imager.git] / conv.c
CommitLineData
02d1d628
AMH
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
14void
15i_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* 0x%x,coeff 0x%x,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++) res[ch]+=(float)(rcolor.channel[ch])*coeff[c];
36 pc+=coeff[c];
37 }
38 for(ch=0;ch<im->channels;ch++) rcolor.channel[ch]=(unsigned char)(((res[ch]/pc>255.0)?255.0:res[ch]/pc));
39 i_ppix(&timg,i,l,&rcolor);
40 }
41 }
42
43 for(l=0;l<im->xsize;l++)
44 {
45 for(i=0;i<im->ysize;i++)
46 {
47 pc=0.0;
48 for(ch=0;ch<im->channels;ch++) res[ch]=0;
49 for(c=0;c<len;c++)
50 if (i_gpix(&timg,l,i+c-center,&rcolor)!=-1)
51 {
52 for(ch=0;ch<im->channels;ch++) res[ch]+=(float)(rcolor.channel[ch])*coeff[c];
53 pc+=coeff[c];
54 }
55 for(ch=0;ch<im->channels;ch++) rcolor.channel[ch]=(unsigned char)(((res[ch]/(float)(pc)>255.0)?255.0:res[ch]/(float)(pc)));
56 i_ppix(im,l,i,&rcolor);
57 }
58 }
59}
60
61
62
63
64
65
66
67
68