代码MATLAB基于形态学的目标检测一

文章《（来点有用的）MATLAB基于形态学的目标检测（一）简单图形统计》的演示代码。在一些例如裂纹测量、细胞统计的工作中，由于目标可能连接、重叠等情况。直接使用区域增长和连通分量计算，往往结果不好。而通过基本形态学的操作能在一定程度上解决这个问题。

% 基于形态学的目标检测（一）已知尺寸规则图形统计
% 测试用
% by HPC_ZY 20190618

clear; close all; clc

%% 生成测试图
% 生成背景
rows = 500;
cols = 500;
bg1 = zeros（rowscols）;
bg2 = zeros（rowscols）;
% 生成图形
N = 8;
R = 45;
rh = 30;
rw = 40;
xy = zeros（N2）; % 用于记录目标真实坐标
for k = 1:N
    x = round（rows*（rand（）*0.8+0.1））;
    y = round（cols*（rand（）*0.8+0.1））;
    bg1 = CreateCircle（bg1 x y R）;
    bg2 = CreateRectangle（bg2 x y rhrw）;
    xy（k:） = [yx];
end
% 对图形进行从左到右，从上到下排序（为了方便后期对比）
idx = sum（xy（:2）+xy（:1）*rows2）;
[idxorder] = sort（idx）;

%% 用形态学计算个数
% 腐蚀
se1 = CircleStrel（round（R*0.95））;
se2 = RectangleStrel（round（rh*0.95）round（rw*0.95））;
ob1 = imerode（bg1se1）;
ob2 = imerode（bg2se2）;
% 连通标记
[l1num1] = bwlabel（ob1）;
[l2num2] = bwlabel（ob2）;
status1 = regionprops（l1‘BoundingBox‘）;
centroid1 = regionprops（l1‘Centroid‘）;
status2 = regionprops（l2‘BoundingBox‘）;
centroid2 = regionprops（l2‘Centroid‘）;

%% 绘制
figure
subplot（231）imshow（bg1）title（‘\fontsize{14}原图‘）
hold on
for i=1:num1
     text（xy（order（i）1）-1xy（order（i）2）-1 num2str（i）‘Color‘ ‘r‘）
end
subplot（232）imshow（ob1）title（‘\fontsize{14}腐蚀结果‘）
subplot（233）imshow（bg1）title（‘\fontsize{14}标记结果‘）
hold on;
for i=1:n