dibview.cpp

清华大学出版社出版的《数字图像处理实训教程》一书以提高动手能力为目标

	int buf,buf1,buf2;
	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方向
			buf1=(int)image_in[i-1][j+1]+2*(int)image_in[i][j+1]+(int)image_in[i+1][j+1];
			buf1=buf1-(int)image_in[i-1][j-1]-2*(int)image_in[i][j-1]-(int)image_in[i+1][j-1];
			//y方向
			buf2=(int)image_in[i-1][j-1]+2*(int)image_in[i-1][j]+(int)image_in[i-1][j+1];
			buf2=buf2-(int)image_in[i+1][j-1]-2*(int)image_in[i+1][j]-(int)image_in[i+1][j+1];
			buf=abs(buf1)+abs(buf2);
			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_sobel.raw");
}

void lapla1(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=4*(int)image_in[i][j]-(int)image_in[i-1][j]-(int)image_in[i][j-1]
				-(int)image_in[i][j+1]-(int)image_in[i+1][j];
			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_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_INCOMP, OnIncomp)
	//}}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::OnIncomp() 
{
 CDibDoc *pDoc=GetDocument();
  HDIB m_hDIB=pDoc->GetHDIB();
  if(m_hDIB==NULL)

  {
	  AfxMessageBox("先打开一幅灰度图");
	  return;
  }

  BYTE *lp; 			
  lp=(BYTE *)::GlobalLock(m_hDIB);
  lp=lp+40;
  int rad=128;//
  long int hist[256];
  long int newhist[256];
  AfxMessageBox("rad为128，可以修改源代码改变rad值");
  GetHist(m_hDIB,hist);//求图像直方图
  qiu(hist,newhist,rad); //求目标直方图newhist
  double old[256];
  double now[256];
  long int totalo;
  long int totaln;
  long int opo,opn;
	  
  totalo=0;
  int i,j;
  for(i=0;i<256;i++)
	 totalo+=newhist[i];
  opo=0;
  for(i=0;i<256;i++)
	   {
		   opo+=newhist[i];
		   old[i]=opo*1.0/totalo;  
     //old[i]中保存目标图像灰度级从0到i的像素比例
	   }  
  totaln=0;
  for(i=0;i<256;i++)
	totaln+=hist[i];
  opn=0;
  for(i=0;i<256;i++)
	   {
		   opn+=hist[i];
		   now[i]=opn*1.0/totaln;
       //now[i]中保存源图像灰度级从0到i的像素比例
	   }
	   
  for( i=0;i<256;i++)
	   {
	   for( j=0;j<256;j++)
	   {
		   if(now[i]<=old[j])
		   {
		   
			   lp[i*4]=lp[i*4+1]=lp[i*4+2]   //更新调色板
				   =j;break;			  

		   }
	   }
	   
	   }

 Invalidate(FALSE);	
}