?? testview.cpp
字號:
double ratio=10 * log10((double)(temp/difference));
::GlobalUnlock (hDIBSource);
*/
// 解除鎖定
::GlobalUnlock((HGLOBAL) pDoc->GetHDIB());
}
void CTestView::ADLwavelet_transform_overlap(int **dst, const int width, const int height, const int levels)
{
int l;
int* tmp_lastrow=(int*)malloc((width+height+height)*sizeof(int));
int* tmp_column_in=tmp_lastrow+width;
int* tmp_column_out=tmp_lastrow+width+height;
for(l=0;l<levels;l++)
{ int w=width>>l;
int h=height>>l;
int i;
//首先是行變換
for(i=0;i<h;i++)
{
wavelet_row(tmp_lastrow,dst[i],w);
memcpy(dst[i],tmp_lastrow,w*sizeof(int));
}
//然后是列變換
for(i=0;i<w;i++)
{ int j;
for(j=0;j<h;j++)
tmp_column_in[j]=dst[j][i];
wavelet_row(tmp_column_out,tmp_column_in,h);
for(j=0;j<h;j++)
dst[j][i]=tmp_column_out[j];
}
}
free(tmp_lastrow);
}
void CTestView::inverse_Horizontal_Update(int levels)
{
int i,j,k;
int temp,temp1,temp2;
int x,y;
int *ttt;
ttt=new int[MyHeight];
for(i=0;i<MyWidth*4;i++)
{
for(j=0;j<MyHeight;j++)
{ttt[j]=ori_QuartPixel_interp[j][i];}
for(j=0;j<MyHeight/2;j++)
{
ori_QuartPixel_interp[2*j][i]=ttt[j];
ori_QuartPixel_interp[2*j+1][i]=ttt[j+MyHeight/2];
}
}
for(j=0;j<MyWidth;j++) //第1行更新
{
k=total_preDrection[levels-1][int(j/N)];
//k=Predict_direc[0][int(j/N)];
temp=ori_QuartPixel_interp[0][4*j];
x=4*j+k;
y=4*j-k;
if(x<0)
{x=-x;}
if(y<0)
{y=-y;}
if(x>4*MyWidth-1)
{x=8*(MyWidth-1)-x;}
if(y>4*MyWidth-1)
{y=8*(MyWidth-1)-y;}
temp1=ori_QuartPixel_interp[1][x];
temp2=ori_QuartPixel_interp[1][y];
ori_QuartPixel_interp[0][4*j]=temp-(temp1+temp2)/4;
}
for(j=0;j<MyWidth;j++)//后MyHeight-1行
{
for(i=2;i<MyHeight;i+=2)
{
k=total_preDrection[levels-1][((int)i/M)*(MyWidth/N)+int(j/N)];
//k=Predict_direc[(int)i/M][int(j/N)];
temp=ori_QuartPixel_interp[i][4*j];
x=4*j+k;
y=4*j-k;
if(x<0)
{x=-x;}
if(y<0)
{y=-y;}
if(x>4*MyWidth-1)
{x=8*(MyWidth-1)-x;}
if(y>4*MyWidth-1)
{y=8*(MyWidth-1)-y;}
temp1=ori_QuartPixel_interp[i-1][x];
temp2=ori_QuartPixel_interp[i+1][y];
ori_QuartPixel_interp[i][4*j]=temp-(temp1+temp2)/4;
}
}
for(i=0;i<MyHeight;i++)
for(j=0;j<MyWidth;j++)
{ origin[i][j]=ori_QuartPixel_interp[i][4*j] ;}
delete[] ttt;
ttt=NULL;
}
void CTestView::inverse_Vertical_Update(int levels)
{
int i,j,k;
int temp,temp1,temp2;
int x,y;
int *ttt;
ttt=new int[MyWidth];
for(i=0;i<MyHeight*4;i++)
{
for(j=0;j<MyWidth;j++)
{ttt[j]=ori_QuartPixel_interp1[i][j];}
for(j=0;j<MyWidth/2;j++)
{
ori_QuartPixel_interp1[i][2*j]=ttt[j];
ori_QuartPixel_interp1[i][2*j+1]=ttt[j+MyWidth/2];
}
}
for(j=0;j<MyHeight;j++) //第1列逆更新
{
k=total_preDrection1[levels-1][(int(j/M))*(MyWidth/N)];
//k=Predict_direc1[int(j/M)][0];
temp=ori_QuartPixel_interp1[4*j][0];
x=4*j+k;
y=4*j-k;
if(x<0)
{x=-x;}
if(y<0)
{y=-y;}
if(x>4*MyHeight-1)
{x=8*(MyHeight-1)-x;}
if(y>4*MyHeight-1)
{y=8*(MyHeight-1)-y;}
temp1=ori_QuartPixel_interp1[x][1];
temp2=ori_QuartPixel_interp1[y][1];
ori_QuartPixel_interp1[4*j][0]=temp-(temp1+temp2)/4;
}
for(j=0;j<MyHeight;j++)//后MyWidth-1列
{
for(i=2;i<MyWidth;i+=2)
{
k=total_preDrection1[levels-1][(int(j/M))*(MyWidth/N)+int(i/N)];
//k=Predict_direc1[(int)j/M][int(i/N)];
temp=ori_QuartPixel_interp1[4*j][i];
x=4*j+k;
y=4*j-k;
if(x<0)
{x=-x;}
if(y<0)
{y=-y;}
if(x>4*MyHeight-1)
{x=8*(MyHeight-1)-x;}
if(y>4*MyHeight-1)
{y=8*(MyHeight-1)-y;}
temp1=ori_QuartPixel_interp1[x][i-1];
temp2=ori_QuartPixel_interp1[y][i+1];
ori_QuartPixel_interp1[4*j][i]=temp-(temp1+temp2)/4;
}
}
for(i=0;i<MyHeight;i++)
for(j=0;j<MyWidth;j++)
{ origin[i][j]=ori_QuartPixel_interp1[4*i][j] ;}
delete[] ttt;
ttt=NULL;
}
void CTestView::inverse_Vertical_Prediction(int levels)
{
int i,j,k;
int temp,temp1,temp2;
int x,y;
for(j=1;j<MyWidth-1;j+=2)//前MyWidth-1列
{
for(i=0;i<MyHeight;i++)
{
k=total_preDrection1[levels-1][(int(i/M))*(MyWidth/N)+int(j/N)];
//k=Predict_direc1[(int)i/M][int(j/N)];
temp=ori_QuartPixel_interp1[4*i][j];
x=4*i+k;
y=4*i-k;
if(x<0)
{x=-x;}
if(y<0)
{y=-y;}
if(x>4*MyHeight-1)
{x=8*(MyHeight-1)-x;}
if(y>4*MyHeight-1)
{y=8*(MyHeight-1)-y;}
temp1=ori_QuartPixel_interp1[x][j-1];
temp2=ori_QuartPixel_interp1[y][j+1];
ori_QuartPixel_interp1[4*i][j]=temp+(temp1+temp2)/2;
}
}
for(j=0;j<MyHeight;j++) //第MyWidth-1列
{
k=total_preDrection1[levels-1][(int(j/M))*(MyWidth/N)+MyWidth/N-1];
//k=Predict_direc1[int(j/M)][MyWidth/N-1];
temp=ori_QuartPixel_interp1[4*j][MyWidth-1];
x=4*j+k;
y=4*j-k;
if(x<0)
{x=-x;}
if(y<0)
{y=-y;}
if(x>4*MyHeight-1)
{x=8*(MyHeight-1)-x;}
if(y>4*MyHeight-1)
{y=8*(MyHeight-1)-y;}
temp1=ori_QuartPixel_interp1[x][MyWidth-2];
temp2=ori_QuartPixel_interp1[y][MyWidth-2];
ori_QuartPixel_interp1[4*j][MyWidth-1]=temp+(temp1+temp2)/2;
}
for(i=0;i<MyHeight;i++)
for(j=0;j<MyWidth;j++)
{origin[i][j]=ori_QuartPixel_interp1[4*i][j];}
}
void CTestView::inverse_Horizontal_Prediction(int levels)
{
int i,j,k;
int x,y;
int temp,temp1,temp2;
for(j=0;j<MyWidth;j++)
{
for(i=1;i<MyHeight-1;i+=2)//前MyHeight-1行
{
k=total_preDrection[levels-1][(int(i/M))*(MyWidth/N)+int(j/N)];
//k=Predict_direc[(int)i/M][int(j/N)];
temp=ori_QuartPixel_interp[i][4*j];
x=4*j+k;
y=4*j-k;
if(x<0)
{x=-x;}
if(y<0)
{y=-y;}
if(x>4*MyWidth-1)
{x=8*(MyWidth-1)-x;}
if(y>4*MyWidth-1)
{y=8*(MyWidth-1)-y;}
temp1=ori_QuartPixel_interp[i-1][x];
temp2=ori_QuartPixel_interp[i+1][y];
ori_QuartPixel_interp[i][4*j]=temp+(temp1+temp2)/2;
}
}
for(j=0;j<MyWidth;j++) //第MyHeight-1行
{
k=total_preDrection[levels-1][(MyHeight/M-1)*(MyWidth/N)+int(j/N)];
//k=Predict_direc[MyHeight/M-1][int(j/N)];
temp=ori_QuartPixel_interp[MyHeight-1][4*j];
x=4*j+k;
y=4*j-k;
if(x<0)
{x=-x;}
if(y<0)
{y=-y;}
if(x>4*MyWidth-1)
{x=8*(MyWidth-1)-x;}
if(y>4*MyWidth-1)
{y=8*(MyWidth-1)-y;}
temp1=ori_QuartPixel_interp[MyHeight-2][x];
temp2=ori_QuartPixel_interp[MyHeight-2][y];
ori_QuartPixel_interp[MyHeight-1][4*j]=temp+(temp1+temp2)/2;
}
for(i=0;i<MyHeight;i++)
for(j=0;j<MyWidth;j++)
{origin[i][j]=ori_QuartPixel_interp[i][4*j];}
}
void CTestView::Sinc_coefficitns()
{
int i,j;
int l;
Sinc_interpolation=new double* [3];
for(i=0;i<3;i++)
Sinc_interpolation[i]=new double [M_sample*2];
for(j=0;j<M_sample*2;j++)
{
l=j-(M_sample-1);
Sinc_interpolation[0][j]=sin(PI*(0.25-l))/(PI*(0.25-l));
Sinc_interpolation[1][j]=sin(PI*(0.5-l))/(PI*(0.5-l));
Sinc_interpolation[2][j]=sin(PI*(0.75-l))/(PI*(0.75-l));
}
}
void CTestView::OnINTQuincunxWAVELETFORWARD()
{
// 獲取文檔
CTestDoc* pDoc = GetDocument();
// 指向DIB的指針
LPSTR lpDIB;
// 指向DIB象素指針
LPSTR lpDIBBits;
// 鎖定DIB
lpDIB = (LPSTR) ::GlobalLock((HGLOBAL) pDoc->GetHDIB());
// 找到DIB圖像象素起始位置
lpDIBBits = ::FindDIBBits(lpDIB);
// 判斷是否是8-bpp位圖（這里為了方便，只處理8-bpp位圖，其它的可以類推）
if (::DIBNumColors(lpDIB) != 256)
{
// 提示用戶
MessageBox("目前只支持查看256色位圖灰度直方圖！", "系統提示" , MB_ICONINFORMATION | MB_OK);
// 解除鎖定
::GlobalUnlock((HGLOBAL) pDoc->GetHDIB());
// 返回
return;
}
int lHeight,lWidth;
int i,j,m;
lHeight=::DIBHeight(lpDIB);
lWidth=::DIBWidth(lpDIB);
DWORD dwLineBytes=WIDTHBYTES(lWidth*8);
MyHeight=lHeight;
MyWidth=lWidth;
M=8;
N=8;
Dialog1 tt;
if(tt.DoModal()==IDOK)
{ levels=tt.m_Levels;
}
free(tt);
//申請空間，初始化原始圖像數組
origin=new int* [MyHeight];
for(i=0;i<MyHeight;i++)
origin[i]=new int [MyWidth];
for(i=0;i<MyHeight;i++)
for(j=0;j<MyWidth;j++)
{
//LPSTR lptemp=lpDIBBits+(lHeight-i-1)*dwLineBytes+j; //-128---+127
//origin[i][j]=(int)(*lptemp);
BYTE *lptemp=(BYTE *)lpDIBBits+(lHeight-1-i)*dwLineBytes+j; //0--255
origin[i][j]=(int)(*lptemp);
}
total_preDrection=new int* [levels];
total_preDrection1=new int* [levels];
for(m=0;m<levels;m++)
{
//預測
quincunx_Prediction();
//更新
quincunx_Update();
//在一次對L和H預測
quincunx_Prediction1();
//在一次對L和H更新
quincunx_Update1();
//保存預測的方向信息
total_preDrection[m]=new int [(MyHeight/M)*(MyWidth/N)];
for(i=0;i<MyHeight/M;i++)
for(j=0;j<MyWidth/N;j++)
{
total_preDrection[m][i*(MyWidth/N)+j]=Predict_direc[i][j];
}
total_preDrection1[m]=new int [(MyHeight/M)*(MyWidth/N)];
for(i=0;i<MyHeight/M;i++)
for(j=0;j<MyWidth/N;j++)
{
total_preDrection1[m][i*(MyWidth/N)+j]=Predict_direc1[i][j];
}
//下一級分解初始化
MyHeight=MyHeight/2;
MyWidth=MyWidth/2;
if(Predict_direc!=NULL)
{delete Predict_direc;
Predict_direc=NULL;}
if(Block_energy!=NULL)
{delete Block_energy;
Block_energy=NULL;}
if(Predict_direc1!=NULL)
{delete Predict_direc1;
Predict_direc1=NULL;}
if(Block_energy1!=NULL)
{delete Block_energy1;
Block_energy1=NULL;}
}
MyHeight=MyHeight<<levels;
MyWidth=MyWidth<<levels;
// 更新視圖
pDoc->UpdateAllViews(NULL);
// 解除鎖定
::GlobalUnlock((HGLOBAL) pDoc->GetHDIB());
//下面是將得到的系數寫入文件，保存起來，文件保存在工程所在的目錄下
CFile CoeMatrix;
CFileException fe;
CoeMatrix.Open("coefficients.txt",CFile::modeCreate | CFile::modeReadWrite | CFile::shareExclusive,&fe);
CString str;
int E=0;
int max=-10000;
int min=100000;
for(i=MyHeight/2;i<MyHeight;i++)
for(j=0;j<MyWidth;j++)
{
if(max<origin[i][j])
{max=origin[i][j];}
if(min>origin[i][j])
{min=origin[i][j];}
E+=abs(origin[i][j]);
}
str.Format("Max:%d,Min:%d,Energy:%d",max,min,E);
// AfxMessageBox(str);
for(i=0;i<lHeight;i++)
for(j=0;j<lWidth;j++)
{
BYTE* lptemp=(BYTE*)lpDIBBits+(lHeight-i-1)*dwLineBytes+j;
(*lptemp)=(BYTE)origin[i][j];
}
}
void CTestView::OnINTQuincunxWAVELETBACKWARD()
{
// 獲取文檔
CTestDoc* pDoc = GetDocument();
// 指向DIB的指針
LPSTR lpDIB;
// 指向DIB象素指針
LPSTR lpDIBBits;
// 鎖定DIB
lpDIB = (LPSTR) ::GlobalLock((HGLOBAL) pDoc->GetHDIB());
// 找到DIB圖像象素起始位置
lpDIBBits = ::FindDIBBits(lpDIB);
// 判斷是否是8-bpp位圖（這里為了方便，只處理8-bpp位圖，其它的可以類推）
if (::DIBNumColors(lpDIB) != 256)
{
// 提示用戶
MessageBox("目前只支持查看256色位圖灰度直方圖！", "系統提示" , MB_ICONINFORMATION | MB_OK);
// 解除鎖定
::GlobalUnlock((HGLOBAL) pDoc->GetHDIB());
// 返回
return;
}
int lHeight,lWidth;
int i,j,m;
lHeight=::DIBHeight(lpDIB);
lWidth=::DIBWidth(lpDIB);
DWORD dwLineBytes=WIDTHBYTES(lWidth*8);
MyHeight=lHeight;
MyWidth=lWidth;
for(m=levels;m>0;m--)
{
MyHeight=MyHeight>>(m-1);
MyWidth=MyWidth>>(m-1);
/* //預測
quincunx_Prediction();
//更新
quincunx_Update();
//在一次對L和H預測
quincunx_Prediction1();
//在一次對L和H更新
quincunx_Update1(); */
//對L和H逆更新
inverse_quincunx_Update1(m);
//對L和H逆預測
inverse_quincunx_Prediction1(m);
//逆更新
inverse_quincunx_Update(m);
//逆預測
inverse_quincunx_Prediction(m);
//下一級逆變換初始化
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -