📄 dibview.cpp

📁 清华大学出版社出版的《数字图像处理实训教程》一书以提高动手能力为目标
💻 CPP
📖 第 1 页 / 共 2 页
字号:
上一页 12
		for(r=0;r<rmax;r++)
			hough[k][r]=0;
	for(i=0;i<h;i++)                  //对原图作hough变换，同算法一
	for(j=0;j<w;j++)
	{
		if(image_in[i][j]==255)
		{
			for(k=0;k<180;k++)
			{
			r=(int)(j*cos(k*t)+i*sin(k*t));
			if(r<0)r+=w+h+1;
			hough[k][r]++;
			}
		}
	} 

	int max=0,max_k,max_r;                //寻找最长的一条直线
	for(k=0;k<180;k++)
	for(r=0;r<rmax;r++)
	{
		if(hough[k][r]>max)
		{
			max=hough[k][r];
			max_k=k;
			max_r=r;
		}
	}
	if(max<houghthresh)break;                //当最大值小于设定的临界值时，循环结束
//	makeline1(image_in,w,h,t,max_k,max_r,image_out);
	makeline2(image_in,w,h,t,max_k,max_r,image_out);
	}
	SaveRawFile(image_out,w,h,"E:\\peacock_hough.raw");
}


//绘制直线算法一
void makeline1(BYTE *image_in[2000],int w,int h,double t,int k,int r,BYTE *image_out[2000])
{
	int i,j,temp;
	for(i=0;i<h;i++)   //重新对数组上的白点作hough变换，如果发现某点落在(k,r)的区域内，在原图与新图中均做标记。
		for(j=0;j<w;j++)
		{
			if(image_in[i][j]==255)
			{
				temp=(int)(j*cos(k*t)+i*sin(k*t));
				if(temp<0)temp+=w+h+1;
				if(temp==r)
				{image_in[i][j]=250;image_out[i][j]=255;}
			}
		}
}


//绘制图像算法二
void makeline2(BYTE *image_in[2000],int w,int h,double t,int k,int r,BYTE *image_out[2000])
{
	int i,j,temp;
	for(i=0;i<h;i++)   //对原图重新作hough变换，标记已经提取到的直线(k,r)上的白点
		for(j=0;j<w;j++)
		{
			if(image_in[i][j]==255)
			{
				temp=(int)(j*cos(k*t)+i*sin(k*t));
				if(temp<0)temp+=w+h+1;
				if(temp==r)
				{image_in[i][j]=250;}
			}
		}
	for(j=0;j<w;j++)   //在输出数组image_out[][]中通过逆函数绘制直线
	{
		if(k!=0)
		{
		if(k>90&&k<180&&r>h)r-=w+h+1;
		i=(int)(r/sin(k*t)-j/tan(k*t));	
		if(i>=0&&i<h)image_out[i][j]=255;
		}
	}
}

//以下为二值化，梯度锐化，中值滤波等函数，不做注释！
void diedaiyz(BYTE *image[2000],int w,int h)
{
	int i,j;
	int count[256];
	for(i=0;i<256;i++)
		count[i]=0;
	for(i=0;i<h;i++)
		for(j=0;j<w;j++)
			count[image[i][j]]++;
	int T=100;
	while(1)
	{
		int temp=T;
		int t1=0,t2=0,t3=0,t4=0;
		for(i=0;i<=temp;i++)
		{t1+=count[i];t3+=count[i]*i;}
		for(i=temp+1;i<256;i++)
		{t2+=count[i];t4+=count[i]*i;}
		T=(int)((t3/t1+t4/t2)/2);
		if(T==temp)break;
	}
	for(i=0;i<h;i++)
		for(j=0;j<w;j++)
			if(image[i][j]>=T)image[i][j]=255;
			else image[i][j]=0;
	SaveRawFile(image,w,h,"E:\\peacock_diedaiyz.raw");
}

void sanfenyz(BYTE *image[2000],int w,int h)
{
	int i,j;
	int fmax=0,fmin=255;
	for(i=0;i<h;i++)
		for(j=0;j<w;j++)
		{
			if(image[i][j]>fmax)fmax=image[i][j];
			if(image[i][j]<fmin)fmin=image[i][j];
		}
	int T=fmax-(fmax-fmin)/3;
	for(i=0;i<h;i++)
		for(j=0;j<w;j++)
			if(image[i][j]>=T)image[i][j]=255;
			else image[i][j]=0;
	SaveRawFile(image,w,h,"E:\\peacock_sanfenyz.raw");
}

void sobel(BYTE *image_in[2000],int w,int h,BYTE *image_out[2000])
{
	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_90, On90)
	ON_COMMAND(ID_0, On0)
	ON_COMMAND(ID_HOR, OnHor)
	ON_COMMAND(ID_180, On180)
	ON_COMMAND(ID_VERTI, OnVerti)
	ON_COMMAND(ID_270, On270)
	ON_COMMAND(ID_RE, OnRe)
	//}}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::On90() 
{
GetDocument()->ReplaceHDIB(
		Rotate(GetDocument()->GetHDIB(),90,SLOW));
Invalidate(TRUE);
}

void CDibView::On0() 
{
GetDocument()->ReplaceHDIB(
		Rotate(GetDocument()->GetHDIB(),0,SLOW));
Invalidate(TRUE);
}

void CDibView::OnHor() 
{
GetDocument()->ReplaceHDIB(
		Rotate(GetDocument()->GetHDIB(),HOR,SLOW));
Invalidate(TRUE);
}

void CDibView::On180() 
{
GetDocument()->ReplaceHDIB(
		Rotate(GetDocument()->GetHDIB(),180,SLOW));
Invalidate(TRUE);
}

void CDibView::OnVerti() 
{
GetDocument()->ReplaceHDIB(
		Rotate(GetDocument()->GetHDIB(),VERTI,SLOW));
Invalidate(TRUE);
}

void CDibView::On270() 
{
GetDocument()->ReplaceHDIB(
		Rotate(GetDocument()->GetHDIB(),270,SLOW));
Invalidate(TRUE);
}

void CDibView::OnRe() 
{
GetDocument()->ReplaceHDIB(
		Rotate(GetDocument()->GetHDIB(),27,MIDFAST));
Invalidate(TRUE);
}
上一页 12
⌨️ 快捷键说明

复制代码 Ctrl + C
搜索代码 Ctrl + F
全屏模式 F11
切换主题 Ctrl + Shift + D
显示快捷键 ?
增大字号 Ctrl + =
减小字号 Ctrl + -