二维变换在计算机图形学中的应用

**************************************************************************
*  文件名：FrequencyTrans.cpp
*
*  正交变换API函数库：
*
*  FFT()    - 快速付立叶变换
*  IFFT()    - 快速付立叶反变换
*  DCT()    - 离散余弦变换
*  IDCT                - 离散余弦反变换
*  WALSH()    - 沃尔什－哈达玛变换
*  IWALSH()            - 沃尔什－哈达玛反变换
*
*  Fourier()   - 图像的付立叶变换
*  InFourier()         - 图象的傅立叶逆变换
*  DIBDct()   - 图像的离散余弦变换
*  DIBIDct()   - 图像的离散余弦反变换
*  DIBWalsh()   - 图像的沃尔什－哈达玛变换
*  DIBIWalsh()   - 图像的沃尔什－哈达玛反变换
*  WT()                - 小波变换
*  IWT()               - 小波反变换
*
*************************************************************************/
#include "FrequencyTrans.h"
#include "stdafx.h"
#include "DIBAPI.h"
#include <math.h>
#include <direct.h>
#include <complex>
using namespace std;
// 常数π

#define PI 3.1415926535
//小波变换所用部分
#define DD 23
float h[23]={-0.002,-0.003,0.006,0.006,-0.013,
        -0.012, 0.030, 0.023,-0.078,-0.035,
        0.307,0.542,0.307,-0.035,-0.078,
        0.023,0.030,-0.012,-0.013,0.006,
        0.006,-0.003,-0.002};
float g[23]={0.002,-0.003,-0.006,0.006,0.013,-0.012,
  -0.030,0.023,0.078,-0.035,-0.307,0.542,
  -0.307,-0.035,0.078,0.023,-0.030,-0.012,
  0.013,0.006,-0.006,-0.003,0.002};
//........................（中间部分略）

/***************************************************************
* *   二维正交小波变换及其逆变换（有待验证）
* ***************************************************************/
int edge_function(int x,int xsize)
{     if(x<0)
       x=-x;
    if(x>=xsize)
      x=xsize*2-x-2;
    return(x);
} ///////////////////////////////////////////////////////////////
/*小波变换
//////////////////////////////////////////////////////////*/
VOID WINAPI WT( LONG xs, LONG ys,LONG xsize,LONG ysize,/*FLOAT *image)unsigned */float * image)
{    int i,j,k,n;
   float temp1,temp2;
   float *bufferx,*buffery;