小波变换 C++

				版权声明：本文为博主原创文章，转载时请注明来源。有问题请e-mail:mr.baishu@gmail.com					https://blog.youkuaiyun.com/u011177305/article/details/51313718				</div>
							<link rel="stylesheet" href="https://csdnimg.cn/release/phoenix/template/css/ck_htmledit_views-f57960eb32.css">
							            <link rel="stylesheet" href="https://csdnimg.cn/release/phoenix/template/css/ck_htmledit_views-f57960eb32.css">
					<div class="htmledit_views" id="content_views">

</pre><pre name="code" class="cpp"><p><span style="font-family:KaiTi_GB2312;font-size:14px;">看完小波变换的理论后想实现一线，同时发现opencv中并没有小波变换的函数（还是新版的有了，我不知道啊，please give Me a feedback if U guy know）</span></p><p>...............省略线。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。</p><p>	int N = 16;</p>	//该部分首先将图划分为16*16的block区域
	//在每一个循环当中，在每一次循环的所有的操作中，都表示在一个block的操作
	//其中包括一次小波变换得到四个区域
	//一次空间Spatial Energy计算，将每个Block的能量块保存到EnergyBlock（[8,8]大小）中
	//并且将这些能量块的数值保存为vector作为SVM的训练样本
	for (int x=0;x<inputImage.rows;x+=N)
	{
		Mat ROIBlock = inputImage(Range(x,x+N),Range(x,x+N));
 
		//小波变换：
		//waveletblockLH为LH列变换后的结果
		//waveletblock为waveletblock的LH行变换的结果
		Mat WaveletBlockLH=Mat::zeros(16, 16, ROIBlock.type());
		Mat WaveletBlock=Mat::zeros(16, 16, ROIBlock.type());
		Mat EnergyBlock = Mat::zeros(8,8,ROIBlock.type());
		
		//计算
		for (int i=0,j=0;j<ROIBlock.cols;i+=2,j++)
		{
			//判断是否为前八列，如果成立则
			//		j=.i+.(i+1)
			//否则：j=.(i-8)-.(i+1-8)
			if (floor(j / 8) == 0)
			{
				
				WaveletBlockLH.col(j) = ROIBlock.col(i) + ROIBlock.col(i + 1);
			}
			else
 
				WaveletBlockLH.col(j)=ROIBlock.col(i-N)-ROIBlock.col(i-N+1);
			
		}
 
		for (int i = 0,j = 0; j<WaveletBlockLH.rows; i += 2, j++)
		{
			//判断是否为前八列，如果成立则
			//		j=.i+.(i+1)
			//否则：j=.(i-8)-.(i+1-8)
			if (floor(j / 8) == 0)
			{
 
				WaveletBlock.row(j) = WaveletBlockLH.row(i) + WaveletBlockLH.row(i + 1);
			}
			else
 
				WaveletBlock.row(j) = WaveletBlockLH.row(i - N) - WaveletBlockLH.row(i - N+1);
 
		}