目录
一、理论基础
算法分为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
原文地址：http://www.cnblogs.com/matlabfpga/p/16901384.html
