细化

细化(thinning)算法有很多,我们在这里介绍的是一种简单而且效果很好的算法,用它就能够实现从文本抽取骨架的功能。我们的对象是白纸黑字的文本,但在程序中为了处理的方便,还是采用256级灰度图,不过只用到了调色板中0和255两项。

所谓细化,就是从原来的图中去掉一些点,但仍要保持原来的形状。实际上,是保持原图的骨架。所谓骨架,可以理解为图象的中轴,例如一个长方形的骨架是它的长方向上的中轴线;正方形的骨架是它的中心点;圆的骨架是它的圆心,直线的骨架是它自身,孤立点的骨架也是自身。文本的骨架嘛,前言中的例子显示的很明白。那么怎样判断一个点是否能去掉呢?显然,要根据它的八个相邻点的情况来判断,我们给几个例子(如图6.22所示)。

图6.22   根据某点的八个相邻点的情况来判断该点是否能删除

图6.22中,(1)不能删,因为它是个内部点,我们要求的是骨架,如果连内部点也删了,骨架也会被掏空的;(2)不能删,和(1)是同样的道理;(3)可以删,这样的点不是骨架;(4)不能删,因为删掉后,原来相连的部分断开了;(5)可以删,这样的点不是骨架;(6)不能删,因为它是直线的端点,如果这样的点删了,那么最后整个直线也被删了,剩不下什么;(7)不能删,因为孤立点的骨架就是它自身。

总结一下,有如下的判据:(1)内部点不能删除;(2)孤立点不能删除;(3)直线端点不能删除;(4)如果P是边界点,去掉P后,如果连通分量不增加,则P可以删除。

我们可以根据上述的判据,事先做出一张表,从0到255共有256个元素,每个元素要么是0,要么是1。我们根据某点(当然是要处理的黑色点了)的八个相邻点的情况查表,若表中的元素是1,则表示该点可删,否则保留。

查表的方法是,设白点为1,黑点为0;左上方点对应一个8位数的第一位(最低位),正上方点对应第二位,右上方点对应的第三位,左邻点对应第四位,右邻点对应第五位,左下方点对应第六位,正下方点对应第七位,右下方点对应的第八位,按这样组成的8位数去查表即可。例如上面的例子中(1)对应表中的第0项,该项应该为0;(2)对应37,该项应该为0;(3)对应173,该项应该为1;(4)对应231,该项应该为0;(5)对应237,该项应该为1;(6)对应254,该项应该为0;(7)对应255,该项应该为0。

这张表我已经替大家做好了,可花了我不少时间呢!

static int erasetable[256]={

0,0,1,1,0,0,1,1,          1,1,0,1,1,1,0,1,

1,1,0,0,1,1,1,1,             0,0,0,0,0,0,0,1,

0,0,1,1,0,0,1,1,             1,1,0,1,1,1,0,1,

1,1,0,0,1,1,1,1,             0,0,0,0,0,0,0,1,

1,1,0,0,1,1,0,0,             0,0,0,0,0,0,0,0,

0,0,0,0,0,0,0,0,             0,0,0,0,0,0,0,0,

1,1,0,0,1,1,0,0,             1,1,0,1,1,1,0,1,

0,0,0,0,0,0,0,0,             0,0,0,0,0,0,0,0,

0,0,1,1,0,0,1,1,             1,1,0,1,1,1,0,1,

1,1,0,0,1,1,1,1,             0,0,0,0,0,0,0,1,

0,0,1,1,0,0,1,1,             1,1,0,1,1,1,0,1,

1,1,0,0,1,1,1,1,             0,0,0,0,0,0,0,0,

1,1,0,0,1,1,0,0,             0,0,0,0,0,0,0,0,

1,1,0,0,1,1,1,1,             0,0,0,0,0,0,0,0,

1,1,0,0,1,1,0,0,             1,1,0,1,1,1,0,0,

1,1,0,0,1,1,1,0,             1,1,0,0,1,0,0,0

};

有了这张表,算法就很简单了,每次对一行一行的将整个图象扫描一遍,对于每个点(不包括边界点),计算它在表中对应的索引,若为0,则保留,否则删除该点。如果这次扫描没有一个点被删除,则循环结束,剩下的点就是骨架点,如果有点被删除,则进行新的一轮扫描,如此反复,直到没有点被删除为止。

实际上,该算法有一些缺陷。举个简单的例子,有一个黑色矩形,如图6.23所示。

图6.23经过细化后,我们预期的结果是一条水平直线,且位于该黑色矩形的中心。实际的结果确实是一条水平直线,但不是位于黑色矩形的中心,而是最下面的一条边。

为什么会这样,我们来分析一下:在从上到下,从左到右的扫描过程中,我们遇到的第一个黑点就是黑色矩形的左上角点,经查表,该点可以删。下一个点是它右边的点,经查表,该点也可以删,如此下去,整个一行被删了。每一行都是同样的情况,所以都被删除了。到了最后一行时,黑色矩形已经变成了一条直线,最左边的黑点不能删,因为它是直线的端点,它右边的点也不能删,因为如果删除,直线就断了,如此下去,直到最右边的点,也不能删,因为它是直线的右端点。所以最下面的一条边保住了,但这并不是我们希望的结果。

解决的办法是,在每一行水平扫描的过程中,先判断每一点的左右邻居,如果都是黑点,则该点不做处理。另外,如果某个黑点被删除了,那么跳过它的右邻居,处理下一个点。这样就避免了上述的问题。

图6.23  黑色矩形

图6.24  图6.23细化后的结果

解决了上面的问题,我们来看看处理后的结果,如图6.24所示。这次变成一小段竖线了,还是不对,是不是很沮丧?别着急,让我们再来分析一下:在上面的算法中,我们遇到的第一个能删除的点就是黑色矩形的左上角点;第二个是第一行的最右边的点,即黑色矩形的右上角点;第三个是第二行的最左边的点;第四个是第二行的最右边的点;……;整个图象处理这样一次后,宽度减少2。每次都是如此,直到剩最中间一列,就不能再删了。为什么会这样呢?原因是这样的处理过程只实现了水平细化,如果在每一次水平细化后,再进行一次垂直方向的细化(只要把上述过程的行列换一下),就可以了。

这样一来,每处理一次,删除点的顺序变成:(先是水平方向扫描)第一行最左边的点;第一行最右边的点;第二行最左边的点;第二行最右边的点;……最后一行最左边的点;最后一行最右边的点;(然后是垂直方向扫描)第二列最上边的点(因为第一列最上边的点已被删除);第二列最下边的点;第三列最上边的点;第三列最下边的点;……倒数第二列最上边的点(因为倒数第一列最上边的点已被删除);倒数第二列最下边的点。我们发现,刚好剥掉了一圈,这也正是细化要做的事。实际的结果也验证了我们的想法。

细化(thinning)相关推荐

  1. 图像的膨胀与腐蚀、细化

    转自:http://www.opencvchina.com/forum.php?mod=viewthread&tid=1124 原理:在特殊领域运算形式--结构元素(Sturcture Ele ...

  2. 只腐蚀毛刺 腐蚀算法_图像的腐蚀 膨胀及细化

    转载:http://blog.sina.com.cn/s/blog_4bdb170b0100sgtj.html 今天所讲的内容属于一门新兴的学科:数学形态学(Mathematical Morpholo ...

  3. 腐蚀、膨胀、细化算法

    今天所讲的内容属于一门新兴的学科:数学形态学(Mathematical Morphology).说起来很有意思,它是法国和德国的科学家在研究岩石结构时建立的一门学科.形态学的用途主要是获取物体拓扑和结 ...

  4. 图像模式识别 (五)

    图像腐蚀.膨胀.细化 这三个内容属于数学形态学(Mathematical Morphology),它是法国和德国的科学家在研究岩石结构时建立的一门学科.形态学的用途主要是获取物体拓扑和结构信息,它通过 ...

  5. 数字图像处理(6)——形态学图像处理

    数字图像处理(6)--形态学图像处理 文章目录 数字图像处理(6)--形态学图像处理 1 基本概念 2 二值形态学基本运算 2.1 膨胀(dilation) 2.2 腐蚀(erosion) 2.3 开 ...

  6. OpenCV4 快速入门 (学习笔记 全)

    第1章 基础知识 1.1 基础结构介绍 作者博客https://blog.csdn.net/shuiyixin?type=blog https://blog.csdn.net/shuiyixin/ar ...

  7. 【文献阅读】 Sorghum segmentation by skeleton extraction

    目录 A 个人 B 文献 摘要 1 介绍introduction 2 相关工作related work 3 概述overview 4 骨架化 4.1 thinning 4.2骨架修剪 5 分割 6 B ...

  8. 冈萨雷斯《数字图像处理》学习笔记(九)形态学图像处理

    目录 序言 一.集合论基础 1.集合的子集与相等 2.集合的基本运算符 3.集合的反射和平移 二.二值形态学的基本运算 1.腐蚀 2.膨胀 3.膨胀与腐蚀运算的对偶性 3.开运算与闭运算 4.一些基本 ...

  9. 图像的膨胀与腐蚀--确认这个大牛的才是正确的

    网上找到很多个版本,经过自己分析对比概念,发现有些说自己方法是正确的,但是我觉得这个大牛的分析才是正确的. 图像的膨胀与腐蚀.细化 原理:在特殊领域运算形式--结构元素(Sturcture Eleme ...

最新文章

  1. IIS7中的站点、应用程序和虚拟目录详细介绍
  2. python工具之myql数据库操作
  3. python调用电脑蜂鸣器一直响_电脑开机蜂鸣器一直响,9声滴~,什么问题?怎么处理?...
  4. python二维散点分布图_深入理解皮尔逊相关系数amp;python代码
  5. 云服务器mqtt协议,云服务器mqtt协议
  6. 文件字符输入流 java
  7. ADO.net中临时表DataTable的简单应用
  8. Android TV开发总结(一)构建一个TV app前要知道的事儿
  9. SecurityError: Blocked a frame with origin from accessing a cross-origin frame
  10. ISO14443、15693、18000体系分析
  11. (黑科技)超链接记录贴
  12. Python RPM包制作
  13. Fluent:如何用Custom Field Function Calculater替代UDMI提高收敛性
  14. 「数据架构」数据模型,数据字典,数据库模式 和ERD的比较
  15. java零基础Ⅰ-- 1.java 概述
  16. pc控制android源码,pc 通过vnc控制android
  17. 寻找你的热情(1)——检查自己所处的位置
  18. 主流视频编码压缩技术基本概念(一)
  19. H5实现无插件视频监控按需直播
  20. 国内唯一!腾讯iOA被权威机构报告列入竞争者能力象限

热门文章

  1. VMBox CentOS安装记录
  2. 哔哩哔哩2020校招原题 - 复数乘法
  3. PHP WebShell源码下载(1)
  4. 梦想在远方,理想在路上
  5. java.lang.IllegalArgumentException: interface UserMapper is not visible from class loader
  6. 【20200422】编译原理课程课业打卡十七之求解文法FirstVTLastVT构造文法算符优先关系表
  7. win10锁屏后默认1分钟进入睡眠状态的解决办法
  8. 注册表写入二进制数据
  9. python123英文字符的鲁棒_Robust Locally Weighted Regression 鲁棒局部加权回归 -R实现
  10. php取雅加达时间,2018年雅加达亚运会电竞赛程表 8月电子竞技比赛时间一览