东南大学系统实验(DSP)实验报告.docxVIP

系统实验(DSP)实验报告

实验题目

=1\*GB2⑴图像的锐化处理〔高通滤波处理〕

处理模板如下：

对应数学表达式：

=2\*GB2⑵图像的边缘检测

的梯度分别为：

总梯度幅度

实验原理

为了更加通俗地描述图像处理方法，这里引入模板的概念。模板就是一个矩阵，模板大小通常为2?2、3?3，模板也好比一个窗口，将模板窗口覆盖在图像上，窗口覆盖到的象素就是将要处理的象素，而象素所对应的模板矩阵元素的值就是加权值，模板运算就是将模板矩阵元素与对应的象素值相乘并求和。如下列图所示，模板为3?3大小，对应到图像上窗口亦为3?3大小，窗口中有九个象素与模板相对应。模板第一行与、与、与

图像模板示意图

相对应，模板第二、三行依此类推。

根据模板运行规那么，3?3模板对应的表达式数学为：

=1\*GB2⑴、图像的锐化处理

锐化式可以用模板表示为：

对应的数学表达式是：

=2\*GB2⑵、图像的边界检测

边界特点：沿边界走向特性变化比拟缓慢，垂直于边界走向特性变化比拟剧烈。

间接检测：特性不连续性，检测梯度。

二维图象的总梯度定义为：

为了处理方便，也可采用如下定义：

或

总梯度简化为。

常用检测方法如下：

=1\*GB3①梯度法

其梯度幅度为

=2\*GB3②Roberts梯度

=3\*GB3③平面拟合算子

=4\*GB3④算子

实验内容

=1\*GB2⑴图像的锐化处理〔高通滤波处理〕

C语言方法

程序

int*buffer_red_org,*buffer_green_org,*buffer_blue_org;

int*buffer_red,*buffer_green,*buffer_blue;

staticRead_RGB_Image();

voidmain()

{

buffer_red_org=(int*)0

buffer_green_org=(int*)0

buffer_blue_org=(int*)0

buffer_red=(int*)0

buffer_green=(int*)0

buffer_blue=(int*)0

ImageWidth=250;

ImageHeight=200;

Read_RGB_Image(Fish250_200_RGB_CCS);

ImageAcute(buffer_red,buffer_red_org);

ImageAcute(buffer_green,buffer_green_org);

ImageAcute(buffer_blue,buffer_blue_org);

//ImageAcuteSub(ImageWidth,ImageHeight,buffer_red,buffer_red_org);

//ImageAcuteSub(ImageWidth,ImageHeight,buffer_green,buffer_green_org);

//ImageAcuteSub(ImageWidth,ImageHeight,buffer_blue,buffer_blue_org);

}

ImageAcute(int*buffer,int*buffer_org)

{inti;

intx,y;

intr;

floata=0.25;

i=ImageWidth+1;

for(y=1;yImageHeight;y++)

for(x=1;xImageWidth;x++)

{r=(1+4*a)*buffer_org[i]-a*(buffer_org[i-1]+buffer_org[i-ImageWidth]+buffer_org[i+1]+buffer_org[i+ImageWidth]);

if(r255)r=255;

if(r0)r=0;

buffer[i++]=r;

}

}

staticRead_RGB_Image(char*FileName)

{inti;

intx,y;

unsignedcharr,g,b,temp;

FILE*fp;