%%%%%%%%%%%%%%擴(kuò)散濾波函數(shù)%%%%%%%%%%%%%%
function?[uttn]?=?CoherenceFilter（uOptions）
%COHERENCEFILTER函數(shù)用來進(jìn)行2D灰度圖像的一致性擴(kuò)散增強(qiáng)，該方法可以在保護(hù)邊緣的同時(shí)降低圖像噪聲，并平滑圖像由于噪聲導(dǎo)致的邊緣斷裂
%?Iout?=?CoherenceFilter（Iin?Options）
%輸入：
%???Iin?????:?2D灰度圖像
%???Options?:?濾波的各個(gè)參數(shù)
%輸出：
%???Iout????:?處理后的圖像
%?Options
%???Options.T??????:??總的擴(kuò)散時(shí)間
%???Options.dt?????:??擴(kuò)散時(shí)間步長?五種擴(kuò)散方案默認(rèn)值為0.2
%???Options.sigma??:??高斯濾波的Sigma值，默認(rèn)為1??????????????????
%???Options.rho????:??計(jì)算結(jié)構(gòu)張量時(shí)高斯平滑的Sigma值，默認(rèn)為1
%???Options.Scheme?:??用數(shù)值表示各擴(kuò)散方案
%?????????????????????‘S‘標(biāo)準(zhǔn)離散化
%?????????????????????‘N‘非負(fù)離散化
%?????????????????????‘I‘隱式離散化
%?????????????????????‘R‘旋轉(zhuǎn)不變方案，標(biāo)準(zhǔn)離散化（隱式）5x5內(nèi)核模板（默認(rèn)）
%?????????????????????‘O‘優(yōu)化的求導(dǎo)內(nèi)核
%???Weickert方程中決定擴(kuò)散幅度的常量
%???Options.C??????:??默認(rèn)值為?1e-10
%???Options.m??????:??默認(rèn)值為?1
%???Options.alpha??:??默認(rèn)值為?0.001
%算法分析：
%???1對圖像進(jìn)行高斯預(yù)處理；
%???2計(jì)算結(jié)構(gòu)張量、特征值、特征向量；
%?????（圖像邊緣給出特征值?特征向量指向最大特征值的方向，用來描述邊緣方向）
%???3?計(jì)算擴(kuò)散張量；
%???4?用一個(gè)有限差分方案來進(jìn)行擴(kuò)散；
%???5?回到步驟1，直到達(dá)到規(guī)定的擴(kuò)散時(shí)間.
%
%?注釋
%?標(biāo)準(zhǔn)和非負(fù)的擴(kuò)散僅允許小的擴(kuò)散時(shí)間步長，不然離散會(huì)變得不穩(wěn)定；隱式、旋轉(zhuǎn)不變和優(yōu)化擴(kuò)散方案允許大的擴(kuò)散時(shí)間步長。?
%?舉例2D
%???I?=?im2double（imread（‘LUNG2‘））;
%???I?=?rgb2gray（I）;
%???JS?=?CoherenceFilter（Istruct（‘T‘1‘rho‘1‘Scheme‘‘S‘））;
%???JN?=?CoherenceFilter（Istruct（‘T‘1‘rho‘1‘Scheme‘‘N‘））;
%???JI?=?CoherenceFilter（Istruct（‘T‘1‘rho‘1‘Scheme‘‘I‘））;
%???JR?=?CoherenceFilter（Istruct（‘T‘1‘rho‘1‘Scheme‘‘R‘））;
%???JO?=?CoherenceFilter（Istruct（‘T‘1‘rho‘1‘Scheme‘‘O‘））;
%???figure?
%???subplot（231）?imshow（I）?title（‘原始圖像‘）;
%???subplot（232）?imshow（JS）?title（‘標(biāo)準(zhǔn)方案‘）;
%???subplot（233）?imshow（JN）?title（‘非負(fù)方案‘）;
%???subplot（234）?imshow（JI）?title（‘隱式方案‘）;
%???subplot（235）?imshow（JR）?title（‘旋轉(zhuǎn)不變方案‘）;
%???subplot（236）?imshow（JO）?title（‘優(yōu)化方案‘）;
%
%增加所有所需函數(shù)的路徑

try?%用到try?catch?函數(shù)：如果執(zhí)行try語段出錯(cuò)，轉(zhuǎn)而執(zhí)行catch語段
????functionname=‘CoherenceFilter.m‘;
????functiondir=which（functionname）;
????functiondir=functiondir（1:end-length（functionname））;
????addpath（[functiondir?‘/functions2D‘]）%將functions2D文件夾添加到matlab搜索路徑
????addpath（[functiondir?‘/functions‘]）%將functions文件夾添加到matlab搜索路徑
catch?me
????disp（me.message）;
end

%設(shè)置默認(rèn)參數(shù)
defaultoptions=struct（‘T‘2‘dt‘[]‘sigma‘?1?‘rho‘?1‘Scheme‘‘S‘??‘C‘?1e-10?‘m‘1‘a(chǎn)lpha‘0.001）;

if（~exist（‘Options‘‘var‘））%如果不存在輸入?yún)?shù)，則使用默認(rèn)參數(shù)
????Options=defaultoptions;
else?%否則的話，使用用戶輸入的參數(shù)
????tags?=?fieldnames（defaultoptions）;%fieldnames：將defaultoptions結(jié)構(gòu)體中的變量名賦值給tags
????for?i=1:length（tags）
????????if（~isfield（Optionstags{i}））??Options.（tags{i}）=defaultoptions.（tags{i}）;?end?%isfield：判斷輸入是否是結(jié)構(gòu)體數(shù)組的成員
????end
????if（length（tags）~=length（fieldnames（Options）））
????????warning（‘CoherenceFilter:unknownoption‘‘unknown?options?found‘）;
????end
end

if（isempty（Options.dt））%如果用戶沒有輸入dt值，則取以下默認(rèn)值
????switch?lower（Options.Scheme）?
??????case?‘s‘?Options.dt=0.2;
??????case?‘n‘?Options.dt=0.2;
??????case?‘i‘?Options.dt=0.2;
??????case?‘r‘?Options.dt=0.2;
??????case?‘o‘?Options.dt=0.2;??
??????otherwise
????????error（‘CoherenceFilter:unknownoption‘‘unknown?scheme‘）;
????end
end
????
%%%%%%%%各向異性擴(kuò)散的主循環(huán)程序%%%%

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件????????429??2013-05-09?16:45??二維擴(kuò)散濾波用于圖像增強(qiáng)\functions2D\ConstructDiffusionTensor2D.m

?????文件????????591??2013-05-23?11:04??二維擴(kuò)散濾波用于圖像增強(qiáng)\functions2D\diffusion_scheme_2D_rotation_invariant.m

?????文件???????1766??2013-05-23?11:04??二維擴(kuò)散濾波用于圖像增強(qiáng)\functions2D\diffusion_scheme_2D_implicit.m

?????文件???????2112??2013-05-23?11:04??二維擴(kuò)散濾波用于圖像增強(qiáng)\functions2D\diffusion_scheme_2D_optimized.m

?????文件???????1221??2013-05-23?11:05??二維擴(kuò)散濾波用于圖像增強(qiáng)\functions2D\diffusion_scheme_2D_standard.m

?????文件???????1016??2013-05-23?11:05??二維擴(kuò)散濾波用于圖像增強(qiáng)\functions2D\EigenVectors2D.m

?????文件???????1636??2013-05-10?14:06??二維擴(kuò)散濾波用于圖像增強(qiáng)\functions2D\diffusion_scheme_2D_non_negativity.m

?????文件????????580??2013-05-10?16:51??二維擴(kuò)散濾波用于圖像增強(qiáng)\functions2D\StructureTensor2D.m

?????文件???????2257??2013-05-23?10:59??二維擴(kuò)散濾波用于圖像增強(qiáng)\functions\imgaussian.m

?????文件????????912??2013-05-23?11:02??二維擴(kuò)散濾波用于圖像增強(qiáng)\functions\derivatives.m

?????文件???????1141??2013-05-23?14:40??二維擴(kuò)散濾波用于圖像增強(qiáng)\test_shiyan1.m

?????文件???????5247??2013-05-23?11:03??二維擴(kuò)散濾波用于圖像增強(qiáng)\CoherenceFilter.m

?????文件???????1511??2013-05-23?14:48??二維擴(kuò)散濾波用于圖像增強(qiáng)\test_shiyan2.m

?????文件???????1256??2013-05-23?11:03??二維擴(kuò)散濾波用于圖像增強(qiáng)\compare.m

?????文件?????177742??2013-04-17?09:29??二維擴(kuò)散濾波用于圖像增強(qiáng)\LUNG2.png

?????目錄??????????0??2013-05-29?16:11??二維擴(kuò)散濾波用于圖像增強(qiáng)\functions2D

?????目錄??????????0??2013-05-29?16:11??二維擴(kuò)散濾波用于圖像增強(qiáng)\functions

?????目錄??????????0??2013-05-29?16:11??二維擴(kuò)散濾波用于圖像增強(qiáng)

-----------?---------??----------?-----??----

???????????????199417????????????????????18