%?BFILTER2?Two?dimensional?bilateral?filtering.
%?This?function?implements?2-D?bilateral?filtering?using?B?=?bfilter2（AWSIGMA）?performs?2-D?bilateral?filtering
%?for?the?grayscale?A.?A?should?be?a?double?precision?matrix?of?size?NxMx1?or?NxMx3?（i.e.?grayscale?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）.

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%?Pre-process?input?and?select?appropriate?filter.
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.
if?~exist（‘w‘‘var‘）?||?isempty（w）?||?...
??????numel（w）?~=?1?||?w????w?=?5;
end
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.
?????????F?=?H.*G（（iMin:iMax）-i+w+1（jMin:jMax）-j+w+1）;
?????????B（ij）?=?sum（F（:）.*I（:））/sum（F（:））;
???????????????
???end
???waitbar（i/dim（1））;
end

%?Close?waitbar.
close（h）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件??????10328??2015-11-07?18:25??79419116Bilateral-Filtering-advanced\Bilateral?Filtering\475.tif

?????文件???????2733??2015-11-28?18:02??79419116Bilateral-Filtering-advanced\Bilateral?Filtering\bfilter2.m

?????文件??????42447??1997-04-29?01:19??79419116Bilateral-Filtering-advanced\Bilateral?Filtering\einstein.jpg

?????文件????????835??2018-01-09?16:45??79419116Bilateral-Filtering-advanced\Bilateral?Filtering\run.m

????..A.SH.?????17920??2015-11-30?09:31??79419116Bilateral-Filtering-advanced\Bilateral?Filtering\Thumbs.db

?????文件??????10328??2015-11-07?18:25??79419116Bilateral-Filtering-advanced\Bilateral?Filtering?advanced\475.tif

?????文件??????42447??1997-04-29?01:19??79419116Bilateral-Filtering-advanced\Bilateral?Filtering?advanced\einstein.jpg

?????文件???????3162??2015-11-28?18:21??79419116Bilateral-Filtering-advanced\Bilateral?Filtering?advanced\jbfilter2.m

?????文件????????841??2018-01-09?16:47??79419116Bilateral-Filtering-advanced\Bilateral?Filtering?advanced\run1.m

????..A.SH.?????14336??2015-11-28?16:50??79419116Bilateral-Filtering-advanced\Bilateral?Filtering?advanced\Thumbs.db

?????目錄??????????0??2017-02-07?14:02??79419116Bilateral-Filtering-advanced\Bilateral?Filtering

?????目錄??????????0??2018-01-23?20:56??79419116Bilateral-Filtering-advanced\Bilateral?Filtering?advanced

?????目錄??????????0??2018-01-09?16:39??79419116Bilateral-Filtering-advanced

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

???????????????145377????????????????????13