适用集料粘连清除的沥青混凝土断层图像建模前处理方法技术

本发明专利技术公开了一种适用集料粘连清除的沥青混凝土断层图像建模前处理方法，基于Gauss滤波算法对断层灰度图像进行像素灰度值平滑处理，可提高断层灰度图像集料‑沥青胶浆边界对比度和集料区域灰度值的稳定度；基于CT扫描曝光度环状分布规律对每张断层灰度图像进行环状分区处理再重组，可提高算法程序生成的阈值算子的针对性和适用性；基于相邻断层图像相似化拟合处理技术对断层二值图像沿采集方向构建层间联系，此方法可清除断层图像扫描过程和二值化处理过程中出现的错误粘连、错误孔洞和独立面片，后利用体绘制方法对具备良好层间联系的断层二值图像建立沥青混凝土三维模型，该模型可适用于后续细观结构特征分析与力学行为数值模拟。

【技术实现步骤摘要】
适用集料粘连清除的沥青混凝土断层图像建模前处理方法
本专利技术涉及断层图像三维重建
，确切地说涉及一种适用集料粘连清除的沥青混凝土断层图像建模前处理方法。
技术介绍
断层图像是利用断层扫描设备(ComputedTomography,CT)探头绕沥青混凝土试件体表旋转采集，或采用磁共振设备(MagneticResonance,MR)在沥青混凝土试件四周同时进行三维信息采集，经处理后沿采集顺序排列而成的层状图。断层图像能准确地显示试件内部集料和沥青胶浆的形态学特征、密度和轮廓边界，处理结果经三维重建后可用于后续细观结构特征分析与力学行为数值模拟。沥青混凝土三维重建技术可通过二值化处理技术对扫描获得连续断层图像进行图像二值化，获得二值化断层图像，后利用体绘制技术对获得的二值化断层图像建立试件的三维微细观结构模型。断层二值图像中图像的像素信息，仅由二进制码“0”和“1”表示，像素信息“0”表示该像素点显示黑暗，像素信息“1”则表示该像素点显示白亮。由多个黑暗像素点平面排列构成的黑暗区域表示沥青胶浆，由多个白亮像素点平面排列构成的白亮区域表示集料。体绘制技术是将层状图像中的像素二值信息沿初始采集方向z扩展为体素二值信息再将体像素信息按照三维坐标组合生成三维模型的技术。Gauss滤波处理技术可利用已设定的Gauss算子通过算法处理图片噪声，并平滑图像。现有三维重建技术，一般可通过对扫描获得的断层图像经滤波处理、形态学处理和图像二值化处理获得二值化图像，该步骤的目的在于筛除断层图像中与三维重建所需数据无关的杂质信息，以达到提高建模针对性和简化三维重建运算量的效果，后对二值化图像采用体绘制或面绘制技术进行三维重建。但由于断层图像受到扫描设备精度限制、扫描角度误差、图像曝光度不均匀和沥青混合料集料密度不均匀的影响，现有的断层成像技术生成的层状图像通常存在大量噪点和集料粘连。在曝光度不足区域和细小集料集中区域易出现集料-沥青胶浆边界模糊现象，造成相邻层状图片间不连续，降低层状图片二值化处理的集料分离效果。该影响将导致经体绘制获得的三维重建模型存在大量的集料错误粘连、集料错误孔洞和片状集料。申请号为CN201210172375.7的专利文献公开了一种基于沥青路面试件X-rayCT图像的集料细观实体模型重构方法，该文献基于目标识别算法对沥青路面试件断层图像中的独立集料的轮廓进行识别和建立集料边缘包围盒数组确定集料粒径的上下限阈值，根据阈值筛选出符合粒径要求的有效边缘，后建立沥青混合料细观实体模型。该方法在基于循环程序筛除无效边缘的同时，会造成生成的模型失去大量的形态学特征，其表现为生成的集料表面具有大量的规则性棱和角，且无法筛除具有有效粒径的层状独立面片。申请号为CN201610553499.8的专利文献公开了一种焦炭显微结构图的三维重建方法，该文献基于体绘制技术利用3DSlicer软件对焦炭断层序列二维显微结构图进行三维重建。采用这种直接对不具有良好层间关系的断层图像三维重建的建模方法生成的建模结果虽保留了良好的形态学特征，但建模结果中存在有大量的错误结构粘连、孔洞和独立面片等问题。以上问题可后期人工通过三维重建软件清除处理，但会大量消耗人力、物力和时间。
技术实现思路
本专利技术的针对现有技术的不足，克服现有技术的缺点，提出了一种适用集料粘连清除的沥青混凝土断层图像建模前处理方法。该方法基于Gauss滤波处理技术、图片形态学处理技术和相邻断层图像相似化拟合技术对断层图像进行预处理，构建断层图像沿采集方向的层间联系，可清除断层图像扫描过程中和二值化处理过程中出现的错误粘连、错误孔洞和独立面片，然后结合体绘制技术完成沥青混凝土试件模型三维重建。为实现上述目的，本专利技术的技术方案如下：一种适用集料粘连清除的沥青混凝土断层图像建模前处理方法，包括断层图像扫描、高像素化插值处理、灰度化处理、滤波平滑处理、形态学膨胀腐蚀处理、二值化处理、像素还原处理、相邻断层图像相似化拟合处理、体绘制建模和断层图像重生成处理等处理步骤，所有步骤均可利用Matlab循环程序或相关建模软件完成；具体步骤如下：第一步，断层图像扫描：沥青混凝土试件CT扫描，对试件进行初始状态的断层图像扫描，获得断层图像，该采集方向设为z轴正向。第二步，高像素化插值处理：基于插值方法对得到的断层图像进行高像素化插值处理，得到长宽像素值均大于原长宽像素值两倍以上的高像素断层图像。第三步，灰度化处理和滤波平滑处理：按采集顺序，对高像素断层图像进行灰度化处理得到高像素断层灰度图像；灰度化完成后，采用Gauss滤波算法对高像素断层灰度图像进行平滑处理，得到平滑后的高像素断层灰度图像。第四步，形态学膨胀腐蚀处理和二值化处理：由于单张断层图像曝光度呈同心环状不均匀分布，因此按图片大小采用环状分割方式将单张断层图像由内向外分成1个圆形和若干个圆环区域，若最外圈图像不是圆环形状，则等同于圆环形状进行分割处理，直到层状图像的所有有效部分被分割，分区域进行形态学处理再组装还原，完成单张断层图像的形态学处理；然后按采集顺序批量完成所有断层图像的二值化处理，得到高像素断层二值图像。第五步，像素还原处理：对得到的高像素断层二值图像进行像素还原处理，得到断层二值图像。第六步，相邻断层图像相似化拟合处理：采用相邻断层图像相似化拟合技术对得到的断层二值图像构建沿采集方向的层间联系，步骤如下；(1)取交膨胀处理或取并腐蚀处理，根据实际需求选用取交膨胀处理：将a张连续的断层二值图像中相同位置的像素信息值相加，相加结果按照像素位置导入到与原断层图像等像素大小的空白新图中，该空白新图的标号为原相邻两断层图像的较小值；新图像素信息值相加结果为“0”的位置像素信息不变，显示黑暗；新图像素信息值相加结果为“1”的位置像素信息修正为“0”，显示黑暗；新图像素信息值相加结果为“≥2”的位置像素信息修正为“1”，显示白亮；修正完成后对新图白亮集料区域进行膨胀处理，并删除集料内部的沥青胶浆孔洞。取并腐蚀处理：将b张连续的断层二值图像中相同位置的像素信息值相加，相加结果按照像素位置导入到与原断层图像等像素大小的空白新图中，该空白新图的标号为原相邻两断层图像的较小值；新图像素信息值相加结果为“0”的位置像素信息不变，显示黑暗；新图像素信息值相加结果为“1”的位置像素信息不变，显示白亮；新图像素信息值相加结果为“≥2”的像素信息修正为“1”，显示白亮；修正完成后对新图白亮集料区域进行腐蚀处理，并删除集料内部的沥青胶浆孔洞。(2)按照顺序完成后续相邻断层图像相似化拟合处理，并依次生成新图，新图的顺序集合构成新断层二值图像图集。第七步，体绘制建模：基于体绘制技术对第六步得到的新断层二值图像图集进行三维重建获得沥青混凝土模型。第八步，对第七步得到的模型沿任意径向进行断层图像重生成，经二值化处理得到断层二值图像。第九步，断层图像重生成处理：对第八步重生成得到的断层二值图像图集重复第六步-第八本文档来自技高网...

【技术保护点】
1.一种适用集料粘连清除的沥青混凝土断层图像建模前处理方法，所有步骤均可利用Matlab循环程序或相关建模软件完成；其特征在于，具体步骤如下：/n第一步，断层图像扫描：沥青混凝土试件CT扫描，对试件进行初始状态的断层图像扫描，获得断层图像，该采集方向设为z轴正向；/n第二步，高像素化插值处理：基于插值方法对得到的断层图像进行高像素化插值处理，得到长宽像素值均大于原长宽像素值两倍以上的高像素断层图像；/n第三步，灰度化处理和滤波平滑处理：按采集顺序，对高像素断层图像进行灰度化处理得到高像素断层灰度图像；灰度化完成后，采用Gauss滤波算法对高像素断层灰度图像进行平滑处理，得到平滑后的高像素断层灰度图像；/n第四步，形态学膨胀腐蚀处理和二值化处理：由于单张断层图像曝光度呈同心环状不均匀分布，因此按图片大小采用环状分割方式将单张断层图像由内向外分成1个圆形和若干个圆环区域，若最外圈图像不是圆环形状，则等同于圆环形状进行分割处理，直到层状图像的所有有效部分被分割，分区域进行形态学处理再组装还原，完成单张断层图像的形态学处理；然后按采集顺序批量完成所有断层图像的二值化处理，得到高像素断层二值图像；/n第五步，像素还原处理：对得到的高像素断层二值图像进行像素还原处理，得到断层二值图像；/n第六步，相邻断层图像相似化拟合处理：采用相邻断层图像相似化拟合技术对得到的断层二值图像构建沿采集方向的层间联系，步骤如下；/n(1)取交膨胀处理或取并腐蚀处理，根据实际需求选用/n取交膨胀处理：将a张连续的断层二值图像中相同位置的像素信息值相加，相加结果按照像素位置导入到与原断层图像等像素大小的空白新图中，该空白新图的标号为原相邻两断层图像的较小值；新图像素信息值相加结果为“0”的位置像素信息不变，显示黑暗；新图像素信息值相加结果为“1”的位置像素信息修正为“0”，显示黑暗；新图像素信息值相加结果为“≥2”的位置像素信息修正为“1”，显示白亮；修正完成后对新图白亮集料区域进行膨胀处理，并删除集料内部的沥青胶浆孔洞；/n取并腐蚀处理：将b张连续的断层二值图像中相同位置的像素信息值相加，相加结果按照像素位置导入到与原断层图像等像素大小的空白新图中，该空白新图的标号为原相邻两断层图像的较小值；新图像素信息值相加结果为“0”的位置像素信息不变，显示黑暗；新图像素信息值相加结果为“1”的位置像素信息不变，显示白亮；新图像素信息值相加结果为“≥2”的像素信息修正为“1”，显示白亮；修正完成后对新图白亮集料区域进行腐蚀处理，并删除集料内部的沥青胶浆孔洞；/n(2)按照顺序完成后续相邻断层图像相似化拟合处理，并依次生成新图，新图的顺序集合构成新断层二值图像图集；/n第七步，体绘制建模：基于体绘制技术对第六步得到的新断层二值图像图集进行三维重建获得沥青混凝土模型；/n第八步，对第七步得到的模型沿任意径向进行断层图像重生成，经二值化处理得到断层二值图像；/n第九步，断层图像重生成处理：对第八步重生成得到的断层二值图像图集重复第六步-第八步，生成新断层二值图像集；/n第十步，采用体绘制技术对第九步得到的新断层二值图像集进行沥青混凝土试件模型三维重建。/n...

【技术特征摘要】
1.一种适用集料粘连清除的沥青混凝土断层图像建模前处理方法，所有步骤均可利用Matlab循环程序或相关建模软件完成；其特征在于，具体步骤如下：
第一步，断层图像扫描：沥青混凝土试件CT扫描，对试件进行初始状态的断层图像扫描，获得断层图像，该采集方向设为z轴正向；
第二步，高像素化插值处理：基于插值方法对得到的断层图像进行高像素化插值处理，得到长宽像素值均大于原长宽像素值两倍以上的高像素断层图像；
第三步，灰度化处理和滤波平滑处理：按采集顺序，对高像素断层图像进行灰度化处理得到高像素断层灰度图像；灰度化完成后，采用Gauss滤波算法对高像素断层灰度图像进行平滑处理，得到平滑后的高像素断层灰度图像；
第四步，形态学膨胀腐蚀处理和二值化处理：由于单张断层图像曝光度呈同心环状不均匀分布，因此按图片大小采用环状分割方式将单张断层图像由内向外分成1个圆形和若干个圆环区域，若最外圈图像不是圆环形状，则等同于圆环形状进行分割处理，直到层状图像的所有有效部分被分割，分区域进行形态学处理再组装还原，完成单张断层图像的形态学处理；然后按采集顺序批量完成所有断层图像的二值化处理，得到高像素断层二值图像；
第五步，像素还原处理：对得到的高像素断层二值图像进行像素还原处理，得到断层二值图像；
第六步，相邻断层图像相似化拟合处理：采用相邻断层图像相似化拟合技术对得到的断层二值图像构建沿采集方向的层间联系，步骤如下；
(1)取交膨胀处理或取并腐蚀处理，根据实际需求选用
取交膨胀处理：将a张连续的断层二值图像中相同位置的像素信息值相加，相加结果按照像素位置导入到与原断层图像等像素大小的空白新图中，该空白新图的标号为原相邻两断层图像的较小值；新图像素信息值相加结果为“0”的位置像素信息不变，显示黑暗；新图像素信息值相加结果为“1”的位置像素信息修正为“0”，显示黑暗；新图像素信息值相加结果为“≥2”的位置像素信息修正为“1”，显示白亮；修正完成后对新图白亮集料区域进行膨胀处理，并删除集料内部的沥青胶浆孔洞；
取并腐蚀处理：将b张连续的断层二值图像中相同位置的像素信息值相加，相加结果按照像素位置导入到与原断层图像等像素大小的空白新图中，该空白新图的标号为原相邻两断层图像的较小值；新图像素信息值相加结果为“0”的位置像素信息不变，显示黑暗；新图像素信息值相加结果为“1”的位置像素信息不变，显示白亮；新图像素信息值相加结果为“≥2”...

【专利技术属性】
技术研发人员：孙依人，王超，张壮，
申请(专利权)人：大连理工大学，
类型：发明
国别省市：辽宁;21