img=imread（‘D:\My?Documents\桌面\紅外小胖哈哈舞\600.jpg‘）;
img_gray=rgb2gray（img）;
figure?（1）;imshow（img_gray）;
my_hist=imhist（img_gray）;

figure（3）;plot（my_hist（1:255）‘）;????%直方圖

L=0;???%統(tǒng)計灰度級數(shù)
A1=90;A2=90;???%%%分段?低?常?高溫區(qū)域
height=240;
width?=320;
s_num=height?*?width;
T=60;?????%5閾值大小
max_gray=256;??%圖像最高灰度等級


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%灰度標記
ff_x=zeros（1256）;
?for?i=1:256
??if（my_hist（i1）>T?|?my_hist（i1）==T）
????ff_x（1i）=1;
????L=L+1;
??end
?end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%計算s1_can、s2_can、s3_can
s1=0;s2=0;s3=0;
for?xx=1:A1
????s1=s1+my_hist（xx）;
end

for?yy=（A1+1）:A2
????s2=s2+my_hist（yy）;
end

for?zz=（A2+1）:max_gray
????s3=s3+my_hist（zz）;
end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%
%模式判斷
%%%%%%%%%%%%%%%%%%%%%%%%%%?
s1_can=s1/s_num;
s2_can=s2/s_num;
s3_can=s3/s_num;

%?if（s1_can>0.8）
%???為低溫目標
%?end
%?
%?if（s2_can>0.8）
%???為常溫目標
%?end
%?
%?if（s3_can>0.8）
%???為高溫目標
%?end
%?
%?if（s1_can+s2_can>0.8）
%???為低溫和常溫目標
%?end
%?
%?if（s2_can+s3_can>0.8）
%???為常溫和高溫目標
%?end
%?
%?if（s3_can+s1_can>0.8）
%???為高溫和低溫目標
%?end
%s1_can+s2_can+s3_can=1;
%%%%%%%%%%%%%%%%%%%%%%%%%

ref_k=zeros（1max_gray+1）;
for?j=1:256
??for?l=1:j
?????ref_k（1j）=ff_x（1l）+ref_k（1j）;
??end
end


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%灰度的等距離排列
his_out=zeros（1256）;
for?m=1:256
???his_out（1m）=255*ref_k（1m）/L;
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
I=zeros（240320）;???%%初始化

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

for?iii=1:256
???I_dengji（iii）=iii;
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

for?ii=1:height
???for?jj=1:width

??????img_gray（iijj）=1+img_gray（iijj）;???%%表示等級


?????if（ref_k（1I_dengji（img_gray（iijj）））<（A1）?|?ref_k（1I_dengji（img_gray（iijj）））==A1）
????????I（iijj）=s1_can*255*ref_k（1I_dengji（img_gray（iijj）））/A1;
????????
?????elseif（ref_k（1I_dengji（img_gray（iijj）））>A1?&?ref_k（1I_dengji（img_gray（iijj）））==A2）
????????I（iijj）=（s2_can*255*（ref_k（1I_dengji（img_gray（iijj）））-A1）/（A2-A1））+s1_can*255;
????????
?????elseif（ref_k（1I_dengji（img_gray（iijj）））>A2）
????????I（iijj）=（s3_can*255*（I_dengji（img_gray（iijj））-A2）/（ref_k（1max_gray）-A2））+（1-s3_can）*255;
?????end
?????
???end
end?
?I1=round（I）;
?figure（2）;
?imshow（uint8（I1））;
??
?my_hist_1=imhist（uint8（I1））;
figure（4）;
plot（my_hist_1（1:255））;????%直方圖
%imshow（my_hist_1）;????

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

?????文件???????2526??2013-07-10?17:31??yizhonghongwaizishiyingduibiduzengqiang.m

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

?????????????????2526????????????????????1