f=imread（‘01_test.tif‘）;
fr=f（::1）;?fg=f（::2）;?fb=f（::3）;
mr=im2double（fr）;?mg=im2double（fg）;?mb=im2double（fb）;
n=141;%定義模板大小。?kid=141；
n1=floor（（n+1）/2）;%確定中心
a1=60;?%定義標準差（尺度）?kid=60；
for?i=1:n
for?j=1:n
b（ij）?=exp（-（（i-n1）^2+（j-n1）^2）/（a1*a1））/（pi*a1*a1）;?%高斯函數。
end
end
nr1=imfilter（mrb‘conv‘‘replicate‘）;
ng1=imfilter（mgb‘conv‘‘replicate‘）;
nb1=imfilter（mbb‘conv‘‘replicate‘）;%卷積濾波。
ur1=log（nr1）;?ug1=log（ng1）;?ub1=log（nb1）;
tr1=log（mr+eps）;tg1=log（mg+eps）;tb1=log（mb+eps）;
yr1=（tr1-ur1）/3;yg1=（tg1-ug1）/3;yb1=（tb1-ub1）/3;
a2=10;?%定義標準差（尺度）
for?i=1:n
33
for?j=1:n
a（ij）?=exp（-（（i-n1）^2+（j-n1）^2）/（a2*a2））/（pi*a2*a2）;?%高斯函數。
end
end
nr2=imfilter（mra‘conv‘‘replicate‘）;
ng2=imfilter（mga‘conv‘‘replicate‘）;
nb2=imfilter（mba‘conv‘‘replicate‘）;%卷積濾波。
ur2=log（nr2）;ug2=log（ng2）;ub2=log（nb2）;
tr2=log（mr+eps）;tg2=log（mg+eps）;tb2=log（mb+eps）;
yr2=（tr2-ur2）/3;yg2=（tg2-ug2）/3;yb2=（tb2-ub2）/3;
a3=150;?%定義標準差（尺度）kid=150；
for?i=1:n
for?j=1:n
e（ij）?=exp（-（（i-n1）^2+（j-n1）^2）/（a3*a3））/（pi*a3*a3）;?%高斯函數。
end
end
nr3=imfilter（mre‘conv‘‘replicate‘）;
ng3=imfilter（mge‘conv‘‘replicate‘）;
nb3=imfilter（mbe‘conv‘‘replicate‘）;%卷積濾波。
ur3=log（nr3）;ug3=log（ng3）;ub3=log（nb3）;
tr3=log（mr+eps）;tg3=log（mg+eps）;tb3=log（mb+eps）;
yr3=（tr3-ur3）/3;yg3=（tg3-ug3）/3;yb3=（tb3-ub3）/3;
dr=yr1+yr2+yr3;dg=yg1+yg2+yg3;db=yb1+yb2+yb3;
cr=im2uint8（dr）;?cg=im2uint8（dg）;?cb=im2uint8（db）;
z=cat（3crcgcb）;
subplot（121）;imshow（f）;title（‘原圖‘）;
subplot（122）;imshow（z）;title（‘MSR?處理后‘）;%顯示處理后的圖像

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????1627??2019-04-14?23:03??MSR.m
?????文件????????2301??2019-04-13?23:07??MSRCR.m
?????文件????????3823??2019-04-13?02:55??SSR.m