dct变换源代码
来源:互联网 发布:网易广州 知乎 编辑:程序博客网 时间:2024/06/05 19:07
#define EXTRA_NAME "@dcttransform."
#include "loadbmp.h"
#define Point(x,y) lpPoints[(x)+(y)*nWidth]
void FFT(COMPLEX * TD, COMPLEX * FD, int power)
{
int count;
int i,j,k,bfsize,p;
double angle;
COMPLEX *W,*X1,*X2,*X;
count=1<<power;
W=(COMPLEX *)malloc(sizeof(COMPLEX)*count/2);
X1=(COMPLEX *)malloc(sizeof(COMPLEX)*count);
X2=(COMPLEX *)malloc(sizeof(COMPLEX)*count);
for(i=0;i<count/2;i++)
{
angle=-i*pi*2/count;
W[i].re=cos(angle);
W[i].im=sin(angle);
}
memcpy(X1,TD,sizeof(COMPLEX)*count);
for(k=0;k<power;k++)
{
for(j=0;j<1<<k;j++)
{
bfsize=1<<(power-k);
for(i=0;i<bfsize/2;i++)
{
p=j*bfsize;
X2[i+p]=Add(X1[i+p],X1[i+p+bfsize/2]);
X2[i+p+bfsize/2]=Mul(Sub(X1[i+p],X1[i+p+bfsize/2]),W[i*(1<<k)]);
}
}
X=X1;
X1=X2;
X2=X;
}
for(j=0;j<count;j++)
{
p=0;
for(i=0;i<power;i++)
{
if (j&(1<<i)) p+=1<<(power-i-1);
}
FD[j]=X1[p];
}
free(W);
free(X1);
free(X2);
}
void DCT_Pre(double *f, double *F, int power)
{
int i,count;
COMPLEX *X;
double s;
count=1<<power;
X=(COMPLEX *)malloc(sizeof(COMPLEX)*count*2);
memset(X,0,sizeof(COMPLEX)*count*2);
for(i=0;i<count;i++)
{
X[i].re=f[i];
}
FFT(X,X,power+1);
s=1/sqrt(count);
F[0]=X[0].re*s;
s*=sqrt(2);
for(i=1;i<count;i++)
{
F[i]=(X[i].re*cos(i*pi/(count*2))+X[i].im*sin(i*pi/(count*2)))*s;
}
free(X);
}
void Dct()
{
int w=1,h=1,wp=0,hp=0;
while(w*2<=nWidth)
{
w*=2;
wp++;
}
while(h*2<=nHeight)
{
h*=2;
hp++;
}
int x,y;
BYTE *lpPoints=new BYTE[nWidth*nHeight];
GetPoints(lpPoints);
double *f=new double[w*h];
double *W=new double[w*h];
for(y=0;y<h;y++)
{
for(x=0;x<w;x++)
{
f[x+y*w]=Point(x,y);
}
}
for(y=0;y<h;y++)
{
DCT_Pre(&f[w*y],&W[w*y],wp);
}
for(y=0;y<h;y++)
{
for(x=0;x<w;x++)
{
f[x*h+y]=W[x+w*y];
}
}
for(x=0;x<w;x++)
{
DCT_Pre(&f[x*h],&W[x*h],hp);
}
double a;
memset(lpPoints,0,nWidth*nHeight);
for(y=0;y<h;y++)
{
for(x=0;x<w;x++)
{
a=fabs(W[x*h+y]);
if (a>255) a=255;
Point(x,nHeight-y-1)=(BYTE)(a);
}
}
delete f;
delete W;
PutPoints(lpPoints);
delete lpPoints;
}
void main(int argc, char *argv[])
{
if(argc==2)
FileName=argv[1];
else
return;
OpenFile();
Dct();
SaveAs();
}
- dct变换源代码
- dct变换源代码
- DCT变换
- DCT变换
- DCT变换
- DCT变换
- dct变换
- DCT变换
- DCT变换
- dct变换
- DCT变换
- DCT变换
- DCT变换
- DCT变换、DCT反变换、分块DCT变换
- DCT变换、DCT反变换、分块DCT变换
- DCT变换、DCT反变换、分块DCT变换
- DCT变换、DCT反变换、分块DCT变换
- [转] DCT变换、DCT反变换、分块DCT变换
- 面向对象之封装与多态
- 高斯模糊处理(图像复原)源代码
- 人脸检测与识别的站点
- 水平镜像处理源代码
- MongoDB C++应用 简单例子
- dct变换源代码
- xtrabackup 1.6.3安装出现cc1: error: unrecognized command line option "-fdiagnostics-show-option"
- 第二章(编译和链接)
- 我学到的一些关于编程的事儿(翻译)
- 硬盘数据恢复
- 对spring框架的理解
- 带指针的通信结构体32位与64位兼容
- php隐藏域
- Unable to compile class for JSP(不能编译class)