小波變換C語言實(shí)現(xiàn)代碼

1、#include<stdio.h>#include<math.h>#defineheight256#definewidth256voidtwoDDWT(doubleorg口width,doubledwt口width,intNumofBand);voidoneDDWT(intflag,doubleorgwidth,doubledwtwidth);voidtwoDIDWT(doubledwt口width,doublerec口width,intNumofBand);voidoneDIDWT(intflag,doubledwt口width,doublerec口width);in

2、tmain()inti,j;doubleorgheightwidth,dwtheightwidth,recheightwidth;FILE*lena;FILE*lena3;FILE*lena4;lena=fopen("c:lena.raw","rb");lena3=fopen("c:lena3.raw","wb");lena4=fopen("c:lena4.raw","wb");for(i=0;i<height;i+)for(j=0;j<width;j+)orgi

3、j=fgetc(lena);twoDDWT(org,dwt,4);for(i=0;i<height/2;i+)for(j=0;j<width/2;j+)orgij=dwtij;twoDDWT(org,dwt,7);for(i=0;i<height;i+)for(j=0;j<width;j+)fputc(dwtij,lena3);twoDIDWT(dwt,rec,7);for(i=0;i<height/2;i+)(for(j=0;j<width/2;j+)(dwtij=recij;)twoDIDWT(dwt,rec,4);for(i=0;i<height

4、;i+)(for(j=0;j<width;j+)(fputc(recij,lena4);)return0;)voidtwoDDWT(doubleorgwidth,doubledwtwidth,intNumofBand)(if(NumofBand=4)(oneDDWT(0,org,dwt);oneDDWT(1,org,dwt);)if(NumofBand=7)(oneDDWT(2,org,dwt);oneDDWT(3,org,dwt);)voidoneDDWT(intflag,doubleorgwidth,doubledwtwidth)(inti,j;doubletempheightwid

5、th,temp1heightwidth;doublef_LPF口=-0.125,0.25,0.75,0.25,-0.125;doublef_HPF口=-0.5,1,-0.5;if(flag=0)/verticaldirectionfor(i=0;i<height;i+)for(j=0;j<width;j+)if(j=0)tempij=f_LPF0*orgij+2+f_LPF1*orgij+1+f_LPF2*orgij+f_LPF3*orgij+1+f_LPF4*orgij+2;elseif(j=1)tempij=f_LPF0*orgij+f_LPF1*orgij-1+f_LPF2*

6、orgij+f_LPF3*orgij+1+f_LPF4*orgij+2;elseif(j>1&&j<height-2)tempij=f_LPF0*orgij-2+f_LPF1*orgij-1+f_LPF2*orgij+f_LPF3*orgij+1+f_LPF4*orgij+2;elseif(j=height-2)tempij=f_LPF0*orgij-2+f_LPF1*orgij-1+f_LPF2*orgij+f_LPF3*orgij+1+f_LPF4*orgij;elseif(j=height-1)tempij=f_LPF0*orgij-2+f_LPF1*orgi

7、j-1+f_LPF2*orgij+f_LPF3*orgij-1+f_LPF4*orgij-2;)for(i=0;i<height;i+)(for(j=0;j<width/2;j+)(dwtij=tempij*2+1;)for(i=0;i<height;i+)(for(j=0;j<width;j+)(if(j=0)tempij=f_HPF0*orgij+1+f_HPF1*orgij+f_HPF2*orgij+1;elseif(j>0&&j<height-1)tempij=f_HPF0*orgij-1+f_HPF1*orgij+f_HPF2*or

8、gij+1;elseif(j=height-1)tempij=f_HPF0*orgij-1+f_HPF1*orgij+f_HPF2*orgij-1;)for(i=0;i<height;i+)(for(j=0;j<width/2;j+)(dwtij+width/2=tempij*2;)if(flag=1)/horizontaldirection(for(i=0;i<height;i+)(for(j=0;j<width/2;j+)(if(i=0)tempij=f_LPF0*dwti+2j+f_LPF1*dwti+1j+f_LPF2*dwtij+f_LPF3*dwti+1j+

9、f_LPF4*dwti+2j;elseif(i=1)tempij=f_LPF0*dwtij+f_LPF1*dwti-1j+f_LPF2*dwtij+f_LPF3*dwti+1j+f_LPF4*dwti+2j;elseif(i>1&&i<width-2)tempij=f_LPF0*dwti-2j+f_LPF1*dwti-1j+f_LPF2*dwtij+f_LPF3*dwti+1j+f_LPF4*dwti+2j;elseif(i=width-2)tempij=f_LPF0*dwti-2j+f_LPF1*dwti-1j+f_LPF2*dwtij+f_LPF3*dwti+1

10、j+f_LPF4*dwtij;elseif(i=width-1)tempij=f_LPF0*dwti-2j+f_LPF1*dwti-1j+f_LPF2*dwtij+f_LPF3*dwti-1j+f_LPF4*dwti-2j;)for(i=0;i<height;i+)(for(j=0;j<width;j+)(temp1ij=tempi*2+1j;)for(i=0;i<height;i+)(for(j=0;j<width/2;j+)(if(i=0)tempij=f_HPF0*dwti+1j+f_HPF1*dwtij+f_HPF2*dwti+1j;elseif(i>0&

11、amp;&i<width-1)tempij=f_HPF0*dwti-1j+f_HPF1*dwtij+f_HPF2*dwti+1j;elseif(i=width-1)tempij=f_HPF0*dwti-1j+f_HPF1*dwtij+f_HPF2*dwti-1j;)for(i=0;i<height;i+)(for(j=0;j<width/2;j+)(temp1ij+width/2=tempi*2j;)for(i=0;i<height;i+)(for(j=width/2;j<width;j+)(if(i=0)tempij=f_LPF0*dwti+2j+f_L

12、PF1*dwti+1j+f_LPF2*dwtij+f_LPF3*dwti+1j+f_LPF4*dwti+2j;elseif(i=1)tempij=f_LPF0*dwtij+f_LPF1*dwti-1j+f_LPF2*dwtij+f_LPF3*dwti+1j+f_LPF4*dwti+2j;elseif(i>1&&i<height-2)tempij=f_LPF0*dwti-2j+f_LPF1*dwti-1j+f_LPF2*dwtij+f_LPF3* dwti+1j+f_LPF4*dwti+2j;elseif(i=height-2)tempij=f_LPF0*dwti-2

13、j+f_LPF1*dwti-1j+f_LPF2*dwtij+f_LPF3*dwti+1j+f_LPF4*dwtij;elseif(i=height-1)tempij=f_LPF0*dwti-2j+f_LPF1*dwti-1j+f_LPF2*dwtij+f_LPF3* dwti-1j+f_LPF4*dwti-2j;)for(i=0;i<height;i+)(for(j=width/2;j<width;j+)(temp1i+height/2j-width/2=tempi*2+1j;for(i=0;i<height;i+)(for(j=width/2;j<width;j+)(

14、if(i=0)tempij=f_HPF0*dwti+1j+f_HPF1*dwtij+f_HPF2*dwti+1j;elseif(i>0&&i<height-1)tempij=f_HPF0*dwti-1j+f_HPF1*dwtij+f_HPF2*dwti+1j;elseif(i=height-1)tempij=f_HPF0*dwti-1j+f_HPF1*dwtij+f_HPF2*dwti-1j;)for(i=0;i<height;i+)(for(j=width/2;j<width;j+)(temp1i+height/2j=tempi*2j;)for(i=0

15、;i<height;i+)(for(j=0;j<width;j+)(dwtij=temp1ij;)if(flag=2)/verticaldirection(for(i=0;i<height/2;i+)(for(j=0;j<width;j+)(if(j=0)tempij=f_LPF0*orgij+2+f_LPF1*orgij+1+f_LPF2*orgij+f_LPF3*orgij+1+f_LPF4*orgij+2;elseif(j=1)tempij=f_LPF0*orgij+f_LPF1*orgij-1+f_LPF2*orgij+f_LPF3*orgij+1+f_LPF4

16、*orgij+2;elseif(j>1&&j<height/2-2)tempij=f_LPF0*orgij-2+f_LPF1*orgij-1+f_LPF2*orgij+f_LPF3* orgij+1+f_LPF4*orgij+2;elseif(j=height/2-2)tempij=f_LPF0*orgij-2+f_LPF1*orgij-1+f_LPF2*orgij+f_LPF3* orgij+1+f_LPF4*orgij;elseif(j=height/2-1)tempij=f_LPF0*orgij-2+f_LPF1*orgij-1+f_LPF2*orgij+f_

17、LPF3* orgij-1+f_LPF4*orgij-2;)for(i=0;i<height/2;i+)()for(i=0;i<height/2;i+)(for(j=0;j<width/2;j+)(if(j=0)tempij=f_HPF0*orgij+1+f_HPF1*orgij+f_HPF2*orgij+1;elseif(j>0&&j<height-1)tempij=f_HPF0*orgij-1+f_HPF1*orgij+f_HPF2*orgij+1;elseif(j=height-1)tempij=f_HPF0*orgij-1+f_HPF1*o

18、rgij+f_HPF2*orgij-1;)for(i=0;i<height/2;i+)(for(j=0;j<width/4;j+)(dwtij+width/4=tempij*2;)if(flag=3)/horizontaldirection(for(i=0;i<height/2;i+)(for(j=0;j<width/4;j+)(if(i=0)tempij=f_LPF0*dwti+2j+f_LPF1*dwti+1j+f_LPF2*dwtij+f_LPF3*dwti+1j+f_LPF4*dwti+2j;elseif(i=1)tempij=f_LPF0*dwtij+f_LP

19、F1*dwti-1j+f_LPF2*dwtij+f_LPF3*dwti+1j+f_LPF4*dwti+2j;elseif(i>1&&i<height/2-2)tempij=f_LPF0*dwti-2j+f_LPF1*dwti-1j+f_LPF2*dwtij+f_LPF3* dwti+1j+f_LPF4*dwti+2j;elseif(i=height/2-2)tempij=f_LPF0*dwti-2j+f_LPF1*dwti-1j+f_LPF2*dwtij+f_LPF3* dwti+1j+f_LPF4*dwtij;elseif(i=height/2-1)tempij=

20、f_LPF0*dwti-2j+f_LPF1*dwti-1j+f_LPF2*dwtij+f_LPF3* dwti-1j+f_LPF4*dwti-2j;)for(i=0;i<height/4;i+)(for(j=0;j<width/4;j+)(temp1ij=tempi*2+1j;)for(j=0;j<width/4;j+)(if(i=0)tempij=f_HPF0*dwti+1j+f_HPF1*dwtij+f_HPF2*dwti+1j;elseif(i>0&&i<height/2-1)tempij=f_HPF0*dwti-1j+f_HPF1*dwti

21、j+f_HPF2*dwti+1j;elseif(i=height/2-1)tempij=f_HPF0*dwti-1j+f_HPF1*dwtij+f_HPF2*dwti-1j;)for(i=0;i<height/4;i+)(for(j=0;j<width/4;j+)(temp1i+height/4j=tempi*2j;)for(i=0;i<height/2;i+)(for(j=width/4;j<width/2;j+)(if(i=0)tempij=f_LPF0*dwti+2j+f_LPF1*dwti+1j+f_LPF2*dwtij+f_LPF3*dwti+1j+f_LPF

22、4*dwti+2j;elseif(i=1)tempij=f_LPF0*dwtij+f_LPF1*dwti-1j+f_LPF2*dwtij+f_LPF3*dwti+1j+f_LPF4*dwti+2j;elseif(i>1&&i<height/2-2)tempij=f_LPF0*dwti-2j+f_LPF1*dwti-1j+f_LPF2*dwtij+f_LPF3* dwti+1j+f_LPF4*dwti+2j;elseif(i=height/2-2)tempij=f_LPF0*dwti-2j+f_LPF1*dwti-1j+f_LPF2*dwtij+f_LPF3*dwti

23、+1j+f_LPF4*dwtij;elseif(i=height/2-1)tempij=f_LPF0*dwti-2j+f_LPF1*dwti-1j+f_LPF2*dwtij+f_LPF3*dwti-1j+f_LPF4*dwti-2j;)for(i=0;i<height/4;i+)(for(j=width/4;j<width/2;j+)(temp1ij=tempi*2+1j;)for(i=0;i<height/2;i+)(for(j=width/4;j<width/2;j+)(if(i=0)tempij=f_HPF0*dwti+1j+f_HPF1*dwtij+f_HPF2

24、*dwti+1j;elseif(i>0&&i<width/2-1)tempij=f_HPF0*dwti-1j+f_HPF1*dwtij+f_HPF2*dwti+1j;elseif(i=width/2)tempij=f_HPF0*dwti-1j+f_HPF1*dwtij+f_HPF2*dwti-1j;)for(i=0;i<height/4;i+)(for(j=width/4;j<width/2;j+)(temp1i+height/4j=tempi*2j;)for(i=0;i<height/2;i+)(for(j=0;j<width/2;j+)(

25、dwtij=temp1ij;)voidtwoDIDWT(doubledwt口width,doublerec口width,intNumofBand)(if(NumofBand=7)(oneDIDWT(2,dwt,rec);oneDIDWT(3,dwt,rec);)if(NumofBand=4)(oneDIDWT(1,dwt,rec);oneDIDWT(0,dwt,rec);)voidoneDIDWT(intflag,doubledwtwidth,doublerecwidth)(inti,j;doubletempheightwidth,temp1heightwidth;doublei_LPF=0,

26、5,1,0.5;doublei_HPF=-0.125,-0.25,0.75,-0.25,-0.125;if(flag=2)/horizontaldirectionfor(i=0;i<height/2;i+)for(j=0;j<width/4;j+)if(i%2=1)tempij=dwti/2j;elsetempij=0;for(i=0;i<height/2;i+)for(j=0;j<width/4;j+)if(i=0)recij=i_LPF0*tempi+1j+i_LPF1*tempij+i_LPF2*tempi+1j;elseif(i>0&&i&

27、lt;height/2-1)recij=i_LPF0*tempi-1j+i_LPF1*tempij+i_LPF2*tempi+1j;elseif(i=height/2-1)recij=i_LPF0*tempi-1j+i_LPF1*tempij+i_LPF2*tempi-1j;for(i=height/4;i<height/2;i+)for(j=0;j<width/4;j+)dwti-height/4j=dwtij;for(i=0;i<height/2;i+)for(j=0;j<width/4;j+)if(i%2=0)tempij=dwti/2j;elsetempij=0

28、;for(i=0;i<height/2;i+)for(j=0;j<width/4;j+)if(i=0)temp1ij=i_HPF0*tempi+2j+i_HPF1*tempi+1j+i_HPF2*tempij+i_HPF3*tempi+1j+i_HPF4*tempi+2j;elseif(i=1)temp1ij=i_HPF0*tempij+i_HPF1*tempi-1j+i_HPF2*tempij+i_HPF3*tempi+1j+i_HPF4*tempi+2j;elseif(i>1&&i<height/2-2)temp1ij=i_HPF0*tempi-2j

29、+i_HPF1*tempi-1j+i_HPF2*tempij+i_HPF3*tempi+1j+i_HPF4*tempi+2j;elseif(i=height/2-2)temp1ij=i_HPF0*tempi-2j+i_HPF1*tempi-1j+i_HPF2*tempij+i_HPF3*tempi+1j+i_HPF4*tempij;elseif(i=height/2-1)temp1ij=i_HPF0*tempi-2j+i_HPF1*tempi-1j+i_HPF2*tempij+i_HPF3*tempi-1j+i_HPF4*tempi-2j;for(i=0;i<height/2;i+)fo

30、r(j=0;j<width/4;j+)recij=recij+temp1ij;for(j=width/4;j<width/2;j+)(dwtij-height/4=dwtij;)for(i=0;i<height/2;i+)(for(j=0;j<width/4;j+)(if(i%2=1)tempij=dwti/2j;elsetempij=0;)for(i=0;i<height/2;i+)(for(j=0;j<width/4;j+)(if(i=0)temp1ij=i_LPF0*tempi+1j+i_LPF1*tempij+i_LPF2*tempi+1j;else

31、if(i>0&&i<height/2-1)temp1ij=i_LPF0*tempi-1j+i_LPF1*tempij+i_LPF2*tempi+1j;elseif(i=height/2-1)temp1ij=i_LPF0*tempi-1j+i_LPF1*tempij+i_LPF2*tempi-1j;)for(i=height/4;i<height/2;i+)(for(j=width/4;j<width/2;j+)(dwti-height/4j-width/4=dwtij;)for(i=0;i<height/2;i+)(for(j=0;j<wid

32、th/4;j+)(if(i%2=0)tempij=dwti/2j;elsetempij=0;)for(i=0;i<height/2;i+)(for(j=0;j<width/4;j+)(if(i=0)dwtij=i_HPF0*tempi+2j+i_HPF1*tempi+1j+i_HPF2*tempij+i_HPF3*tempi+1j+i_HPF4*tempi+2j;elseif(i=1)dwtij=i_HPF0*tempij+i_HPF1*tempi-1j+i_HPF2*tempij+i_HPF3*tempi+1j+i_HPF4*tempi+2j;elseif(i>1&

33、&i<height/2-2)dwtij=i_HPF0*tempi-2j+i_HPF1*tempi-1j+i_HPF2*tempij+i_HPF3*tempi+1j+i_HPF4*tempi+2j;elseif(i=height/2-2)dwtij=i_HPF0*tempi-2j+i_HPF1*tempi-1j+i_HPF2*tempij+i_HPF3*tempi+1j+i_HPF4*tempij;elseif(i=height/2-1)dwtij=i_HPF0*tempi-2j+i_HPF1*tempi-1j+i_HPF2*tempij+i_HPF3*tempi-1j+i_HPF

34、4*tempi-2j;)for(i=0;i<height/2;i+)(for(j=0;j<width/4;j+)(recij+width/4=temp1ij+dwtij;)if(flag=3)/verticaldirection(for(i=0;i<height/2;i+)(for(j=0;j<width/2;j+)(if(j%2=1)tempij=recij/2;elsetempij=0;)for(i=0;i<height/2;i+)(for(j=0;j<width/2;j+)(if(j=0)temp1ij=i_LPF0*tempij+1+i_LPF1*t

35、empij+i_LPF2*tempij+1;elseif(j>0&&j<height/2-1)temp1ij=i_LPF0*tempij-1+i_LPF1*tempij+i_LPF2*tempij+1;elseif(j=height/2-1)temp1ij=i_LPF0*tempij-1+i_LPF1*tempij+i_LPF2*tempij-1;)for(i=0;i<height/2;i+)(for(j=width/4;j<width/2;j+)(recij-width/4=recij;)for(j=0;j<width/2;j+)(if(j%2=

36、0)tempij=recij/2;elsetempij=0;for(i=0;i<height/2;i+)(for(j=0;j<width/2;j+)(if(j=0)recij=i_HPF0*tempij+2+i_HPF1*tempij+1+i_HPF2*tempij+i_HPF3*tempij+1+i_HPF4*tempij+2;elseif(j=1)recij=i_HPF0*tempij+i_HPF1*tempij-1+i_HPF2*tempij+i_HPF3*tempij+1+i_HPF4*tempij+2;elseif(j>1&&j<height/

37、2-2)recij=i_HPF0*tempij-2+i_HPF1*tempij-1+i_HPF2*tempij+i_HPF3*tempij+1+i_HPF4*tempij+2;elseif(j=height/2-2)recij=i_HPF0*tempij-2+i_HPF1*tempij-1+i_HPF2*tempij+i_HPF3*tempij+1+i_HPF4*tempij;elseif(j=height/2-1)recij=i_HPF0*tempij-2+i_HPF1*tempij-1+i_HPF2*tempij+i_HPF3*tempij-1+i_HPF4*tempij-2;for(i=

38、0;i<height/2;i+)(for(j=0;j<width/2;j+)(recij=recij+temp1ij;if(flag=1)/horizontaldirection(for(i=0;i<height;i+)(for(j=0;j<width/2;j+)(if(i%2=1)tempij=dwti/2j;elsetempij=0;for(i=0;i<height;i+)(for(j=0;j<width/2;j+)(if(i=0)recij=i_LPF0*tempi+1j+i_LPF1*tempij+i_LPF2*tempi+1j;elseif(i&g

39、t;0&&i<height-1)recij=i_LPF0*tempi-1j+i_LPF1*tempij+i_LPF2*tempi+1j;elseif(i=height-1)recij=i_LPF0*tempi-1j+i_LPF1*tempij+i_LPF2*tempi-1j;for(i=0;i<height/2;i+)for(j=width/2;j<width;j+)dwtij-height/2=dwtij;for(i=0;i<height;i+)for(j=0;j<width/2;j+)if(i%2=0)tempij=dwti/2j;elsete

40、mpij=0;for(i=0;i<height;i+)for(j=0;j<width/2;j+)if(i=0)temp1ij=i_HPF0*tempi+2j+i_HPF1*tempi+1j+i_HPF2*tempij+i_HPF3*tempi+1j+i_HPF4*tempi+2j;elseif(i=1)temp1ij=i_HPF0*tempij+i_HPF1*tempi-1j+i_HPF2*tempij+i_HPF3*tempi+1j+i_HPF4*tempi+2j;elseif(i>1&&i<height-2)temp1ij=i_HPF0*tempi-

41、2j+i_HPF1*tempi-1j+i_HPF2*tempij+i_HPF3*tempi+1j+i_HPF4*tempi+2j;elseif(i=height-2)temp1ij=i_HPF0*tempi-2j+i_HPF1*tempi-1j+i_HPF2*tempij+i_HPF3*tempi+1j+i_HPF4*tempij;elseif(i=height-1)temp1ij=i_HPF0*tempi-2j+i_HPF1*tempi-1j+i_HPF2*tempij+i_HPF3*tempi-1j+i_HPF4*tempi-2j;for(i=0;i<height;i+)for(j=

42、0;j<width/2;j+)recij=recij+temp1ij;for(i=height/2;i<height;i+)(for(j=0;j<width/2;j+)(dwti-width/2j=dwtij;)for(i=0;i<height;i+)(for(j=0;j<width/2;j+)(if(i%2=1)tempij=dwti/2j;elsetempij=0;)for(i=0;i<height;i+)(for(j=0;j<width/2;j+)(if(i=0)temp1ij=i_LPF0*tempi+1j+i_LPF1*tempij+i_LP

43、F2*tempi+1j;elseif(i>0&&i<height-1)temp1ij=i_LPF0*tempi-1j+i_LPF1*tempij+i_LPF2*tempi+1j;elseif(i=height-1)temp1ij=i_LPF0*tempi-1j+i_LPF1*tempij+i_LPF2*tempi-1j;)for(i=height/2;i<height;i+)(for(j=width/2;j<width;j+)(dwti-height/2j-width/2=dwtij;)for(i=0;i<height;i+)(for(j=0;j&

44、lt;width/2;j+)(if(i%2=0)tempij=dwti/2j;elsetempij=0;)for(i=0;i<height;i+)(for(j=0;j<width/2;j+)(if(i=0)dwtij=i_HPF0*tempi+2j+i_HPF1*tempi+1j+i_HPF2*tempij+i_HPF3*tempi+1j+i_HPF4*tempi+2j;elseif(i=1)dwtij=i_HPF0*tempij+i_HPF1*tempi-1j+i_HPF2*tempij+i_HPF3*tempi+1j+i_HPF4*tempi+2j;elseif(i>1&&i<height-2)dwtij=i_HPF0*tempi-2j+i_HPF1*tempi-1j+i_HPF2*tempij+i_