- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
#include stdio.h
#include math.h
#define height 256
#define width 256
void twoDDWT( double org[][width], double dwt[][width], int NumofBand);
void oneDDWT( int flag, double org[][width], double dwt[][width]);
void twoDIDWT( double dwt[][width], double rec[][width], int NumofBand);
void oneDIDWT( int flag, double dwt[][width], double rec[][width]);
int main()
{
int i, j;
double org[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);
}
}
return 0;
}
void twoDDWT( double org[][width], double dwt[][width], int NumofBand)
{
if (
原创力文档


文档评论(0)