开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > 图形图像 > 波变换 > 查看源码
WFilter.cpp
资源名称：imagefusionwavelettransform.rar　[点击查看]
上传用户：aqingfeng
上传日期：2014-03-25
资源大小：1839k
文件大小：6k
源码类别：

波变换

开发平台：

Visual C++

WFilter.cpp：源码内容
  1. // WFilter.cpp: implementation of the CWFilter class.
  2. //
  3. //////////////////////////////////////////////////////////////////////
  5. #include "stdafx.h"
  6. #include "小波变换.h"
  7. #include "WFilter.h"
  8. #include "WvltTrans.h"
  9. #include <math.h>
  11. #ifdef _DEBUG
  12. #undef THIS_FILE
  13. static char THIS_FILE[]=__FILE__;
  14. #define new DEBUG_NEW
  15. #endif
  17. //////////////////////////////////////////////////////////////////////
  18. // Construction/Destruction
  19. //////////////////////////////////////////////////////////////////////
  21. CWFilter::CWFilter()
  22. {
  24. }
  26. CWFilter::~CWFilter()
  27. {
  29. }
  30. /**********************************************************************************
  31. *函数描述： LPass_Filter实现小波变换的低通滤波，减少图像的边缘信息，使图像更加平滑*
  32. *函数参数： unsigned char **pData：图像小波变换后的小波系数矩阵   *
  33. * int nHeight  ：图像属性参数，数值为原始图像的高度值   *
  34. * int nWidth  ：图像属性参数，数值为原始图像的宽度值   *
  35. * int nLayer  ：低通滤波器的滤波阶数，数值为小波变换的层数   *
  36. *函数返回值：函数无返回值，小波系数参数是利用指针进行调用   *
  37. **********************************************************************************/
  39. void CWFilter::LPass_Filter(short **pData, int nHeight, int nWidth, int nLayer, int scale)
  40. {
  41. short **spOriginData, **spTransData0, **spTransData1;
  42. int iHeight = nHeight, iWidth = nWidth;
  43. spOriginData = pData;
  44. spTransData0 = new short *[nHeight];
  45. spTransData1 = new short *[nHeight];
  46. m_FilterData = new short * [nHeight];
  47. for(int i = 0; i < nHeight; i ++)
  48. {
  49. spTransData0[i] = new short [nWidth];
  50. spTransData1[i] = new short [nWidth];
  51. m_FilterData[i] = new short [nWidth];
  52. for(int j = 0; j < nWidth; j ++)
  53. {
  54. m_FilterData[i][j] = 0;
  55. }
  56. }
  57. CWvltTrans *Trans;
  58. Trans->DWT_TriLayers(spOriginData, spTransData0, spTransData1, nHeight, nHeight/2, nWidth, nWidth/2, nLayer, 1.414);
  59. for(i = 0; i < (int) nHeight / (4 * scale); i ++)
  60. {
  61. for(int j = 0; j < (int) nWidth / (4 * scale); j ++)
  62. {
  63. m_FilterData[i][j] = spTransData1[i][j];
  64. }
  65. }
  66. delete spTransData0;
  67. delete spTransData1;
  68. DWT_Inverse(iHeight, iWidth, 3);
  69. }
  71. /****************************************************************************
  72. *函数描述： HPass_Filter完成图像的高通滤波处理，得到图像的边缘信息 *
  73. *函数参数： short **pData ：二维指针，存放正则化后的小波系数 *
  74. * int   nHeight ：图像属性参数，数值为原始图像的高度值 *
  75. * int   nHeight_H ：图像属性参数，数值为原始图像高度值的一半 *
  76. * int   nWidth ：图像属性参数，数值为原始图像的宽度值 *
  77. * int   nLayer ：小波变换的层数 *
  78. *函数返回值：函数无返回值，小波系数参数是利用指针进行调用 *
  79. *****************************************************************************/
  80. void CWFilter::HPass_Filter(short **pData, int nHeight, int nWidth, int nLayer)
  81. {
  82. short **spOriginData, **spTransData0, **spTransData1;
  83. int iHeight = nHeight, iWidth = nWidth;
  84. spOriginData = pData;
  85. spTransData0 = new short *[nHeight];
  86. spTransData1 = new short *[nHeight];
  87. m_FilterData = new short * [nHeight];
  88. for(int i = 0; i < nHeight; i ++)
  89. {
  90. spTransData0[i] = new short [nWidth];
  91. spTransData1[i] = new short [nWidth];
  92. m_FilterData[i] = new short [nWidth];
  93. for(int j = 0; j < nWidth; j ++)
  94. {
  95. m_FilterData[i][j] = 0;
  96. }
  97. }
  98. CWvltTrans *Trans;
  99. Trans->DWT_TriLayers(spOriginData, spTransData0, spTransData1, nHeight, nHeight/2, nWidth, nWidth/2, nLayer, 1.414);
  100. for(i = (int) nHeight / 8; i <  nHeight; i ++)
  101. {
  102. for(int j = (int) nWidth / 8; j <  nWidth ; j ++)
  103. {
  104. m_FilterData[i][j] = spTransData1[i][j];
  105. }
  106. }
  107. delete spTransData0;
  108. delete spTransData1;
  109. DWT_Inverse(iHeight, iWidth, 3);
  110. }
  112. /************************************************************************
  113. *函数描述： DWT_Inverse完成若干层小波变换的逆变换操作 *
  114. *函数参数： int   nHeight ：图像属性参数，数值为原始图像的高度值 *
  115. * int   nWidth ：图像属性参数，数值为原始图像的宽度值 *
  116. * int   nLayer ：小波变换的层数 *
  117. *************************************************************************/
  118. void CWFilter::DWT_Inverse(int nHeight, int nWidth, int nLayer)
  119. {
  120. int iHeight = nHeight / 4, iWidth = nWidth / 4;
  121. int iHeight_H = iHeight / 2, iWidth_H = iWidth / 2;
  122. short **spTransData0, **spTransData1;
  123. spTransData0 = new short * [nHeight];
  124. spTransData1 = new short * [nHeight];
  126. for(int i = 0; i < nHeight; i++)
  127. {
  128. spTransData0[i] = new short [nWidth];
  129. spTransData1[i] = new short [nWidth];
  130. }
  132. CWvltTrans WTrans;
  133. for( i = nLayer; i >= 1; i--)
  134. {
  135. WTrans.DWTi_Once(spTransData0, spTransData1, m_FilterData, iHeight, iHeight_H, iWidth, iWidth_H, i, 1.414);
  136. iHeight <<= 1; iWidth <<= 1;
  137. iHeight_H = iHeight / 2; iWidth_H = iWidth / 2;
  138. }
  139. for(int y = 0; y < nHeight; y++)
  140. {
  141. for(int x = 0; x < nWidth; x++)
  142. {
  143. m_FilterData[y][x] = spTransData0[y][x];
  144. }
  145. }
  148. }
  150. /****************************************************************************
  151. *函数描述： HPass_Filter完成图像的高通滤波处理，得到图像的边缘信息 *
  152. *函数参数： short **pData ：二维指针，存放正则化后的小波系数 *
  153. * int   nHeight ：图像属性参数，数值为原始图像的高度值 *
  154. * int   nHeight_H ：图像属性参数，数值为原始图像高度值的一半 *
  155. * int   nWidth ：图像属性参数，数值为原始图像的宽度值 *
  156. * int   nLayer ：小波变换的层数 *
  157. *函数返回值：函数无返回值，小波系数参数是利用指针进行调用 *
  158. *****************************************************************************/
  160. void CWFilter::HPass_Filter2(short **pData, int nHeight, int nWidth, int nLayer)
  161. {
  162. short **spOriginData, **spTransData0, **spTransData1;
  163. int iHeight = nHeight, iWidth = nWidth;
  164. spOriginData = pData;
  165. spTransData0 = new short *[nHeight];
  166. spTransData1 = new short *[nHeight];
  167. m_FilterData = new short * [nHeight];
  168. for(int i = 0; i < nHeight; i ++)
  169. {
  170. spTransData0[i] = new short [nWidth];
  171. spTransData1[i] = new short [nWidth];
  172. m_FilterData[i] = new short [nWidth];
  173. }
  174. CWvltTrans *Trans;
  175. Trans->DWT_TriLayers(spOriginData, spTransData0, spTransData1, nHeight, nHeight/2, nWidth, nWidth/2, nLayer, 1.414);
  176. for(i = 0; i <  nHeight; i ++)
  177. {
  178. for(int j = 0; j <  nWidth ; j ++)
  179. {
  180. m_FilterData[i][j] = spTransData1[i][j];
  181. if(i < nHeight / 8 && j < nWidth / 8)
  182. m_FilterData[i][j] /= 2;
  183. }
  184. }
  185. delete spTransData0;
  186. delete spTransData1;
  187. DWT_Inverse(iHeight, iWidth, 3);
  188. }