%?_ANIU?
%?本程序?yàn)殪o態(tài)參數(shù)測量系統(tǒng)m程序?
%?分段調(diào)試用	??
%??
%???
%?本程序不調(diào)用子函數(shù)??
%?2007年3月27日?
???
tic???????????%?Stopwatch?timer?
clc?
clear
?
?
I=imread（‘f:\guang.jpg‘）;?%?讀入圖像?
I=rgb2gray（I）;???????%?轉(zhuǎn)換為灰度圖像?
[RowCol]=size（I）;???%?行數(shù)Row＝576，列數(shù)Col＝768?
Row?
Col?
?
[pixel_numbergray_scale]=imhist（I255）;????%?獲得圖像直方圖，gray_scale為灰度值，pixel_number為該灰度值上的總像素數(shù)?
prob=pixel_number/（Row*Col）;????????????????%?為灰度為某一值的概率?
?
?
ave_gray=0;?????????????????????????%?圖像平均灰度初值為0?
for?i=1:255?
????ave_gray=ave_gray+i*prob（i）;????%?圖像平均灰度?
end?
?
?
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%?
?
%%%%%%?以下幾個值要在實(shí)際環(huán)境中最后確定?
?
?
R0=350;??????%?光靶有效區(qū)域的最大半徑，大于該值為需要灰度重建部分（需要先去掉的背景區(qū)域）?
a=0.5*Row-10;???%?圓心的初始位置?
b=0.5*Col-25;???%?這兩個偏移值，在最終程序中要去掉?
k0=1.0;??%?去掉黑色區(qū)域的灰度調(diào)整因子（用于去掉黑邊＆背景區(qū)域）?
k=0.75;??%?分區(qū)域自適應(yīng)雙閾值分割因子?
k1=0.9;??%?全局自適應(yīng)雙閾值小因子（用于提取同心圓，配合4D）?
k2=1.6;??%?全局自適應(yīng)雙閾值大因子（用于提取激光的光斑）?
D=10;????%?十字絲寬度、圓寬度范圍（用于提取兩個直線標(biāo)尺，圓的寬度）?
?
?
%?經(jīng)驗(yàn)值?
%?5個同心圓的估計初值?
r（1）=135;?
r（2）=180;?
r（3）=225;?
r（4）=270;?
r（5）=310;?
?
?
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%?
?
figure（1）?
imshow（I）;?
title（‘原始圖像‘）;?
?
?
%?圖像預(yù)處理:中值濾波（為了加快速度可將這一步去掉，作用效果也不明顯）?
I2=medfilt2（I）;?
?
figure（2）?
imshow（I2）;?
title（‘中值濾波‘）;?
?
?
I3=I2;?
%?去掉采集圖像時可能形成的黑邊?（去掉多少由具體情況定，一般為3個象素一邊）?
I3（1:）=ave_gray*k0;?
I3（Row:）=ave_gray*k0;?
I3（:1）=ave_gray*k0;?
I3（:Col）=ave_gray*k0;?
?
?
%圖像預(yù)處理:去掉黑色背景區(qū)域?
for?i=1:Row?
????for?j=1:Col??????????????????????
?????????
????????%?去掉黑色背景區(qū)域?
????????if（（i-a）*（i-a）+（j-b）*（j-b））>R0*R0?
????????????I3（ij）=ave_gray*k0;?
????????end?
????end?
end?
?
%?測試去掉黑色背景區(qū)域圖像?
figure（3）?
imshow（I3）;?
title（‘去掉黑色區(qū)域‘）;?
?
?
%%%?圓擬合部分？？？？？?
N=5;??????????????????????%?本測試圖像共有5個圓??
point_number=zeros（1N）;??%?圓圖像上點(diǎn)數(shù)統(tǒng)計值（數(shù)組用于記錄點(diǎn)數(shù)）?
ave_moment=zeros（N9）;????%?平均矩?
?
?
?
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%?
%?全局閾值分割方法獲得同心圓標(biāo)尺圖像（亮度相對小）、光斑圖像（亮度相對大）?
?
I4=I3;??%%%%%%%?獲取同心圓?
I5=I3;??%%%%%%%?獲取激光的光斑?
for?i=1:Row?
????for?j=1:Col?????????
?????????
????????%?獲得同心圓標(biāo)尺圖像時，必須?
????????%?去掉兩個直線標(biāo)尺?
?????????
????????if（I3（ij）4*D）&&（abs（j-b）>4*D）??%%%大于4D去掉兩個直線標(biāo)尺?
????????????I4（ij）=I3（ij）;?
????????else?
????????????I4（ij）=255;?
????????end?
?????????
%?????????%?驗(yàn)證半徑初值?
%?????????l=（i-a）*（i-a）+（j-b）*（j-b）;?
%?????????l=sqrt（l）;?
%?????????if（abs（l-r（5））<0.1*D）?
%?????????????I4（ij）=0;?
%?????????end?
?????????
?
????????%?初次分割的光斑圖像?
?
????????if（I3（ij）>ave_gray*k2）?
????????????I5（ij）=I3（ij）;?
????????else?
????????????I5（ij）=0;?
???????end??
???????%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%?
?
?
?
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%?
?
?????????
????????%?同心圓圖像分區(qū)域擬合的有關(guān)計算?
????????if?（I4（ij）????????????L=（i-b）*（i-b）+（j-a）*（j-a）;?
????????????L=sqrt（L）;????%?計算圓上一點(diǎn)到圓心初值的距離?
?????????????
????????????x1=j;?????????%?計算邊沿xy坐標(biāo)的各階矩，共9個值?
????????????x2=j^2;?
????????????x3=j^3;?
????????????

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件???????15075??2014-09-18?14:54??zhixin.m