通过matlab实现水产养殖鱼类成熟度自动分析系统

Posted 2022-12-08 fpga和matlab

欢迎订阅《FPGA学习入门100例教程》、《MATLAB学习入门100例教程》

目录

一、理论基础

二、核心程序

三、测试结果

---

一、理论基础

算法分为2个部分：

第一部分是鱼提取：

        首先更加像素判断，将鱼的背景因素去除，这个主要是像素判断，包括背景中的白色背景以及蓝色字体背景。

       这些背景去除之后，将得到一个鱼的大致轮廓。

       然后将提取的鱼进行颜色分解，即RGB转换为HSV，然后使用S通道，进行二值化处理，即，即im2bw这个函数，二值化处理之后，去掉图像中的小的噪声，即bwareaopen这个函数，

然后进行填充处理，即将鱼区域中的小的空洞区域进行填充，

      最后，通过腐蚀操作，将鱼边缘部分一些不平稳的消除。实现鱼的抠图操作。

     对鱼边缘提取和记录计算：

     通过edge这个函数，进行sobel边缘提取。可以提取鱼的边缘轮廓。

有了轮廓之后，进行边缘坐标点的定位。

这里，首先确定鱼嘴巴上的点，然后根据鱼身体的各个部分比例，依次提取各个点的位置。

二、核心程序

 

clc;
clear;
close all;
warning off;
addpath 'func\\'



%读取图片
I0      = imread('images\\test.jpg');


I1      = 1-im2bw(rgb2gray(I0),0.9999);
%去除小的连通域
I2      = bwareaopen(I1,100000);
%腐蚀
se1     = strel('disk',5);
I3      = imerode(I2,se1);
%填充
I3      = imfill(I3,'holes');

figure;
subplot(211);imshow(I0);
subplot(212);imshow(I3);
%边缘提取
I4 = edge(I3,'sobel');


%端点提取和连线
%第80个点基本就是B点
[R,C] = size(I4);
[Y,X] = find(I4==1);
Xb    = X(80);
Yb    = Y(80);
%D
for i = 1:length(X)
    if sqrt((Xb - X(i))^2 + (Yb - Y(i))^2) >= 600 & sqrt((Xb - X(i))^2 + (Yb - Y(i))^2) <= 601 & i<3000
       Xd = X(i);
       Yd = Y(i);
    end
end
%A
for i = 1:length(X)
    if sqrt((Xb - X(i))^2 + (Yb - Y(i))^2) >= 700 & sqrt((Xb - X(i))^2 + (Yb - Y(i))^2) <= 710 & i<3000 & Y(i)>1.2*Yd
       Xa = X(i)+20;
       Ya = Y(i)-150;
    end
end
%F
for i = 1:length(X)
    if sqrt((Xd - X(i))^2 + (Yd - Y(i))^2) >= 1100 & sqrt((Xd - X(i))^2 + (Yd - Y(i))^2) <= 1110 & i > 1000
       Xf = X(i);
       Yf = Y(i);
    end
end
%H
for i = 1:length(X)
    if sqrt((Xf - X(i))^2 + (Yf - Y(i))^2) >= 580 & sqrt((Xf - X(i))^2 + (Yf - Y(i))^2) <= 585 & Y(i)>1.05*Yf & i > 2000
       Xh = X(i);
       Yh = Y(i);
    end
end
%J
for i = 1:length(X)
    if sqrt((Xh - X(i))^2 + (Yh - Y(i))^2) >= 880 & sqrt((Xh - X(i))^2 + (Yh - Y(i))^2) <= 890 & i > 2000
       Xj = X(i);
       Yj = Y(i);
    end
end
%I
for i = 1:length(X)
    if sqrt((Xj - X(i))^2 + (Yj - Y(i))^2) >= 325 & sqrt((Xj - X(i))^2 + (Yj - Y(i))^2) <= 355 & Y(i)>1.2*Yj & X(i) < 1.05*Xj &  X(i) > 0.95*Xj
       Xi = X(i);
       Yi = Y(i);
    end
end

%G
for i = 1:length(X)
    if sqrt((Xi - X(i))^2 + (Yi - Y(i))^2) >= 200 & sqrt((Xi - X(i))^2 + (Yi - Y(i))^2) <= 210 & X(i) < Xj & Y(i)>Yi
       Xg = X(i);
       Yg = Y(i);
    end
end


%E
for i = 1:length(X)
    if sqrt((Xg - X(i))^2 + (Yg - Y(i))^2) >=830 & sqrt((Xg - X(i))^2 + (Yg - Y(i))^2) <= 840 & X(i) < Xg & Y(i)>Yi
       Xe = X(i);
       Ye = Y(i);
    end
end



%C
for i = 1:length(X)
    if sqrt((Xe - X(i))^2 + (Ye - Y(i))^2) >=700 & sqrt((Xe - X(i))^2 + (Ye - Y(i))^2) <= 710 & X(i) < Xe
       Xc = X(i);
       Yc = Y(i)-50;
    end
end

三、测试结果

 

 

A16-62