clear;
%打開任意一幅圖片
[filepath]=uigetfile（{‘*.jpg‘‘（*.jpg）?jpg圖片‘;‘*.bmp‘‘（*.bmp）?BMP圖片‘;‘*.jpeg‘‘（*.jpeg）?jpeg圖片‘;‘*.*‘‘（*.*）?所有文件‘}[‘請選擇圖片！（^_^）‘]）;???
I=imread（[pathfile]）;
%subplot（121）;
figureimshow（I）;?????????????%figure4待處理圖像
%[lcd]=size（I）;
x=I（::1）;
xx=I（::2）;
xxx=I（::3）;
I1=rgb2hsi（I）;???%將RGB中的原圖轉(zhuǎn)換到HSV空間
f=I1（::1）;
ff=I1（::2）;
fff=I1（::3）;
%%%%%%%%%%%%%構(gòu)造高斯核%%%%%%%%%%%%%%%%%%%%%%%
k1=8;
k2=8;
%r=161;
r=100
alf=1600;
nn=floor（（r+1）/2）;
for?i=1:r
????for?j=1:r
????????b（ij）?=exp（-（（i-nn）^2+（j-nn）^2）/（k1*alf））/（k2*pi*alf*10000）;?
?????????%b（ij）?=exp（-（（i-nn）^2+（j-nn）^2）/（k1*alf））;?%高斯函數(shù)
???end
end
%%%%%%%%%%%對R分量的處理%%%%%%%%%%%%%
%Img?=?double（f）;
%K=imfilter（Imgb）;??
%G=log（Img+1）-log（K+1）;
%G=Img.*K;
%mi=min（min（G））;
%ma=max（max（G））;
%???????L=（G-mi）*255/（ma-mi）;
%%%%%%%%%%%%%%對G分量的處理%%%%%%%%%%%%%
%Img?=?double（ff）;
%K=imfilter（Imgb）;??
%G=log（Img+1）-log（K+1）;
%G=Img.*K;
%mi=min（min（G））;
%ma=max（max（G））;
???????%LL=（G-mi）*255/（ma-mi）;
%%%%%%%%%%%%%隨B分量的處理%%%%%%%%%