CN109060813A - 基于图像优化分类的土石料级配自动检测方法 - Google Patents

基于图像优化分类的土石料级配自动检测方法 Download PDF

Info

Publication number
CN109060813A
CN109060813A CN201810570185.8A CN201810570185A CN109060813A CN 109060813 A CN109060813 A CN 109060813A CN 201810570185 A CN201810570185 A CN 201810570185A CN 109060813 A CN109060813 A CN 109060813A
Authority
CN
China
Prior art keywords
image
earth
stone material
stone
building stones
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.)
Pending
Application number
CN201810570185.8A
Other languages
English (en)
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.)
China University of Mining and Technology Beijing CUMTB
Original Assignee
China University of Mining and Technology Beijing CUMTB
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 China University of Mining and Technology Beijing CUMTB filed Critical China University of Mining and Technology Beijing CUMTB
Priority to CN201810570185.8A priority Critical patent/CN109060813A/zh
Publication of CN109060813A publication Critical patent/CN109060813A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/84Systems specially adapted for particular applications
    • G01N21/85Investigating moving fluids or granular solids
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0004Industrial image inspection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/11Region-based segmentation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/136Segmentation; Edge detection involving thresholding
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/155Segmentation; Edge detection involving morphological operators
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/84Systems specially adapted for particular applications
    • G01N21/85Investigating moving fluids or granular solids
    • G01N2021/8592Grain or other flowing solid samples
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10004Still image; Photographic image
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30108Industrial image inspection
    • G06T2207/30132Masonry; Concrete

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Health & Medical Sciences (AREA)
  • Quality & Reliability (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Image Analysis (AREA)

Abstract

本发明公开了一种基于图像优化分类的土石料级配自动检测方法,首先采用γ变换处理不同光照条件下采集的粒级混配的土石料图像,自适应处理各种明暗光照条件下的图像;采用动态otsu优化技术对土石料图像进行自适应阈值法处理,区分出土石两类目标;然后对二值化图像进行密度变换,保留完整的石料区域,弱化土石部分的噪声影响;将二值图像进行形态学变换,增强密度变换过程削弱的矿石边缘信息,提取图像中的矿石目标边界并进行几何形态统计;对图像中目标进行伪目标去除,保留的矿石目标进行粒级、体积、质量计算,获取土石料的累级级配曲线分布。适用于不同粒级混配的土石料自动检测,也适用于各类石料粒级的自动检测。

Description

基于图像优化分类的土石料级配自动检测方法
技术领域
本发明涉及一种不同粒级石料、不同材质石料与沙土混配而成的土石料级配的检测技术,尤其涉及一种基于图像优化分类的土石料级配自动检测方法。
背景技术
土石料由块度各异的石料与沙土混配而成,用于混凝土建筑施工、土石路基施工、桥梁及水利坝体。土石料的粒级分布检测,可提高土石料的材料使用效率,加快上述施工过程的建设速度及增强施工目标的坚固程度,对于上述工程建设具有重要的意义。施工现场的土石料中,土石材料各异,粒级及混配随机出现,目前没有自动化精确检测技术。
用于施工过程中的土石料由于光照的不均匀使得采集图像亮度相差很大,其中存在若干阴影。土石材料呈堆积状,石料内部颜色复杂多样,图像中目标众多、大小差异大,大石块覆盖些许沙土而使得内部噪声严重,小石块与沙土混叠。混合土石料图像中分割石料目标及进行粒级分析的难度很大。
目前,土石料粒级主要是通过人工筛分的方式获得,这种方法费时费力,效率很低。现有技术中可以对土石料粒级自动分析的方法少之又少。
发明内容
本发明的目的是提供一种准确、有效的基于图像优化分类的土石料级配自动检测方法。
本发明的目的是通过以下技术方案实现的:
本发明的基于图像优化分类的土石料级配自动检测方法,包括步骤:
图像中的石料及沙土目标区域分离:采用高分辨率工业相机采集土石图像,并增强所述土石料图像的对比度,对增强后的图像采用动态otsu优化二值分割;
对石料区域精确提取:采用密度变换技术去弱化沙土区域的噪声信息,去除噪声目标,增强矿石目标的边缘信息;
土石料粒级自动统计:进行矿石目标的几何特征、粒级、体积及级配自动统计。
由上述本发明提供的技术方案可以看出,本发明实施例提供的基于图像优化分类的土石料级配自动检测方法,能准确、有效的进行土石料图像分割,并实现混配土石料的粒级自动检测和自动提取,适用于不同粒级混配的土石料自动检测,也适用于各类石料及粒级的自动检测。
附图说明
图1为本发明的具体实施例中混合粒级土石料的灰度图;
图2为本发明的具体实施例中的伽马变换效果图;
图3为本发明的具体实施例中混合粒级土石料的原图;
图4为本发明的具体实施例中的动态自适应Otsu优化阈值图;
图5为本发明的具体实施例中的二值化结果图;
图6为本发明的具体实施例中的密度变换效果图;
图7为本发明的具体实施例中的膨胀示意图;
图8为本发明的具体实施例中的膨胀效果图;
图9为本发明实施例提供的基于图像优化分类的土石料级配自动检测方法的流程示意图。
具体实施方式
下面将对本发明实施例作进一步地详细描述。本发明实施例中未作详细描述的内容属于本领域专业技术人员公知的现有技术。
本发明的基于图像优化分类的土石料级配自动检测方法,其较佳的具体实施方式是:
包括步骤:
图像中的石料及沙土目标区域分离:采用高分辨率工业相机采集土石图像,并增强所述土石料图像的对比度,对增强后的图像采用动态otsu优化二值分割;
对石料区域精确提取:采用密度变换技术去弱化沙土区域的噪声信息,去除噪声目标,增强矿石目标的边缘信息;
土石料粒级自动统计:进行矿石目标的几何特征、粒级、体积及级配自动统计。
所述图像中的石料及沙土目标区域分离具体包括:
A、增强所述土石料图像的对比度,将所采集的图像进行γ变换处理,将图像中土石、不同石料像素的灰度值进行增强,扩大土石像素间、石料内部及边缘像素间的灰度差异性,减少石料内部像素间的差异;
B、将步骤A处理后的土石图像采用动态otsu优化方法进行二值化,对于每幅图片的每个像素都通过最小阈值进行绑定,将它的两个窗口由图像和对象大小进行自适应。
所述对石料区域进行精确提取具体包括:
C、二值化图像中,沙土部分噪声严重,采用像素邻域密度变换方法,保留完整的石料区域,弱化土石部分的噪声像素;
D、将图像进行形态学膨胀变换处理,增强密度变换过程削弱的矿石边缘信息,提取图像中的矿石目标边界并进行几何形态统计;
所述进行土石料粒级自动统计包括:
E、对图像中目标进行边界提取,计算每个目标的体积、质量与粒级,根据土石料每个粒级的质量百分比,获取土石料的累级级配曲线分布。
本发明的基于图像优化分类的土石料级配自动检测方法,首先采用γ变换处理不同光照条件下采集的粒级混配的土石料图像,自适应处理各种明暗光照条件下的图像;之后采用动态otsu优化技术对土石料图像进行自适应阈值法处理,区分出土石两类目标;然后对二值化图像进行密度变换,保留完整的石料信息,弱化土石部分的噪声影响;将图像进行形态学变换,增强密度变换过程削弱的矿石边缘信息,提取图像中的矿石目标边界并进行几何形态统计;对图像中目标进行面积过滤,保留的矿石目标进行粒径计算,获取土石料的粒级分布。能准确、有效的进行土石料图像分割,并实现混配土石料的粒级自动检测和自动提取,适用于不同粒级混配的土石料自动检测,也适用于各类石料及粒级的自动检测。
具体实施例,如图1至图9所示:
首先对土石料通过工业摄像头采集图像,其主要包括:
1、伽马变换增强土石料图像
由于光照或者灯光变化不均匀的原因会造成图像上的阴影,而γ曲线矫正就是通过一定的方法来矫正图像的这种偏差的方法。一般情况下,当用于γ矫正的值大于1时,图像的高光部分被压缩而暗调部分被扩展,当γ矫正的值小于1时,图像的高光部分被扩展而暗调部分被压缩,γ矫正一般用于平滑的扩展暗调的细节。公式如下:
f(x,y)′=f(x,y)γ (1)
原图像表示为f(x,y),指数值表示为γ。
2、动态自适应Otsu优化分离土、石目标
采用动态自适应Otsu优化,将Otsu的基于最大类间方差的最优阈值与Simphiwe提出的双窗自适应阈值法相结合,根据图像中矿石图像及目标的相对大小,针对同一中心像素确定双窗尺寸,再对双窗应用Otsu法分别求得各自局部最优阈值,取两者中的较小者为最佳阈值,双边滤波图中灰度值大于最佳阈值的为矿石,小于最佳阈值的为背景,从而生成二值图。复杂矿石图像中目标大小不一,光照不均,该方法通过小尺度的窗检测小粒度的沙土目标,较大尺度的窗检测大粒度的矿石目标,具有较强的适应性,二值化效果优于传统方法。
γ矫正后图像表达为f(x,y),首先分别求得每个像素p的两个邻域窗口,尺寸为n×n和m×m(即Npn和Npm),窗口的最优阈值分别为Totsu(Npn,t)和Totsu(Npm,t)。
式中,t为灰度级范围取值[0,1];l为最大灰度级(255);σB2t为最大类间方差,方差越大代表着错分率越小。
其次,根据最优阈值将图像f(x,y)二值化,fb(x,y)为二值化结果。
土石料图像(图3)二值化是后续种子区域提取的基础,二值化结果见图4。
3、密度变换去除背景噪声
前述图像背景存在大量的颗粒性噪声,采用密度过滤器去除噪声,前述二值化图像表述为f(x,y),经过密度过滤器去除噪声后的图片像素值用fb(x,y)表示,其中(x,y)为像素点坐标。过滤器的尺寸选取为a×a。输入图像中以每个白色像素点为中心,取一个a×a邻域,统计邻域范围中存在的白色像素点数目,达到下面公式中的阈值要求,是图像中此像素值是作为白色目标的条件,达不到阈值此白色像素点设置为黑色背景。
4、面积统计去除噪声目标
上述密度变换后的图像噪声目标信息被弱化,提取土石料图像中石块目标区域的面积、周长、圆形度等特征,去除噪声性伪目标。
(1)土石方图像中石块边界用链码记为a1a2…an,Δli为每码段ai的边界线段长度,周长P为:
其中,ne为石块边界链码序列中偶数码个数;n为边界链码中码的总个数。
(2)石块区域面积S
其中Δy=yi-yi-1,Δx=xi-xi-1,xi和yi分别为土石方图像中石块边界横纵坐标,x0和y0分别为石块边界起始点横纵坐标。
(3)石块区域圆形度R0
R0=4πS/P2 0<R0≤1 (7)
其中,S和P分别为土石方图像中石块区域的面积和周长。
面积非常小的区域作为噪声,直接去除。面积小于阈值且圆形度严重偏离1的目标为伪目标,直接去除,剩余部分为石料目标。
5、形态学膨胀变换恢复石料边界
密度变换损失掉石料目标的边缘信息,采用形态学膨胀变换恢复图像中的石料目标区域,并可以填充石料内部的噪声。每个白色像素点为中心,在图像中取一个结构元素B,结构B范围中的像素点,不论原像素值是白色还是黑色,均设置像素值作为白色。
膨胀(dilation)运算的定义是,把结构元素B平移目标a后得到Ba,若Ba击中X,则记下a点。所有满足上述条件的a点组成的集合称做X被B膨胀的结果,如图3所示,用公式表示为:
图7中X是被处理的对象,B是结构元素,不难知道,对于任意一个在阴影部分的点a,Ba击中X,X被B膨胀的结果就是那个阴影部分。
保留下的石料目标进行粒级计算,并生成粒级分布统计结果。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明披露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求书的保护范围为准。

Claims (4)

1.一种基于图像优化分类的土石料级配自动检测方法,其特征在于,包括步骤:
图像中的石料及沙土目标区域分离:采用高分辨率工业相机采集土石图像,并增强所述土石料图像的对比度,对增强后的图像采用动态otsu优化二值分割;
对石料区域精确提取:采用密度变换技术去弱化沙土区域的噪声信息,去除噪声目标,增强矿石目标的边缘信息;
土石料粒级自动统计:进行矿石目标的几何特征、粒级、体积及级配自动统计。
2.根据权利要求1所述的基于图像优化分类的土石料级配自动检测方法,其特征在于,所述图像中的石料及沙土目标区域分离具体包括:
A、增强所述土石料图像的对比度,将所采集的图像进行γ变换处理,将图像中土石、不同石料像素的灰度值进行增强,扩大土石像素间、石料内部及边缘像素间的灰度差异性,减少石料内部像素间的差异;
B、将步骤A处理后的土石图像采用动态otsu优化方法进行二值化,对于每幅图片的每个像素都通过最小阈值进行绑定,将它的两个窗口由图像和对象大小进行自适应。
3.根据权利要求2所述的基于图像优化分类的土石料级配自动检测方法,其特征在于,所述对石料区域进行精确提取具体包括:
C、二值化图像中,沙土部分噪声严重,采用像素邻域密度变换方法,保留完整的石料区域,弱化土石部分的噪声像素;
D、将图像进行形态学膨胀变换处理,增强密度变换过程削弱的矿石边缘信息,提取图像中的矿石目标边界并进行几何形态统计。
4.根据权利要求1、2或3所述的基于图像优化分类的土石料级配自动检测方法,其特征在于,所述进行土石料粒级自动统计包括:
E、对图像中目标进行边界提取,计算每个目标的体积、质量与粒级,根据土石料每个粒级的质量百分比,获取土石料的累级级配曲线分布。
CN201810570185.8A 2018-06-05 2018-06-05 基于图像优化分类的土石料级配自动检测方法 Pending CN109060813A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810570185.8A CN109060813A (zh) 2018-06-05 2018-06-05 基于图像优化分类的土石料级配自动检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810570185.8A CN109060813A (zh) 2018-06-05 2018-06-05 基于图像优化分类的土石料级配自动检测方法

Publications (1)

Publication Number Publication Date
CN109060813A true CN109060813A (zh) 2018-12-21

Family

ID=64820406

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810570185.8A Pending CN109060813A (zh) 2018-06-05 2018-06-05 基于图像优化分类的土石料级配自动检测方法

Country Status (1)

Country Link
CN (1) CN109060813A (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108375529A (zh) * 2018-02-12 2018-08-07 中国水利水电科学研究院 土石料级配自动检测方法及检测系统
CN112017164A (zh) * 2020-08-18 2020-12-01 中国水利水电科学研究院 基于深度阈值卷积模型的土石料级配检测方法
CN113610771A (zh) * 2021-07-15 2021-11-05 三峡大学 基于b样条曲线和粒子群算法的堆石料图像检测级配曲线矫正方法
CN114356249A (zh) * 2022-01-10 2022-04-15 苏州久道信息科技有限公司 一种自动检测和优化图像质量的处理系统及其处理方法
WO2023279556A1 (zh) * 2021-07-09 2023-01-12 福建南方路面机械股份有限公司 一种骨料粒径抽样监测方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104063866A (zh) * 2014-06-26 2014-09-24 中国矿业大学(北京) 一种矿石传送过程中的粒度检测方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104063866A (zh) * 2014-06-26 2014-09-24 中国矿业大学(北京) 一种矿石传送过程中的粒度检测方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
YSC6688: "《https://blog.csdn.net/ysc6688/article/details/50772382/》", 18 July 2017, CSDN博客 *
宓超等: "《装卸机器视觉及其应用》", 31 January 2016, 上海科学技术出版社 *
贺添等: "《ERDAS遥感影像处理基础实验教程》", 31 March 2014, 黄河水利出版社 *
靳晓颖等: "基于标记分水岭分割算法在矿石异物检测应用研究", 《计算机科学与应用》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108375529A (zh) * 2018-02-12 2018-08-07 中国水利水电科学研究院 土石料级配自动检测方法及检测系统
CN112017164A (zh) * 2020-08-18 2020-12-01 中国水利水电科学研究院 基于深度阈值卷积模型的土石料级配检测方法
WO2023279556A1 (zh) * 2021-07-09 2023-01-12 福建南方路面机械股份有限公司 一种骨料粒径抽样监测方法
CN113610771A (zh) * 2021-07-15 2021-11-05 三峡大学 基于b样条曲线和粒子群算法的堆石料图像检测级配曲线矫正方法
CN113610771B (zh) * 2021-07-15 2024-05-28 三峡大学 基于b样条曲线和粒子群算法的堆石料图像检测级配曲线矫正方法
CN114356249A (zh) * 2022-01-10 2022-04-15 苏州久道信息科技有限公司 一种自动检测和优化图像质量的处理系统及其处理方法

Similar Documents

Publication Publication Date Title
CN109060813A (zh) 基于图像优化分类的土石料级配自动检测方法
CN112419250B (zh) 路面裂缝数字图像提取、裂缝修补及裂缝参数计算方法
Singh et al. Shadow detection and removal from remote sensing images using NDI and morphological operators
CN104792792B (zh) 一种逐步求精的路面裂缝检测方法
CN109242791B (zh) 一种针对破损植物叶片的批量修复方法
CN102663726B (zh) 对物料堆积图像进行物料分割的方法和装置
CN114677601B (zh) 一种基于无人机巡检的结合深度学习的大坝裂缝检测方法
CN110956183A (zh) 一种沥青路面裂缝形态提取方法
CN108038883A (zh) 一种应用于公路路面视频图像的裂缝检测及识别方法
CN103942803A (zh) 基于sar图像的水域自动检测方法
CN112861654A (zh) 一种基于机器视觉的名优茶采摘点位置信息获取方法
CN104008528A (zh) 基于阈值分割的非均匀光场水下目标探测图像增强方法
CN112288726B (zh) 一种井下带式输送机带面异物检测方法
CN116342635B (zh) 一种地质测绘中裂缝轮廓提取方法
CN106023134A (zh) 一种钢材晶粒的自动晶界提取方法
CN108764132B (zh) 一种湖泊湿地遥感图像误差检测方法
CN102262778A (zh) 基于改进的分数阶微分掩模的图像增强方法
CN104167005A (zh) 一种基于相似度函数的自适应加权椒盐噪声滤波方法
CN105092597A (zh) 一种硬塑材料表面的裂纹检测方法
CN101655977B (zh) 基于差分图像检测、多窗口滤波去除图像脉冲噪声的方法
CN103455985B (zh) 基于Hessian结构分析的道路裂缝增强方法
CN107545557A (zh) 粪便图像中的虫卵检测方法和装置
JP2013105245A (ja) 画像処理方法、装置及びプログラム
Han et al. An enhanced image binarization method incorporating with Monte-Carlo simulation
CN116630321A (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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20181221