opencv小波变换（OpenCV小波变换）.docVIP

opencv小波变换（OpenCV小波变换） #包括数据\ highgui \ highgui HPP 。 #包括数学。” #包含iostream 使用命名空间； / /二维离散变换小波源代码 无效的DWT（IplImage * pimage，int nlayer） { / /执行条件 如果（pimage） { 如果（pimage - nchannels = = 1 pimage -深度= = ipl_depth_32f （（pimage -宽度 nlayer） nlayer）= = pimage -宽度 （（pimage -高 nlayer） nlayer）= = pimage -高度） { int i，x，y，n； 浮：= 0； 浮fradius = SQRT（2.0f）； int nWidth = pimage -宽度； 国际nheight = pimage -高度； 国际nhalfw = nWidth / 2； 国际nhalfh = nheight / 2； 浮** pdata =新浮* [ pimage -高度]； 浮*船头=新浮pimage - [宽度]； 浮* pcolumn =新浮pimage - [高度]； 为（i = 0；i pimage -高度；i++）数据[我] =（float *）（pimage -数据+ pimage - widthstep *我）； / /多层小波变换 对于（n = 0；n＜nlayer；n + +，nWidth / = 2，nheight / = 2，nhalfw / = 2，nhalfh / = 2） { / /水平变换 为（y = 0；Y＜nheight；Y + +） { / /奇偶分离 memcpy（船头，pdata [Y]，sizeof（float）* nWidth）； 为（i = 0；i nhalfw；i++） { X =我* 2； 数据[年] [我] =船头[X]； 数据[年] [ nhalfw +我=船头] [ x + 1 ]； } / /提升小波变换 为（i = 0；i nhalfw - 1；i++） { F值=（pdata [年] [我] +数据[年] [我] + 1）/ 2； 数据[年] [我]：nhalfw + =； } F值=（pdata [年] [ 1 ] nhalfw +数据[年] [ 2 ] nhalfw）/ 2； 数据[年] [ 1 ] =：nWidth； F值=（pdata [ ] [ ] [ ] nhalfw +数据+ 1 ] [ nhalfw）/ 4； 数据[年] [ 0 ] + =； 为（i = 1；i nhalfw；i++） { F值=（pdata [ ] [ ] nhalfw +我+数据[年] [我] nhalfw + 1）/ 4； 数据[年] [我] + =； } / /频带系数 为（i = 0；i nhalfw；i++） { 数据[年] [我] * = fradius； 数据[年] [ nhalfw +我] / = fradius； } } / /垂直变换 为（x = 0；x＜nWidth；X + +） { / / 奇偶分离 for (i = 0; i nhalfh; i + +) { y = i * 2; pcolumn [in] = pdata [y] [x]; pcolumn [nhalfh +] = pdata [y + 1] [x]; } for (i = 0; i nheight; i + +) { pdata [in] [x] = pcolumn [in]; } / / 提升小波变换 for (i = 0; i nhalfh - 1; i + +) { fvalue = (pdata [in] [x] + pdata [i + 1] [x]) / 2; pdata [nhalfh +] [x] = fvalue; } fvalue = (pdata [nhalfh - 1] [x] + pdata [nhalfh - 2] [x]) / 2; pdata [nheight - 1] [x] = fvalue; fvalue = (pdata [nhalfh] [x] + pdata [nhalfh + 1] [x]) / 4; pdata [0] [x] + = fvalue; for (i = 1; in nhalfh; i + +) { fvalue = (pdata [nhalfh +] [x] + pdata [nhalfh + - 1] [x]) / 4; pdata [in] [x] + = fvalue; } / / 频带系数 for (i = 0; i nhalfh; i