halcon基础之区域特征分析和区域提取

这是网上写得比较好的一篇贴子,勇哥转载一下以方便大家。

讲的都是基础的region操作。


1、区域的运算

在学习特征分析之前,先了解一下几个区域的基本运算方式。

区域的运算主要包含:区域作差、区域合并、区域集合、区域的选择、区域的填充、区域的骨架等。在项目中,根据实际需要,灵活的选择相应的计算方式。

在Halcon创建两个矩形,作为后面算子测试用:

*生成两个矩形区域
gen_rectangle1 (Rectangle1, 123, 196, 186, 356)
gen_rectangle1 (Rectangle2, 123+20, 196+20, 186-20, 356-20)

区域作差:将一个区域减去另一个区域,便得到区域的差值。

算子:difference 

如:

*两个矩形区域做差
difference (Rectangle1, Rectangle2, RegionDifference)

区域合并(求和):将一个区域和另外一个区域求和,便可以将两个区域合并成一个区域。

算子:union1、union2 

如:

*区域的合并,union1是对集合的合并,union2是对区域的加
union2 (RegionDifference, Rectangle2, RegionUnion)

区域集合:将两个或者两个以上的区域放入一个集合中,这个集合称为区域的集合。特别注意的是:在集合中,每个区域都是孤立的,并没有进行合并,要与区域的和进行区分,每一个区域称为一个obj

算子:concat_obj 

如:

*两个区域放到一个集合concat_obj (Rectangle1, Rectangle2, ObjectsConcat)

区域的选择:从区域集合中进行区域的选择。

算子:select_obj。参数3:要选择的区域ID号。

如:

*区域的选择:参数3表示要选择的第几个区域
select_obj (ObjectsConcat, ObjectSelected, 2)

区域的填充:对区域中的孔洞进行封闭,称为孔洞的填充。

算子:fill_up 

如:

*区域的填充
fill_up (RegionDifference, RegionFillUp)

区域骨架的求取:区域的骨架可以理解为图像的中轴,骨架求取就是求取区域中轴的过程。

算子:skeleton、fill_up_shape

如:

*区域骨架的求取
skeleton (RegionFillUp1, Skeleton)
*指定面积范围的区域的填充:参数4和参数5表示区域面积范围。不在面积范围的不填充fill_up_shape (RegionDifference, RegionFillUp1, 'area', 1, 10000)

 

2、区域的特征分析

在学习区域特征分析前,先了自我问一下,为什么要学习区域特征分析。根据前面所讲的,已经能够进行区域的分割和区域的提取了。这是因为,如果我们只想从分割出来的区域中提取指定要求的一个区域,那该如何选择呢?那这就需要给定提取区域的特征来进行提取。所以要对区域的特征分析,要根据每一个区域的特征,选择出我们需要的区域。

常用的区域区域特征:区域的面积、区域的最小外接矩形(长、宽跟角度)、区域的角度、区域的凸性、区域的圆度、区域的最小外接圆大小等。下面来对其进行具体讲解:

首先绘制一个矩形区域,变量名Region。

区域的面积:区域内像素点数量总。

算子:area_center

area_center (Region, Area, Row4, Column4)

区域的最小外接矩形:外接矩形长、宽的一半长、宽半轴,对应的矩形长轴方向称为区域的方向。区域的最小外接矩形,分为两种外接方式,如下图所示

算子:

smallest_rectangle2:最小外接矩形2

smallest_rectangle1:最小外接矩形1

区域的角度:区域的角度就是外接矩形的角度。

区域的凸性:每个区域都存在一个最小的凸性区域,这个区域内的任意位置两点连线都在该区域内,这样的区域称为区域的凸性区域。其中原区域与凸性区域的比值称为区域的凸性,范围0-1.

算子:

shape_trans:将区域转换为凸性区域(凸包)。

convexity:获取区域的凸性值

如:

*凸集转换,'convex'生成凸集,也可以选择生成外接圆
shape_trans (Region, RegionTrans, 'convex')
*获取区域的凸性值,值保存在Convexity中。表示的是该区域面积和凸集的比值,范围0-1

区域的圆度:指区域边界是否接近圆形的程度,范围0-1。

算子:circularity

*获取区域的圆度circularity (Region, Circularity)

区域的外接圆:能够将区域完全包含的圆区域,称为区域的最小外接圆。

算子:

smallest_circle:获取区域的最小外接圆

region_features :获得区域的任意特征

*获取区域外接圆
smallest_circle (Region, Row3, Column3, Radius)
gen_circle (Circle, Row3, Column3, Radius)
*获取区域的特征,比如面积:'area'region_features (Circle, 'area', Value)

3、案例:遥控器上的字符特征提取

经过前面的阈值分割、形态学膨胀腐蚀、区域特征分析、区域运算等学习,利用这些知识可以对一幅图像中指定的区域进行提取出来。比如对一幅遥控器背面图像的操作说明的文字进行提取。

需求:提取遥控器上的字符区域,并计算每一个字符汉字区域的面积、外接矩形的长宽,然后将计算的值存储到txt文本文件中。

遥控器图片

分析:先将图片转换为灰度图,根据灰度值将遥控器区域提取出来,然后再去求取每个字符的区域,最后再对每一个字符串进行计算区域面积、外接矩形长宽。

 

程序实现讲解:

读取图片并灰度转换

read_image (Image, '素材.jpg')*灰度转换rgb1_to_gray (Image, GrayImage)

通过静态全局阈值分割(根据灰度值不同来分割),将图片分割成不同区域。遥控器的图像区域比较亮,灰度值基本在180以上。灰度值查看方法前面已经说过,在图片窗口上,按住Ctrl键,鼠标跳出来的一个是鼠标坐标值、另一个就是灰度值。

*静态阈值分割threshold (GrayImage, Region, 150, 255)

阈值分割后,每个区域的边缘比较凹凸,形态学开操作(腐蚀)。

*开操作-腐蚀操作opening_circle (Region, RegionOpening, 3.5)

连通域标记操作

*连通域标记,将几个区域分隔开connection (RegionOpening, ConnectedRegions)

将遥控器的区域选择出来,遥控器的区域面积最大,因此可以用面积来进行筛选。

*方法1:遥控器区域面积大,可以用选择面积大的方法。将遥控的区域单独选择出来。
*参数3:最大面积,参数4:百分比。
select_shape_std (ConnectedRegions, SelectedRegions, 'max_area', 70)
*方法2:根据凸性、或者遥控的外接矩形凸性来选择,

注意的是,提取出来的遥控器区域,字符所占的位置都是空的,因为字符是黑色的,在阈值分割的时候都已经分割开来了,因此要获得整个遥控图片所占区域,要将字符串的也位置填充起来。

*填充:将字符串的填充起来fill_up (SelectedRegions, RegionFillUp)

将遥控器图像,从灰度图中裁剪下来

*裁剪,将遥控器的图片区域从原图中裁剪下来reduce_domain (GrayImage, RegionFillUp, ImageReduced)

此时就可以对裁剪下来的遥控器图像进行特征分析了。字符是黑色的,遥控面板是白色的,灰度值区分还是很明显的。使用阈值分割即可

*全局阈值:前两个数值:选择分割的图像灰度范围,最后一个是要选择的分割的面积(可以将一些点去掉)fast_threshold (ImageReduced, Region1, 0, 100, 20)

可以看到,此时提取出来的字符,是斜着的,并且很多字的笔画是分开的并没有连接起来,这样子的话程序是会将其当成多个字符区域。

因此先将这些字符转到水平。

*求取区域的方向
orientation_region (Region1, Phi)
*求取区域的中心
area_center (Region1, Area, Row, Column)
*求取区域的角度,参数6:新的旋转之后的角度
vector_angle_to_rigid (Row, Column, Phi, Row, Column, -3.1415*0.5, HomMat2D)
*对区域旋转操作:转到水平
affine_trans_region (Region1, RegionAffineTrans, HomMat2D, 'nearest_neighbor')

然后,再膨胀操作。先形态学闭运算、再竖直方向膨胀

*膨胀操作-闭运算:
*同一个字符相邻的区域先浓缩成一团:将左右连接在一起
closing_circle (RegionAffineTrans, RegionClosing, 3.5)
*然后再竖直方向膨胀,将竖直方向的都连接起来。参数3是宽度表示水平方向的,则值要设置的小dilation_rectangle1 (RegionClosing, RegionDilation, 1, 11)

此时连通域操作,将每个字标记分开

*连通域操作:connection (RegionDilation, ConnectedRegions1)

观察发现,有一个逗号,还有“约3秒”、“按键”连接在了一起成了一个区域。

先去掉逗号,通过面积选择来操作

*选择指定的区域:去除逗号,,逗号面积小,大概80.参数5和参数6表示的是筛选的区域面积范围select_shape (ConnectedRegions1, SelectedRegions1, 'area', 'and', 150, 99999)

分开连接在一起的字符比较麻烦,按下面的方法来实现。

将Region1旋转至水平和SelectedRegions1求交集才能求出完整的字符

*旋转至水平
affine_trans_region (Region1, RegionAffineTrans1, HomMat2D, 'nearest_neighbor')
*两个区域求交集,得出遥控器上的字符串
intersection (SelectedRegions1, RegionAffineTrans1, RegionIntersection)

然后对每个字符区域进行排序

*区域的排序,参数5:选择排序方式以行排序
sort_region (RegionIntersection, SortedRegions, 'first_point', 'true', 'row')

*“约3秒”、“按键”没有分割开,需要分割开。用区域外接矩形1,长度法提取出来

select_shape (SortedRegions, SelectedRegions2, 'rect2_len1', 'and', 30, 99999)

然后再对其进行分开操作

*提取出来后再闭运算-膨胀
closing_circle (SelectedRegions2, RegionClosing1, 1.5)
*再竖直方向膨胀(使每个字符笔画都是连接的)
dilation_rectangle1 (RegionClosing1, RegionDilation1, 1, 11)
*连通域connection (RegionDilation1, ConnectedRegions2)
*再求交集,获得字符
intersection ( ConnectedRegions2,SelectedRegions2, RegionIntersection1)

*生成空区域
gen_empty_obj (EmptyObject)
*RegionIntersection1求交集EmptyObject,放到生成的空区域中
concat_obj (EmptyObject, RegionIntersection1, EmptyObject)

约3秒”、“按键”连起来的字符区域去掉,再加上“约3秒”、“按键”重新分割开的,形成完整的每个汉字都独立的区域

***“约3秒”、“按键”连起来的字符区域去掉,再加上“约3秒”、“按键”重新分割开的。
*区域做差。
difference (RegionIntersection, RegionIntersection1, RegionDifference)
concat_obj (EmptyObject, RegionDifference, EmptyObject)

至此,遥控器上的每个独立的字符都已经提取出来了。下面要做的就是,每一个字符汉字区域的面积、外接矩形的长宽,然后将计算的值存储到txt文本文件中。

可以先统计总共有多少个字符,然后将字符进行排序,再写一个循环遍历程序对字符区域进行求面积、计算外接矩形的长宽,将所求的结果保存到数组中。最后将保存区域面积、外接矩形长款的数组数组写到txt文本中。

功能比较简单,直接贴出代码程序了,代码都有对应注释。

*统计总共有多少个字符区域
count_obj (EmptyObject, Number)
*按行排序
sort_region (EmptyObject, SortedRegions1, 'first_point', 'true', 'row')
 
*创建一个数组,用来保存区域面积
a:=[]
*创建一个数组,用来保存外接矩形1边长L1
L1:=[]
*创建一个数组,用来保存外接矩形1边长L2
L2:=[]
 
*遍历
for i:=1 to Number-2 by 1
    *
    select_obj (SortedRegions1, ObjectSelected, i)
    *区域面积
    area_center (ObjectSelected, Area1, Row1, Column1)
    *遍历的结果添加到数组中,并且添加tab符
    a:=[a,Area1,'\t']
    *外接矩形1的两个边长
    smallest_rectangle2 (ObjectSelected, Row2, Column2, Phi1, Length1, Length2)
    *遍历的结果添加到数组中,并且添加tab符
    L1:=[L1,Length1,'\t']
    *遍历的结果添加到数组中,并且添加tab符
    L2:=[L2,Length2,'\t']
endfor
 
*添加换行符
a:=[a,'\r\n']
*添加换行符
L1:=[L1,'\r\n']
*添加换行符
L2:=[L2,'\r\n']
 
****将37个字符区域的面积和外接矩形1的两个边长参数写到txt文本中
 
*新建txt。文件保存在程序根目录下
open_file ('1.txt', 'output', FileHandle)
fnew_line (FileHandle)
*写到文本中.面积a
fwrite_string (FileHandle, a)
*另起一行再写。L1
fnew_line (FileHandle)
fwrite_string (FileHandle, L1)
*另起一行再写。L2
fnew_line (FileHandle)
fwrite_string (FileHandle, L2)
*关闭文件
close_file (FileHandle)


最终保存到txt文本的数据

工程完整程序下载链接:

https://download.csdn.net/download/panjinliang066333/12206334


本文出自勇哥的网站《少有人走的路》wwww.skcircle.com,转载请注明出处!讨论可扫码加群:

发表评论:

◎欢迎参与讨论,请在这里发表您的看法、交流您的观点。

会员中心
搜索
«    2024年3月    »
123
45678910
11121314151617
18192021222324
25262728293031
网站分类
标签列表
最新留言
    热门文章 | 热评文章 | 随机文章
文章归档
友情链接
  • 订阅本站的 RSS 2.0 新闻聚合
  • 扫描加本站机器视觉QQ群,验证答案为:halcon勇哥的机器视觉
  • 点击查阅微信群二维码
  • 扫描加勇哥的非标自动化群,验证答案:C#/C++/VB勇哥的非标自动化群
  • 扫描加站长微信:站长微信:abc496103864
  • 扫描加站长QQ:
  • 扫描赞赏本站:
  • 留言板:

Powered By Z-BlogPHP 1.7.2

Copyright Your skcircle.com Rights Reserved.

鄂ICP备18008319号


站长QQ:496103864 微信:abc496103864