mhm_090505_1.m

来自「利用重心法找到图像中光斑的中心」· M 代码 · 共 11 行

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);