CN103593831A - 一种自动去除水泥浆体背散射电子图像制样缺陷的方法 - Google Patents
一种自动去除水泥浆体背散射电子图像制样缺陷的方法 Download PDFInfo
- Publication number
- CN103593831A CN103593831A CN201310513980.0A CN201310513980A CN103593831A CN 103593831 A CN103593831 A CN 103593831A CN 201310513980 A CN201310513980 A CN 201310513980A CN 103593831 A CN103593831 A CN 103593831A
- Authority
- CN
- China
- Prior art keywords
- image
- sample preparation
- backscattered electron
- images
- cement
- 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
Links
Images
Landscapes
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
本发明涉及一种自动去除水泥浆体背散射电子图像制样缺陷的方法,该方法包括以下步骤:(1)图像读取;(2)灰度图像二值化处理;(3)用滤波方法自动去除制样缺陷;(4)计算未水化水泥熟料的面积。当图像数量较多时,采用批量处理的方法。首先,将全部图像文件进行规范化命名,然后,按照步骤(1)和(2)将原始灰度图像逐一转化为二值图像,最后,循环执行步骤(3)和(4),自动去除全部图像的制样缺陷,同时输出每张图中未水化熟料颗粒面积的计算结果。本发明所述的图像处理方法可以为背散射电子图像的后期处理节约时间,避免人工操作标准不一致所产生的误差。
Description
技术领域
本发明涉及材料分析领域,尤其涉及一种自动去除水泥浆体背散射电子图像制样缺陷的方法。
背景技术
在水泥浆体微观结构研究中,背散射电子图像是一种非常有效的微观测试手段。大量研究表明,用它进行反应进程、界面特征、水化产物的物相组成及其形态与分布、孔结构特征等多项指标的观测都是卓有成效的。目前,这项技术最主要的应用在于测定水泥浆体中熟料的水化程度。
然而,背散射电子图像以样品剖面为观测对象,因此,所用的样品需要经过环氧树脂镶嵌、不同粒级的砂纸和抛光液的研磨、抛光等多道工序加以处理,制样过程十分繁琐。由于水泥浆体是多相、非均质材料,内部还含有为数不少的孔隙,加之水胶比、龄期、辅助胶凝材料的种类及其掺量等多种因素都会影响浆体的硬度以及微观结构的致密程度,因此,其制样工艺需要根据样品特点进行相应调整,不能一成不变。但是,磨抛机的压力、转速、研磨时间、研磨介质的用法、用量等诸多因素都会对磨抛效果产生影响,因此,操作人员往往在样品制备环节耗费了大量时间之后,得到的图像仍然是带有一些瑕疵的。其中,比较常见的是在熟料颗粒上存在一些细小的凹坑或者划痕。
如果用这种带有制样缺陷的图像计算水化程度,会造成计算结果偏高。解决这一问题的传统做法是,将背散射电子图像二值化之后,由操作人员借助专业的图像处理软件,根据材料学的经验,人工去除熟料颗粒上残留的瑕疵。但是,这种做法不仅处理速度慢,而且还会因为操作人员的判断标准不一致而带来误差。
发明内容
本发明的目的是提供一种自动去除水泥浆体背散射电子图像制样缺陷的方法,它能够对水泥浆体背散射电子图像的二值化结果进行图像分析,自动识别并去除其中存在的划痕、点蚀、凹坑等制样缺陷。采用图像法计算水泥水化程度时,将该方法纳入其中,可以为背散射电子图像的后期处理节约时间,避免人工操作标准不一致所产生的误差,而且还可以摆脱图像修复阶段对专业图像处理软件的依赖。
本发明提出的自动去除水泥浆体背散射电子图像制样缺陷的方法,具体步骤如下:
(1) 图像读取:采用Photoshop、Matlab或Visual C++等软件,读取水泥浆体原始背散射电子图像的有效观测区域;
(2) 灰度图像二值化处理:计算出步骤(1)所得有效观测区域的灰度直方图,
根据其灰度分布特点,结合原始背散射电子图像特征,人工设定灰度阈值,对有效观测区域进行二值化处理,形成黑白图像,其中白色区域为未水化的水泥熟料颗粒,黑色区域为环氧树脂、水化产物、孔隙和制样缺陷等;
(3) 用滤波方法自动去除制样缺陷:根据黑白图像的放大倍数,选择滤波模板的像素区域,在待处理的黑白图像中逐点地移动滤波模板,将每一点(x,y)处的取值转变为滤波模板覆盖像素区域内所占比例最高的物相,由此去除大面积白色区域中残留的黑色孤立像素,即熟料中残留的制样缺陷,以及大片黑色区域中残留的白色孤立像素,即通常是灰度图像二值化时由于灰度阈值设定偏低而产生的熟料噪声信号;
(4) 计算未水化水泥熟料的面积:对步骤(3)所得的经过滤波后的图像进行统计分析,计算其中白色区域的像素总数,即为未水化水泥熟料颗粒的总面积。
本发明中,步骤(3)中所述滤波模板的像素区域通常是3×3、5×5或7×7等尺寸的像素区域。
在本发明中,当需要处理的原始背散射电子图像数量较多时,采用批量处理的方法对多张图像进行统一处理,进一步提高图像处理效率。在这种情况下,本发明还可以通过如下技术措施来实现:首先,将全部图像文件进行规范化命名,然后,按照步骤(1)和(2)将原始灰度图像逐一转化为二值化的黑白图像,最后,循环调用步骤(3)和(4),自动去除全部图像的制样缺陷,同时输出每张图像中未水化水泥熟料颗粒面积的计算结果。
与传统的背散射电子图像缺陷去除方法相比,本发明具有以下特点:
(1) 去除图像缺陷的途径:通过滤波程序,对二值图像中熟料颗粒上残存的制样缺陷进行自动滤波,这可以提高图像处理速度,并且避免人工操作中因缺陷判断标准不一致而产生的误差,此外还可以摆脱对专业图像处理软件的依赖;
(2) 批量处理多张图像:当需要处理的原始背散射电子图像数量较多时,采用批量处理的方法对多张二值化的黑白图像进行统一处理,并将未水化熟料颗粒残核的面积汇总、输出,这可以显著提高图像处理效率。
附图说明
图1是实施例1中硅酸盐水泥水化1d的背散射电子图像(w/c=0.3,图像放大倍数为500倍)。
图2 是实施例1中图1所示背散射电子图像的灰度直方图。
图3 是实施例1中图1所示背散射电子图像的二值化处理结果。
图4 是实施例1中图3所示二值化图像的滤波结果。
图5是实施例2中硅酸盐水泥水化90d的背散射电子图像(w/c=0.3,图像放大倍数为4000倍)。
图6 是实施例2中图5所示背散射电子图像的灰度直方图。
图7 是实施例2中图5所示背散射电子图像的二值化处理结果。
图8 是实施例2中图7所示二值化图像的滤波结果。
图9是实施例3中粉煤灰水泥水化1d的背散射电子图像(w/c=0.3,图像放大倍数为500倍)。
图10 是实施例3中图9所示背散射电子图像的灰度直方图。
图11 是实施例3中图9所示背散射电子图像的二值化处理结果。
图12 是实施例3中图11所示二值化图像的滤波结果。
图13是滤波方法的流程图。
其中:f(x,y):黑白图像中每个像素的原始灰度值,取值为0或1;
fmax:模板覆盖区域内出现几率最高的物相,取值为0或1;
g(x,y):滤波后图像中每一像素(x,y)的数值,取值为0或1。
在上述图像中,0代表黑色,可能是环氧树脂、水化产物、孔隙、制样缺陷等;1代表白色,是未水化的水泥熟料颗粒。
具体实施方式
下面通过实施例进一步说明本发明。
实施例1
用背散射电子图像研究一种硅酸盐水泥在0.3水灰比下水化1d的水化程度,图像放大倍数为500倍。在个别图像中,熟料颗粒上存在一些微小的凹坑以及褶皱状的缺陷。对于带有此类缺陷的图像,可以按照以下步骤进行处理:
(1) 图像读取:用Matlab软件读取水泥浆体原始背散射电子图像的有效观测区域(1024×884像素);
(2) 灰度图像二值化处理:计算出步骤(1)所得有效观测区域的灰度直方图,根据其灰度分布特点,结合原始背散射电子图像特征,将灰度阈值设为130,然后用这一阈值对有效观测区域进行二值化处理,形成黑白图像,其中白色区域为未水化的水泥熟料颗粒;
(3) 用滤波方法自动去除制样缺陷:读取步骤(2)所得的黑白图像,以5×5的像素区域为滤波模板,在待处理的黑白图像中逐点地移动模板,将每一点(x,y)处的取值转变为模板覆盖的25个像素中所占比例最高的物相,由此去除熟料颗粒上残留的制样缺陷;
(4) 计算未水化水泥熟料的面积:对步骤(3)所得的图像进行统计分析,计算其中白色区域的像素总数,即为未水化水泥熟料颗粒的总面积。
实施例2
用背散射电子图像研究一种硅酸盐水泥在0.3水灰比下水化90d后浆体的水化程度以及孔结构特征,图像放大倍数为4000倍。拍摄过程中发现,在少数图像中熟料颗粒上存在明显的划痕,因此,按照以下步骤处理这些缺陷:
(1) 图像读取:用Matlab软件读取水泥浆体原始背散射电子图像的有效观测区域(1024×884像素);
(2) 灰度图像二值化处理:计算出步骤(1)所得有效观测区域的灰度直方图,根据其灰度分布特点,结合原始背散射电子图像特征,将灰度阈值设为115,然后用这一阈值对有效观测区域进行二值化处理,形成黑白图像,其中白色区域为未水化的水泥熟料颗粒;
(3) 用滤波方法自动去除制样缺陷:读取步骤(2)所得的黑白图像,以7×7的像素区域为滤波模板,在待处理的黑白图像中逐点地移动模板,将每一点(x,y)处的取值转变为模板覆盖的49个像素中所占比例最高的物相,由此去除熟料颗粒上残留的制样缺陷;
(4) 计算未水化水泥熟料的面积:对步骤(3)所得的图像进行统计分析,计算其中白色区域的像素总数,即为未水化水泥熟料颗粒的总面积。
实施例3
用背散射电子图像研究一种粉煤灰水泥在0.3水胶比下水化1d后熟料的水化程度,图像放大倍数为500倍。拍摄过程中发现,部分图像中熟料颗粒上存在微小凹坑以及褶皱状的缺陷,对于这些图像,可以按照以下步骤进行处理:
(1) 图像读取:用Matlab软件读取水泥浆体原始背散射电子图像的有效观测区域(1024×884像素);
(2) 灰度图像二值化处理:计算出步骤(1)所得有效观测区域的灰度直方图,根据其灰度分布特点,结合原始背散射电子图像特征,将灰度阈值设为136,然后用这一阈值对有效观测区域进行二值化处理,形成黑白图像,其中白色区域为未水化的水泥熟料颗粒;
(3) 用滤波方法自动去除制样缺陷:读取步骤(2)所得的黑白图像,以5×5的像素区域为滤波模板,在待处理的黑白图像中逐点地移动模板,将每一点(x,y)处的取值转变为模板覆盖的25个像素中所占比例最高的物相,由此去除熟料颗粒上残留的制样缺陷;
(4) 计算未水化水泥熟料的面积:对步骤(3)所得的图像进行统计分析,计算其中白色区域的像素总数,即为未水化水泥熟料颗粒的总面积。
图1和图9分别是硅酸盐水泥和粉煤灰水泥水化1d的背散射电子图像有效观测区域。这两张图像的问题主要是熟料颗粒上存在微小的凹坑。产生这种缺陷的原因是,水化早期硅酸盐水泥和粉煤灰水泥的浆体结构都比较疏松,在制备电镜样品时,多采用相对缓和的研磨抛光工艺,这样可以防止粘结不牢固的未水化熟料以及粉煤灰颗粒发生脱落,但是,往往也会因为磨抛不够彻底而在熟料颗粒表面留下一些微小的凹坑,或是高低不平、类似“褶皱”的缺陷。因此,这类试样在制备电镜样品时,按照相应的灰度直方图(图2和图10)对灰度图像(图1和图9)进行二值化处理之后,上述缺陷多以黑色小点的形式散布于白色的熟料颗粒中(图3和图11)。采用5×5像素的滤波模板,对其进行自动滤波。从图4和图12所示的滤波结果来看,主要缺陷均已去除,两张图像中未水化的熟料颗粒分别为193299像素和110125像素。
图5所示的水泥浆体水化龄期长,而且没有掺加辅助胶凝材料,因此浆体结构致密。在研磨抛光过程中,这类样品可以承受较高的转速和压力,因此,将其打磨平整相对容易。但是,正因为磨抛机的转速高、压力大,磨抛过程中未被及时冲走的磨料和碎屑也更容易划伤样品,在熟料颗粒上留下难以消除的划痕。通过分析灰度直方图(图6),设定阈值将灰度图像二值化处理之后,划痕就以黑线的形式存在于白色熟料颗粒上(图7)。采用7×7像素的滤波模板,对其进行自动滤波。从图8所示的滤波结果来看,主要缺陷均已去除,图中未水化的熟料颗粒面积为95198像素。
从以上三个实例可以看出,本发明所述的图像处理方法能够有效消除水泥浆体背散射电子图像二值化结果中的制样缺陷,提高图像处理效率,减少操作人员的主观因素对处理结果的干扰,而且,滤波过程无需借助专业的图像处理软件就能完成。此外,本发明所述的图像处理方法对于硅酸盐水泥和复合水泥,早龄期样品和长龄期样品均适用。
Claims (2)
1.一种自动去除水泥浆体背散射电子图像制样缺陷的方法,其特征在于具体步骤如下:
(1)图像读取:采用Photoshop、Matlab或Visual C++软件,读取水泥浆体原始背散射电子图像的有效观测区域;
(2)灰度图像二值化处理:计算出步骤(1)所得有效观测区域的灰度直方图,根据其灰度分布特点,结合原始背散射电子图像特征,人工设定灰度阈值,对有效观测区域进行二值化处理,形成黑白图像,其中白色区域为未水化的水泥熟料颗粒,黑色区域为环氧树脂、水化产物、孔隙和制样缺陷;
(3)用滤波方法自动去除制样缺陷:根据黑白图像的放大倍数,选择滤波模板的像素区域,在待处理的黑白图像中逐点地移动滤波模板,将每一点(x,y)处的取值转变为滤波模板覆盖像素区域内所占比例最高的物相,由此去除大面积白色区域中残留的黑色孤立像素,即熟料中残留的制样缺陷,以及大片黑色区域中残留的白色孤立像素,即通常是灰度图像二值化时由于灰度阈值设定偏低而产生的熟料噪声信号;
(4)计算未水化水泥熟料的面积:对步骤(3)所得的经过滤波后的图像进行统计分析,计算其中白色区域的像素总数,即为未水化水泥熟料颗粒的总面积。
2.根据权利要求1所述的自动去除水泥浆体背散射电子图像制样缺陷的方法,其特征在于当需要处理的原始背散射电子图像数量较多时,采用批量处理的方法对多张图像进行统一处理,进一步提高图像处理效率,具体步骤如下:首先,将全部待处理的原始背散射电子图像文件进行规范化命名,然后,按照步骤(1)和(2)将原始灰度图像逐一转化为二值图像,最后,循环执行步骤(3)和(4),自动去除全部图像的制样缺陷,同时输出每张图像中未水化熟料颗粒面积的计算结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310513980.0A CN103593831A (zh) | 2013-10-25 | 2013-10-25 | 一种自动去除水泥浆体背散射电子图像制样缺陷的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310513980.0A CN103593831A (zh) | 2013-10-25 | 2013-10-25 | 一种自动去除水泥浆体背散射电子图像制样缺陷的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN103593831A true CN103593831A (zh) | 2014-02-19 |
Family
ID=50083956
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310513980.0A Pending CN103593831A (zh) | 2013-10-25 | 2013-10-25 | 一种自动去除水泥浆体背散射电子图像制样缺陷的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103593831A (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107615042A (zh) * | 2015-06-26 | 2018-01-19 | 栗田工业株式会社 | 树脂分析方法以及树脂处理方法 |
CN108961252A (zh) * | 2018-07-27 | 2018-12-07 | Oppo(重庆)智能科技有限公司 | 检测logo粘贴不良的方法、电子装置及计算机可读存储介质 |
CN111060540A (zh) * | 2019-12-18 | 2020-04-24 | 汇鸿智能科技(辽宁)有限公司 | 一种夹杂物自动识别系统中夹杂物颗粒自动识别方法 |
CN111402206A (zh) * | 2020-02-27 | 2020-07-10 | 红云红河烟草(集团)有限责任公司 | 一种烟支刮破视觉检测方法及系统 |
CN113762266A (zh) * | 2021-09-01 | 2021-12-07 | 北京中星天视科技有限公司 | 目标检测方法、装置、电子设备和计算机可读介质 |
CN114295660A (zh) * | 2021-12-29 | 2022-04-08 | 中建西部建设湖南有限公司 | 一种水泥基材料初始堆积密度的计算方法 |
CN115035303A (zh) * | 2022-06-17 | 2022-09-09 | 郑州磨料磨具磨削研究所有限公司 | 一种电镀有色cBN砂轮的磨料浓度检测方法 |
CN116297581A (zh) * | 2023-03-06 | 2023-06-23 | 中国长江三峡集团有限公司 | 一种混凝土中方镁石反应程度的定量分析方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0500267A2 (en) * | 1991-02-22 | 1992-08-26 | AT&T Corp. | Image processing system |
-
2013
- 2013-10-25 CN CN201310513980.0A patent/CN103593831A/zh active Pending
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0500267A2 (en) * | 1991-02-22 | 1992-08-26 | AT&T Corp. | Image processing system |
Non-Patent Citations (3)
Title |
---|
吴浩,姚燕等: "数字图像处理技术在水泥混凝土研究中的应用", 《混凝土与水泥制品》, no. 4, 31 August 2007 (2007-08-31) * |
王培铭,丰曙霞,刘贤萍: "背散射电子图像分析在水泥基材料微观结构研究中的应用", 《硅酸盐学报》, no. 10, 31 October 2011 (2011-10-31) * |
雷海林: "水泥水化过程的计算机模拟算法分析", 《中国优秀硕士学位论文全文数据库(工程科技1辑)》, no. 4, 15 April 2006 (2006-04-15) * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107615042A (zh) * | 2015-06-26 | 2018-01-19 | 栗田工业株式会社 | 树脂分析方法以及树脂处理方法 |
CN108961252A (zh) * | 2018-07-27 | 2018-12-07 | Oppo(重庆)智能科技有限公司 | 检测logo粘贴不良的方法、电子装置及计算机可读存储介质 |
CN111060540A (zh) * | 2019-12-18 | 2020-04-24 | 汇鸿智能科技(辽宁)有限公司 | 一种夹杂物自动识别系统中夹杂物颗粒自动识别方法 |
CN111402206A (zh) * | 2020-02-27 | 2020-07-10 | 红云红河烟草(集团)有限责任公司 | 一种烟支刮破视觉检测方法及系统 |
CN113762266A (zh) * | 2021-09-01 | 2021-12-07 | 北京中星天视科技有限公司 | 目标检测方法、装置、电子设备和计算机可读介质 |
CN113762266B (zh) * | 2021-09-01 | 2024-04-26 | 北京中星天视科技有限公司 | 目标检测方法、装置、电子设备和计算机可读介质 |
CN114295660A (zh) * | 2021-12-29 | 2022-04-08 | 中建西部建设湖南有限公司 | 一种水泥基材料初始堆积密度的计算方法 |
CN115035303A (zh) * | 2022-06-17 | 2022-09-09 | 郑州磨料磨具磨削研究所有限公司 | 一种电镀有色cBN砂轮的磨料浓度检测方法 |
CN115035303B (zh) * | 2022-06-17 | 2024-04-26 | 郑州磨料磨具磨削研究所有限公司 | 一种电镀有色cBN砂轮的磨料浓度检测方法 |
CN116297581A (zh) * | 2023-03-06 | 2023-06-23 | 中国长江三峡集团有限公司 | 一种混凝土中方镁石反应程度的定量分析方法 |
CN116297581B (zh) * | 2023-03-06 | 2024-04-26 | 中国长江三峡集团有限公司 | 一种混凝土中方镁石反应程度的定量分析方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103593831A (zh) | 一种自动去除水泥浆体背散射电子图像制样缺陷的方法 | |
Tong et al. | A new image-based method for concrete bridge bottom crack detection | |
CN101648092B (zh) | 蓝藻水脱水处理方法 | |
CN110135275A (zh) | 一种集料颗粒识别方法 | |
CN105092593A (zh) | 一种基于高强度照明的钢板缺陷检测方法 | |
CN110261391A (zh) | 一种led芯片外观检测系统及方法 | |
CN110969620A (zh) | 一种磁瓦波纹类缺陷的检测方法及装置 | |
CN112723756A (zh) | 用于熔石英元件表面激光辐照损伤缺陷的组合修复方法 | |
CN109697717A (zh) | 一种基于图像自动搜索的衬砌裂缝识别方法 | |
CN101059450A (zh) | 冶金矿相智能识别方法 | |
US10365207B2 (en) | Method of automatically modifying imaging range in biological sample microscopic imaging | |
CN117392066A (zh) | 缺陷检测方法、装置、设备及存储介质 | |
CN110728286A (zh) | 一种基于火花图像的砂带磨削材料去除率识别方法 | |
CN115272362A (zh) | 一种数字病理全场图像有效区域分割方法、装置 | |
CN109604044A (zh) | 一种赤泥选铁工艺 | |
CN111986802A (zh) | 一种肺腺癌病理分化等级的辅助判定系统及方法 | |
CN111127450B (zh) | 一种基于图像的桥梁裂缝检测方法及系统 | |
CN108181324A (zh) | 一种基于机器视觉的木板表面检测方法 | |
CN111307070A (zh) | 基于数字图像处理对混凝土粗集料的棱角度测量方法 | |
CN110487789B (zh) | 基于磨片装置的岩石细观结构三维重建系统及方法 | |
CN107909567B (zh) | 数字图像的细长型连通区域提取方法 | |
CN115683962A (zh) | 一种用于选煤厂矿浆跑粗监测的矿浆粒度分析方法 | |
CN110309966B (zh) | 基于偏最小二乘法的疏浚能耗预测方法 | |
CN112906639B (zh) | 一种铬合金钢材中铁素体的图像识别方法及装置 | |
Diógenes et al. | Grain size measurement by image analysis: An application in the ceramic and in the metallic industries |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | 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: 20140219 |