function?g=bilateral_color（f）
%?f為待處理的彩色圖像，必須歸一化到【01】范圍內(nèi)
%?對(duì)于彩色圖像，像素值的接近與否不能用RGB空間值，建議使用CIE顏色空間
%applycform用于顏色空間轉(zhuǎn)換，格式為：
%?cform=makecform（‘srgb2lab‘）;?
%?f=applycform（fcform）;

figureimshow（f[]）;

f=double（f）/255;

f=applycform（fmakecform（‘srgb2lab‘））;

[mnc]=size（f）;
r=5;
fr=padarray（f[rr]‘symmetric‘）;
sigma_s=input（‘input?spatial?similarity:‘）;?%?2?4?8……
sigma_r=input（‘input?range?similarity:‘）;?%?0.1?0.2?0.4……
sigma_r=sigma_r*100;

w1=fspecial（‘gaussian‘[2*r+12*r+1]sigma_s）;

h=waitbar（0‘wait……‘）;

for?i=r+1:r+m
????for?j=r+1:r+n
????????diffpix（::1）=fr（i-r:i+rj-r:j+r1）-fr（ij1）;?%?計(jì)算像素差值
????????diffpix（::2）=fr（i-r:i+rj-r:j+r2）-fr（ij2）;
????????diffpix（::3）=fr（i-r:i+rj-r:j+r3）-fr（ij3）;
?
????????w2=exp（-（diffpix（::1）.^2+diffpix（::2）.^2+diffpix（::3）.^2）/（2*sigma_r.^2））;
????????
????????w=w1.*w2;
?????????
?????????fr（ij1）=sum（fr（i-r:i+rj-r:j+r1）.*w）/sum（w）;
?????????fr（ij2）=sum（fr（i-r:i+rj-r:j+r2）.*w）/sum（w）;
?????????fr（ij3）=sum（fr（i-r:i+rj-r:j+r3）.*w）/sum（w）;
????end
????waitbar（i/m）;
end
close（h）;

g（::1）=fr（r+1:r+mr+1:r+n1）;
g（::2）=fr（r+1:r+mr+1:r+n2）;
g（::3）=fr（r+1:r+mr+1:r+n3）;

g=applycform（gmakecform（‘lab2srgb‘））;
figureimshow（g[]）;

end

????????????
????????????
????????????
????????????
????????????
????????????
????????????
????????????
????????????
????????????

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

?????文件???????1499??2014-08-10?08:36??bilateral_color.m

?????文件????????295??2014-08-10?08:47??bilateral_filter.m

?????文件???????1000??2014-08-10?09:05??bilateral_gray.m

?????文件??????16439??2014-08-09?07:44??hui1.jpg

?????文件??????16937??2014-08-09?07:45??yellowrose.jpg

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

????????????????36170????????????????????5