clear
clc

for?n=1:36
k=sprintf（‘%d.tif‘n）;
image=imread（k）;


rgb_r=image（::1）;
rgb_g=image（::2）;
rgb_b=image（::3）;
[widthheight]=size（rgb_r）;
x1=width*height;

xx_R_m=medfilt2（rgb_r[35]）;%中值濾波
xx_G_m=medfilt2（rgb_g[35]）;
xx_B_m=medfilt2（rgb_b[35]）;
?
sum_R=sum（sum（xx_R_m））;%求和RGB??
sum_G=sum（sum（xx_G_m））;
sum_B=sum（sum（xx_B_m））;
?
mean_R=sum_R/x1;%RGB平均值
mean_G=sum_G/x1;
mean_B=sum_B/x1;

hsi=zeros（widthheight3）;?
hsi_h=hsi（::1）;?
hsi_s=hsi（::2）;?
hsi_i=hsi（::3）;?
PI=3.1415956;?
for?i=1:width?
?????for?j=1:height?
?????????R=double（rgb_r（ij））;G=double（rgb_g（ij））;B=double（rgb_b（ij））;?
????????min?=?R;?
if?（G????min?=?G;?
?????end?
if?（B?