mhm_090505_1.m

利用重心法找到图像中光斑的中心

clear
I=imread('E:\image\20远-中01.bmp');%读取原始图像
figure,imshow(I);
BW1=edge(I,'sobel',0.035);%边界提取
figure,imshow(BW1)
n=8;
[xdata,ydata]=ginput(n);%取点
x0=[0,0,100];
f=@(x)sqrt((xdata-x(1)).^2+(ydata-x(2)).^2)-x(3);
x=lsqnonlin(f,x0);