数字图像处理集中实习报告封面和书写格式

《数字图像处理》

集中实习报告

（ 2010 - 2011学年第 2 学期）

---------------------------------------------- 实习成绩：

教师评语：

教师签名：

年月日

实验项目：图片的边缘检测

一、 实验内容

1 边缘检测的基本思想是先利用边缘增强算子，突出图像中的局部边缘。

2 对bmp 格式的图片进处理，并用边缘检测算子，检测图像的边缘。

3 算子的原理就是利用边缘灰度变化的一阶或二阶导数特点，可以将边缘点检测出来。

4 用Sobel 算子，Prewitt 算子，拉普拉斯算子，检测图像的边缘。

Sobel 算子的模板

二、算法实现的主要内容

(说明：填写每个算子实现的核心代码)

void CDib::RoberPrewiSobel(BYTE *&m_pBitmap, DWORD m_dwHeight, DWORD m_dwWidth, WORD m_flag,int kk)

{

if(m_flag==24)

{

m_dwWidth *= 3;

}

else if(m_flag != 8)

{

AfxMessageBox("只处理8位和24位图像！");

return;

}

//为输出图像申请空间

BYTE *pData = new BYTE[m_dwWidth*m_dwHeight];

BYTE *pDataOut = new BYTE[m_dwWidth*m_dwHeight];

BYTE *pOld = pData;

memcpy(pData, m_pBitmap, m_dwWidth*m_dwHeight);

memcpy(pDataOut, m_pBitmap, m_dwWidth*m_dwHeight);

int byte=m_flag/8;

int a[9]={0,0,0,0,0,0,0,0,0},b[9]={0,0,0,0,0,0,0,0,0};

if(kk==1)//Prewitt算子

{

a[0]=-1; a[2]=1; a[3]=-1; a[5]=1; a[6]=-1; a[8]=1;

b[0]=-1; b[1]=-1; b[2]=-1; b[6]=1; b[7]=1; b[8]=1;

}

else if(kk==2)//Sobel算子

{

a[0]=-1; a[2]=1; a[3]=-2; a[5]=2;a[6]=-1; a[8]=1;

b[0]=-1;b[1]=-2; b[2]=-1; b[6]=1;b[7]=2; b[8]=1;

}

DWORD i,j,k,t;

int x,y;

for(k=0;k

{for(i=1;i

{for(j=1;j<(m_dwWidth-1)/byte;j++)

{ x=a[0]**(pData+(i-1)*m_dwWidth+(j-1)*byte+k)

+a[1]**(pData+(i-1)*m_dwWidth+j*byte+k)

+a[2]**(pData+(i-1)*m_dwWidth+(j+1)*byte+k)

+a[3]**(pData+i*m_dwWidth+(j-1)*byte+k)

+a[4]**(pData+i*m_dwWidth+j*byte+k)

+a[5]**(pData+i*m_dwWidth+(j+1)*byte+k)

+a[6]**(pData+(i+1)*m_dwWidth+(j-1)*byte+k)

+a[7]**(pData+(i+1)*m_dwWidth+j*byte+k)

+a[8]**(pData+(i+1)*m_dwWidth+(j+1)*byte+k);

y=b[0]**(pData+(i-1)*m_dwWidth+(j-1)*byte+k) +b[1]**(pData+(i-1)*m_dwWidth+j*byte+k)

+b[2]**(pData+(i-1)*m_dwWidth+(j+1)*byte+k)

+b[3]**(pData+i*m_dwWidth+(j-1)*byte+k)

+b[4]**(pData+i*m_dwWidth+j*byte+k)

+b[5]**(pData+i*m_dwWidth+(j+1)*byte+k)

+b[6]**(pData+(i+1)*m_dwWidth+(j-1)*byte+k)

+b[7]**(pData+(i+1)*m_dwWidth+j*byte+k)

+b[8]**(pData+(i+1)*m_dwWidth+(j+1)*byte+k);

t=int(sqrt(x*x+y*y)+0.5);

if(t>255)

t=255;

*(pDataOut+i*m_dwWidth+j*byte+k)=t;

}

}

}

pData = pOld;

memcpy(m_pBitmap, pDataOut, m_dwWidth*m_dwHeight);

if(pData) delete []pData;

if(pDataOut) delete []pDataOut;

}

//Laplacian算子,Laplacian拓展算子

void CDib::Suanzi(BYTE *&m_pBitmap, DWORD m_dwHeight, DWORD m_dwWidth, WORD m_flag,int kk) {

if(m_flag==24)

{m_dwWidth *= 3; }

else if(m_flag != 8)

{

AfxMessageBox("只处理8位和24位图像！");

return;

}

BYTE *pData = new BYTE[m_dwWidth*m_dwHeight];

BYTE *pDataOut = new BYTE[m_dwWidth*m_dwHeight];

BYTE *pOld = pData;

memcpy(pData, m_pBitmap, m_dwWidth*m_dwHeight);

memcpy(pDataOut, m_pBitmap, m_dwWidth*m_dwHeight);

int byte=m_flag/8;

int a[9]={0,0,0,0,0,0,0,0,0};

if(kk==0) //Laplacian算子

{

a[1]=-1; a[3]=-1; a[4]=4; a[5]=-1; a[7]=-1;

}

else if(kk==1)// Laplacian增强算子

{

a[1]=-1; a[3]=-1; a[4]=8; a[5]=-1; a[7]=-1;

}

DWORD i,j,k,t;

for(k=0;k

{

for(i=1;i

{

for(j=1;j<(m_dwWidth-1)/byte;j++)

{

//第一行

t=a[0]**(pData+(i-1)*m_dwWidth+(j-1)*byte+k)

+a[1]**(pData+(i-1)*m_dwWidth+j*byte+k)

+a[2]**(pData+(i-1)*m_dwWidth+(j+1)*byte+k) //第二行

+a[3]**(pData+i*m_dwWidth+(j-1)*byte+k)

+a[4]**(pData+i*m_dwWidth+j*byte+k)

+a[5]**(pData+i*m_dwWidth+(j+1)*byte+k)

//第三行

+a[6]**(pData+(i+1)*m_dwWidth+(j-1)*byte+k)

+a[7]**(pData+(i+1)*m_dwWidth+j*byte+k)

+a[8]**(pData+(i+1)*m_dwWidth+(j+1)*byte+k);

t=int(abs(t)+0.5);

if(t>255)

t=255;

*(pDataOut+i*m_dwWidth+j*byte+k)=t;

}

}

}

pData = pOld;

memcpy(m_pBitmap, pDataOut, m_dwWidth*m_dwHeight);

if(pData) delete []pData;

if(pDataOut) delete []pDataOut;

}

三、处理结果

Prewitt算子处理结果

Soble算子

Laplacian算子运算模板

四、实习体会

这次实习让我对MFC编程的概念有了更深一步的认识，比如说封装性。我还发现用CDib很方便，很好地利用了类的封装性，可以用把程序写到一个主程序Dib.cpp里面,然后在view里面调用。

刚开始做的时候什么都不懂，实在是不知从哪里做起，只好把书上的代码写进去，不断改正，慢慢地也看懂了代码的意思。后面的就是根据前面的该一下。但是算法还是基本上和书上差不多。不过自己编的时候还是有很多细节的部分没有注意到，比如说，强制数据类型转换，数组定义赋值等问题，我自己编的时候没有注意这个问题，结果出了很多错，经过不断改正终于使程序正常运行。

正如老师说的那样，C++的学习很大程度上靠自学，平时要多练多问，多和同学交流学习，这样才能不断进步，还有就是要不怕麻烦，多多练习！