opencv小波文档.docVIP

OpenCV里面没有提供小波变换函数确实是一大遗憾，网上找点小波变换的代码吧，大多是针对BMP直接操作的，还没有发现非常适合OpenCV的，没办法，自己动手手写了一个，意在抛砖引玉，欢迎大家批评，拍砖头，砸鸡蛋。提供函数DWT()和IDWT()，前者完成任意层次的小波变换，后者完成任意层次的小波逆变换。输入图像要求必须是单通道浮点图像，对图像大小也有要求（1层变换：w,h必须是2的倍数；2层变换：w,h必须是4的倍数；3层变换：w,h必须是8的倍数......），变换后的结果直接保存在输入图像中。1、函数参数简单，图像指针pImage和变换层数nLayer。2、一个函数直接完成多层次二维小波变换，尽量减少下标运算，避免不必要的函数调用，以提高执行效率。3、变换过程中，使用了一个指针数组pData用于保存每行数据的起始位置，pRow和pColumn用于保存一行和一列临时数据，用于奇偶分离或合并，内存消耗较少。// 二维离散小波变换（单通道浮点图像）void DWT(IplImage *pImage, int nLayer){? ?// 执行条件? ?if (pImage)? ?{? ?? 爄f (pImage-nChannels == 1 ?牋 牋 爌Image-depth == IPL_DEPTH_32F ?牋 牋 ?(pImage-width nLayer) nLayer) == pImage-width ?牋 牋 ?(pImage-height nLayer) nLayer) == pImage-height)?牋 爗?牋 牋 爄nt??爄, x, y, n;?牋 牋 爁loat?爁Value?? 0;?牋 牋 爁loat?爁Radius?= sqrt(2.0f);?牋 牋 爄nt??爊Width?? pImage-width;?牋 牋 爄nt??爊Height?= pImage-height;?牋 牋 爄nt??爊HalfW?? nWidth / 2;?牋 牋 爄nt??爊HalfH?? nHeight / 2;?牋 牋 爁loat **pData??= new float*[pImage-height];?牋 牋 爁loat?*pRow??? new float[pImage-width];?牋 牋 爁loat?*pColumn?= new float[pImage-height];?牋 牋 爁or (i = 0; i pImage-height; i++)?牋 牋 爗?牋 牋 牋 爌Data[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[y][nHalfW + i] = pRow[x + 1];? ?? ?? ?? ?? ?}? ?? ?? ?? ?? ?// 提升小波变换? ?? ?? ?? ?? 爁or (i = 0; i nHalfW - 1; i++)?牋 牋 牋 牋 爗?牋 牋 牋 牋 牋 爁Value = (pData[y][i] + pData[y][i + 1]) / 2;?牋 牋 牋 牋 牋 爌Data[y][nHalfW + i] -= fValue;?牋 牋 牋 牋 爙?牋 牋 牋 牋 爁Value = (pData[y][nHalfW