CN1766588A - 修复混凝土裂缝表面图像的智能化处理方法 - Google Patents
修复混凝土裂缝表面图像的智能化处理方法 Download PDFInfo
- Publication number
- CN1766588A CN1766588A CN 200510094637 CN200510094637A CN1766588A CN 1766588 A CN1766588 A CN 1766588A CN 200510094637 CN200510094637 CN 200510094637 CN 200510094637 A CN200510094637 A CN 200510094637A CN 1766588 A CN1766588 A CN 1766588A
- Authority
- CN
- China
- Prior art keywords
- image
- pixel
- gray
- scale value
- sigma
- 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.)
- Granted
Links
Images
Landscapes
- Investigating Materials By The Use Of Optical Means Adapted For Particular Applications (AREA)
Abstract
本发明涉及的是修复混凝土裂缝表面图像的智能化处理方法,它的操作步骤分:图像采集、图像分割、表面覆盖率的计算;优点:由于最大类间方差阈值分割法根据图像自动确定阈值,是一种面向均匀性的算法,依据均匀性度量得到最佳性能,所以图像分割精度较高,计算稳定可靠,自动化程度高、计算速度也较快,可有效的来测定表面覆盖率,评价电沉积效果,并有望建立表面覆盖率与混凝土耐久性提高之间的关系,以此来控制电沉积处理的时间。本发明还可以用于电化学沉积方法修复混凝土裂缝中裂缝处横断面覆盖率、纳米镀层中纳米颗粒含量的确定以及类似图像的处理。
Description
技术领域
本发明涉及的是修复混凝土裂缝表面图像的智能化处理方法,属于土木建筑工程与计算机科学与工程相结合的技术领域。
背景技术
对于电沉积效果的评价指标——表面覆盖率的测定,由于沉积物分布极不均匀,采用一般的方法将带来较大的误差,如利用Photoshop实现图像分割,其实现分割的操作繁复,自动化程度不高,“扩大选取”“选取相似”尺度不易把握以及色彩范围无法准确定位,带有很强的主观性,导致分割的精度很难准确控制,从而影响表面覆盖率计算的精确性;而对于自定义阈值的二值化图像分割法,如MATLAB软件中其二值化阈值参数是人为设定的0~1之间的数值,很难找到理想的阈值,也带有很强的主观性,影响图像分割的准确性,从而使表面覆盖率的计算精度很难准确控制。
发明内容
本发明的目的在于针对上述存在的缺陷,提出一种对修复混凝土裂缝表面图像的智能化处理方法。它具有图像分割精度高,计算稳定可靠,自动化程度高、计算速度快等特点。
本发明的目的是通过下列措施实现的:
修复混凝土裂缝表面图像的智能化处理方法,其特征是它的操作步骤分:
一、图像采集
拍摄经电化学沉积方法修复后带裂缝混凝土表面的数码照片,相机像数在300万以上,所述的修复混凝土裂缝的电化学沉积方法是采用ZnSO4或MgSO4、浓度0.05mol/L~0.25mol/L的溶液,浸泡带裂缝的混凝土,施加电流,混凝土表面的电流密度为0.25A/m2~1A/m2,对带裂缝的混凝土进行修复;
二、图像分割
调整图像的对比度,将对比度控制在20%~30%,拉大沉积物图像与混凝土表面图像的反差;采用最大类间方差阈值分割法,将直方图的某一阈值处理分割成两组,当被分成的两组的方差最大时,决定阈值,
设图像的灰度值在m1至m2(-1<m1<m2<256,m1,m2为整数)之间,则这幅图像具有(m2-m1+1)个灰级,图像象素总数为:N=图片长×图片宽(单位:象素),灰度值为i(m1-1<i<m2+1)的象素个数为H(i),
用k(m1-1<k<m2+1,k为整数)将图像象素分成两组:A={m1~k},表示灰度值在m1~k之间的象素群和B={k+1~m2},表示灰度值在k+1~m2之间的象素群,则各组中的象素个数与图像总象素个数比值分别为:
两组象素的数学期望值(均值)分别为:
两组间的方差值为:
σ2(k)=PA·PB(EA-EB)2
从m1至m2之间变换k值,分别计算以上方差值,取使得方差值最大时的k值作为图像二值化分割阈值T,
阈值确定以后,便可对图像进行二值化处理,判断图像上每个象素的灰度值,灰度值不小于阈值T的象素被赋予白色,灰度值小于阈值T的象素被赋为黑色,形成一幅二值化图像;
三、表面覆盖率的计算
设置背景色(灰度值i),统计分割后图像中灰度为背景色(即灰度值为i)的象素个数H(i),得出其分布概率为
(N为图像象素总数,N=图片长×图片宽(单位:象素),则沉积物表面覆盖率为1-F(i),
对于二值化分割后的黑白图像,黑色为背景色,灰度值为0,表面覆盖率按下式计算:
本发明的有益效果:
由于最大类间方差阈值分割法根据图像自动确定阈值,是一种面向均匀性的算法,依据均匀性度量得到最佳性能,所以图像分割精度较高,计算稳定可靠,可有效的来测定表面覆盖率,评价电沉积效果,并有望建立表面覆盖率与混凝土耐久性提高之间的关系,以此来控制电沉积处理的时间。本发明提供的表面覆盖率测定方法还可以用于电化学沉积方法修复混凝土裂缝中裂缝处横断面覆盖率、纳米镀层中纳米颗粒含量的确定以及类似图像的处理。
本发明提供的表面覆盖率测定方法图像分割精度较高,计算稳定可靠,自动化程度高、计算速度也较快。
附图说明
图1是经电化学沉积方法修复后带裂缝混凝土表面的图像
图2是图像分割过程中的对原始图像剪裁后的图像
图3是图像对比度调整后的图像
图4是经二值化后的图像
图5是本发明的流程图
具体实施方式
实施例
一、图像采集
电化学沉积方法采用MgSO4、浓度0.05mol/L的溶液,浸泡带裂缝的混凝土,施加电流,混凝土表面的电流密度为1A/m2,对混凝土裂缝进行修复,用像数在300万以上的数码相机对修复后带裂缝混凝土的表面拍摄,如图1所示;
二、图像分割
对原始图像进行剪裁,除去混凝土表面不用于计算表面覆盖率的部分,如图2所示,然后调整图像的对比度,将对比度控制在20%~30%,拉大沉积物图像与混凝土表面图像的反差,如图3所示,采用最大类间方差阈值分割法,将直方图的某一阈值处理分割成两组,当被分成的两组的方差最大时,决定阈值,
设图像的灰度值在m1至m2(-1<m1<m2<256,m1,m2为整数)之间,则这幅图像具有(m2-m1+1)个灰级,图像象素总数为:N=图片长×图片宽(单位:象素),灰度值为i(m1-1<i<m2+1)的象素个数为H(i),
用k(m1-1<k<m2+1,k为整数)将图像象素分成两组:A={m1~k},表示灰度值在m1~k之间的象素群和B={k+1~m2},表示灰度值在k+1~m2之间的象素群,则各组中的象素个数与图像总象素个数比值分别为:
两组象素的数学期望值(均值)分别为:
两组间的方差值为:
σ2(k)=PA·PB(EA-EB)2
从m1至m2之间变换k值,分别计算以上方差值,取方差值最大的k值168,作为图像二值化分割阈值T,
对图像进行二值化处理,判断图像上每个象素的灰度值,灰度值不小于阈值168的象素被赋予白色,灰度值小于阈值168的象素被赋为黑色,形成一幅二值化图像,如图4所示;
三、表面覆盖率的计算
设置背景色(灰度值i),统计分割后图像中灰度为背景色(即灰度值为i)的象素个数H(i),得出其分布概率为
(N为图像象素总数,N=图片长×图片宽(单位:象素),则沉积物表面覆盖率为1-F(i),
对于二值化分割后的黑白图像,黑色为背景色,灰度值为0,表面覆盖率按下式计算:
最后选中沉积物计算出带裂缝混凝土表面覆盖率为64.57%。
Claims (1)
1、修复混凝土裂缝表面图像的智能化处理方法,其特征是它的操作步骤分:
一、图像采集
拍摄经电化学沉积方法修复后带裂缝混凝土表面的数码照片,相机像数在300万以上,所述的修复混凝土裂缝的电化学沉积方法是采用ZnSO4或MgSO4、浓度0.05mol/L~0.25mol/L的溶液,浸泡带裂缝的混凝土,施加电流,混凝土表面的电流密度为0.25A/m2~1A/m2,对带裂缝的混凝土进行修复;
二、图像分割
调整图像的对比度,将对比度控制在20%~30%,拉大沉积物图像与混凝土表面图像的反差;采用最大类间方差阈值分割法,将直方图的某一阈值处理分割成两组,当被分成的两组的方差最大时,决定阈值,
设图像的灰度值在m1至m2(-1<m1<m2<256,m1,m2为整数)之间,则这幅图像具有(m2-m1+1)个灰级,图像象素总数为:N=图片长×图片宽(单位:象素),灰度值为i(m1-1<i<m2+1)的象素个数为H(i),
用k(m1-1<k<m2+1,k为整数)将图像象素分成两组:A={m1~k},表示灰度值在m1~k之间的象素群和B={k+1~m2},表示灰度值在k+1~m2之间的象素群,则各组中的象素个数与图像总象素个数比值分别为:
两组象素的数学期望值(均值)分别为:
两组间的方差值为:
σ2(k)=PA·PB(EA-EB)2
从m1至m2之间变换k值,分别计算以上方差值,取使得方差值最大时的k值作为图像二值化分割阈值T,
阈值确定以后,便可对图像进行二值化处理,判断图像上每个象素的灰度值,灰度值不小于阈值T的象素被赋予白色,灰度值小于阈值T的象素被赋为黑色,形成一幅二值化图像;
三、表面覆盖率的计算
设置背景色(灰度值i),统计分割后图像中灰度为背景色(即灰度值为i)的象素个数H(i),得出其分布概率为
(N为图像象素总数,N=图片长×图片宽(单位:象素),则沉积物表面覆盖率为1-F(i),
对于二值化分割后的黑白图像,黑色为背景色,灰度值为0,表面覆盖率按下式计算:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2005100946372A CN100478681C (zh) | 2005-09-30 | 2005-09-30 | 对混凝土裂缝修复后呈现的表面图像的处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2005100946372A CN100478681C (zh) | 2005-09-30 | 2005-09-30 | 对混凝土裂缝修复后呈现的表面图像的处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1766588A true CN1766588A (zh) | 2006-05-03 |
CN100478681C CN100478681C (zh) | 2009-04-15 |
Family
ID=36742598
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2005100946372A Expired - Fee Related CN100478681C (zh) | 2005-09-30 | 2005-09-30 | 对混凝土裂缝修复后呈现的表面图像的处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100478681C (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101726266B (zh) * | 2009-12-07 | 2011-08-31 | 东华大学 | 利用计算机图像处理测定热轧非织造布轧点面积比率的方法 |
CN102496006A (zh) * | 2011-12-05 | 2012-06-13 | 北方工业大学 | 一种空白纸张检测方法 |
CN104915652A (zh) * | 2015-06-09 | 2015-09-16 | 迪安杰科技无锡有限公司 | 指纹图像的修复方法和装置 |
-
2005
- 2005-09-30 CN CNB2005100946372A patent/CN100478681C/zh not_active Expired - Fee Related
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101726266B (zh) * | 2009-12-07 | 2011-08-31 | 东华大学 | 利用计算机图像处理测定热轧非织造布轧点面积比率的方法 |
CN102496006A (zh) * | 2011-12-05 | 2012-06-13 | 北方工业大学 | 一种空白纸张检测方法 |
CN102496006B (zh) * | 2011-12-05 | 2013-07-03 | 北方工业大学 | 一种空白纸张检测方法 |
CN104915652A (zh) * | 2015-06-09 | 2015-09-16 | 迪安杰科技无锡有限公司 | 指纹图像的修复方法和装置 |
Also Published As
Publication number | Publication date |
---|---|
CN100478681C (zh) | 2009-04-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106290388B (zh) | 一种绝缘子故障自动检测方法 | |
CN105510195B (zh) | 一种堆叠骨料的粒度粒形在线检测方法 | |
CN110458119A (zh) | 一种非接触式测量的混凝土骨料级配快速识别方法 | |
CN109035232A (zh) | 一种基于深度学习和机器视觉融合的夹片缺陷检测方法 | |
CN108346135B (zh) | 一种改进的Criminisi图像修复方法 | |
CN116630813B (zh) | 一种公路路面施工质量智能检测系统 | |
CN112964712A (zh) | 一种快速检测沥青路面状态的方法 | |
CN100351853C (zh) | 高噪声图像特征点自动提取方法 | |
CN115496692B (zh) | 一种润滑油磨粒图像增强方法 | |
CN1766588A (zh) | 修复混凝土裂缝表面图像的智能化处理方法 | |
CN112258529A (zh) | 一种基于实例分割算法的路面裂缝像素级别检测方法 | |
CN110956183B (zh) | 一种沥青路面裂缝形态提取方法 | |
CN111291814A (zh) | 基于卷积神经网络与信息熵数据融合策略的裂纹识别算法 | |
CN109120230B (zh) | 一种太阳能电池片el图像检测和缺陷识别方法 | |
CN111968079B (zh) | 基于断面局部极值及分段稀疏的三维路面裂缝提取方法 | |
CN113916894A (zh) | 基于量子鸽群峰谷区间分割的铜冲压件表面裂缝检测方法 | |
CN1821755A (zh) | 水中絮体分形脉动影像实时检测方法 | |
CN113570651B (zh) | 基于sem图像的碳酸盐岩储层孔隙半径分布定量方法 | |
CN111738131A (zh) | 一种合金两相微结构参数特征提取方法 | |
CN110987993A (zh) | 一种原位有机质定量统计分析方法 | |
CN111724331A (zh) | 一种基于生成网络的多孔介质图像重建方法 | |
CN108629765B (zh) | 基于序列阈值差的精子显微视频序列滤波质量客观评价方法 | |
CN116862871A (zh) | 一种基于混合特征的木材计数方法 | |
CN103761525A (zh) | 基于显微ct数字图像处理铝电解用阴极炭块组分检测方法 | |
CN1845175A (zh) | 基于小波和共生矩阵的纹理表面缺损检测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20090415 Termination date: 20091030 |