OpenCV实现图像小波变换与逆变换代码分享

/********************************************************************************
*函数描述： DWT 实现小波反变换 单通道浮点图像    
*函数参数： int nLayer 小波变换的层数 1层变换：w,h必须是2的倍数；2层变换：w,h必须是4的倍数；3层变换：w,h必须是8的倍数......
*函数返回值：函数无返回值
*********************************************************************************/
void DWT(IplImage *pImage,int nLayer)
{
if (pImage->nChannels == 1 &&
pImage->depth == IPL_DEPTH_32F && 
((pImage->width >> nLayer) << nLayer) == pImage->width &&
((pImage->height >> nLayer) << nLayer) == pImage->height) //1层变换wh是2的倍数；2层变换wh是4的倍数；3层变换wh是8的倍数...
{//仅处理单通道浮点图像
int     i, x, y, n;
float   fValue   = 0;
float   fRadius  = sqrt(3.0f);  //设置频带系数
int     nWidth   = pImage->width;
int     nHeight  = pImage->height;
int     nHalfW   = nWidth / 2;
int     nHalfH   = nHeight / 2;
float **pData    = new float*[pImage->height];  //指针数组pData用于保存每行数据的起始位置
float  *pRow     = new float[pImage->width];  //pRow用于保存一行临时数据
float  *pColumn  = new float[pImage->height]; //pColumn用于保存一列临时数据
for (i = 0; i < pImage->height; i++)
{
pData[i] = (float*) (pImage->imageData + pImage->widthStep * i);  //保存每行数据的起始位置
}

// 多层小波变换
for (n = 0; n < nLayer; n++, nWidth /= 2, nHeight /= 2, nHalfW /= 2, nHalfH /= 2)
{

// 水平变换
for (y = 0; y < nHeight; y++)
{
// 奇偶分离
memcpy(pRow, pData[y], sizeof(float) * nWidth);
for (i = 0; i < nHalfW; i++)
{
x = i * 2;
pData[y][i] = pRow[x]; //pData左半部分--原图偶
pData[y][nHalfW + i] = pRow[x + 1]; //pData右半部分--原图奇
}
// 提升小波变换
for (i = 0; i < nHalfW - 1; i++)
{
fValue = (pData[y][i] + pData[y][i + 1]) / 2;
pData[y][nHalfW + i] -= fValue; //pData行后半部分
}
fValue = (pData[y][nHalfW - 1] + pData[y][nHalfW - 2]) / 2;
pData[y][nWidth - 1] -= fValue; //pData行最后一个元素的处理
fValue = (pData[y][nHalfW] + pData[y][nHalfW + 1]) / 4;
pData[y][0] += fValue;  //pData行第一个元素的处理
for (i = 1; i < nHalfW; i++)
{
fValue = (pData[y][nHalfW + i] + pData[y][nHalfW + i - 1]) / 4;
pData[y][i] += fValue;  //pData行前半部分
}
// 频带系数
for (i = 0; i < nHalfW; i++)
{
pData[y][i] *= fRadius;  //pData最终行前半部分