基于各向异性扩散的数字图像处理 作业

《基于各向异性扩散的数字图像处理 作业》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、基于各向异性扩散的数字图像处理1Perona-Malik算法概述各向异性扩散作为现行的一种非常流行的偏微分方程数字图像处理技术，是由传统的Gaussian滤波发展而来的，有着强大的理论基础，并有着传统的数字图像方法无法企及的良好特性，其特点是可以在平滑的同时保持边缘特征。由于这种优良的特性，使其在图像的平滑、去噪、恢复、增强和分割等方面得到了广泛的应用。Koenderink和Witkin两位学者把尺度空间的严格理论引入到了数字图像处理之中，而尺度空间理论正是现代偏微分方程数字图像处理的理论基础，他们的工作主要是将多尺度

2、图像表示为Gaussian滤波器处理的结果，相当于将原图像输入热传导方程进而得到的序列时间图像。多尺度滤波得到的图像序列，可以看作各向同性热扩散方程的解，这就是用扩散方程的观点来看待高斯滤波：(1.1)不仅是热扩散方程可以产生尺度空间，其他的抛物线方程也可以产生尺度空间，进而满足极大值原理的演化方程也能定义一类尺度空间。各向异性扩散方程的提出是这个领域中一个里程碑式的标志，它开辟了数字图像处理中偏微分方程理论和应用的一个新领域，提出了各向异性扩散（AnisotropicDiffusion）方程，一般称为Perona-M

3、alik（简称P-M）扩散模型。(1.2)其中是散度算子，是图像的梯度，是扩散系数。各向异性扩散是一个能量散发的过程，这个过程与能量曲面的形状有关，从数学上来看，其等价于一个能量最小化问题的求解。从数学意义上讲Perona-Malik模型是改进的热传导偏微分方程。它是通过函数自适应地控制扩散速度。理论上希望边缘内部的区域，c远离0，图像可以平滑，而在边缘附近，，图像不再平滑。铃铛形的径向函数g(w)可取(1.3)或者(1.4)当时形成了各向异性的扩散模型：(1.5)可以利用Taylor级数展开进行线性近似(1.6)由于

4、(1.7)所以可以得出P-K模型简化推导：(1.8)这样把较长的时间t分割为，始终从步计算到步即可(1.9)把改写为，改写为，其中(1.10)(1.11)则(1.12)其中(1.13)则得到Perona-Malik算法的下述迭代形式(1.14)其中，以Perona-Malik算法为代表的各向异性扩散滤波模型既能有效地去除噪声，又能较好地保持图像的边缘和重要的细节信息。但是考虑到PM模型对小尺度噪声敏感以及部分边缘和细节失真问题的存在，对模型的扩散系数的改进仍然需要进一步的改进。2算法实现本算法的实验在VC++6.0环境

5、下实现。先建一个类CDib，用来实现位图的载入以及一些参数的获取。新建一个类CPeronaMalik，类中函数PM(CDib*pDib)用来实现PM算法。voidCPeronaMalik::PM(CDib*pDib)｛inti=0,j=0;intnWidth=pDib->GetWidth();intnHeight=pDib->GetHeight();intnWidthBytes=pDib->GetWidth();unsignedchar*pDibBits=pDib->GetpDibBits();//每行字节数intli

6、neByte=(bmpWidth*biBitCount/8+3)/4*4;//循环变量，针对彩色图像，遍历每像素的三个分量intk;intn;doubletemp1=0;doubletemp2=0;intnn;//迭代次数doubledeltat;doubledelta;//步长，deltat小于0.25//c1、c2表示扩散函数intp;//范数doublesigma;//平均灰度doubleK;if(biBitCount==8){//对于灰度图像for(i=0;i

7、idth;j++){temp1=double(pDibBits[i*nWidthBytes+j])+temp1;}}//添加图像处理的代码////////////定参数////////////////////////nn=200;//迭代次数deltat=0.02;delta=deltat/2;//步长，deltat小于0.25//c1、c2表示扩散函数p=2;//范数sigma=0.4*temp1/(double(bmpHeight)*double(bmpWidth));//平均灰度inta=1;//参数,a>0//

8、////////////////////////////////////////for(i=0;i