主要是用來平滑圖像的，克服了高斯模糊的缺陷，各向異性擴(kuò)散在平滑圖像時是保留圖像邊緣的（和雙邊濾波很像）。

通常我們有將圖像看作矩陣的，看作圖的，看作隨機(jī)過程的，記得過去還有看作力場的。

這次新鮮，將圖像看作熱量場了。每個像素看作熱流，根據(jù)當(dāng)前像素和周圍像素的關(guān)系，來確定是否要向周圍擴(kuò)散。比如某個鄰域像素和當(dāng)前像素差別較大，則代表這個鄰域像素很可能是個邊界，那么當(dāng)前像素就不向這個方向擴(kuò)散了，這個邊界也就得到保留了。

先看下效果吧：

具體的推導(dǎo)公式都是熱學(xué)上的，自己也不太熟悉，感興趣的可以去看原論文，引用量超7000吶。

我這里只介紹一下最終結(jié)論用到的公式。

主要迭代方程如下：

?I就是圖像了，因?yàn)槭莻€迭代公式，所以有迭代次數(shù)t。

四個散度公式是在四個方向上對當(dāng)前像素求偏導(dǎo)，news就是東南西北嘛，公式如下：

?而cN/cS/cE/cW則代表四個方向上的導(dǎo)熱系數(shù)，邊界的導(dǎo)熱系數(shù)都是小的。公式如下：

?

?最后整個公式需要先前設(shè)置的參數(shù)主要有三個，迭代次數(shù)t，根據(jù)情況設(shè)置；導(dǎo)熱系數(shù)相關(guān)的k，取值越大越平滑，越不易保留邊緣；lambda同樣也是取值越大越平滑。

最后是matlab代碼：

clear all; close all; clc;k=15; %導(dǎo)熱系數(shù),控制平滑 lambda=0.15; %控制平滑 N=20; %迭代次數(shù) img=double(imread('lena.jpg')); imshow(img,[]); [m n]=size(img);imgn=zeros(m,n); for i=1:Nfor p=2:m-1for q=2:n-1%當(dāng)前像素的散度，對四個方向分別求偏導(dǎo)，局部不同方向上的變化量，%如果變化較多，就證明是邊界，想方法保留邊界NI=img(p-1,q)-img(p,q);SI=img(p+1,q)-img(p,q);EI=img(p,q-1)-img(p,q);WI=img(p,q+1)-img(p,q);%四個方向上的導(dǎo)熱系數(shù)，該方向變化越大，求得的值越小，從而達(dá)到保留邊界的目的cN=exp(-NI^2/(k*k));cS=exp(-SI^2/(k*k));cE=exp(-EI^2/(k*k));cW=exp(-WI^2/(k*k));imgn(p,q)=img(p,q)+lambda*(cN*NI+cS*SI+cE*EI+cW*WI); %擴(kuò)散后的新值 endendimg=imgn; %整個圖像擴(kuò)散完畢，用已擴(kuò)散圖像的重新擴(kuò)散。 endfigure; imshow(imgn,[]);

參考：

《特征提取與圖像處理（第二版）》

轉(zhuǎn)載于:https://www.cnblogs.com/tiandsp/archive/2013/04/18/3029468.html

總結(jié)

以上是生活随笔為你收集整理的matlab练习程序（各向异性扩散）的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。