数字图像处理实验模块.doc

数字图像处理实验模块/ jihanfeiDoc.cpp : implementation of the CJihanfeiDoc class/#include stdafx.h#include jihanfei.h#include math.h#include jihanfeiDoc.h#ifdef _DEBUG#define new DEBUG_NEW#undef THIS_FILEstatic char THIS_FILE = _FILE_;#endif/ CJihanfeiDocIMPLEMENT_DYNCREATE(CJihanfeiDoc, CDocument)BEGIN_MESSAGE_MAP(CJihanfeiDoc, CDocument)/AFX_MSG_MAP(CJihanfeiDoc)ON_COMMAND(ID_fei, Onfei)ON_COMMAND(ID_niu, Onniu)ON_COMMAND(ID_ji, Onji)ON_COMMAND(ID_a, Ona)ON_COMMAND(ID_b, Onb)ON_COMMAND(ID_zq, Onzq)ON_COMMAND(ID_count, Oncount)ON_COMMAND(ID_neibor, Onneibor)ON_COMMAND(ID_neibor8, Onneibor8)ON_COMMAND(ID_mneibor, Onmneibor)ON_COMMAND(ID_gaosi, Ongaosi)ON_COMMAND(ID_center, Oncenter)ON_COMMAND(ID_layer, Onlayer)ON_COMMAND(ID_one, Onone)ON_COMMAND(ID_two, Ontwo)ON_COMMAND(ID_Robert, OnRobert)ON_COMMAND(ID_Sobel, OnSobel)ON_COMMAND(ID_zhijie, Onzhijie)ON_COMMAND(ID_erjie, Onerjie)ON_COMMAND(ID_average, Onaverage)ON_COMMAND(ID_Hough, OnHough)ON_COMMAND(ID_gaositidu, Ongaositidu)ON_COMMAND(ID_trail, Ontrail)/AFX_MSG_MAPEND_MESSAGE_MAP()/ CJihanfeiDoc construction/destructionCJihanfeiDoc:CJihanfeiDoc()/ TODO: add one-time construction code hereCJihanfeiDoc:CJihanfeiDoc()BOOL CJihanfeiDoc:OnNewDocument()if (!CDocument:OnNewDocument()return FALSE;/ TODO: add reinitialization code here/ (SDI documents will reuse this document)return TRUE;/ CJihanfeiDoc serializationvoid CJihanfeiDoc:Serialize(CArchive& ar)if (ar.IsStoring()/ TODO: add storing code hereelse/ TODO: add loading code here/ CJihanfeiDoc diagnostics#ifdef _DEBUGvoid CJihanfeiDoc:AssertValid() constCDocument:AssertValid();void CJihanfeiDoc:Dump(CDumpContext& dc) constCDocument:Dump(dc);#endif /_DEBUG/ CJihanfeiDoc commandsBOOL CJihanfeiDoc:OnOpenDocument(LPCTSTR lpszPathName) if (!CDocument:OnOpenDocument(lpszPathName)return FALSE;if (!m_Dib.Read_8BitsGray_Image(lpszPathName)/ 判断读取文件是否成功CString strMsg;/ 失败，可能非BMP格式strMsg = 读取图像时出错！可能是不支持该类型的图像文件！;/ 提示出错MessageBox(NULL, strMsg, 系统提示, MB_ICONINFORMATION | MB_OK);return FALSE;/ 返回FALSEelsem_Dib.MakePalette();SetPathName(lpszPathName);/ 设置文件名称SetModifiedFlag(FALSE);/ 初始化脏标记为FALSEreturn TRUE;/ 返回TRUEreturn TRUE;BOOL CJihanfeiDoc:CanCloseFrame(CFrameWnd* pFrame) / TODO: Add your specialized code here and/or call the base classreturn CDocument:CanCloseFrame(pFrame);void CJihanfeiDoc:Onfei() / TODO: Add your command handler code hereLPSTR lpStartBit = m_Dib.GetPixelBit();/ 象素起始位置指针if (lpStartBit = NULL)/ 找不到象素起始位置指针return ;/ 返回FALSEint iWidth = m_Dib.GetWidth();/ 图象宽度int iHeight = m_Dib.GetHeight();/ 图象高度int iBitsWidth = m_Dib.GetBitsWidth();/ 图象中每行的字节数float temp=0.0f;/ 临时变量double dtemp=0;/ 临时变量unsigned char * lpPixel;/ 指向源DIB图像象素位置的指针int i;/ 循环变量int j;float fSlope=1.4f;float fIntercept=20.0f;for (i=0;iiHeight;i+)for (j=0;j128)*lpPixel=255;else*lpPixel=0;UpdateAllViews(NULL);return;void CJihanfeiDoc:Onniu() / TODO: Add your command handler code hereLPSTR lpStartBit = m_Dib.GetPixelBit();/ 象素起始位置指针if (lpStartBit = NULL)/ 找不到象素起始位置指针return ;/ 返回FALSEint iWidth = m_Dib.GetWidth();/ 图象宽度int iHeight = m_Dib.GetHeight();/ 图象高度int iBitsWidth = m_Dib.GetBitsWidth();/ 图象中每行的字节数float temp=0.0f;/ 临时变量double dtemp=0;/ 临时变量unsigned char * lpPixel;/ 指向源DIB图像象素位置的指针int i;/ 循环变量int j;float fSlope=1.4f;float fIntercept=20.0f;for (i=0;iiHeight;i+)for (j=0;jiWidth;j+)/读取每象素存储位置lpPixel=(unsigned char*)lpStartBit + (iHeight-1-i) * iBitsWidth + j;/ 以下位置填写代码实现线性变换*lpPixel=255-*lpPixel;UpdateAllViews(NULL);return;void CJihanfeiDoc:Onji() / TODO: Add your command handler code hereLPSTR lpStartBit = m_Dib.GetPixelBit();/ 象素起始位置指针if (lpStartBit = NULL)/ 找不到象素起始位置指针return ;/ 返回FALSEint iWidth = m_Dib.GetWidth();/ 图象宽度int iHeight = m_Dib.GetHeight();/ 图象高度int iBitsWidth = m_Dib.GetBitsWidth();/ 图象中每行的字节数float temp=0.0f;/ 临时变量double dtemp=0;/ 临时变量unsigned char * lpPixel;/ 指向源DIB图像象素位置的指针int i;/ 循环变量int j;float fSlope=1.4f;float fIntercept=20.0f;long plcount256;for (i=0;iiHeight;i+)for (j=0;jiWidth;j+)/读取每象素存储位置lpPixel=(unsigned char*)lpStartBit + (iHeight-1-i) * iBitsWidth + j;plcount*(lpPixel)+;/ 以下位置填写代码实现线性变换UpdateAllViews(NULL);return;void CJihanfeiDoc:Ona() / TODO: Add your command handler code hereLPSTR lpStartBit = m_Dib.GetPixelBit();/ 象素起始位置指针if (lpStartBit = NULL)/ 找不到象素起始位置指针return ;/ 返回FALSEint iWidth = m_Dib.GetWidth();/ 图象宽度int iHeight = m_Dib.GetHeight();/ 图象高度int iBitsWidth = m_Dib.GetBitsWidth();/ 图象中每行的字节数float temp=0.0f;/ 临时变量float dtemp=0;/ 临时变量unsigned char * lpPixel;/ 指向源DIB图像象素位置的指针int i;/ 循环变量int j;double fSlope=1.4;double fIntercept=20.0;for (i=0;iiHeight;i+)for (j=0;j255.0)*lpPixel = 255;else if (temp0.0)*lpPixel = 0;else*lpPixel = (unsigned char) (temp+0.5);UpdateAllViews(NULL);return;void CJihanfeiDoc:Onb() / TODO: Add your command handler code hereLPSTR lpStartBit = m_Dib.GetPixelBit();/ 象素起始位置指针if (lpStartBit = NULL)/ 找不到象素起始位置指针return ;/ 返回FALSEint iWidth = m_Dib.GetWidth();/ 图象宽度int iHeight = m_Dib.GetHeight();/ 图象高度int iBitsWidth = m_Dib.GetBitsWidth();/ 图象中每行的字节数float temp=0.0f;/ 临时变量double dtemp=0;/ 临时变量unsigned char * lpPixel;/ 指向源DIB图像象素位置的指针int i;/ 循环变量int j;double fSlope=0.4;double fIntercept=20.0;for (i=0;iiHeight;i+)for (j=0;j255.0)*lpPixel = 255;else if (temp0.0)*lpPixel = 0;else*lpPixel = (unsigned char) (temp+0.5);UpdateAllViews(NULL);return;void CJihanfeiDoc:Onzq() / TODO: Add your command handler code hereLPSTR lpStartBit = m_Dib.GetPixelBit();/ 象素起始位置指针if (lpStartBit = NULL)/ 找不到象素起始位置指针return ;/ 返回FALSEint iWidth = m_Dib.GetWidth();/ 图象宽度int iHeight = m_Dib.GetHeight();/ 图象高度int iBitsWidth = m_Dib.GetBitsWidth();/ 图象中每行的字节数/ 临时变量unsigned char * lpPixel;/ 指向源DIB图像象素位置的指针int i;/ 循环变量int j;long plcount256;long laycount256;for (i= 0; i256; i+)/ 计数清空为0plcounti=0;for (i=0;iiHeight;i+) for (j=0;jiWidth;j+)lpPixel=(unsigned char*)lpStartBit + (iHeight-1-i) * iBitsWidth + j;plcount*lpPixel+; double temp=0;for (i= 0; i255) laycounti = 255;for (i=0;iiHeight;i+)for (j=0;jiWidth;j+)/读取每象素存储位置lpPixel=(unsigned char*)lpStartBit + (iHeight-1-i) * iBitsWidth + j;*(lpPixel)=unsigned char(laycount*(lpPixel);UpdateAllViews(NULL);return;/统计量算法增强void CJihanfeiDoc:Oncount() / TODO: Add your command handler code hereLPSTR lpStartBit = m_Dib.GetPixelBit();/ 象素起始位置指针if (lpStartBit = NULL)/ 找不到象素起始位置指针return ;/ 返回FALSEint iWidth = m_Dib.GetWidth();/ 图象宽度int iHeight = m_Dib.GetHeight();/ 图象高度int iBitsWidth = m_Dib.GetBitsWidth();/ 图象中每行的字节数/ 临时变量unsigned char * lpPixel;/ 指向源DIB图像象素位置的指针int i;/ 循环变量int j;float oave;/ 源图像灰度均值，依据实际图像赋值float ostd;/ 源图像灰度标准差，依据实际图像赋值float nave=150;/ 增强后图像的灰度均值，依据实际要求赋值float nstd=64;/ 增强后图像的灰度标准差，依据实际要求赋值float temp;float a;float b;float sum1=0,sum2=0;for (i=0;iiHeight;i+)for (j=0;jiWidth;j+)/读取每象素存储位置lpPixel=(unsigned char*)lpStartBit + (iHeight-1-i) * iBitsWidth + j;sum1+=*lpPixel;oave=sum1/(iHeight*iWidth);for (i=0;iiHeight;i+)for (j=0;jiWidth;j+)/读取每象素存储位置lpPixel=(unsigned char*)lpStartBit + (iHeight-1-i) * iBitsWidth + j;sum2+=float(fabs(*lpPixel-oave)*(fabs(*lpPixel-oave);ostd=sum2/(iHeight*iWidth);a=float(sqrt(nstd/ostd);b=nave-a*oave;for (i=0;iiHeight;i+)for (j=0;j255.0)*lpPixel = 255;else if (temp0.0)*lpPixel = 0;else*lpPixel = (unsigned char) (temp+0.5);UpdateAllViews(NULL);return;/四邻域滤波void CJihanfeiDoc:Onneibor() / TODO: Add your command handler code hereLPSTR lpStartBit = m_Dib.GetPixelBit();/ 象素起始位置指针if (lpStartBit = NULL)/ 找不到象素起始位置指针return ;/ 返回FALSEint iWidth = m_Dib.GetWidth();/ 图象宽度int iHeight = m_Dib.GetHeight();/ 图象高度int iBitsWidth = m_Dib.GetBitsWidth();/ 图象中每行的字节数/ 临时变量double dtemp=0;/ 临时变量unsigned char * lpPixel;/ 指向源DIB图像象素位置的指针int i;/ 循环变量int j;for (i=1;iiHeight-1;i+)for (j=1;jiWidth-1;j+)/读取每象素存储位置lpPixel=(unsigned char*)lpStartBit + (iHeight-1-i) * iBitsWidth + j;*lpPixel=(*(lpPixel+iBitsWidth)+*(lpPixel-iBitsWidth)+*(lpPixel+1)+*(lpPixel-1)/4;UpdateAllViews(NULL);return;/8邻域滤波void CJihanfeiDoc:Onneibor8() / TODO: Add your command handler code hereLPSTR lpStartBit = m_Dib.GetPixelBit();/ 象素起始位置指针if (lpStartBit = NULL)/ 找不到象素起始位置指针return ;/ 返回FALSEint iWidth = m_Dib.GetWidth();/ 图象宽度int iHeight = m_Dib.GetHeight();/ 图象高度int iBitsWidth = m_Dib.GetBitsWidth();/ 图象中每行的字节数/ 临时变量double dtemp=0;/ 临时变量unsigned char * lpPixel;/ 指向源DIB图像象素位置的指针int i;/ 循环变量int j;for (i=1;iiHeight-1;i+)for (j=1;jiWidth-1;j+)/读取每象素存储位置lpPixel=(unsigned char*)lpStartBit + (iHeight-1-i) * iBitsWidth + j;*lpPixel=(*(lpPixel-iBitsWidth-1)+*(lpPixel-iBitsWidth+1)+*(lpPixel+iBitsWidth-1)+*(lpPixel+iBitsWidth+1)+*(lpPixel+iBitsWidth)+*(lpPixel-iBitsWidth)+*(lpPixel+1)+*(lpPixel-1)/8;UpdateAllViews(NULL);return;/多邻域滤波void CJihanfeiDoc:Onmneibor() / TODO: Add your command handler code hereLPSTR lpStartBit = m_Dib.GetPixelBit();/ 象素起始位置指针if (lpStartBit = NULL)/ 找不到象素起始位置指针return ;/ 返回FALSEint iWidth = m_Dib.GetWidth();/ 图象宽度int iHeight = m_Dib.GetHeight();/ 图象高度int iBitsWidth = m_Dib.GetBitsWidth();/ 图象中每行的字节数/ 临时变量double dtemp=0;/ 临时变量unsigned char * lpPixel;/ 指向源DIB图像象素位置的指针unsigned char * lpPixel2;int i;/ 循环变量int j;int p,q;double sum=0;intn=3;for (i=n/2;iiHeight-n/2;i+)for (j=n/2;jiWidth-n/2;j+)/读取每象素存储位置lpPixel=(unsigned char*)lpStartBit + (iHeight-1-i) * iBitsWidth + j;for(p=i-n/2;p=i+n/2;p+)for(q=j-n/2;q=j+n/2;q+)lpPixel2=(unsigned char*)(lpStartBit + (iHeight-1-p) * iBitsWidth + q);sum+=*lpPixel2;*lpPixel=unsigned char(sum/(n*n);sum=0;UpdateAllViews(NULL);return;/高斯滤波void CJihanfeiDoc:Ongaosi() / TODO: Add your command handler code hereLPSTR lpStartBit = m_Dib.GetPixelBit();/ 象素起始位置指针if (lpStartBit = NULL)/ 找不到象素起始位置指针return ;/ 返回FALSEint iWidth = m_Dib.GetWidth();/ 图象宽度int iHeight = m_Dib.GetHeight();/ 图象高度int iBitsWidth = m_Dib.GetBitsWidth();/ 图象中每行的字节数/ 临时变量double dtemp=0;/ 临时变量unsigned char * lpPixel;/ 指向源DIB图像象素位置的指针int i;/ 循环变量int j;for (i=1;iiHeight-1;i+)for (j=1;jiWidth-1;j+)/读取每象素存储位置lpPixel=(unsigned char*)lpStartBit + (iHeight-1-i) * iBitsWidth + j;*lpPixel=(*(lpPixel-iBitsWidth-1)+*(lpPixel-iBitsWidth+1)+*(lpPixel+iBitsWidth-1)+*(lpPixel+iBitsWidth+1)+*(lpPixel+iBitsWidth)*3+*(lpPixel-iBitsWidth)*3+*(lpPixel+1)*3+*(lpPixel-1)*3+*lpPixel*8)/24;UpdateAllViews(NULL);return;/中值滤波void CJihanfeiDoc:Oncenter() / TODO: Add your command handler code hereLPSTR lpStartBit = m_Dib.GetPixelBit();/ 象素起始位置指针if (lpStartBit = NULL)/ 找不到象素起始位置指针return ;/ 返回FALSEint iWidth = m_Dib.GetWidth();/ 图象宽度int iHeight = m_Dib.GetHeight();/ 图象高度int iBitsWidth = m_Dib.GetBitsWidth();/ 图象中每行的字节数/ 临时变量double dtemp=0;/ 临时变量unsigned char * lpPixel;/ 指向源DIB图像象素位置的指针unsigned char * lpPixel2;int i;/ 循环变量int j;int p,q;char temp;double sum=0;intn=3;for (i=n/2;iiHeight-n/2;i+)for (j=n/2;jiWidth-n/2;j+)/读取每象素存储位置lpPixel=(unsigned char*)lpStartBit + (iHeight-1-i) * iBitsWidth + j;for(p=i-n/2;pi+n/2;p+)for(q=j-n/2;q(*(lpPixel2+1)temp=*lpPixel2;*lpPixel2=*(lpPixel2+1);*(lpPixel2+1)=temp;if(*lpPixel5)*lpPixel=*(lpPixel2-n/2*iBitsWidth-n/2);else*lpPixel=*lpPixel;UpdateAllViews(NULL);return;/梯度法锐化void CJihanfeiDoc:Onlayer() / TODO: Add your command handler code hereLPSTR lpStartBit = m_Dib.GetPixelBit();/ 象素起始位置指针if (lpStartBit = NULL)/ 找不到象素起始位置指针return ;/ 返回FALSEint iWidth = m_Dib.GetWidth();/ 图象宽度int iHeight = m_Dib.GetHeight();/ 图象高度int iBitsWidth = m_Dib.GetBitsWidth();/ 图象中每行的字节数/ 临时变量double dtemp=0;/ 临时变量unsigned char * lpPixel;/ 指向源DIB图像象素位置的指针int i;/ 循环变量int j;char x;char y;char p;for (i=0;iiHeight-1;i+)for (j=0;j80)*lpPixel=255;else *lpPixel=0;UpdateAllViews(NULL);return;/水平方向微分算子void CJihanfeiDoc:Onone() / TODO: Add your command handler code hereLPSTR lpStartBit = m_Dib.GetPixelBit();/ 象素起始位置指针if (lpStartBit = NULL)/ 找不到象素起始位置指针return ;/ 返回FALSEint iWidth = m_Dib.GetWidth();/ 图象宽度int iHeight = m_Dib.GetHeight();/ 图象高度int iBitsWidth = m_Dib.GetBitsWidth();/ 图象中每行的字节数/ 临时变量double dtemp=0;/ 临时变量unsigned char * lpPixel;/ 指向源DIB图像象素位置的指针int i;/ 循环变量int j;for (i=1;iiHeight-1;i+)for (j=1;j130)*lpPixel=255;else*lpPixel=0;UpdateAllViews(NULL);return;/垂直方向微分算子void CJihanfeiDoc:Ontwo() / TODO: Add your command handler code hereLPSTR lpStartBit = m_Dib.GetPixelBit();/ 象素起始位置指针if (lpStartBit = NULL)/ 找不到象素起始位置指针return ;/ 返回FALSEint iWidth = m_Dib.GetWidth();/ 图象宽度int iHeight = m_Dib.GetHeight();/ 图象高度int iBitsWidth = m_Dib.GetBitsWidth();/ 图象中每行的字节数/ 临时变量double dtemp=0;/ 临时变量 unsigned char * lpPixel;/ 指向源DIB图像象素位置的指针int i;/ 循环变量int j;for (i=1;iiHeight-1;i+)for (j=1;j130)*lpPixel=255;else *lpPixel=0;UpdateAllViews(NULL);return;/Robert算子锐化void CJihanfeiDoc:OnRobert() / TODO: