?? wvlttransdib.cpp
字號:
// WvltTrans.cpp: implementation of the CWvltTrans class.
//
//////////////////////////////////////////////////////////////////////
#include "stdafx.h"
#include "windowsx.h"
#include "math.h"
#include "MainFrm.h"
#include "DynSplitView2.h"
#include "WvltTransDib.h"
#ifdef _DEBUG
#undef THIS_FILE
static char THIS_FILE[]=__FILE__;
#define new DEBUG_NEW
#endif
//////////////////////////////////////////////////////////////////////
// Construction/Destruction
//////////////////////////////////////////////////////////////////////
CWvltTransDib::CWvltTransDib()
{
}
CWvltTransDib::~CWvltTransDib()
{
}
/********************************************************************************
*函數(shù)描述： Hangbianhuan實現(xiàn)小波行變換
*函數(shù)參數(shù)： 無
*函數(shù)返回值：函數(shù)無返回值
*********************************************************************************/
void CWvltTransDib::Hangbianhuan()
{
int i,j;
LONG wide,height;
LPBYTE temp1,m_pData2;
wide=this->GetWidth();
height=this->GetHeight();
m_pData2=this->GetData2();
int nWide=wide/2;
//分配臨時數(shù)據(jù)空間
temp1 = new BYTE[height*wide];
//從設(shè)備緩存中獲得原始圖像數(shù)據(jù)
for(j = 0; j < height; j ++)
{
for(i = 0; i < nWide; i ++)
{
int w = i *2;
temp1[j*wide+i] = m_pData2[j*wide+w]; //偶
temp1[j*wide+nWide+i] = m_pData2[j*wide+w+1]; //奇
}
}
//通過圖像的差分，完成小波變換
for(j=0; j<height; j++)
{
for(i=0; i<nWide-1; i++)
{
temp1[j*wide+nWide-1+i] =temp1[j*wide+(nWide-1)+i] - temp1[j*wide+i]+128;
}
}
//小波經(jīng)過處理后，放入顯示緩存中
for(j=0; j<(int)height; j++)
{
for(i=0; i<(int)wide; i++)
{
m_pData2[j*wide+i] = temp1[j*wide+i];
}
}
//刪除臨時的數(shù)據(jù)空間
delete temp1;
}
/********************************************************************************
*函數(shù)描述： Liebianhuan實現(xiàn)小波列變換
*函數(shù)參數(shù)： 無
*函數(shù)返回值：函數(shù)無返回值
*********************************************************************************/
void CWvltTransDib::Liebianhuan()
{
int i,j;
LONG wide,height;
LPBYTE temp,m_pData2;
wide=this->GetWidth();
height=this->GetHeight();
m_pData2=this->GetData2();
int nHeight=height/2;
//分配臨時數(shù)據(jù)空間
temp = new BYTE[height*wide];
for(i = 0; i < wide; i ++)
{
for(j = 0; j < nHeight; j ++)
{
int h = j *2 ;
temp[j*wide+i] = m_pData2[h*wide+i]; //even
temp[(nHeight + j)*wide+i] = m_pData2[(h+1)*wide+i]; //odd
}
}
//通過圖像的差分，完成小波變換
for(i=0; i<wide; i++)
{
for(j=0; j<nHeight-1; j++)
{
temp[ j*wide+i]=temp[(nHeight + j)*wide+i]- temp[j*wide+i]+128;
}
}
for(j=0; j<(int)height; j++)
{
for(i=0; i<(int)wide; i++)
{
m_pData2[j*wide+i] = temp[j*wide+i];
}
}
//刪除臨時的數(shù)據(jù)空間
delete temp;
}
/********************************************************************************
*函數(shù)描述： Once實現(xiàn)小波變換
*函數(shù)參數(shù)： int n 小波變換的層數(shù)
*函數(shù)返回值：函數(shù)無返回值
*********************************************************************************/
void CWvltTransDib::Once(int n)
{
int i,j;
LONG wide,height,nWide,nHeight;
LPBYTE temp1, temp2,m_pData2;
for(int k=1;k<=n;k++)
{
wide=this->GetWidth();
height=this->GetHeight();
m_pData2=this->GetData2();
//分配臨時數(shù)據(jù)空間
temp1 = new BYTE[height*wide];
temp2 = new BYTE[height*wide];
nWide=wide/pow(2,k);
nHeight=height/pow(2,k);
//完成行變換
for(j = height*(1-1/pow(2,(k-1))); j < height; j ++)
{
for(i = 0; i < nWide; i ++)
{
int w = i *2;
temp1[j*wide+i] = m_pData2[j*wide+w]; //偶
temp1[j*wide+nWide+i] = m_pData2[j*wide+w+1]; //奇
}
}
//通過圖像的差分，完成小波變換
for(j= height*(1-1/pow(2,(k-1))); j<height; j++)
{
for(i=0; i<nWide-1; i++)
{
temp1[j*wide+nWide-1+i] =temp1[j*wide+nWide-1+i] - temp1[j*wide+i]+128;
}
}
//完成列變換
for(i = 0; i < wide/pow(2,k-1); i ++)
{
for(j = height*(1-1/pow(2,(k-1))); j < height*(1-1/pow(2,(k-1)))+nHeight; j ++)
{
int m, h;
m=height*(1-1/pow(2,(k-1)));
h= (j-m)*2 ;
temp2[j*wide+i] = temp1[(m+h)*wide+i]; //even
temp2[(nHeight + j)*wide+i] = temp1[(m+h+1)*wide+i]; //odd
}
}
//通過圖像的差分，完成小波變換
for(i=0; i<wide/pow(2,k-1); i++)
{
for(j = height*(1-1/pow(2,(k-1))); j < height*(1-1/pow(2,(k-1)))+nHeight; j ++)
{
temp2[ j*wide+i] = temp2[j*wide+i]-temp2[(nHeight + j)*wide+i]+128;
}
}
//小波經(jīng)過處理后，放入顯示緩存中
for(j = height*(1-1/pow(2,(k-1))); j < height; j ++)
{
for(i=0; i<wide/pow(2,k-1); i++)
{
m_pData2[j*wide+i] = temp2[j*wide+i];
}
}
//刪除臨時的數(shù)據(jù)空間
delete temp1;
delete temp2;
}
}
/********************************************************************************
*函數(shù)描述： IDWT 實現(xiàn)小波反變換
*函數(shù)參數(shù)： int n 小波變換的層數(shù)
*函數(shù)返回值：函數(shù)無返回值
*********************************************************************************/
void CWvltTransDib::IDWT(int n)
{
int i,j;
LONG wide,height,nWide,nHeight;
LPBYTE temp1, temp2,temp3,m_pData2;
for(int k=n;k>=1;k--)
{
wide=(this->GetWidth());
height=(this->GetHeight());
m_pData2=this->GetData2();
//分配臨時數(shù)據(jù)空間
temp1 = new BYTE[height*wide];
temp2 = new BYTE[height*wide];
temp3 = new BYTE[height*wide];
nWide=wide/pow(2,k);
nHeight=height/pow(2,k);
memcpy(temp1,m_pData2,height*wide);
for(i=0; i<wide/pow(2,k-1); i++)
{
for(j = height*(1-1/pow(2,(k-1))); j < height*(1-1/pow(2,(k-1)))+nHeight; j ++)
{
temp1[ j*wide+i]=temp1[(nHeight + j)*wide+i]+ temp1[j*wide+i]-128;
}
}
for(i = 0; i < wide/pow(2,k-1); i ++)
{
for(j = height*(1-1/pow(2,(k-1))); j < height*(1-1/pow(2,(k-1)))+nHeight; j ++)
{
int m, h;
m=height*(1-1/pow(2,(k-1)));
h= (j-m)*2 ;
temp2[(m+h)*wide+i]=temp1[j*wide+i];
temp2[(m+h+1)*wide+i]=temp1[(nHeight + j)*wide+i];
}
}
for(j= height*(1-1/pow(2,(k-1))); j<height; j++)
{
for(i=0; i<nWide-1; i++)
{
temp2[j*wide+nWide+i] += temp2[j*wide+i]-128;
}
}
for(j = height*(1-1/pow(2,(k-1))); j < height; j ++)
{
for(i = 0; i < nWide; i ++)
{
int w = i *2;
temp3[j*wide+w]=temp2[j*wide+i];
temp3[j*wide+w+1]=temp2[j*wide+nWide+i];
}
}
//小波經(jīng)過處理后，放入顯示緩存中
for(j = height*(1-1/pow(2,(k-1))); j < height; j ++)
{
for(i=0; i<wide/pow(2,k-1); i++)
{
m_pData2[j*wide+i] = temp3[j*wide+i];
}
}
//刪除臨時的數(shù)據(jù)空間
delete temp1;
delete temp2;
}
}
/********************************************************************************
*函數(shù)描述： LowFilter實現(xiàn)小波低通濾波效果
*函數(shù)參數(shù)： int n 小波變換層數(shù)
*函數(shù)返回值：函數(shù)無返回值
*********************************************************************************/
void CWvltTransDib::LowFilter(int n)
{
int i,j;
LONG wide,height,nWide,nHeight;
LPBYTE m_pData2;
wide=this->GetWidth();
height=this->GetHeight();
m_pData2=this->GetData2();
nWide=wide/2;
nHeight=height/2;
Once(n);
for(j=0;j<height;j++)
{
for(i=0;i<wide;i++)
{
if((j>=height*(1-1/pow(2,(n-1))))&&(i<wide/pow(2,n)))
continue;
else
m_pData2[j*wide+i] =(BYTE)128;
}
}
IDWT(n);
}
/********************************************************************************
*函數(shù)描述： HighFilter 實現(xiàn)高通濾波效果
*函數(shù)參數(shù)： int n 小波變換層數(shù)
*函數(shù)返回值：函數(shù)無返回值
*********************************************************************************/
void CWvltTransDib::HighFilter(int n)
{
int i,j;
LONG wide,height,nWide,nHeight;
LPBYTE m_pData2;
wide=this->GetWidth();
height=this->GetHeight();
m_pData2=this->GetData2();
nWide=wide/2;
nHeight=height/2;
LONG lLineBytes = (wide*3+3)/4 *4;
Once(n);
for(j=height*(1-1/pow(2,(n-1)));j<height;j++)
{
for(i=0;i<wide/pow(2,n);i++)
{
m_pData2[j*wide+i] =(BYTE)128;
}
}
IDWT(n);
}
?? 快捷鍵說明
復(fù)制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -