数学形态学应用——腐蚀、膨胀、击中（二）

最新推荐文章于 2022-05-17 19:46:01 发布

Rjex

最新推荐文章于 2022-05-17 19:46:01 发布

阅读量478

点赞数

分类专栏：数字图像处理

本文链接：https://blog.csdn.net/qq_42746999/article/details/102896916

版权

数字图像处理专栏收录该内容

12 篇文章 7 订阅

订阅专栏

击中——在多个多边形图中找到所要求的图形。

目标物如下图所示，找到较大正方形D的位置，令每种形状的原点位于其重心处。
在这里插入图片描述

过程

(1).假设其中长方形为C，小正方形为E，则A为C，D，E的并集，Ac为A的补集。设D被一个小窗口包围。
(2).D对A的腐蚀如下左图所示,(W-D)对Ac的腐蚀如下右图所示：

在这里插入图片描述
以上两图的交集如下所示，该交集显示了D的原点的位置。

将上图中D的位置放大，如下图所示：

代码

%write by 长安 Rjex
I = zeros(300,400);
I(21:200,21:140) = 1;  %矩形的宽为120，高为180
I(151:251,201:301) = 1;  %大正方形的边长为100
I(61:121,321:381) = 1;  %小正方形的边长为60
figure(1),imshow(I);
[m n]=size(I);

k=101;
t = (k-1)/2;
D = zeros(k);
K = zeros(k);
I2 =zeros(m,n);
for i = t+1 : m-t
    for j = t+1 : n-t
        K = D+I(i-t:i+t,j-t:j+t);
        M = find(K == 0);
        if isempty(M) == 1
            I2(i,j) = I(i,j);
        end
    end
end
figure(2);
imshow(I2);

I3 = 1-I;
I4 = zeros(m,n);
T = zeros(k+20);
T(11:k+10,11:k+10)=-2;
k = 121;
t = (k-1)/2;
for i = t+1 : m-t
    for j = t+1 : n-t
        K =T+I3(i-t:i+t,j-t:j+t);
        M = find(K == 0);
        if isempty(M) == 1
            I4(i,j) = I(i,j);
        end
    end
end
figure(3);
imshow(I4);
figure(4);
I5 = I2&I4;
imshow(I5);