CN109116002B - 一种路用粗集料棱角性量化方法 - Google Patents

一种路用粗集料棱角性量化方法 Download PDF

Info

Publication number
CN109116002B
CN109116002B CN201810667800.7A CN201810667800A CN109116002B CN 109116002 B CN109116002 B CN 109116002B CN 201810667800 A CN201810667800 A CN 201810667800A CN 109116002 B CN109116002 B CN 109116002B
Authority
CN
China
Prior art keywords
aggregate
image
skeleton
points
edge angle
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201810667800.7A
Other languages
English (en)
Other versions
CN109116002A (zh
Inventor
沙爱民
李伟
孙朝云
岳宝峰
李菲菲
郝雪丽
李滢滢
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Changan University
Original Assignee
Changan University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Changan University filed Critical Changan University
Priority to CN201810667800.7A priority Critical patent/CN109116002B/zh
Publication of CN109116002A publication Critical patent/CN109116002A/zh
Application granted granted Critical
Publication of CN109116002B publication Critical patent/CN109116002B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/42Road-making materials

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Food Science & Technology (AREA)
  • Medicinal Chemistry (AREA)
  • Physics & Mathematics (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Image Analysis (AREA)
  • Image Processing (AREA)

Abstract

本发明公开了一种路用粗集料棱角性量化方法,包括输入集料棱角性检测系统得到的集料图像,集料图像灰度化等一系列处理,读取集料图像,采用M结构元素进行形态学集料骨架提取,得到集料边缘棱角的全部信息;初始化各个变量,形态学运算获取集料颗粒骨架图像;读取集料颗粒骨架图像,形态学闭运算进行集料骨架图像像素连接操作,并求骨架图像的大小利用形态学操作移除目标边界像素,组成图像的骨架,将集料骨架图像二值化,对每幅集料骨架图像棱角点数求和;以平均棱角点表征单个集料的棱角性,用平均棱角点和面积权值共同量化整批集料的棱角性。可实现对集料棱角性高效,精确的检测和定量分析,具有检测效率高,客观性强,检测成本低的优点。

Description

一种路用粗集料棱角性量化方法
技术领域
本发明属于道路工程领域,涉及集料特性量化,具体涉及一种路用粗集料棱角性量化方法。
背景技术
在新增高速道路网建设中,沥青混凝土路面占到了90%以上,为了延长路面的寿命、提高使用性能,需要改善混合料的质量,人们经常只从混合料的级配组成和沥青作为切入点,对集料特性量化研究相对较少。集料形态特性是影响混合料性能的重要因素之一,主要包括尺寸、形状、棱角性和表面纹理等。集料尺寸的分布影响沥青混合料耐久度与强度特性;集料的形状影响混合料的稳定性与工作性;集料棱角性影响路面结构的嵌挤作用、抗滑、抗车辙特性;集料表面纹理影响沥青路面结构的抗剪特性。粗集料的棱角性体现了集料表面的突出特性,棱角性较好的集料可以提高颗粒与颗粒之间的嵌挤结构进而提高混凝土路面的路用性能。我国现行规范中,对于粗集料的棱角性还没有统一的评价方法和定量指标,并且其相关研究较少,目前粗集料棱角性量化分析技术还处于研究阶段,现有方法主要存在以下几方面问题:
1)检测技术落后
由于大部分现有集料棱角性检测方法采用全人工或半人工手动经验式检测,检测设备庞大且机动性较差,检测过程繁杂、耗时、费力。
2)检测精度差
现有棱角性检测方法大多数是定性分析而非定量分析,检测结果受人员水平和操作习惯影响较大,采集的数据精度不高并且记录过程繁琐,工程应用局限性高。
3)规范性不一
美国和欧洲对粗集料的棱角性有各自的评价方法但大部分是定性指标,目前,我国尚无明确的集料棱角性规范和定量指标,因人而异,按照不同的标准评价出的集料棱角性指标千差万别,直接导致沥青混合料选料标准混乱、质量不高,沥青路面出现早期路用性能下降或破坏现象。
发明内容
针对上述现有技术存在的缺陷或不足,本发明的目的在于,提供一种路用粗集料棱角性量化方法,从而实现对集料棱角性高效,精确的检测和定量分析。
为了实现上述任务,本发明采取如下的技术解决方案予以实现:
一种路用粗集料棱角性量化方法,其特征在于,按以下步骤实施:
步骤一,输入集料棱角性检测系统得到的集料图像;
步骤二,将集料图像灰度化、中值滤波去噪、高亮细节增强、图像熵分割、闭运算目标提取、坑洞填充等一系列处理;
步骤三,读取步骤二中的集料图像,采用M结构元素进行形态学集料骨架提取,得到集料边缘棱角的全部信息;
步骤四,初始化各个变量,按照式
Figure BDA0001708198950000021
对目标物体进行形态学运算,判断当前Corr_TXi是否为空,若未空,则循环停止,并获取集料颗粒骨架图像;
步骤五,读取集料颗粒骨架图像,形态学闭运算进行集料骨架图像像素连接操作,并求骨架图像的大小x,y,即为行数和列数;从第二行逐个扫描集料骨架点的像素值大小,判断八邻域像素是否符合第n列的一个点,周围邻域中,若它的m+1行,n+1列和m-1行,n-1列所在点的像素值为1,或者它的m+1行,n-1列和m-1行,n+1列所在点的像素值为1,若符合,则不需要断点修补,同时去除闭运算操作后的一些干扰点;否则需要断裂点的修补,将像素点0替换成1进行骨架连通改善;扫描完集料图像像素的所有点,循环终止;
步骤六,读取子文件夹里面所有的修补完的集料骨架图像,利用形态学操作移除目标边界像素,但不允许目标分隔开,保留下来的像素组成图像的骨架;将集料骨架图像二值化,在四邻域提取边界点,对集料骨架二值化图像进行形态学操作,提取骨架棱角点,并去除图像中的孤立点;
步骤七,填充集料骨架二值化图像中的空洞区域,并对每幅集料骨架图像棱角点数求和;
步骤八,以平均棱角点表征单个集料的棱角性,用平均棱角点和面积权值共同量化整批集料的棱角性,将骨架棱角点法棱角性量化指标定义如下:
Figure BDA0001708198950000031
Figure BDA0001708198950000032
Figure BDA0001708198950000033
由于每个集料是在三种不同状态下得到的二维图形,式中,i代表单个集料的不同三种状态,si为集料在其中某一个状态下的棱角点,P为单个集料的平均棱角点数;ai为集料在其中某一个状态下的集料面积,A为单个集料的平均面积;Ia为整批集料棱角性量化指标,N为总集料个数。
本发明的路用粗集料棱角性量化方法,与现有技术相比,具有以下技术效果:
(1)采用的骨架棱角点量化方法可以准确、快捷、高效地量化集料的棱角性,从传统的定性分析转为定量分析。
(2)将每个集料的3个特征面分别进行颗粒棱角特征参数表征并将每个面的面积引入进行加权,得到加权后的集料棱角指标更加合理、准确,更能突出集料的形态宏观特性。
(3)降低了检测算法的复杂度、运行时间短,可用于集料棱角性的实时检测,有利于集料质量好坏的监测和过程控制。
(4)克服了传统检测方法受人员水平和操作习惯影响的缺陷,同时基于骨架棱角性量化方法开发的集料棱角性检测系统,具有检测效率高,客观性强,检测成本低的优点。
附图说明
图1为结构元素示意图,其中(a)图表示L结构元素,(b)图表示M结构元素;
图2为本发明的路用粗集料棱角性量化方法得到集料骨架棱角点提取图。
下面结合附图和具体实施例对本发明作进一步详细地解释和说明。
具体实施方式
参见图1所示,申请人研究发现,L结构元素图1的(a图)进行细化骨架提取,得到的集料骨架图像边缘信息不足,表征集料棱角性不准确,而图1的(b图)M结构元素进行形态学集料骨架提取,得到的骨架图像边缘信息丰富,能够代表集料边缘棱角的全部信息。因此,本申请采用M结构元素进行形态学集料骨架提取。
参见图2,本实施例给出一种路用粗集料棱角性量化方法,具体包括如下步骤:
步骤一,输入集料棱角性检测系统得到的集料图像。
步骤二,将集料图像灰度化、中值滤波去噪、高亮细节增强、图像熵分割、闭运算目标提取、坑洞填充等一系列处理。
步骤三、读取步骤二中的集料图像,选择如图1(b)图的M结构元素(se)进行形态学集料骨架提取,得到集料边缘棱角的全部信息;
步骤四,初始化各个变量,按照式
Figure BDA0001708198950000041
对目标物体进行形态学运算,判断当前Corr_TXi是否为空,若未空,则循环停止,并获取骨架图像。
步骤五,读取集料颗粒骨架图像,形态学闭运算进行集料骨架图像像素连接操作,并求骨架图像的大小x,y,即为行数和列数。从第二行逐个扫描集料骨架点的像素值大小,判断八邻域像素是否符合第n列的一个点,周围邻域中,若它的m+1行,n+1列和m-1行,n-1列所在点的像素值为1,或者它的m+1行,n-1列和m-1行,n+1列所在点的像素值为1,若符合,则不需要断点修补,同时去除闭运算操作后的一些干扰点;否则需要断裂点的修补,将像素点0替换成1进行骨架连通改善。扫描完集料图像像素的所有点,循环终止。
步骤六,读取子文件夹里面所有的修补完的集料骨架图像。利用形态学操作移除目标边界像素,但不允许目标分隔开,保留下来的像素组成图像的骨架。集料骨架图像二值化,在四邻域提取边界点。对集料骨架二值化图像进行形态学操作,提取骨架棱角点,并去除图像中的孤立点。
步骤七,填充集料骨架二值化图像中的空洞区域,并对每幅集料骨架图像棱角点数求和,本实施例中,对图2集料骨架棱角点提取图中的集料总棱角点求和为si=51。
步骤八,以平均棱角点表征单个集料的棱角性,用平均棱角点和面积权值共同量化整批集料的棱角性,将骨架棱角点法棱角性量化指标定义如下:
Figure BDA0001708198950000051
Figure BDA0001708198950000052
Figure BDA0001708198950000053
由于每个集料是在三种不同状态下得到的二维图形,式中,i代表单个集料的不同三种状态,si为集料在其中某一个状态下的棱角点,P为单个集料的平均棱角点数;ai为集料在其中某一个状态下的集料面积,A为单个集料的平均面积;Ia为整批集料棱角性量化指标,N为总集料个数。
本实施例中,以
Figure BDA0001708198950000061
平均棱角点表征单个集料的棱角性,求得相对应集料的面积
Figure BDA0001708198950000062
用平均棱角点P和面积A权值
Figure BDA0001708198950000063
共同量化整批集料的棱角性,例如用5颗集料代表整批集料的棱角性,棱角指标值
Figure BDA0001708198950000064
其中面积A用像素值表示。

Claims (1)

1.一种路用粗集料棱角性量化方法,其特征在于,按以下步骤实施:
步骤一,输入集料棱角性检测系统得到的集料图像;
步骤二,将集料图像灰度化、中值滤波去噪、高亮细节增强、图像熵分割、闭运算目标提取、坑洞填充一系列处理;
步骤三,读取步骤二中的集料图像,采用M结构元素进行形态学集料骨架提取,得到集料边缘棱角的全部信息;
步骤四,初始化各个变量,按照式
Figure FDA0002911089080000011
对目标物体进行形态学运算,判断当前Corr_TXi是否为空,若为空,则循环停止,并获取集料颗粒骨架图像;
步骤五,读取集料颗粒骨架图像,形态学闭运算进行集料骨架图像像素连接操作,并求骨架图像的大小x,y,即为行数和列数;从第二行逐个扫描集料骨架点的像素值大小,判断八邻域像素是否符合第n列的一个点,周围邻域中,若它的m+1行,n+1列和m-1行,n-1列所在点的像素值为1,或者它的m+1行,n-1列和m-1行,n+1列所在点的像素值为1,则不需要断点修补,同时去除闭运算操作后的一些干扰点;否则需要断裂点的修补,将像素点0替换成1进行骨架连通改善;扫描完集料图像像素的所有点,循环终止;
步骤六,读取子文件夹里面所有的修补完的集料骨架图像,利用形态学操作移除目标边界像素,但不允许目标分隔开,保留下来的像素组成图像的骨架;将集料骨架图像二值化,在四邻域提取边界点,对集料骨架二值化图像进行形态学操作,提取骨架棱角点,并去除图像中的孤立点;
步骤七,填充集料骨架二值化图像中的空洞区域,并对每幅集料骨架图像棱角点数求和;
步骤八,以平均棱角点表征单个集料的棱角性,用平均棱角点和面积权值共同量化整批集料的棱角性,将骨架棱角点法棱角性量化指标定义如下:
Figure FDA0002911089080000021
Figure FDA0002911089080000022
Figure FDA0002911089080000023
由于每个集料是在三种不同状态下得到的二维图形,式中,i代表单个集料的不同三种状态,si为集料在其中某一个状态下的棱角点,P为单个集料的平均棱角点数;ai为集料在其中某一个状态下的集料面积,A为单个集料的平均面积;Ia为整批集料棱角性量化指标,N为总集料个数。
CN201810667800.7A 2018-06-26 2018-06-26 一种路用粗集料棱角性量化方法 Active CN109116002B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810667800.7A CN109116002B (zh) 2018-06-26 2018-06-26 一种路用粗集料棱角性量化方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810667800.7A CN109116002B (zh) 2018-06-26 2018-06-26 一种路用粗集料棱角性量化方法

Publications (2)

Publication Number Publication Date
CN109116002A CN109116002A (zh) 2019-01-01
CN109116002B true CN109116002B (zh) 2021-04-20

Family

ID=64821939

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810667800.7A Active CN109116002B (zh) 2018-06-26 2018-06-26 一种路用粗集料棱角性量化方法

Country Status (1)

Country Link
CN (1) CN109116002B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110658117A (zh) * 2019-09-30 2020-01-07 中冶南方工程技术有限公司 一种粗集料棱角性评价装置及方法
CN111307070B (zh) * 2019-11-05 2021-06-15 长安大学 基于数字图像处理对混凝土粗集料的棱角度测量方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU614380A1 (ru) * 1977-03-01 1978-07-05 Коммунарский горно-металлургический институт Устройство дл определени угловатости зерен сыпучих материалов
CN101354241A (zh) * 2008-07-11 2009-01-28 长安大学 集料数字图像评价系统及评价方法
WO2015191779A1 (en) * 2014-06-10 2015-12-17 Advanced Concrete Technologies, Llc Methods for increasing aggregate hardness, hardened aggregate, and structures including the hardened aggregate
CN106442231A (zh) * 2016-08-31 2017-02-22 北京建筑大学 一种基于数字图像分析技术的粗集料棱角性评价方法
CN106918553A (zh) * 2017-03-26 2017-07-04 长安大学 基于aims集料成像系统的粗集料抗滑性能评价方法
CN206431062U (zh) * 2017-01-13 2017-08-22 长安大学 一种用于冷料供给系统粗集料棱角性实时监测装置

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU614380A1 (ru) * 1977-03-01 1978-07-05 Коммунарский горно-металлургический институт Устройство дл определени угловатости зерен сыпучих материалов
CN101354241A (zh) * 2008-07-11 2009-01-28 长安大学 集料数字图像评价系统及评价方法
WO2015191779A1 (en) * 2014-06-10 2015-12-17 Advanced Concrete Technologies, Llc Methods for increasing aggregate hardness, hardened aggregate, and structures including the hardened aggregate
CN106442231A (zh) * 2016-08-31 2017-02-22 北京建筑大学 一种基于数字图像分析技术的粗集料棱角性评价方法
CN206431062U (zh) * 2017-01-13 2017-08-22 长安大学 一种用于冷料供给系统粗集料棱角性实时监测装置
CN106918553A (zh) * 2017-03-26 2017-07-04 长安大学 基于aims集料成像系统的粗集料抗滑性能评价方法

Non-Patent Citations (9)

* Cited by examiner, † Cited by third party
Title
Evaluation of Fine Aggregate Morphology by Image Method and Its Effect on Skid-Resistance of Micro-Surfacing;Yue Xiao etal;《Materials》;20180529;第11卷;全文 *
Image-Based Coarse-Aggregate Angularity Analysis and Evaluation;Ju Huyan etal;《J.Mater.Giv.Eng》;20200630;第32卷(第6期);全文 *
Three-dimensional characterization of morphology and abrasion decay laws for coarse aggregates;Haitao Ge etal;《Construction and Building Materials》;20180823;第188卷;全文 *
基于数字图像处理的粗集料棱角性量化研究;李嘉等;《公路交通科技》;20080731;第25卷(第7期);全文 *
基于距离场细化的骨架提取算法;李川等;《微电子学与计算机》;20111031;第28卷(第10期);全文 *
沥青混凝土再生集料的棱角性分析;唐伯明等;《中国公路学报》;20150131(第1期);全文 *
粗集料棱角性的图像评价方法;汪海年等;《东南大学学报(自然科学版)》;20080731;第38卷(第4期);全文 *
粗集料棱角性试验方法研究;王火明等;《公路交通技术》;20121231(第6期);全文 *
细集料棱角性试验方法及标准研究;颜家学;《黑龙江交通科技》;20141231(第8期);全文 *

Also Published As

Publication number Publication date
CN109116002A (zh) 2019-01-01

Similar Documents

Publication Publication Date Title
CN107527067B (zh) 一种基于探地雷达的铁路路基病害智能识别方法
Zhang et al. Accurate detection and evaluation method for aggregate distribution uniformity of asphalt pavement
AU2015392660B2 (en) Stepwise-refinement pavement crack detection method
CN109584286B (zh) 一种基于广义回归神经网络的沥青路面构造深度计算方法
CN105445160B (zh) 一种沥青混合料的空隙特征及其提取方法
CN109087396A (zh) 一种基于混凝土ct图像像素特征的细观结构重构方法
CN101153850A (zh) 一种沥青混合料的检测方法及系统
Wang et al. Pavement crack segmentation algorithm based on local optimal threshold of cracks density distribution
CN112561903B (zh) 一种适用于寒冷地区沥青路面抗温缩裂缝的方法
CN112964712A (zh) 一种快速检测沥青路面状态的方法
CN109116002B (zh) 一种路用粗集料棱角性量化方法
Yu et al. Intelligent Crack Detection and Quantification in the Concrete Bridge: A Deep Learning‐Assisted Image Processing Approach
CN106529593A (zh) 路面病害检测方法和系统
Mubashshira et al. An unsupervised approach for road surface crack detection
CN113256805A (zh) 一种基于三维点云重构的路面线性裂缝信息快速计算方法
Shao et al. A kind of pavement crack detection method based on digital image processing
Safaei et al. A novel adaptive pixels segmentation algorithm for pavement crack detection
CN112085699A (zh) 一种基于二维图像的路面裂缝提取方法
Hu et al. A two-dimensional entropy-based method for detecting the degree of segregation in asphalt mixture
CN111307070B (zh) 基于数字图像处理对混凝土粗集料的棱角度测量方法
CN108986078A (zh) 一种基于凸包几何特征的粗集料棱角性量化方法
CN106568925A (zh) 水泥乳化沥青砂浆厚度分析方法
Safaei Pixel and region-based image processing algorithms for detection and classification of pavement cracks
CN110517220A (zh) 一种基于激光三维数据的集料表面数量检测方法
CN107524065B (zh) 一种多空路面再生时骨架恢复方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant