im?=?imread（‘myFig.jpg‘）;?%?讀圖‘myFig.jpg‘改成你的圖片名.

figGray?=?double（rgb2gray（im））;?%?轉(zhuǎn)灰度圖

%?假設(shè)高斯型直方圖?設(shè)置亮斑閾值

mR?=?mean（figGray（:））;?

varR?=?var（figGray（:））;

thr?=?mR+5*sqrt（varR）;?%?5?sigma

%?返回亮斑平均中心坐標(biāo)

[row?col]?=?find（figGray>thr）;

mr?=?mean（row）;?%?行坐標(biāo)?

mc?=?mean（col）;??%?列坐標(biāo)

%text（1.50.3‘cos（x）‘）?

%

figure（10）;imshow（im）;?

hold?on

plot（mc?mr?‘+g‘?‘markersize‘?15?‘linewidth‘2）;%?畫圖看一下
legend（‘中心點(diǎn)坐標(biāo)‘）;
xlabel（mr）;
ylabel（mc）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件?????????539??2018-12-18?11:42??rec.m
?????文件?????????142??2018-12-18?11:42??README.md