一种基于Sobel边缘检测的水泥基材料CT图像中熟料相识别方法技术

一种基于Sobel边缘检测的水泥基材料CT图像中熟料相识别方法，用于水泥基材料CT图像的识别与处理。所述方法包括：对于水泥基材料的CT图像进行切割，选取代表性区域；对于待识别区域使用Sobel算子在水平和竖直方向求梯度，计算整体梯度；根据图像情况选取阈值，进行图像二值化，从而得到水泥基材料CT图像中的熟料相。本发明专利技术操作简便，对于计算消耗小，且能达到较好的CT图像识别效果，具有较强的实用价值。

A clinker phase recognition method based on Sobel edge detection in CT images of cement-based materials

A clinker phase recognition method based on Sobel edge detection in cement-based material CT image is used to recognize and process cement-based material CT image. The method includes: cutting the CT image of cement-based materials and selecting representative areas; calculating the overall gradient by using Sobel operator in horizontal and vertical directions for the areas to be identified; selecting threshold according to the image situation, binarizing the image, and obtaining clinker phase in the CT image of cement-based materials. The invention has simple operation, low calculation consumption, good CT image recognition effect and strong practical value.

【技术实现步骤摘要】
一种基于Sobel边缘检测的水泥基材料CT图像中熟料相识别方法
本专利技术属于图像处理
，特别涉及一种基于Sobel边缘检测的水泥基材料CT图像中熟料相识别方法。
技术介绍
CT扫描技术是利用X射线等对于物体进行断层扫描的技术，能够无损检测试样内部的物相组成、孔隙结构和裂纹分布，从而生成试样完整的三维模型。CT技术的这一特点使得其在材料微观结构研究方面得到了广泛应用，在水泥基材料中的应用也开始得到关注。水泥基材料微观结构复杂，包含熟料相、水化产物相和孔隙等，熟料相的分布与形貌跟水泥基材料的水化反应密切相关，而水化反应是水泥基材料强度和耐久性的基础。因此，需要准确识别水泥基材料CT图像中熟料相，分析熟料相的含量与分布情况，才能构建水泥基材料的微观结构模型。目前，对于水泥基材料CT图像尚未形成标准的处理方法，不同的图像处理方法在熟料相的识别效果和处理效率上均有或多或少的缺陷，难以有效地达到微观结构识别与重构的水平。因此急需开发一种有效的水泥基材料CT图像熟料相识别方法。
技术实现思路
为了克服上述现有技术的缺点，本专利技术的目的在于提供一种基于Sobel边缘检测的水泥基材料CT图像中熟料相识别方法，能够对水泥基材料CT图像中的熟料相进行准确、快速的识别。为了实现上述目的，本专利技术采用的技术方案是：一种基于Sobel边缘检测的水泥基材料CT图像中熟料相识别方法，其特征在于，计算水泥基材料中熟料相和其他相的灰度梯度，设定阈值，利用灰度梯度的差异进行物相识别。在计算梯度之前，先进行图像滤波，图像滤波过程可采用curvature算法。所述CT图像为灰度图像，主要滤除图像中的背景噪声，减少图像中异常的灰度值突变情况。所述滤波后的图像采用Sobel算子进行边缘检测，在检测过程中计算得到灰度梯度。所述灰度梯度包括水平梯度和竖直梯度两部分，以水平梯度和竖直梯度之平方和的算术平方根作为整体梯度。所述灰度梯度计算过程中着重考虑目标像素点上下左右的相邻四个点的灰度值贡献。所得灰度梯度需要与给定的阈值作比较，大于给定阈值的判定为物相边缘，小于给定阈值的作判定为物相内部。对于水泥基材料的熟料相识别，阈值的范围在0至1之间为宜，可根据检测结果进行调整。给定阈值之后，对图像的所有像素点计算灰度梯度并与所设阈值进行比较，得到二值化的边缘检测图像，熟料相的边缘为白色，熟料相及其他物相的内部为黑色，从而准确地识别熟料的形貌和含量。与现有技术相比，本专利技术的有益效果是：本专利技术利用水泥基材料中熟料相和其他相的灰度梯度差异进行物相识别，不直接依赖于灰度信息本身，相对于直接采用物相灰度进行区分的方法更为准确，易于计算机编程计算，具有重要的实用价值。附图说明图1是本专利技术实施流程图。图2是本专利技术算法示意图，其中(a)，像素点f(x,y)周边3×3邻域灰度值标记；(b)，水泥基材料CT图像3×3邻域；(c)，水平卷积模板；(d)，竖直卷积模板。具体实施方式下面结合附图和实施例详细说明本专利技术的实施方式。如图1和图2所示，一种基于Sobel边缘检测的水泥基材料CT图像中熟料相识别方法，其完整流程包括图像滤波、计算梯度、设定阈值、图像二值化以及熟料相识别。主要方案在于通过计算水泥基材料中熟料相和其他相的灰度梯度，设定阈值，利用灰度梯度的差异进行物相识别。具体地，图像滤波采用curvature算法进行。CT图像为灰度图像，主要滤除图像中的背景噪声。滤波后的图像采用Sobel算子进行边缘检测。检测过程需要计算灰度梯度，包括水平梯度和竖直梯度两部分。梯度的计算可以通过卷积模板实现，如图2(c)(d)，记图像的灰度分布函数为f(x,y)，其中x和y分别是图像灰度矩阵的行标和列标，函数值为像素点灰度；水平梯度记为Sx，竖直梯度记为Sy。则水平梯度和竖直梯度的计算公式为：Sx＝(a3+2a5+a8)-(a1+2a4+a6)Sy＝(a6+2a7+a8)-(a1+2a2+a3)记图像梯度为G(x,y)，则图像梯度的计算公式为：G(x,y)＝(Sx2+Sy2)1/2灰度梯度G(x,y)需要与给定的阈值T作比较，对于水泥基材料的熟料相识别，阈值的范围在0至1之间为宜，可根据检测结果进行调整。若G(x,y)>T，判定为物相边缘；若G(x,y)<T，判定为物相内部。对于图像的所有像素点计算梯度G(x,y)并与T值进行比较，即可得到二值化的边缘检测图像。熟料相的边缘为白色，熟料相及其他物相的内部为黑色，从而准确地识别熟料的形貌和含量。本文档来自技高网...

【技术保护点】
1.一种基于Sobel边缘检测的水泥基材料CT图像中熟料相识别方法，其特征在于，计算水泥基材料中熟料相和其他相的灰度梯度，设定阈值，利用灰度梯度的差异进行物相识别。

【技术特征摘要】
1.一种基于Sobel边缘检测的水泥基材料CT图像中熟料相识别方法，其特征在于，计算水泥基材料中熟料相和其他相的灰度梯度，设定阈值，利用灰度梯度的差异进行物相识别。2.根据权利要求1所述基于Sobel边缘检测的水泥基材料CT图像中熟料相识别方法，其特征在于，在计算梯度之前，先进行图像滤波，所述CT图像为灰度图像，主要滤除图像中的背景噪声。3.根据权利要求2所述基于Sobel边缘检测的水泥基材料CT图像中熟料相识别方法，其特征在于，所述图像滤波过程采用curvature算法。4.根据权利要求2所述基于Sobel边缘检测的水泥基材料CT图像中熟料相识别方法，其特征在于，所述滤波后的图像采用Sobel算子进行边缘检测，在检测过程中计算得到灰度梯度。5.根据权利要求1或4所述基于Sobel边缘检测的水泥基材料CT图像中熟料相识别方法，其特征在于，所述灰度梯度计算过程中着重考虑目标像素点上下左右的相邻四个点的灰度值贡献。6.根据权利要求1或4所述基于Sobel边缘检测的水泥基材料CT图像中熟料相识别方法，其特征在于，所述灰度梯度包括水平梯度和竖直梯度两部分，以水平梯度和竖直梯度之平方和的算术平方根作为整体梯度。7.根据权利要求...

【专利技术属性】
技术研发人员：魏亚，吴泽弘，梁思明，
申请(专利权)人：清华大学，
类型：发明
国别省市：北京,11