%%%%%%質心求取程序

%
function?[plot_xplot_y]=calculate_zhixin（I）
?I=imread（‘原圖.png‘）;%讀取原圖像??
figure（1）;imshow（I）;%顯示原圖像??
??
I_gray=rgb2gray（I）;%原圖像變?yōu)榛叶葓D像??



I=imcomplement（I_gray）;
T=imtophat（Istrel（‘disk‘10））;
I_bw=~im2bw（Tgraythresh（T））;
figure（2）imshow（I_bw）





%level=graythresh（I_gray）;%計算圖像I_gray的全局閾值，level為標準化灰度值，其范圍為[0?1]??
[heightwidth]=size（I_gray）;%計算灰度圖像的長寬??
%I_bw=im2bw（I_graylevel）;%im2bw使用閾值level將灰度圖像轉換為二值圖像．??
%figure（2）;imshow（I_bw）;%顯示二值圖像（背景為白色）??
??


for?i=1:height?%%循環(huán)中進行反色??
????for?j=1:width?????
????????if?I_bw（ij）==1????????
????????????I_bw（ij）=0;????
????????else?I_bw（ij）=1;???
????????end??
????end??
end??
figure（3）;imshow（I_bw）;%顯示取反后的二值圖像（背景為黑色）??


??
[Lnum]=bwlabel（I_bw8）;%bwlabel標注二值圖像I_bw中的目標物體，返回標識矩陣Ｌ和I_bw中目標物體的數(shù)量num，８表示連通數(shù)．??
plot_x=zeros（num1）;%%zeros（mn）產(chǎn)生m×n的全0矩陣．用于記錄質心位置的橫坐標??
plot_y=zeros（num1）;%zeros（mn）產(chǎn)生m×n的全0矩陣．用于記錄質心位置的縱坐標??
??
for?k=1:num??%%num個區(qū)域依次統(tǒng)計質心位置??????
????sum_x=0;????s

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2015-04-09?14:25??圖像校正\
?????文件????????1876??2015-03-27?17:23??圖像校正\calculate_zhixin.asv
?????文件????????1874??2015-03-31?15:23??圖像校正\calculate_zhixin.m
?????文件????????1782??2015-03-26?21:01??圖像校正\calcultate_zhixin.asv
?????文件????????1300??2015-03-26?21:01??圖像校正\fengezuobiao.asv
?????文件????????1402??2015-03-26?21:01??圖像校正\fengezuobiao.m
?????文件????????2001??2015-03-27?17:23??圖像校正\image_huifu.asv
?????文件????????1961??2015-03-26?22:48??圖像校正\image_huifu.m
?????文件????????4860??2015-03-26?21:01??圖像校正\main.asv
?????文件????????5083??2015-03-31?15:24??圖像校正\main.m
?????文件???????23183??2015-03-26?21:01??圖像校正\原圖.png
?????文件???????26913??2015-03-26?21:01??圖像校正\畸變圖.png