%?BFILTER2?Two?dimensional?bilateral?filtering.
%????This?function?implements?2-D?bilateral?filtering?using
%????the?method?outlined?in:
%
%???????C.?Tomasi?and?R.?Manduchi.?Bilateral?Filtering?for?
%???????Gray?and?Color?Images.?In?Proceedings?of?the?IEEE?
%???????International?Conference?on?Computer?Vision?1998.?
%
%????B?=?bfilter2（AWSIGMA）?performs?2-D?bilateral?filtering
%????for?the?grayscale?or?color?image?A.?A?should?be?a?double
%????precision?matrix?of?size?NxMx1?or?NxMx3?（i.e.?grayscale
%????or?color?images?respectively）?with?normalized?values?in
%????the?closed?interval?[01].?The?half-size?of?the?Gaussian
%????bilateral?filter?window?is?defined?by?W.?The?standard
%????deviations?of?the?bilateral?filter?are?given?by?SIGMA
%????where?the?spatial-domain?standard?deviation?is?given?by
%????SIGMA（1）?and?the?intensity-domain?standard?deviation?is
%????given?by?SIGMA（2）.
%
%?Douglas?R.?Lanman?Brown?University?September?2006.
%?dlanman@brown.edu?http://mesh.brown.edu/dlanman
%?我理解的意思是B=bfilter的第一個參數是待處理的圖片，第二個參數是算子長度，SIGMA（1）是空間的公式里的sigma，SIGMA（2）是公式像素值的SIGMA

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%?Pre-process?input?and?select?appropriate?filter.
%?下面幾個if判斷都是判斷輸入的參數是否正確
function?B?=?bfilter2（Awsigma）

%?Verify?that?the?input?image?exists?and?is?valid.
if?~exist（‘A‘‘var‘）?||?isempty（A）
???error（‘Input?image?A?is?undefined?or?invalid.‘）;
end
if?~isfloat（A）?||?~sum（[13]?==?size（A3））?||?...
??????min（A（:））??1
???error（[‘Input?image?A?must?be?a?double?precision?‘...
??????????‘matrix?of?size?NxMx1?or?NxMx3?on?the?closed?‘...
??????????‘interval?[01].‘]）;??????
end

%?Verify?bilateral?filter?window?size.
%?nume1（w）判斷數組w中元素個數
if?~exist（‘w‘‘var‘）?||?isempty（w）?||?...
??????numel（w）?~=?1?||?w????w?=?5;
end
%?將w朝正四舍五入
w?=?ceil（w）;

%?Verify?bilateral?filter?standard?deviations.
if?~exist（‘sigma‘‘var‘）?||?isempty（sigma）?||?...
??????numel（sigma）?~=?2?||?sigma（1）?<=?0?||?sigma（2）?<=?0
???sigma?=?[3?0.1];
end

%?Apply?either?grayscale?or?color?bilateral?filtering.
if?size（A3）?==?1
???B?=?bfltGray（Awsigma（1）sigma（2））;
else
???B?=?bfltColor（Awsigma（1）sigma（2））;
end


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%?Implements?bilateral?filtering?for?grayscale?images.
function?B?=?bfltGray（Awsigma_dsigma_r）

%?Pre-compute?Gaussian?distance?weights.
[XY]?=?meshgrid（-w:w-w:w）;
G?=?exp（-（X.^2+Y.^2）/（2*sigma_d^2））;
%?Create?waitbar.
h?=?waitbar（0‘Applying?bilateral?filter...‘）;
set（h‘Name‘‘Bilateral?Filter?Progress‘）;

%?Apply?bilateral?filter.
dim?=?size（A）;
B?=?zeros（dim）;
for?i?=?1:dim（1）
???for?j?=?1:dim（2）
??????
?????????%?Extract?local?region.
?????????iMin?=?max（i-w1）;
?????????iMax?=?min（i+wdim（1））;
?????????jMin?=?max（j-w1）;
?????????jMax?=?min（j+wdim（2））;
?????????I?=?A（iMin:iMaxjMin:jMax）;
??????
?????????%?Compute?Gaussian?intensity?weights.
?????????H?=?exp（-（I-A（ij））.^2/（2*sigma_r^2））;
??????
?????????%?Calculate?bilateral?filter?response.
??????

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2019-08-29?20:42??SAR圖像變化檢測方法，包括了對數比和均值比兩種較為經典的，\
?????文件????????4905??2018-01-16?12:11??SAR圖像變化檢測方法，包括了對數比和均值比兩種較為經典的，\bfilter2.m
?????文件???????????0??2019-08-27?17:23??SAR圖像變化檢測方法，包括了對數比和均值比兩種較為經典的，\cols2.m
?????文件????????2863??2019-08-26?20:06??SAR圖像變化檢測方法，包括了對數比和均值比兩種較為經典的，\LOG_D_image.m
?????文件????????3485??2019-08-26?20:03??SAR圖像變化檢測方法，包括了對數比和均值比兩種較為經典的，\Maxture.m
?????文件????????3504??2019-08-26?20:00??SAR圖像變化檢測方法，包括了對數比和均值比兩種較為經典的，\MEAN_D_image.m
?????文件??????305254??2016-04-15?16:26??SAR圖像變化檢測方法，包括了對數比和均值比兩種較為經典的，\sar1gray.bmp
?????文件??????305254??2016-04-15?16:26??SAR圖像變化檢測方法，包括了對數比和均值比兩種較為經典的，\sar2gray.bmp
?????文件??????305254??2016-04-15?16:26??SAR圖像變化檢測方法，包括了對數比和均值比兩種較為經典的，\sarreference.bmp
?????文件???????????0??2019-08-26?17:51??SAR圖像變化檢測方法，包括了對數比和均值比兩種較為經典的，\Xm.m
?????文件??????103278??2019-08-27?17:41??SAR圖像變化檢測方法，包括了對數比和均值比兩種較為經典的，\Xm2.bmp