用matlab检测斑马鱼重心

clc;

clear;

img1=imread('E:\matlab\zebrafish\zebrafish.jpg');

figure('name','显示原始图像');

imshow(img1);

img2 = rgb2gray(img1); %转换成灰度图像

figure('name','显示灰度图像');

imshow(img2);

%thresh = graythresh(img2); %获取转换成二值图像的阈值

img3 = im2bw(img2,0.64); %转换为二值图像,手动调整阈值，保证斑马鱼外轮廓比较准确

figure('name','显示二值图像');

imshow(img3);

img4 = bwareaopen(img3,4300); %删除二值图像img3中面积小于500的对象

figure('name','显示滤除内部噪声后的图像');

imshow(img4);

img5=1-img4;%图像反转

figure('name','显示反转后的图像');

imshow(img5);

img5 =bwareaopen(img5,53000); %去噪，阈值选取很大，其实是要去掉鱼身体外部的痕迹，只留下整个鱼身

figure('name','显示滤除外部噪声及重心标记后的图像');

imshow(img5);

hold on ; %维持图像，保证坐标能被标记在img5上

S=regionprops(img5,'Centroid'); %求出鱼身的重心坐标,坐标信息存在S.Centroid

plot(S.Centroid(1),S.Centroid(2),'+'); %用+标记重心点

result=sprintf(' (%.1f,%.1f)',S.Centroid);

text(S.Centroid(1),S.Centroid(2),result);%将坐标以文本的形式标记在img5上

hold off;

%求出鱼身体部分的像素个数

S=regionprops(img5,'Area');

a=S.Area;

%a=131686