小波變換C語(yǔ)言實(shí)現(xiàn)代碼_第1頁(yè)
小波變換C語(yǔ)言實(shí)現(xiàn)代碼_第2頁(yè)
小波變換C語(yǔ)言實(shí)現(xiàn)代碼_第3頁(yè)
小波變換C語(yǔ)言實(shí)現(xiàn)代碼_第4頁(yè)
小波變換C語(yǔ)言實(shí)現(xiàn)代碼_第5頁(yè)
已閱讀5頁(yè)，還剩5頁(yè)未讀， 繼續(xù)免費(fèi)閱讀

下載本文檔

版權(quán)說(shuō)明：本文檔由用戶(hù)提供并上傳，收益歸屬內(nèi)容提供方，若內(nèi)容存在侵權(quán)，請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)

文檔簡(jiǎn)介

#include<stdio.h>#include<math.h>#defineheight256#definewidth256voidtwoDDWT(doubleorg[][width],doubledwt[][width],intNumofBand);voidoneDDWT(intflag,doubleorg[][width],doubledwt[][width]);voidtwoDIDWT(doubledwt[][width],doublerec[][width],intNumofBand);voidoneDIDWT(intflag,doubledwt[][width],doublerec[][width]);intmain()(

inti,j;

doubleorg[height][width],dwt[height][width],rec[height][width];

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++)(org[i][j]=fgetc(lena);)

)

twoDDWT(org,dwt,4);

for(i=0;i<height/2;i++)

(for(j=0;j<width/2;j++)(org[i][j]=dwt[i][j];)

)

twoDDWT(org,dwt,7);

for(i=0;i<height;i++)

(for(j=0;j<width;j++)(fputc(dwt[i][j],lena3);)

)

twoDIDWT(dwt,rec,7);

for(i=0;i<height/2;i++)

(for(j=0;j<width/2;j++)(dwt[i][j]=rec[i][j];)

)

twoDIDWT(dwt,rec,4);

for(i=0;i<height;i++)

(for(j=0;j<width;j++)(fputc(rec[i][j],lena4);)

)

return0;)voidtwoDDWT(doubleorg[][width],doubledwt[][width],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,doubleorg[][width],doubledwt[][width])(

inti,j;

doubletemp[height][width],temp1[height][width];

doublef_LPF[]=(-0.125,0.25,0.75,0.25,-0.125);

doublef_HPF[]=(-0.5,1,-0.5);

if(flag==0)//verticaldirection

(for(i=0;i<height;i++)(for(j=0;j<width;j++)(if(j==0)temp[i][j]=f_LPF[0]*org[i][j+2]+f_LPF[1]*org[i][j+1]+f_LPF[2]*org[i][j]+f_LPF[3]org[i][j+1]+f_LPF[4]*org[i][j+2];elseif(j==1)temp[i][j]=f_LPF[0]*org[i][j]+f_LPF[1]*org[i][j-1]+f_LPF[2]*org[i][j]+f_LPF[3]*org[i][j+1]+f_LPF[4]*org[i][j+2];eseif(j>1&&j<height-2)temp[i][j]=f_LPF[0]*org[i][j-2]+f_LPF[1]*org[i][j-1]+f_LPF[2]*org[i][j]+f_LPF[3]org[i][j+1]+f_LPF[4]*org[i][j+2];elseif(j==height-2)temp[i][j]=f_LPF[0]*org[i][j-2]+f_LPF[1]*org[i][j-1]+f_LPF[2]*org[i][j]+f_LPF[3]org[i][j+1]+f_LPF[4]*org[i][j];elseif(j==height-1)temp[i][j]=f_LPF[0]*org[i][j-2]+f_LPF[1]*org[i][j-1]+f_LPF[2]*org[i][j]+f_LPF[3]org[i][j-1]+f_LPF[4]*org[i][j-2];))for(i=0;i<height;i++)(for(j=0;j<width/2;j++)(dwt[i][j]=temp[i][j*2+1];))for(i=0;i<height;i++)(for(j=0;j<width;j++)(if(j==0)temp[i][j]=f_HPF[0]*org[i][j+1]+f_HPF[1]*org[i][j]+f_HPF[2]*org[i][j+1];elseif(j>0&&j<height-1)temp[i][j]=f_HPF[0]*org[i][j-1]+f_HPF[1]*org[i][j]+f_HPF[2]*org[i][j+1];elseif(j==height-1)temp[i][j]=f_HPF[0]*org[i][j-1]+f_HPF[1]*org[i][j]+f_HPF[2]*org[i][j-1];))for(i=0;i<height;i++)(for(j=0;j<width/2;j++)(dwt[i][j+width/2]=temp[i][j*2];))

)

if(flag==1)//horizontaldirection

(for(i=0;i<height;i++)(for(j=0;j<width/2;j++)(if(i==0)

temp[i][j]=f_LPF[0]*dwt[i+2][j]+f_LPF[1]*dwt[i+1][j]+f_LPF[2]*dwt[i][j]+f_LPF[3]*dwt[i+1][j]+f_LPF[4]*dwt[i+2][j];elseif(i==1)

temp[i][j]=f_LPF[0]*dwt[i][j]+f_LPF[1]*dwt[i-1][j]+f_LPF[2]*dwt[i][j]+f_LPF[3]*dwt[i+1][j]+f_LPF[4]*dwt[i+2][j];elseif(i>1&&i<width-2)temp[i][j]=f_LPF[0]*dwt[i-2][j]+f_LPF[1]*dwt[i-1][j]+f_LPF[2]*dwt[i][j]+f_LPF[3]dwt[i+1][j]+f_LPF[4]*dwt[i+2][j];elseif(i==width-2)temp[i][j]=f_LPF[0]*dwt[i-2][j]+f_LPF[1]*dwt[i-1][j]+f_LPF[2]*dwt[i][j]+f_LPF[3]dwt[i+1][j]+f_LPF[4]*dwt[i][j];elseif(i==width-1)temp[i][j]=f_LPF[0]*dwt[i-2][j]+f_LPF[1]*dwt[i-1][j]+f_LPF[2]*dwt[i][j]+f_LPF[3]dwt[i-1][j]+f_LPF[4]*dwt[i-2][j];)for(i=0;i<height;i++)(for(j=0;j<width;j++)(temp1[i][j]=temp[i*2+1][j];))for(i=0;i<height;i++)(for(j=0;j<width/2;j++)(if(i==0)temp[i][j]=f_HPF[0]*dwt[i+1][j]+f_HPF[1]*dwt[i][j]+f_HPF[2]*dwt[i+1][j];elseif(i>0&&i<width-1)temp[i][j]=f_HPF[0]*dwt[i-1][j]+f_HPF[1]*dwt[i][j]+f_HPF[2]*dwt[i+1][j];elseif(i==width-1)temp[i][j]=f_HPF[0]*dwt[i-1][j]+f_HPF[1]*dwt[i][j]+f_HPF[2]*dwt[i-1][j];))for(i=0;i<height;i++)(for(j=0;j<width/2;j++)(temp1[i][j+width/2]=temp[i*2][j];))for(i=0;i<height;i++)(for(j=width/2;j<width;j++)(if(i==0)

temp[i][j]=f_LPF[0]*dwt[i+2][j]+f_LPF[1]*dwt[i+1][j]+f_LPF[2]*dwt[i][j]+f_LPF[3]*dwt[i+1][j]+f_LPF[4]*dwt[i+2][j];elseif(i==1)

temp[i][j]=f_LPF[0]*dwt[i][j]+f_LPF[1]*dwt[i-1][j]+f_LPF[2]*dwt[i][j]+f_LPF[3]*dwt[i+1][j]+f_LPF[4]*dwt[i+2][j];elseif(i>1&&i<height-2)temp[i][j]=f_LPF[0]*dwt[i-2][j]+f_LPF[1]*dwt[i-1][j]+f_LPF[2]*dwt[i][j]+f_LPF[3]dwt[i+1][j]+f_LPF[4]*dwt[i+2][j];elseif(i==height-2)temp[i][j]=f_LPF[0]*dwt[i-2][j]+f_LPF[1]*dwt[i-1][j]+f_LPF[2]*dwt[i][j]+f_LPF[3]dwt[i+1][j]+f_LPF[4]*dwt[i][j];elseif(i==height-1)temp[i][j]=f_LPF[0]*dwt[i-2][j]+f_LPF[1]*dwt[i-1][j]+f_LPF[2]*dwt[i][j]+f_LPF[3]dwt[i-1][j]+f_LPF[4]*dwt[i-2][j];))for(i=0;i<height;i++)(for(j=width/2;j<width;j++)(temp1[i+height/2][j-width/2]=temp[i*2+1][j];for(i=0;i<height;i++)(for(j=width/2;j<width;j++)(if(i==0)temp[i][j]=f_HPF[0]*dwt[i+1][j]+f_HPF[1]*dwt[i][j]+f_HPF[2]*dwt[i+1][j];elseif(i>0&&i<height-1)temp[i][j]=f_HPF[0]*dwt[i-1][j]+f_HPF[1]*dwt[i][j]+f_HPF[2]*dwt[i+1][j];elseif(i==height-1)temp[i][j]=f_HPF[0]*dwt[i-1][j]+f_HPF[1]*dwt[i][j]+f_HPF[2]*dwt[i-1][j];))for(i=0;i<height;i++)(for(j=width/2;j<width;j++)(temp1[i+height/2][j]=temp[i*2][j];))for(i=0;i<height;i++)(for(j=0;j<width;j++)(dwt[i][j]=temp1[i][j];))

)

if(flag==2)//verticaldirection

(for(i=0;i<height/2;i++)(for(j=0;j<width;j++)(if(j==0)temp[i][j]=f_LPF[0]*org[i][j+2]+f_LPF[1]*org[i][j+1]+f_LPF[2]*org[i][j]+f_LPF[3]org[i][j+1]+f_LPF[4]*org[i][j+2];elseif(j==1)

temp[i][j]=f_LPF[0]*org[i][j]+f_LPF[1]*org[i][j-1]+f_LPF[2]*org[i][j]+f_LPF[3]*org[i][j+1]+f_LPF[4]*org[i][j+2];eseif(j>1&&j<height/2-2)temp[i][j]=f_LPF[0]*org[i][j-2]+f_LPF[1]*org[i][j-1]+f_LPF[2]*org[i][j]+f_LPF[3]org[i][j+1]+f_LPF[4]*org[i][j+2];elseif(j==height/2-2)temp[i][j]=f_LPF[0]*org[i][j-2]+f_LPF[1]*org[i][j-1]+f_LPF[2]*org[i][j]+f_LPF[3]org[i][j+1]+f_LPF[4]*org[i][j];elseif(j==height/2-1)temp[i][j]=f_LPF[0]*org[i][j-2]+f_LPF[1]*org[i][j-1]+f_LPF[2]*org[i][j]+f_LPF[3]org[i][j-1]+f_LPF[4]*org[i][j-2];))for(i=0;i<height/2;i++)(for(j=0;j<width/4;j++)

(dwt[i][j]=temp[i][j*2+1];))for(i=0;i<height/2;i++)(for(j=0;j<width/2;j++)(if(j==0)temp[i][j]=f_HPF[0]*org[i][j+1]+f_HPF[1]*org[i][j]+f_HPF[2]*org[i][j+1];elseif(j>0&&j<height-1)temp[i][j]=f_HPF[0]*org[i][j-1]+f_HPF[1]*org[i][j]+f_HPF[2]*org[i][j+1];elseif(j==height-1)temp[i][j]=f_HPF[0]*org[i][j-1]+f_HPF[1]*org[i][j]+f_HPF[2]*org[i][j-1];))for(i=0;i<height/2;i++)(for(j=0;j<width/4;j++)(dwt[i][j+width/4]=temp[i][j*2];))

)

if(flag==3)//horizontaldirection

(for(i=0;i<height/2;i++)(for(j=0;j<width/4;j++)(if(i==0)

temp[i][j]=f_LPF[0]*dwt[i+2][j]+f_LPF[1]*dwt[i+1][j]+f_LPF[2]*dwt[i][j]+f_LPF[3]*dwt[i+1][j]+f_LPF[4]*dwt[i+2][j];elseif(i==1)

temp[i][j]=f_LPF[0]*dwt[i][j]+f_LPF[1]*dwt[i-1][j]+f_LPF[2]*dwt[i][j]+f_LPF[3]*dwt[i+1][j]+f_LPF[4]*dwt[i+2][j];elseif(i>1&&i<height/2-2)temp[i][j]=f_LPF[0]*dwt[i-2][j]+f_LPF[1]*dwt[i-1][j]+f_LPF[2]*dwt[i][j]+f_LPF[3]dwt[i+1][j]+f_LPF[4]*dwt[i+2][j];elseif(i==height/2-2)temp[i][j]=f_LPF[0]*dwt[i-2][j]+f_LPF[1]*dwt[i-1][j]+f_LPF[2]*dwt[i][j]+f_LPF[3]dwt[i+1][j]+f_LPF[4]*dwt[i][j];elseif(i==height/2-1)temp[i][j]=f_LPF[0]*dwt[i-2][j]+f_LPF[1]*dwt[i-1][j]+f_LPF[2]*dwt[i][j]+f_LPF[3]dwt[i-1][j]+f_LPF[4]*dwt[i-2][j];))for(i=0;i<height/4;i++)(for(j=0;j<width/4;j++)(temp1[i][j]=temp[i*2+1][j];))

for(j=0;j<width/4;j++)(if(i==0)temp[i][j]=f_HPF[0]*dwt[i+1][j]+f_HPF[1]*dwt[i][j]+f_HPF[2]*dwt[i+1][j];elseif(i>0&&i<height/2-1)temp[i][j]=f_HPF[0]*dwt[i-1][j]+f_HPF[1]*dwt[i][j]+f_HPF[2]*dwt[i+1][j];elseif(i==height/2-1)temp[i][j]=f_HPF[0]*dwt[i-1][j]+f_HPF[1]*dwt[i][j]+f_HPF[2]*dwt[i-1][j];))for(i=0;i<height/4;i++)(for(j=0;j<width/4;j++)(temp1[i+height/4][j]=temp[i*2][j];))for(i=0;i<height/2;i++)(for(j=width/4;j<width/2;j++)(if(i==0)

temp[i][j]=f_LPF[0]*dwt[i+2][j]+f_LPF[1]*dwt[i+1][j]+f_LPF[2]*dwt[i][j]+f_LPF[3]*dwt[i+1][j]+f_LPF[4]*dwt[i+2][j];elseif(i==1)

temp[i][j]=f_LPF[0]*dwt[i][j]+f_LPF[1]*dwt[i-1][j]+f_LPF[2]*dwt[i][j]+f_LPF[3]*dwt[i+1][j]+f_LPF[4]*dwt[i+2][j];elseif(i>1&&i<height/2-2)temp[i][j]=f_LPF[0]*dwt[i-2][j]+f_LPF[1]*dwt[i-1][j]+f_LPF[2]*dwt[i][j]+f_LPF[3]dwt[i+1][j]+f_LPF[4]*dwt[i+2][j];elseif(i==height/2-2)temp[i][j]=f_LPF[0]*dwt[i-2][j]+f_LPF[1]*dwt[i-1][j]+f_LPF[2]*dwt[i][j]+f_LPF[3]dwt[i+1][j]+f_LPF[4]*dwt[i][j];elseif(i==height/2-1)temp[i][j]=f_LPF[0]*dwt[i-2][j]+f_LPF[1]*dwt[i-1][j]+f_LPF[2]*dwt[i][j]+f_LPF[3]dwt[i-1][j]+f_LPF[4]*dwt[i-2][j];))for(i=0;i<height/4;i++)(for(j=width/4;j<width/2;j++)(temp1[i][j]=temp[i*2+1][j];))for(i=0;i<height/2;i++)(for(j=width/4;j<width/2;j++)(if(i==0)temp[i][j]=f_HPF[0]*dwt[i+1][j]+f_HPF[1]*dwt[i][j]+f_HPF[2]*dwt[i+1][j];elseif(i>0&&i<width/2-1)temp[i][j]=f_HPF[0]*dwt[i-1][j]+f_HPF[1]*dwt[i][j]+f_HPF[2]*dwt[i+1][j];elseif(i==width/2)temp[i][j]=f_HPF[0]*dwt[i-1][j]+f_HPF[1]*dwt[i][j]+f_HPF[2]*dwt[i-1][j];)for(i=0;i<height/4;i++)(for(j=width/4;j<width/2;j++)(temp1[i+height/4][j]=temp[i*2][j];))for(i=0;i<height/2;i++)(for(j=0;j<width/2;j++)(dwt[i][j]=temp1[i][j];))

))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,doubledwt[][width],doublerec[][width])(

inti,j;

doubletemp[height][width],temp1[height][width];

doublei_LPF[]=(0.5,1,0.5);

doublei_HPF[]=(-0.125,-0.25,0.75,-0.25,-0.125);

if(flag==2)//horizontaldirection

(for(i=0;i<height/2;i++)(for(j=0;j<width/4;j++)(if(i%2==1)temp[i][j]=dwt[i/2][j];elsetemp[i][j]=0;))for(i=0;i<height/2;i++)(for(j=0;j<width/4;j++)(

if(i==0)

rec[i][j]=i_LPF[0]*temp[i+1][j]+i_LPF[1]*temp[i][j]+i_LPF[2]*temp[i+1][j];elseif(i>0&&i<height/2-1)rec[i][j]=i_LPF[0]*temp[i-1][j]+i_LPF[1]*temp[i][j]+i_LPF[2]*temp[i+1][j];elseif(i==height/2-1)rec[i][j]=i_LPF[0]*temp[i-1][j]+i_LPF[1]*temp[i][j]+i_LPF[2]*temp[i-1][j];))for(i=height/4;i<height/2;i++)(for(j=0;j<width/4;j++)(dwt[i-height/4][j]=dwt[i][j];))for(i=0;i<height/2;i++)(for(j=0;j<width/4;j++)(if(i%2==0)temp[i][j]=dwt[i/2][j];elsetemp[i][j]=0;))for(i=0;i<height/2;i++)(for(j=0;j<width/4;j++)(if(i==0)

temp1[i][j]=i_HPF[0]*temp[i+2][j]+i_HPF[1]*temp[i+1][j]+i_HPF[2]*temp[i][j]+i_HPF[3]*temp[i+1][j]+i_HPF[4]*temp[i+2][j];elseif(i==1)

temp1[i][j]=i_HPF[0]*temp[i][j]+i_HPF[1]*temp[i-1][j]+i_HPF[2]*temp[i][j]+i_HPF[3]*temp[i+1][j]+i_HPF[4]*temp[i+2][j];elseif(i>1&&i<height/2-2)

temp1[i][j]=i_HPF[0]*temp[i-2][j]+i_HPF[1]*temp[i-1][j]+i_HPF[2]*temp[i][j]+i_HPF[3]*temp[i+1][j]+i_HPF[4]*temp[i+2][j];elseif(i==height/2-2)

temp1[i][j]=i_HPF[0]*temp[i-2][j]+i_HPF[1]*temp[i-1][j]+i_HPF[2]*temp[i][j]+i_HPF[3]*temp[i+1][j]+i_HPF[4]*temp[i][j];elseif(i==height/2-1)

temp1[i][j]=i_HPF[0]*temp[i-2][j]+i_HPF[1]*temp[i-1][j]+i_HPF[2]*temp[i][j]+i_HPF[3]*temp[i-1][j]+i_HPF[4]*temp[i-2][j];))for(i=0;i<height/2;i++)(for(j=0;j<width/4;j++)(rec[i][j]=rec[i][j]+temp1[i][j];))

for(j=width/4;j<width/2;j++)(dwt[i][j-height/4]=dwt[i][j];))for(i=0;i<height/2;i++)(for(j=0;j<width/4;j++)(if(i%2==1)temp[i][j]=dwt[i/2][j];elsetemp[i][j]=0;))for(i=0;i<height/2;i++)(for(j=0;j<width/4;j++)(if(i==0)temp1[i][j]=i_LPF[0]*temp[i+1][j]+i_LPF[1]*temp[i][j]+i_LPF[2]*temp[i+1][j];elseif(i>0&&i<height/2-1)temp1[i][j]=i_LPF[0]*temp[i-1][j]+i_LPF[1]*temp[i][j]+i_LPF[2]*temp[i+1][j];elseif(i==height/2-1)temp1[i][j]=i_LPF[0]*temp[i-1][j]+i_LPF[1]*temp[i][j]+i_LPF[2]*temp[i-1][j];))for(i=height/4;i<height/2;i++)(for(j=width/4;j<width/2;j++)(dwt[i-height/4][j-width/4]=dwt[i][j];))for(i=0;i<height/2;i++)(for(j=0;j<width/4;j++)(if(i%2==0)temp[i][j]=dwt[i/2][j];elsetemp[i][j]=0;))for(i=0;i<height/2;i++)(for(j=0;j<width/4;j++)(if(i==0)

dwt[i][j]=i_HPF[0]*temp[i+2][j]+i_HPF[1]*temp[i+1][j]+i_HPF[2]*temp[i][j]+i_HPF[3]*temp[i+1][j]+i_HPF[4]*temp[i+2][j];elseif(i==1)

dwt[i][j]=i_HPF[0]*temp[i][j]+i_HPF[1]*temp[i-1][j]+i_HPF[2]*temp[i][j]+i_HPF[3]*temp[i+1][j]+i_HPF[4]*temp[i+2][j];elseif(i>1&&i<height/2-2)

dwt[i][j]=i_HPF[0]*temp[i-2][j]+i_HPF[1]*temp[i-1][j]+i_HPF[2]*temp[i][j]+i_HPF[3]*temp[i+1][j]+i_HPF[4]*temp[i+2][j];elseif(i==height/2-2)

dwt[i][j]=i_HPF[0]*temp[i-2][j]+i_HPF[1]*temp[i-1][j]+i_HPF[2]*temp[i][j]+i_HPF[3]*temp[i+1][j]+i_HPF[4]*temp[i][j];elseif(i==height/2-1)

dwt[i][j]=i_HPF[0]*temp[i-2][j]+i_HPF[1]*temp[i-1][j]+i_HPF[2]*temp[i][j]+i_HPF[3]*temp[i-1][j]+i_HPF[4]*temp[i-2][j];))for(i=0;i<height/2;i++)(for(j=0;j<width/4;j++)(rec[i][j+width/4]=temp1[i][j]+dwt[i][j];))

)

if(flag==3)//verticaldirection

(for(i=0;i<height/2;i++)(for(j=0;j<width/2;j++)(if(j%2==1)temp[i][j]=rec[i][j/2];elsetemp[i][j]=0;))for(i=0;i<height/2;i++)(for(j=0;j<width/2;j++)(if(j==0)temp1[i][j]=i_LPF[0]*temp[i][j+1]+i_LPF[1]*temp[i][j]+i_LPF[2]*temp[i][j+1];elseif(j>0&&j<height/2-1)temp1[i][j]=i_LPF[0]*temp[i][j-1]+i_LPF[1]*temp[i][j]+i_LPF[2]*temp[i][j+1];elseif(j==height/2-1)temp1[i][j]=i_LPF[0]*temp[i][j-1]+i_LPF[1]*temp[i][j]+i_LPF[2]*temp[i][j-1];))for(i=0;i<height/2;i++)(for(j=width/4;j<width/2;j++)(rec[i][j-width/4]=rec[i][j];))

for(j=0;j<width/2;j++)(if(j%2==0)temp[i][j]=rec[i][j/2];elsetemp[i][j]=0;))for(i=0;i<height/2;i++)(for(j=0;j<width/2;j++)(if(j==0)

rec[i][j]=i_HPF[0]*temp[i][j+2]+i_HPF[1]*temp[i][j+1]+i_HPF[2]*temp[i][j]+i_HPF[3]*temp[i][j+1]+i_HPF[4]*temp[i][j+2];elseif(j==1)rec[i][j]=i_HPF[0]*temp[i][j]+i_HPF[1]*temp[i][j-1]+i_HPF[2]*temp[i][j]+i_HPF[3]*temp[i][j+1]+i_HPF[4]*temp[而+2];elseif(j>1&&j<height/2-2)

rec[i][j]=i_HPF[0]*temp[i][j-2]+i_HPF[1]*temp[i][j-1]+i_HPF[2]*temp[i][j]+i_HPF[3]*temp[i][j+1]+i_HPF[4]*temp[i][j+2];elseif(j==height/2-2)

rec[i][j]=i_HPF[0]*temp[i][j-2]+i_HPF[1]*temp[i][j-1]+i_HPF[2]*temp[i][j]+i_HPF[3]*temp[i][j+1]+i_HPF[4]*temp[i][j];elseif(j==height/2-1)

rec[i][j]=i_HPF[0]*temp[i][j-2]+i_HPF[1]*temp[i][j-1]+i_HPF[2]*temp[i][j]+i_HPF[3]*temp[i][j-1]+i_HPF[4]*temp[i][j-2];))for(i=0;i<height/2;i++)(for(j=0;j<width/2;j++)(rec[i][j]=rec[i][j]+temp1[i][j];))

)

if(flag==1)//horizontaldirection

(for(i=0;i<height;i++)(for(j=0;j<width/2;j++)(if(i%2==1)temp[i][j]=dwt[i/2][j];elsetemp[i][j]=0;))for(i=0;i<height;i++)(for(j=0;j<width/2;j++)(if(i==0)

rec[i][j]=i_LPF[0]*temp[i+1][j]+i_LPF[1]*temp[i][j]+i_LPF[2]*temp[i+1][j];elseif(i>0&&i<height-1)rec[i][j]=i_LPF[0]*temp[i-1][j]+i_LPF[1]*temp[i][j]+i_LPF[2]*temp[i+1][j];elseif(i==height-1)rec[i][j]=i_LPF[0]*temp[i-1][j]+i_LPF[1]*temp[i][j]+i_LPF[2]*temp[i-1][j];))for(i=0;i<height/2;i++)(for(j=width/2;j<width;j++)(dwt[i][j-height/2]=dwt[i][j];))for(i=0;i<height;i++)(for(j=0;j<width/2;j++)(if(i%2==0)temp[i][j]=dwt[i/2][j];elsetemp[i][j]=0;))for(i=0;i<height;i++)(for(j=0;j<width/2;j++)(if(i==0)

temp1[i][j]=i_HPF[0]*temp[i+2][j]+i_HPF[1]*temp[i+1][j]+i_HPF[2]*temp[i][j]+i_HPF[3]*temp[i+1][j]+i_HPF[4]*temp[i+2][j];elseif(i==1)

temp1[i][j]=i_HPF[0]*temp[i][j]+i_HPF[1]*temp[i-1][j]+i_HPF[2]*temp[i][j]+i_HPF[3]*temp[i+1][j]+i_HPF[4]*temp[i+2][j];elseif(i>1&&i<height-2)

temp1[i][j]=i_HPF[0]*temp[i-2][j]+i_HPF[1]*temp[i-1][j]+i_HPF[2]*temp[i][j]+i_HPF[3]*temp[i+1][j]+i_HPF[4]*temp[i+2][j];elseif(i==height-2)

temp1[i][j]=i_HPF[0]*temp[i-2][j]+i_HPF[1]*temp[i-1][j]+i_HPF[2]*temp[i][j]+i_HPF[3]*temp[i+1][j]+i_HPF[4]*temp[i][j];elseif(i==height-1)

temp1[i][j]=i_HPF[0]*temp[i-2][j]+i_HPF[1]*temp[i-1][j]+i_HPF[2]*temp[i][j]+i_HPF[3]*temp[i-1][j]+i_HPF[4]*temp[i-2][j];))for(i=0;i<height;i++)(for(j=0;j<width/2;j++)(rec[i][j]=rec[i][j]+temp1[i][j];

for(i=height/2;i<height;i++)(for(j=0;j<width/2;j++)(dwt[i-width/2][j]=dwt[i][j];))for(i=0;i<height;i++)(for(j=0;j<width/2;j++)(if(i%2==1)temp[i][j]=dwt[i/2][j];elsetemp[i][j]=0;))for(i=0;i<height;i++)(for(j=0;j<width/2;j++)(if(i==0)temp1[i][j]=i_LPF[0]*temp[i+1][j]+i_LPF[1]*temp[i][j]+i_LPF[2]*temp[i+1][j];elseif(i>0&&i<height-1)temp1[i][j]=i_LPF[0]*temp[i-1][j]+i_LPF[1]*temp[i][j]+i_LPF[2]*temp[i+1][j];elseif(i==height-1)temp1[i][j]=i_LPF[0]*temp[i-1][j]+i_LPF[1]*temp[i][j]+i_LPF[2]*temp[i-1][j];))for(i=height/2;i<height;i++)(for(j=width/2;j<width;j++)(dwt[i-height/2][j-width/2]=dwt[i][j];))for(i=0;i<height;i++)(for(j=0;j<width/2;j++)(if(i%2==0)temp[i][j]=dwt[i/2][j];elsetemp[i][j]=0;))for(i=0;i<height;i++)(for(j=0;j<width/2;j++)(if(i==0)

dwt[i][j]=i_HPF[0]*temp[i+2][j]+i_HPF[1]*temp[i+1][j]+i_HPF[2]*temp[i][j]+i_HPF[3]*temp[i+1][j]+i_HPF[4]*temp[i+2][j];elseif(i==1)

dwt[i][j]=i_HPF[0]*temp[i][j]+i_HPF[1]*temp[i-1][j]+i_HPF[2]*temp[i][j]+i_HPF[3]*temp[i+1][j]+i_HPF[4]*temp[i+2][j];elseif(i>1&&i<height-2)

dwt[i][j]=i_HPF[0]*temp[i-2][j]+i_HPF[1]*temp[i-1][j]+i_HPF[2]*temp[i][j]+i_HPF[3]*temp[i+1][j]+i_HPF[4]*temp[i+2][j];el

溫馨提示

  • 1. 本站所有資源如無(wú)特殊說(shuō)明，都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等，如果需要附件，請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶(hù)所有。
  • 3. 本站RAR壓縮包中若帶圖紙，網(wǎng)頁(yè)內(nèi)容里面會(huì)有圖紙預(yù)覽，若沒(méi)有圖紙預(yù)覽就沒(méi)有圖紙。
  • 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
  • 5. 人人文庫(kù)網(wǎng)僅提供信息存儲(chǔ)空間，僅對(duì)用戶(hù)上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理，對(duì)用戶(hù)上傳分享的文檔內(nèi)容本身不做任何修改或編輯，并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
  • 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容，請(qǐng)與我們聯(lián)系，我們立即糾正。
  • 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶(hù)因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。

評(píng)論

0/150

提交評(píng)論