mhm_090505_2.m

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

clear
I=imread('E:\image\20远-中01.bmp');%读取原始图像
figure,imshow(I);
BW1=edge(I,'sobel',0.035);%边界提取
figure,imshow(BW1)
n=8;
[x,y]=ginput(n);%取点
for x0=10:640
    for y0=10:480
        for R=10:200
            Q(x0,y0,R)=sigma((sqrt((x-x0)^2+(y-y0)^2)-R);
        end
    end
end
x0,y0,R