matlab求取光斑中心

clear;
clc;
test_im=imread（‘E:\实验\47cm7月2日\靶-3\Y 5\5 （18）.bmp‘）; %原始图像（黑色背景，白色圆点）
test_im_gray=rgb2gray（test_im）;
[mn]=size（test_im_gray）;
bw（1:m1:n）=0; 
for i=1:m
for j=1:n
if test_im_gray（ij）>10;%二值
bw（ij）=1;
end
end
end
L= bwlabel（bw）;
L=imclearborder（L）;
L=bwareaopen（L10）;%除去多余光点
s = regionprops（L ‘centroid‘‘area‘ ‘boundingbox‘）;%求中心
centroids = cat（1 s.Centroid）
areas = [s.Area];
rects = cat（1  s.BoundingBox）;
figure;
imshow（test_im）;