分水岭算法应用---连接断开的线圈

    前言,帮助好友解决如何讲图片中断开的线圈连接起来的问题。首先的想法试利用线圈断裂处断点的8领域,断点的八领域不为0的像素点应该只有一个,利用这个方法将所有的断点找出来,再依次将相邻的断点连接起来。这个想法有个问题,在复杂一点的情况下,断点不能通过这个方法寻找,比如在线圈断裂处中间有个孤立的像素,这个像素的8领域的像素值全为0,所以这个断点通过8领域是找不出来的;于是我想到通过分水岭算法来拟合断开的线圈;代码如下:

clc;
clear all;
close all;
I=imread(['C:\Users\wlh\Desktop\png\4_5','.png']);
if size(I, 3) == 1
    gray = I;
elseif size(I, 3) == 3
    gray = rgb2gray(I);
else
    gray = I(:, :, 1);
end
figure(1);imshow(gray);title('Original ');
se1 = strel('diamond',1);          %定义腐蚀的操作的核;
bw = gray > 0;  %黑白图片经测试tmp取80左右,红色图片经测试tmp取20~27左右  
%将tmp中像素值大于50的置为1,小于50的点置0;(黑色为0,白色为1)
dist = 255 - bwdist(bw);   %非零元素与其最靠近的零元素之间的距离
mask = imextendedmin(dist,6);
%BW = imextendedmin(I,h)计算扩展最小值变换,它是H-minima变换的区域最小值。
%区域最小值是连通像元的分量,具有恒定的强度值,并且其外部边界像素都具有较高的值。 
%h是一个非负的标量。
dist2 = imimposemin(dist, mask);
%imimposemin(I,BW)使用形态重构来修改强度图像I,所以在BW不为零的地方只有区域最小值。
%BW是一个与I相同大小的二进制图像。
bd = watershed(dist2) == 0;
bd = bwareaopen(bd, 5);
bd = imdilate(bd, se1);
[m,n] = size(bd);       
bd([1 m-1], :) = 0;%将分割后的图像与边界的链接处给删除;
bd(:, [1 n-1]) = 0;  
bd = watershed(bd) == 0;
bd = double(bd);
figure(2);imshow(bd,[]);title('bd');

结果展示:(前后结果图是一样的,上传问题)

原图.jpg
结果.jpg
原图1.jpg
结果1.jpg
原图2.jpg
结果2.jpg

    喜欢就点个赞吧 ! (* ^ _ ^ *) ! ! ! !

©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

推荐阅读更多精彩内容

  • 选择题部分 1.(),只有在发生短路事故时或者在负荷电流较大时,变流器中才会有足够的二次电流作为继电保护跳闸之用。...
    skystarwuwei阅读 13,564评论 0 7
  • 专业考题类型管理运行工作负责人一般作业考题内容选项A选项B选项C选项D选项E选项F正确答案 变电单选GYSZ本规程...
    小白兔去钓鱼阅读 9,070评论 0 13
  • 国家电网公司企业标准(Q/GDW)- 面向对象的用电信息数据交换协议 - 报批稿:20170802 前言: 排版 ...
    庭说阅读 11,242评论 6 13
  • Swift1> Swift和OC的区别1.1> Swift没有地址/指针的概念1.2> 泛型1.3> 类型严谨 对...
    cosWriter阅读 11,161评论 1 32
  • 混沌世界,人魔混战。 人军首领嬴,率人族众将,弑魔族大军,胜利在望。 嬴有爱人婵,至爱至亲。 魔族末将莽濒临绝境,...
    不会飞的熊_cf85阅读 535评论 2 8