📄 dibview.cpp
字号:
image_out[i][j]=(BYTE)0;
}
SaveRawFile(image_out,w,h,"E:\\peacock_lapla1.raw");
}
void lapla2(BYTE *image_in[2000],int w,int h,BYTE *image_out[2000])
{
int buf;
int i,j;
for(i=0;i<h;i++)
for(j=0;j<w;j++)
{
if(i-1>=0&&j-1>=0&&i+1<h&&j+1<w)
{
//x方向
buf=8*(int)image_in[i][j]-(int)image_in[i-1][j-1]-(int)image_in[i-1][j]-(int)image_in[i-1][j+1]
-(int)image_in[i][j-1]-(int)image_in[i][j+1]-(int)image_in[i+1][j-1]
-(int)image_in[i+1][j]-(int)image_in[i+1][j+1];
if(buf>255)buf=255;
if(buf<0)buf=0;
image_out[i][j]=(BYTE)buf;
}
else
image_out[i][j]=(BYTE)0;
}
SaveRawFile(image_out,w,h,"E:\\peacock_lapla2.raw");
}
void zzlb(BYTE *image_in[2000],int w,int h,BYTE *image_out[2000])
{
int value[9];
int i,j;
for(i=1;i<h-1;i++)
for(j=1;j<w-1;j++)
{
value[0]=image_in[i-1][j-1];
value[1]=image_in[i-1][j];
value[2]=image_in[i-1][j+1];
value[3]=image_in[i][j-1];
value[4]=image_in[i][j];
value[5]=image_in[i][j+1];
value[6]=image_in[i+1][j-1];
value[7]=image_in[i+1][j];
value[8]=image_in[i+1][j+1];
image_out[i][j]=medianvalue(value);
}
SaveRawFile(image_out,w,h,"E:\\peacock_zzlb.raw");
}
int medianvalue(int *value)
{
int i,j,temp;
for(i=0;i<=7;i++)
for(j=0;j<=7-i;j++)
{
if(value[j]>value[j+1])
{
temp=value[j+1];
value[j+1]=value[j];
value[j]=temp;
}
}
return value[4];
}
IMPLEMENT_DYNCREATE(CDibView, CScrollView)
BEGIN_MESSAGE_MAP(CDibView, CScrollView)
//{{AFX_MSG_MAP(CDibView)
ON_COMMAND(ID_EDIT_COPY, OnEditCopy)
ON_UPDATE_COMMAND_UI(ID_EDIT_COPY, OnUpdateEditCopy)
ON_COMMAND(ID_EDIT_PASTE, OnEditPaste)
ON_UPDATE_COMMAND_UI(ID_EDIT_PASTE, OnUpdateEditPaste)
ON_MESSAGE(WM_DOREALIZE, OnDoRealize)
ON_COMMAND(ID_EDGE, OnEdge)
ON_COMMAND(ID_EN, OnEn)
ON_COMMAND(ID_FAST_EN, OnFastEn)
//}}AFX_MSG_MAP
// Standard printing commands
ON_COMMAND(ID_FILE_PRINT, CScrollView::OnFilePrint)
ON_COMMAND(ID_FILE_PRINT_PREVIEW, CScrollView::OnFilePrintPreview)
END_MESSAGE_MAP()
/////////////////////////////////////////////////////////////////////////////
// CDibView construction/destruction
CDibView::CDibView()
{
}
CDibView::~CDibView()
{
}
/////////////////////////////////////////////////////////////////////////////
// CDibView drawing
void CDibView::OnDraw(CDC* pDC)
{
CDibDoc* pDoc = GetDocument();
HDIB hDIB = pDoc->GetHDIB();
if (hDIB != NULL)
{
LPSTR lpDIB = (LPSTR) ::GlobalLock((HGLOBAL) hDIB);
int cxDIB = (int) ::DIBWidth(lpDIB); // Size of DIB - x
int cyDIB = (int) ::DIBHeight(lpDIB); // Size of DIB - y
::GlobalUnlock((HGLOBAL) hDIB);
CRect rcDIB;
rcDIB.top = rcDIB.left = 0;
rcDIB.right = cxDIB;
rcDIB.bottom = cyDIB;
CRect rcDest;
if (pDC->IsPrinting()) // printer DC
{
// get size of printer page (in pixels)
int cxPage = pDC->GetDeviceCaps(HORZRES);
int cyPage = pDC->GetDeviceCaps(VERTRES);
// get printer pixels per inch
int cxInch = pDC->GetDeviceCaps(LOGPIXELSX);
int cyInch = pDC->GetDeviceCaps(LOGPIXELSY);
//
// Best Fit case -- create a rectangle which preserves
// the DIB's aspect ratio, and fills the page horizontally.
//
// The formula in the "->bottom" field below calculates the Y
// position of the printed bitmap, based on the size of the
// bitmap, the width of the page, and the relative size of
// a printed pixel (cyInch / cxInch).
//
rcDest.top = rcDest.left = 0;
rcDest.bottom = (int)(((double)cyDIB * cxPage * cyInch)
/ ((double)cxDIB * cxInch));
rcDest.right = cxPage;
}
else // not printer DC
{
rcDest = rcDIB;
}
::PaintDIB(pDC->m_hDC, &rcDest, pDoc->GetHDIB(),
&rcDIB, pDoc->GetDocPalette());
}
}
/////////////////////////////////////////////////////////////////////////////
// CDibView printing
BOOL CDibView::OnPreparePrinting(CPrintInfo* pInfo)
{
// default preparation
return DoPreparePrinting(pInfo);
}
/////////////////////////////////////////////////////////////////////////////
// CDibView commands
LRESULT CDibView::OnDoRealize(WPARAM wParam, LPARAM)
{
ASSERT(wParam != NULL);
CDibDoc* pDoc = GetDocument();
if (pDoc->GetHDIB() == NULL)
return 0L; // must be a new document
CPalette* pPal = pDoc->GetDocPalette();
if (pPal != NULL)
{
CMainFrame* pAppFrame = (CMainFrame*) AfxGetApp()->m_pMainWnd;
ASSERT_KINDOF(CMainFrame, pAppFrame);
CClientDC appDC(pAppFrame);
// All views but one should be a background palette.
// wParam contains a handle to the active view, so the SelectPalette
// bForceBackground flag is FALSE only if wParam == m_hWnd (this view)
CPalette* oldPalette = appDC.SelectPalette(pPal, ((HWND)wParam) != m_hWnd);
if (oldPalette != NULL)
{
UINT nColorsChanged = appDC.RealizePalette();
if (nColorsChanged > 0)
pDoc->UpdateAllViews(NULL);
appDC.SelectPalette(oldPalette, TRUE);
}
else
{
TRACE0("\tSelectPalette failed in CDibView::OnPaletteChanged\n");
}
}
return 0L;
}
void CDibView::OnInitialUpdate()
{
CScrollView::OnInitialUpdate();
ASSERT(GetDocument() != NULL);
SetScrollSizes(MM_TEXT, GetDocument()->GetDocSize());
}
void CDibView::OnActivateView(BOOL bActivate, CView* pActivateView,
CView* pDeactiveView)
{
CScrollView::OnActivateView(bActivate, pActivateView, pDeactiveView);
if (bActivate)
{
ASSERT(pActivateView == this);
OnDoRealize((WPARAM)m_hWnd, 0); // same as SendMessage(WM_DOREALIZE);
}
}
void CDibView::OnEditCopy()
{
CDibDoc* pDoc = GetDocument();
// Clean clipboard of contents, and copy the DIB.
if (OpenClipboard())
{
BeginWaitCursor();
EmptyClipboard();
SetClipboardData (CF_DIB, CopyHandle((HANDLE) pDoc->GetHDIB()) );
CloseClipboard();
EndWaitCursor();
}
}
void CDibView::OnUpdateEditCopy(CCmdUI* pCmdUI)
{
pCmdUI->Enable(GetDocument()->GetHDIB() != NULL);
}
void CDibView::OnEditPaste()
{
HDIB hNewDIB = NULL;
if (OpenClipboard())
{
BeginWaitCursor();
hNewDIB = (HDIB) CopyHandle(::GetClipboardData(CF_DIB));
CloseClipboard();
if (hNewDIB != NULL)
{
CDibDoc* pDoc = GetDocument();
pDoc->ReplaceHDIB(hNewDIB); // and free the old DIB
pDoc->InitDIBData(); // set up new size & palette
pDoc->SetModifiedFlag(TRUE);
SetScrollSizes(MM_TEXT, pDoc->GetDocSize());
OnDoRealize((WPARAM)m_hWnd,0); // realize the new palette
pDoc->UpdateAllViews(NULL);
}
EndWaitCursor();
}
}
void CDibView::OnUpdateEditPaste(CCmdUI* pCmdUI)
{
pCmdUI->Enable(::IsClipboardFormatAvailable(CF_DIB));
}
void MobanEdge(HDIB hDIB, int mm[3][3])
{
BYTE *image[5000];
GetImage(hDIB,image);
HDIB newD=(HDIB)CopyHandle(hDIB);
BYTE *image1[5000];
LPSTR lp=(LPSTR)::GlobalLock(hDIB);
GetImage(newD,image1);
for(int i=0;i<DIBHeight(lp);i++)
for(int j=0;j<DIBWidth(lp);j++)
{
if(i<1||i>DIBHeight(lp)-2||
j<1||j>DIBWidth(lp)-2)
continue;
int temp=0;
for(int ii=0;ii<3;ii++)
for(int jj=0;jj<3;jj++)
temp+=image1[i+ii-1][j+jj-1]*mm[ii][jj];
int yuzhi=100;
image[i][j]=abs(temp)>100?255:0;
}
}
void CDibView::OnEdge()
{
CDibDoc *pDoc=GetDocument();
HDIB hDIB=pDoc->GetHDIB();
int mm[3][3];
for(int i=0;i<3;i++)
for(int j=0;j<3;j++)
mm[i][j]=0;
mm[0][0]=-1;
mm[0][2]=1;
mm[1][0]=-2;
mm[1][2]=2;
mm[2][0]=-1;
mm[2][2]=1;
MobanEdge(hDIB, mm);
Invalidate();
}
void CDibView::OnEn()
{
HDIB hDIB=GetDocument()->GetHDIB(); //获取图像句柄
BYTE *image[2000];
int imageWidth,imageHeight;
BYTE *lp;
LPBITMAPINFOHEADER lh;
if(!hDIB) return;
lp=(BYTE *)::GlobalLock(hDIB);
lh=(LPBITMAPINFOHEADER)lp;
imageWidth=lh->biWidth;
imageHeight=lh->biHeight;
int biwidth;
biwidth=lh->biWidth;
if(biwidth%4!=0)
biwidth+=(4-biwidth%4);
if (lh->biBitCount==8)
{
lp=lp+40+256*4;
int i;
for(i=0;i<lh->biHeight;i++)
image[i]=lp+(lh->biHeight-i-1)*biwidth;
}
else
{
::GlobalUnlock(hDIB);
return ;
}
//以上代码将hDIB转化为image数组
double *imageo[2000];
double *imagen[2000];
for(int i=0;i<imageHeight;i++)
{
imagen[i]=new double[imageWidth];
imageo[i]=new double[imageWidth];
for(int j=0;j<imageWidth;j++)
{
imagen[i][j]=imageo[i][j]=(double)image[i][j]; //初始化
}
}
int adj=3;//should be odd
for ( int num=0;num<2;num++)
{
int K=GetK(imageo,imageWidth,imageHeight); //得到K的值
K=K*4;
for(int ii=(int)(adj*1.5);ii<(int)(imageHeight-adj*1.5);ii++)
for(int jj=(int)(adj*1.5);jj<(int)(imageWidth-adj*1.5);jj++)
{
double doo=GetAdj(imageo,ii,jj,adj,'o');
double dn=GetAdj(imageo,ii,jj,adj,'n')-doo;
double ds=GetAdj(imageo,ii,jj,adj,'s')-doo;
double de=GetAdj(imageo,ii,jj,adj,'e')-doo;
double dw=GetAdj(imageo,ii,jj,adj,'w')-doo;
double dm=GetAdj(imageo,ii,jj,adj,'m')-doo;
double dl=GetAdj(imageo,ii,jj,adj,'l')-doo;
double dr=GetAdj(imageo,ii,jj,adj,'r')-doo;
double dt=GetAdj(imageo,ii,jj,adj,'t')-doo;
double cn=g(dn,K,15);
double cs=g(ds,K,15);
double ce=g(de,K,15);
double cw=g(dw,K,15);
double cm=g(dm,K,15);
double ct=g(dt,K,15);
double cr=g(dr,K,15);
double cl=g(dl,K,15);
imagen[ii][jj]=imageo[ii][jj]-0.1*(cn*dn+cs*ds+ce*de+cw*dw+
cm*dm+cl*dl+cr*dr+ct*dt); //迭代
}
for(int k=0;k<imageHeight;k++)
for(int l=0;l<imageWidth;l++)
{
imageo[k][l]=imagen[k][l];
if(imageo[k][l]<0) imageo[k][l]=0;
if(imageo[k][l]>255)imageo[k][l]=255; //溢界处理
}
}
for(int k=0;k<imageHeight;k++)
for(int l=0;l<imageWidth;l++)
{
double temp=imagen[k][l]+0.5;
if(temp>255)
{
temp=255;
}
if(temp<0)
{
temp=0;}
image[k][l]=(BYTE)temp;
}
Invalidate(FALSE); //通知窗口重绘
}
void CDibView::OnFastEn()
{
int r=75;
AfxMessageBox("r为75，可以在源程序中更改");
BYTE *image[2000];
int imageWidth,imageHeight;
BYTE *lp;
LPBITMAPINFOHEADER lh;
HDIB hDIB=GetDocument()->GetHDIB();
if(!hDIB) return;
lp=(BYTE *)::GlobalLock(hDIB);
lh=(LPBITMAPINFOHEADER)lp;
imageWidth=lh->biWidth;
imageHeight=lh->biHeight;
int biwidth;
biwidth=lh->biWidth;
if(biwidth%4!=0)
biwidth+=(4-biwidth%4);
if (lh->biBitCount==8)
{
lp=(BYTE *)::FindDIBBits((LPSTR)lp);
int i;
for(i=0;i<lh->biHeight;i++)
image[i]=lp+(lh->biHeight-i-1)*biwidth;
}
else
{::GlobalUnlock(hDIB);
return ;
}
//以上实现hDIB到image的转换
r=((int)r)/2*2+1;
//保证r为奇数
BYTE *total[2000];
long int *first;
for(int k=0;k<imageHeight-r;k++)
{
total[k]=(BYTE *)
new BYTE[imageWidth-r];
}
first=(long int *)new long int[imageWidth-r];
long int tt;
tt=0;
for(int m=0;m<r;m++)
for(int n=0;n<r;n++)
tt+=image[m][n]; //求邻域内灰度和
first[0]=tt;
for(int l=1;l<imageWidth-r;l++)
{
for(int pp=0;pp<r;pp++)
{
tt+=(long int)(image[pp][(int)r+l-1])
-( long int)(image[pp][l-1]);
}
first[l]=tt;
}
//由first[0]求first[l], 即书中14.6
for( k=0;k<imageWidth-r;k++)
{
tt=first[k]; total[0][k]=tt/(r*r);
for(int l=1;l<imageHeight-r;l++)
{
for( int jj=0;jj<r;jj++)
{
tt+=(long int)(image[(int)r+l-1][jj+k])
-(long int)(image[l-1][jj+k]);
}
total[l][k]=tt/(r*r);
}
}
int temp;
int K;
K=(int)r/2;
int g,h;
for(k=0;k<imageHeight;k++)
for(l=0;l<imageWidth;l++)
{
g=k;
h=l;
g=g-K;
h=l-K;
if(g<0) g=0;
if(h<0) h=0;
if(g>=imageHeight-r) g=imageHeight-r-1;
if(h>=imageWidth-r) h=imageWidth-r-1;
temp=total[g][h];
temp=image[k][l]-temp+128;
if(temp<0) temp=0;
if(temp>255) temp=255;
image[k][l]=temp;
} //边界处理
delete []first;
for(l=0;l<imageHeight-r;l++)
if(total[l]) delete ((total[l]));
::GlobalUnlock(hDIB);
Invalidate(FALSE);
}
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -