本发明涉及一种用于建筑垃圾回收的再生骨料分级评价方法，属于图像处理技术领域。方法包括以下步骤：根据各像素点对应的各隶属度，得到灰度图像对应的各目标隶属度特征图；根据各目标隶属度特征图上各像素点对应的隶属度，得到各目标隶属度特征图对应的各隶属度延伸游程矩阵；根据隶属度延伸游程矩阵，得到各目标隶属度特征图对应的延伸特征指标；根据延伸特征指标和灰度图像上像素点总数量，得到再生骨料中各类型杂质的占比；根据再生骨料中各类型杂质的占比，得到再生骨料对应的质量评价级别。该方法能够相对准确的得到再生骨料中杂质含量，进而也能够相对准确的得到再生骨料的质量评价级别。
1.一种用于建筑垃圾回收的再生骨料分级评价方法，其特征在于，该方法包括如下步骤：获取再生骨料平铺灰度图像；
根据EM算法和所述灰度图像对应的灰度直方图，得到所述灰度图像上各像素点对应的各隶属度；根据所述各像素点对应的各隶属度，得到所述灰度图像对应的各目标隶属度特征图；所述目标隶属度特征图中一个像素点对应一个隶属度；
根据所述各目标隶属度特征图上各像素点对应的隶属度，得到所述各目标隶属度特征图对应的各隶属度延伸游程矩阵；根据所述隶属度延伸游程矩阵，得到各目标隶属度特征图对应的延伸特征指标；
根据所述延伸特征指标和所述灰度图像上像素点总数量，得到再生骨料中各类型杂质的占比；根据所述再生骨料中各类型杂质的占比，得到再生骨料对应的质量评价级别；
所述根据EM算法和所述灰度图像对应的灰度直方图，得到所述灰度图像上各像素点对应的各隶属度的方法，包括：对再生骨料平铺灰度图像上的灰度值进行统计，构建得到所述灰度图像对应的灰度直方图，所述灰度直方图上的横坐标为灰度值，纵坐标为分布概率；
根据所述灰度直方图和EM算法，拟合得到所述灰度图像对应的高斯混合模型；
将各子高斯模型对应的均值按照从小到大的顺序进行排序，得到各子高斯模型对应的序号；
计算所述灰度直方图各灰度值在各子高斯模型中的取值占比，记为所述灰度图像对应的灰度直方图上的各灰度值对于各子高斯模型的隶属度；
将所述各灰度值对于各子高斯模型的隶属度记为所述灰度图像对应的灰度直方图上的各灰度值对应的各隶属度；
获取所述灰度直方图上的各灰度值对应的各像素点；将所述各灰度值对应的各隶属度记为灰度直方图上对应灰度值对应的各像素点的各隶属度；
所述根据所述各像素点对应的各隶属度，得到所述灰度图像对应的各目标隶属度特征图的方法，包括：根据所述各灰度值对于相同序号子高斯模型的隶属度，构建得到所述灰度图像对应的各隶属度特征图；所述各隶属度特征图包括木片塑料杂质隶属度特征图、钢筋杂质隶属度特征图以及混凝土隶属度特征图；所述混凝土不属于再生骨料中的杂质；所述隶属度特征图的数量为再生骨料平铺灰度图像上的任一像素点对应的隶属度的数量；所述隶属度特征图和所述再生骨料平铺灰度图像中的像素点一一对应；
将所述木片塑料杂质隶属度特征图和钢筋杂质隶属度特征图记为目标隶属度特征图。
2.如权利要求1所述的一种用于建筑垃圾回收的再生骨料分级评价方法，其特征在于，根据如下公式计算所述灰度图像对应的灰度直方图上的各灰度值对于各子高斯模型的隶属度：其中， 为所述灰度直方图上的第 个灰度值对于第k个子高斯模型的隶属度，为所述灰度直方图上的第 个灰度值， 为所述高斯混合模型中第k个子高斯模型的权重值， 所述灰度直方图上的第 个灰度值在第k个子高斯模型中的取值，K为高斯混合模型中的子高斯模型个数。
3.如权利要求1所述的一种用于建筑垃圾回收的再生骨料分级评价方法，其特征在于，所述根据所述各目标隶属度特征图上各像素点对应的隶属度，得到所述各目标隶属度特征图对应的各隶属度延伸游程矩阵；根据所述隶属度延伸游程矩阵，得到各目标隶属度特征图对应的延伸特征指标的方法，包括：所述各目标隶属度特征图是指木片塑料杂质隶属度特征图和钢筋杂质隶属度特征图；
获取隶属度取值范围；将所述隶属度取值范围均分为10个隶属度级别，得到木片塑料杂质隶属度特征图上各像素点对应的隶属度级别和钢筋杂质隶属度特征图上各像素点对应的隶属度级别；
构建得到钢筋杂质隶属度特征图在0度方向下的隶属度延伸游程矩阵、在45度方向下的隶属度延伸游程矩阵、在90度方向下的隶属度延伸游程矩阵以及在135度方向下的隶属度延伸游程矩阵；所述钢筋杂质隶属度特征图对应的四个方向下的隶属度延伸游程矩阵的最大游程长度都为D，将D记为钢筋杂质隶属度特征图对应的最大延伸长度；所述钢筋杂质隶属度特征图对应的四个方向下的隶属度延伸游程矩阵的行数和列数相同；
构建得到木片塑料杂质隶属度特征图在0度方向下的隶属度延伸游程矩阵、在45度方向下的隶属度延伸游程矩阵、在90度方向下的隶属度延伸游程矩阵以及在135度方向下的隶属度延伸游程矩阵；所述木片塑料杂质隶属度特征图对应的四个方向下的隶属度延伸游程矩阵的最大游程长度都为M；所述木片塑料杂质隶属度特征图对应的四个方向下的隶属度延伸游程矩阵的行数和列数相同；
根据钢筋杂质隶属度特征图在0度方向下的隶属度延伸游程矩阵和在90度方向下的隶属度延伸游程矩阵中的参数值和游程长度，得到钢筋杂质隶属度特征图中各隶属度级别对应的第一延伸偏向；
根据钢筋杂质隶属度特征图在45度方向下的隶属度延伸游程矩阵和在135度方向下的隶属度延伸游程矩阵中的参数值和游程长度，得到钢筋杂质隶属度特征图中各隶属度级别对应的第二延伸偏向；
根据钢筋杂质隶属度特征图在0度方向下的隶属度延伸游程矩阵和在90度方向下的隶属度延伸游程矩阵中的参数值和游程长度，得到木片塑料杂质隶属度特征图中各隶属度级别对应的第一延伸偏向；
根据钢筋杂质隶属度特征图在45度方向下的隶属度延伸游程矩阵和在135度方向下的隶属度延伸游程矩阵中的参数值和游程长度，得到木片塑料杂质隶属度特征图中各隶属度级别对应的第二延伸偏向；
计算钢筋杂质隶属度特征图中各像素点隶属度的均值，记为钢筋杂质隶属度特征图对应的平均隶属度；
根据钢筋杂质隶属度特征图对应的平均隶属度、钢筋杂质隶属度特征图中各隶属度级别对应的第一延伸偏向和对应的第二延伸偏向，得到钢筋杂质隶属度特征图对应的延伸特性；
计算木片塑料杂质隶属度特征图中各像素点隶属度的均值，记为木片塑料杂质隶属度特征图对应的平均隶属度；
根据木片塑料杂质隶属度特征图对应的平均隶属度、木片塑料杂质隶属度特征图中各隶属度级别对应的第一延伸偏向和对应的第二延伸偏向，得到木片塑料杂质隶属度特征图对应的延伸特征指标。
4.如权利要求3所述的一种用于建筑垃圾回收的再生骨料分级评价方法，其特征在于，根据如下公式计算钢筋杂质隶属度特征图中各隶属度级别对应的第一延伸偏向：其中， 为钢筋杂质隶属度特征图中第g个隶属度级别对应的第一延伸偏向， 为钢筋杂质隶属度特征图对应的最大游程长度， 为钢筋杂质隶属度特征图在0度方向下的隶属度延伸游程矩阵中第g行、第d列的元素值， 为钢筋杂质隶属度特征图在90度方向下的隶属度延伸游程矩阵中第g行、第d列的元素值。
5.如权利要求3所述的一种用于建筑垃圾回收的再生骨料分级评价方法，其特征在于，根据如下公式计算钢筋杂质隶属度特征图中各隶属度级别对应的第二延伸偏向：其中， 为钢筋杂质隶属度特征图中第g个隶属度级别对应的第二延伸偏向，为钢筋杂质隶属度特征图对应的最大游程长度， 为钢筋杂质隶属度特征图在
45度方向下的隶属度延伸游程矩阵中第g行、第d列的元素值， 为钢筋杂质隶属度特征图在135度方向下的隶属度延伸游程矩阵中第g行、第d列的元素值。
6.如权利要求3所述的一种用于建筑垃圾回收的再生骨料分级评价方法，其特征在于，根据如下公式计算钢筋杂质隶属度特征图对应的延伸特征指标：其中， 为钢筋杂质隶属度特征图对应的延伸特征指标， 为钢筋杂质隶属度特征图对应的平均隶属度， 为钢筋杂质隶属度特征图中第g个隶属度级别对应的第二延伸偏向， 为钢筋杂质隶属度特征图中第g个隶属度级别对应的第一延伸偏向。
7.如权利要求3所述的一种用于建筑垃圾回收的再生骨料分级评价方法，其特征在于，根据如下公式计算木片塑料杂质隶属度特征图对应的延伸特征指标：其中， 为木片塑料杂质隶属度特征图对应的延伸特征指标， 为木片塑料杂质隶属度特征图对应的平均隶属度， 为木片塑料杂质隶属度特征图中第g个隶属度级别对应的第一延伸偏向， 为木片塑料杂质隶属度特征图中第g个隶属度级别对应的第二延伸偏向。
8.如权利要求3所述的一种用于建筑垃圾回收的再生骨料分级评价方法，其特征在于，所述根据所述延伸特征指标和所述灰度图像上像素点总数量，得到再生骨料中各类型杂质的占比的方法，包括：获取再生骨料平铺灰度图像上像素点总数量；
将木片塑料杂质隶属度特征图对应的延伸特征指标与再生骨料平铺灰度图像上像素点总数量的比值记为再生骨料中木片塑料杂质的占比；
将钢筋杂质隶属度特征图对应的延伸特征指标与再生骨料平铺灰度图像上像素点总数量的比值记为再生骨料中钢筋杂质的占比。
一种用于建筑垃圾回收的再生骨料分级评价方法
技术领域
[0001] 本发明涉及图像处理技术领域，具体涉及一种用于建筑垃圾回收的再生骨料分级评价方法。
背景技术
[0002] 随着建筑行业的快速发展，产生了大量的废弃混凝土，废弃混凝土是建筑垃圾的主要组分，并且对产生的废弃混凝土的回收利用是提高资源利用率，实现经济效益、社会效益与环保效益统一的重要手段。当前，将废弃混凝土加工成再生骨料是建筑垃圾回收利用的一种主要方式，并且加工生产出的再生骨料的质量不同，然后这些加工出来的不同质量的再生骨料可以被应用制造不同强度需求的建筑材料。因此再生骨料的质量分级对后续再生骨料的应用有重要的指导作用。
[0003] 现有的再生骨料进行评价分级方法一般是依据采样检测方式得到的再生骨料中杂质含量或者依据图像阈值分割得到的杂质含量来对再生骨料进行质量分级，一般情况下利用采样检测得到再生骨料中杂质含量的方法需要复杂的实验室指标，过程较为繁琐；而利用图像阈值分割得到杂质含量的方法一般是利用人工设置阈值对再生骨料图像进行图像分割，将杂质点与混凝土骨料点分开，统计杂质点数量的占比作为杂质含量，而有些杂质点与骨料点相似，单单以像素的灰度信息无法准确将杂质点与混凝土骨料点很好的区分，因此可能导致统计结果存在较大的误差，进而导致后续对再生骨料质量进行评价分级的准确度较低。
发明内容
[0004] 本发明提供一种用于建筑垃圾回收的再生骨料分级评价方法，用于解决现有方法对再生骨料质量进行评价分级准确度较低的问题，所采用的技术方案具体如下：
[0005] 本发明提供了一种用于建筑垃圾回收的再生骨料分级评价方法包括以下步骤：
[0006] 获取再生骨料平铺灰度图像；
[0007] 根据EM算法和所述灰度图像对应的灰度直方图，得到所述灰度图像上各像素点对应的各隶属度；根据所述各像素点对应的各隶属度，得到所述灰度图像对应的各目标隶属度特征图；所述目标隶属度特征图中一个像素点对应一个隶属度；
[0008] 根据所述各目标隶属度特征图上各像素点对应的隶属度，得到所述各目标隶属度特征图对应的各隶属度延伸游程矩阵；根据所述隶属度延伸游程矩阵，得到各目标隶属度特征图对应的延伸特征指标；
[0009] 根据所述延伸特征指标和所述灰度图像上像素点总数量，得到再生骨料中各类型杂质的占比；根据所述再生骨料中各类型杂质的占比，得到再生骨料对应的质量评价级别。
[0010] 有益效果：本发明依据再生骨料平铺灰度图像对应的灰度直方图和EM算法，得到再生骨料平铺灰度图像上各像素点对应的各隶属度；然后通过对各像素点对应的各隶属度进行分析，得到再生骨料平铺灰度图像对应的各目标隶属度特征图；并且隶属度特征图中一个像素点对应一个隶属度；又根据各目标隶属度特征图上各像素点对应的隶属度，得到各目标隶属度特征图对应的各隶属度延伸游程矩阵；之后根据隶属度延伸游程矩阵，得到各目标隶属度特征图对应的延伸特征指标；最后根据延伸特征指标和再生骨料平铺灰度图像上像素点总数量，得到再生骨料中各类型杂质的占比；根据得到的再生骨料中各类型杂质的占比，得到再生骨料对应的质量级别。本发明提供的用于建筑垃圾回收的再生骨料分级评价方法能够相对准确的得到再生骨料中杂质含量，进而也能够相对准确的得到再生骨料的质量评价级别，并且该方法与需要复杂采样才能检测得到再生骨料中杂质含量的方式相比，具有快速、便捷的特点。
[0011] 优选的，根据EM算法和所述灰度图像对应的灰度直方图，得到所述灰度图像上各像素点对应的各隶属度的方法，包括：
[0012] 对再生骨料平铺灰度图像上的灰度值进行统计，构建得到所述灰度图像对应的灰度直方图，所述灰度直方图上的横坐标为灰度值，纵坐标为分布概率；
[0013] 根据所述灰度直方图和EM算法，拟合得到所述灰度图像对应的高斯混合模型；
[0014] 将子高斯模型对应的均值按照从小到大的顺序进行排序，得到各子高斯模型对应的序号；
[0015] 计算所述灰度直方图各灰度值在各子高斯模型中的取值占比，记为所述灰度图像对应的灰度直方图上的各灰度值对于各子高斯模型的隶属度；
[0016] 将所述各灰度值对于各子高斯模型的隶属度记为所述灰度图像对应的灰度直方图上的各灰度值对应的各隶属度；
[0017] 获取所述灰度直方图上的各灰度值对应的各像素点；将所述各灰度值对应的各隶属度记为灰度直方图上对应灰度值对应的各像素点的各隶属度。
[0018] 优选的，根据如下公式计算所述灰度图像对应的灰度直方图上的各灰度值对于各子高斯模型的隶属度：
[0019]
[0020] 其中， 为所述灰度直方图上的第 个灰度值对于第k个子高斯模型的隶属度，为所述灰度直方图上的第 个灰度值， 为所述高斯混合模型中第 个子高斯模型的权重值， 为所述灰度直方图上的第 个灰度值在第k个子高斯模型中的取值，K为高斯混合模型对应的子高斯模型个数。
[0021] 优选的，根据所述各像素点对应的各隶属度，得到所述灰度图像对应的各目标隶属度特征图的方法，包括：
[0022] 根据所述各灰度值对于相同序号子高斯模型的隶属度，构建得到所述灰度图像对应的各隶属度特征图；所述各隶属度特征图包括木片塑料杂质隶属度特征图、钢筋杂质隶属度特征图以及混凝土隶属度特征图；所述混凝土不属于再生骨料中的杂质；所述隶属度特征图的数量为再生骨料平铺灰度图像上的任一像素点对应的隶属度的数量；所述隶属度特征图和所述再生骨料平铺灰度图像中的像素点一一对应；
[0023] 将所述木片塑料杂质隶属度特征图和钢筋杂质隶属度特征图记为目标隶属度特征图。
[0024] 优选的，根据所述各目标隶属度特征图上各像素点对应的隶属度，得到所述各目标隶属度特征图对应的各隶属度延伸游程矩阵；根据所述隶属度延伸游程矩阵，得到各目标隶属度特征图对应的延伸特征指标的方法，包括：
[0025] 所述各目标隶属度特征图是指木片塑料杂质隶属度特征图和钢筋杂质隶属度特征图；
[0026] 获取隶属度取值范围；将所述隶属度取值范围均分为10个隶属度级别，得到木片塑料杂质隶属度特征图上各像素点对应的隶属度级别和钢筋杂质隶属度特征图上各像素点对应的隶属度级别；
[0027] 构建得到钢筋杂质隶属度特征图在0度方向下的隶属度延伸游程矩阵、在45度方向下的隶属度延伸游程矩阵、在90度方向下的隶属度延伸游程矩阵以及在135度方向下的隶属度延伸游程矩阵；所述钢筋杂质隶属度特征图对应的四个方向下的隶属度延伸游程矩阵的最大游程长度都为D，将D记为钢筋杂质隶属度特征图对应的最大延伸长度；所述钢筋杂质隶属度特征图对应的四个方向下的隶属度延伸游程矩阵的行数和列数相同；
[0028] 构建得到木片塑料杂质隶属度特征图在0度方向下的隶属度延伸游程矩阵、在45度方向下的隶属度延伸游程矩阵、在90度方向下的隶属度延伸游程矩阵以及在135度方向下的隶属度延伸游程矩阵；所述木片塑料杂质隶属度特征图对应的四个方向下的隶属度延伸游程矩阵的最大游程长度都为M；所述木片塑料杂质隶属度特征图对应的四个方向下的隶属度延伸游程矩阵的行数和列数相同；
[0029] 根据钢筋杂质隶属度特征图在0度方向下的隶属度延伸游程矩阵和在90度方向下的隶属度延伸游程矩阵中的参数值和游程长度，得到钢筋杂质隶属度特征图中各隶属度级别对应的第一延伸偏向；
[0030] 根据钢筋杂质隶属度特征图在45度方向下的隶属度延伸游程矩阵和在135度方向下的隶属度延伸游程矩阵中的参数值和游程长度，得到钢筋杂质隶属度特征图中各隶属度级别对应的第二延伸偏向；
[0031] 根据钢筋杂质隶属度特征图在0度方向下的隶属度延伸游程矩阵和在90度方向下的隶属度延伸游程矩阵中的参数值和游程长度，得到木片塑料杂质隶属度特征图中各隶属度级别对应的第一延伸偏向；
[0032] 根据钢筋杂质隶属度特征图在45度方向下的隶属度延伸游程矩阵和在135度方向下的隶属度延伸游程矩阵中的参数值和游程长度，得到木片塑料杂质隶属度特征图中各隶属度级别对应的第二延伸偏向；
[0033] 计算钢筋杂质隶属度特征图中各像素点隶属度的均值，记为钢筋杂质隶属度特征图对应的平均隶属度；
[0034] 根据钢筋杂质隶属度特征图对应的平均隶属度、钢筋杂质隶属度特征图中各隶属度级别对应的第一延伸偏向和对应的第二延伸偏向，得到钢筋杂质隶属度特征图对应的延伸特性；
[0035] 计算木片塑料杂质隶属度特征图中各像素点隶属度的均值，记为木片塑料杂质隶属度特征图对应的平均隶属度；
[0036] 根据木片塑料杂质隶属度特征图对应的平均隶属度、木片塑料杂质隶属度特征图中各隶属度级别对应的第一延伸偏向和对应的第二延伸偏向，得到木片塑料杂质隶属度特征图对应的延伸特征指标。
[0037] 优选的，根据如下公式计算钢筋杂质隶属度特征图中各隶属度级别对应的第一延伸偏向：
[0038]
[0039] 其中， 为钢筋杂质隶属度特征图中第g个隶属度级别对应的第一延伸偏向，为钢筋杂质隶属度特征图对应的最大游程长度， 为钢筋杂质隶属度特征图在0度方向下的隶属度延伸游程矩阵中第g行、第d列的元素值， 为钢筋杂质隶属度特征图在90度方向下的隶属度延伸游程矩阵中第g行、第d列的元素值。
[0040] 优选的，根据如下公式计算钢筋杂质隶属度特征图中各隶属度级别对应的第二延伸偏向：
[0041]
[0042] 其中， 为钢筋杂质隶属度特征图中第g个隶属度级别对应的第二延伸偏向，为钢筋杂质隶属度特征图对应的最大游程长度， 为钢筋杂质隶属度特征图在45度方向下的隶属度延伸游程矩阵中第g行、第d列的元素值， 为钢筋杂质隶属度特征图在135度方向下的隶属度延伸游程矩阵中第g行、第d列的元素值。
[0043] 优选的，根据如下公式计算钢筋杂质隶属度特征图对应的延伸特征指标：
[0044]
[0045] 其中， 为钢筋杂质隶属度特征图对应的延伸特征指标， 为钢筋杂质隶属度特征图对应的平均隶属度， 为钢筋杂质隶属度特征图中第g个隶属度级别对应的第二延伸偏向， 为钢筋杂质隶属度特征图中第g个隶属度级别对应的第一延伸偏向。
[0046] 优选的，根据如下公式计算木片塑料杂质隶属度特征图对应的延伸特征指标：
[0047]
[0048] 其中， 为木片塑料杂质隶属度特征图对应的延伸特征指标， 为木片塑料杂质隶属度特征图对应的平均隶属度， 为木片塑料杂质隶属度特征图中第g个隶属度级别对应的第一延伸偏向， 为木片塑料杂质隶属度特征图中第g个隶属度级别对应的第二延伸偏向。
[0049] 优选的，根据所述延伸特征指标和所述灰度图像上像素点总数量，得到再生骨料中各类型杂质的占比的方法，包括：
[0050] 获取再生骨料平铺灰度图像上像素点总数量；
[0051] 将木片塑料杂质隶属度特征图对应的延伸特征指标与再生骨料平铺灰度图像上像素点总数量的比值记为再生骨料中木片塑料杂质的占比；
[0052] 将钢筋杂质隶属度特征图对应的延伸特征指标与再生骨料平铺灰度图像上像素点总数量的比值记为再生骨料中钢筋杂质的占比。
附图说明
[0053] 为了更清楚地说明本发明实施例或现有技术中的技术方案和优点，下面将对实施例或现有技术描述中所需要使用的附图作简单的介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其它附图。
[0054] 图1为本发明一种用于建筑垃圾回收的再生骨料分级评价方法的流程图。
具体实施方式
[0055] 下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整的描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例，基于本发明实施例，本领域普通技术人员所获得的所有其它实施例，都属于本发明实施例保护的范围。
[0056] 除非另有定义，本文所使用的所有的技术和科学技术语与属于本发明的技术人员通常理解的含义相同。
[0057] 本实施例提供了一种用于建筑垃圾回收的再生骨料分级评价方法，详细说明如下：
[0058] 如图1所示，该一种用于建筑垃圾回收的再生骨料分级评价方法，包括以下步骤：
[0059] 步骤S001，获取再生骨料平铺灰度图像。
[0060] 本实施例主要依据再生骨料中杂质含量或者杂质占比对再生骨料的质量进行评价分级；本实施例依据灰度直方图和EM算法进行分类；然后基于分类结果计算再生骨料平铺图像上各像素点对应的各隶属度；然后基于各像素点对应的各隶属度，得到再生骨料平铺灰度图像对应的各目标隶属度特征图；之后分析得到各目标隶属度特征图的延伸特征；其次根据各目标隶属度特征图对应的延伸特性和再生骨料平铺图像上像素点数量，得到再生骨料中各类型杂质的占比；最后根据再生骨料中各类型杂质的占比，得到再生骨料对应的质量评价级别。本实施例提供的用于建筑垃圾回收的再生骨料分级评价方法能够相对准确的得到再生骨料中杂质占比，进而也能够相对准确的得到再生骨料的质量评价级别，并且该方法与需要复杂采样才能检测得到再生骨料中杂质含量的方式相比，具有快速、便捷的特点。
[0061] 本实施例在废弃混凝土作为原料所生产的再生骨料生产线出料口处将生产出的再生骨料进行平铺，利用工业相机获取再生骨料生产线出料口处将生产出的再生骨料进行平铺后的图像，记为再生骨料平铺图像；对再生骨料平铺图像进行灰度化处理，得到再生骨料平铺灰度图像。
[0062] 步骤S002，根据EM算法和所述灰度图像对应的灰度直方图，得到所述灰度图像上各像素点对应的各隶属度；根据所述各像素点对应的各隶属度，得到所述灰度图像对应的各目标隶属度特征图；所述目标隶属度特征图中一个像素点对应一个隶属度。
[0063] 本实施例利用再生骨料平铺灰度图像对应的灰度直方图和EM算法实现对再生骨料平铺灰度图像上的像素点分类，并得到再生骨料平铺灰度图像上各像素点对于各类别的隶属度，即得到再生骨料平铺灰度图像上各像素点对应的各隶属度；后续基于再生骨料平铺灰度图像上各像素点对于各类别的隶属度，得到各目标隶属度特征图；所述目标隶属度特征图可以反映再生骨料中的杂质含量。具体为：
[0064] 对再生骨料平铺灰度图像上的灰度值进行统计，构建得到再生骨料平铺灰度图像对应的灰度直方图，所述灰度直方图上的横坐标为灰度值，纵坐标为分布概率；以所述灰度值和灰度值对应的分布概率为样本数据，根据样本数据并利用EM算法拟合再生骨料平铺灰度图像对应的高斯混合模型，本实施例用该高斯混合模型描述各灰度值在再生骨料平铺灰度图像上出现的概率，记为目标概率，这个概率是由K个子高斯模型的计算结果乘以其对应的权重得到的，并且高斯混合模型中的子高斯模型个数为K；本实施例将K设置为3，主要是为了区分再生骨料平铺灰度图像中的混凝土点、钢筋杂质点、木料塑料杂质点；然后将K个子高斯模型对应的均值按照从小到大的顺序进行排序，得到序号1对应的子高斯模型、序号2对应的子高斯模型以及序号3对应的子高斯模型，即一个子高斯模型对应一个序号；对于一个新输入的灰度值，即：
[0065]
[0066] 其中， 为再生骨料平铺灰度图像对应的灰度直方图上的第 个灰度值，为再生骨料平铺灰度图像对应的灰度直方图上的第 个灰度值在所述灰度图像上出现的概率， 为再生骨料平铺灰度图像对应的高斯混合模型中第 个子高斯模型的权重值，为再生骨料平铺灰度图像对应的灰度直方图上的第 个灰度值在第 个子高斯模型中的取值。
[0067] 由于钢筋杂质像素点灰度值一般比混凝土像素点灰度值低，而木片塑料杂质像素点灰度值一般比混凝土像素点灰度值高，因此序号为1的子高斯模型对应的是钢筋杂质像素点灰度值的分布，序号为2的子高斯模型对应的是混凝土像素点灰度值的分布，序号为3的子高斯模型对应的是木片塑料杂质像素点灰度值的分布；即序号1，2，3分别对应为：钢筋杂质点、混凝土点、木片塑料杂质点。
[0068] 对于再生骨料平铺灰度图像对应的灰度直方图上的各灰度值，根据如下公式计算再生骨料平铺灰度图像对应的灰度直方图上的各灰度值对于各子高斯模型的隶属度（各灰度值在各子高斯模型中的取值占比）：
[0069]
[0070] 其中， 为再生骨料平铺灰度图像对应的灰度直方图上的第 个灰度值对于第k个子高斯模型的隶属度，即第 个灰度值在第k个子高斯模型中的取值占比，即。
[0071] 因此通过上述方式可以得到再生骨料平铺灰度图像对应的灰度直方图上的各灰度值对于各子高斯模型的隶属度，记为再生骨料平铺灰度图像对应的灰度直方图上的各灰度值对应的各隶属度，一个灰度值对应3个隶属度；获取灰度直方图上的各灰度值对应的各像素点；将各灰度值对应的各隶属度记为灰度直方图上对应灰度值对应的各像素点的各隶属度，即一个像素点对应3个隶属度。
[0072] 本实施例根据再生骨料平铺灰度图像上的各像素点对应的各隶属度，构建得到再生骨料平铺灰度图像对应的各隶属度特征图；所述隶属度特征图的数量为再生骨料平铺灰度图像上的任一像素点对应的隶属度的数量，本实施例隶属度特征图的数量为3。构建得到各隶属度特征图的具体过程为：
[0073] 再生骨料平铺灰度图像上的各像素点对应的各隶属度分别为：再生骨料平铺灰度图像上的第 个像素点对于序号为1的子高斯模型的隶属度（第 个像素点对应的灰度值在序号为1的子高斯模型中的取值占比）、第 个像素点对于序号为2的子高斯模型的隶属度（第 个像素点对应的灰度值在序号为2的子高斯模型中的取值占比）、第 个像素点对于序号为3的子高斯模型的隶属度（第 个像素点对应的灰度值在序号为3的子高斯模型中的取值占比）。然后根据再生骨料平铺灰度图像上的各像素点对于序号为1的子高斯模型的隶属度，构建得到再生骨料平铺灰度图像对应的第一隶属度特征图，将所述第一隶属度特征图记为钢筋杂质隶属度特征图；根据再生骨料平铺灰度图像上的各像素点对于序号为2的子高斯模型的隶属度，构建得到再生骨料平铺灰度图像对应的第二隶属度特征图，将所述第二隶属度特征图记为混凝土隶属度特征图；根据再生骨料平铺灰度图像上的各像素点对于序号为3的子高斯模型的隶属度，构建得到再生骨料平铺灰度图像对应的第三隶属度特征图，将所述第三隶属度特征图记为木片塑料杂质隶属度特征图；所述隶属度特征图和再生骨料平铺灰度图像中的像素点一一对应。
[0074] 由于本实施例主要分析的是再生骨料中的杂质含量，而混凝土不属于杂质，因此本实施例不对混凝土隶属度特征图进行分析，只对木片塑料杂质隶属度特征图和钢筋杂质隶属度特征图进行分析，将所述木片塑料杂质隶属度特征图和钢筋杂质隶属度特征图记为目标隶属度特征图。
[0075] 步骤S003，根据所述各目标隶属度特征图上各像素点对应的隶属度，得到所述各目标隶属度特征图对应的各隶属度延伸游程矩阵；根据所述隶属度延伸游程矩阵，得到各目标隶属度特征图对应的延伸特征指标。
[0076] 不同类别杂质含有不同的延伸特征，例如钢筋杂质本身具有细长的特点，而其被粉碎后作为再生骨料中的杂质其本身也保留着一定的延伸特性，即钢筋杂质总是在一定方向具有短暂的延伸特性；而木片塑料杂质本身不具备延伸特征，其被粉碎后倾向于块状、片状结构；因此以这两种杂质的延伸特征为基础，构建木片塑料杂质隶属度特征图和钢筋杂质隶属度特征图对应的隶属度延伸游程矩阵，后续基于隶属度延伸游程矩阵得到木片塑料杂质隶属度特征图和钢筋杂质隶属度特征图对应的延伸特征指标，后续基于延伸特性得到再生骨料中杂质的占比。具体为：
[0077] 由于隶属度为0到1之间的数值，则按照隶属度取值范围，可以将隶属度均分为10个隶属度级别，级别1对应的隶属度范围为 ，级别2对应的隶属度范围 ，以此类推，级别10对应的隶属度范围为 ；因此可以得到木片塑料杂质隶属度特征图上各像素点对应的隶属度级别和钢筋杂质隶属度特征图上各像素点对应的隶属度级别，并且所述木片塑料杂质隶属度特征图和钢筋杂质隶属度特征图对应的隶属度级别范围都为1到10级。
[0078] 之后构建得到钢筋杂质隶属度特征图在各 方向下的隶属度延伸游程矩阵，所述的取值为0度，45度，90度和135度，即得到钢筋杂质隶属度特征图在0度方向下的隶属度延伸游程矩阵、在45度方向下的隶属度延伸游程矩阵、在90度方向下的隶属度延伸游程矩阵以及在135度方向下的隶属度延伸游程矩阵；所述钢筋杂质隶属度特征图对应的四个方向下的隶属度延伸游程矩阵的最大游程长度都为D，将D记为钢筋杂质隶属度特征图对应的最大延伸长度，并且钢筋杂质隶属度特征图对应的四个方向下的隶属度延伸游程矩阵的行数和列数相同；并且矩阵中第g行、第d列的元素值表示为隶属度级别为g、游程长度为d的游程在钢筋杂质隶属度特征图上出现的次数。同理，得到木片塑料杂质隶属度特征图在0度方向下的隶属度延伸游程矩阵、在45度方向下的隶属度延伸游程矩阵、在90度方向下的隶属度延伸游程矩阵以及在135度方向下的隶属度延伸游程矩阵；同样的木片塑料杂质隶属度特征图对应的四个方向下的隶属度延伸游程矩阵的最大游程长度都为M，并且对应的四个方向下的隶属度延伸游程矩阵的行数和列数相同。
[0079] 对于任一隶属度特征图在四个方向下的隶属度延伸游程矩阵，该隶属度特征图在相互垂直的两个方向下的游程长度差异，代表了其延伸特性，即若两个相互垂直方向的延伸特征差异不同，表示这些像素点集中于某一方向而在与其垂直的方向上的延伸很小，符合钢筋杂质的杆状特性；如果若两个相互垂直方向的延伸特征差异相同，则说明符合木片塑料杂质的片状、块状特征；因此本实施例根据钢筋杂质隶属度特征图在0度方向下的隶属度延伸游程矩阵和在90度方向下的隶属度延伸游程矩阵中的参数值和游程长度，得到钢筋杂质隶属度特征图中各隶属度级别对应的第一延伸偏向；根据如下公式计算钢筋杂质隶属度特征图中各隶属度级别对应的第一延伸偏向：
[0080]
[0081] 其中， 为钢筋杂质隶属度特征图中第g个隶属度级别对应的第一延伸偏向，为钢筋杂质隶属度特征图对应的最大游程长度， 为钢筋杂质隶属度特征图在0度方向下的隶属度延伸游程矩阵中第g行、第d列的元素值， 为钢筋杂质隶属度特征图在90度方向下的隶属度延伸游程矩阵中第g行、第d列的元素值； 的值越大，表明钢筋杂质隶属度特征图在0度和90度两个相互垂直方向的延伸特征差异越大；为钢筋杂质隶属度特征图在0°方向下的隶属度延伸游程矩
阵内依据游程长度加权而得到的延伸程度，即钢筋杂质隶属度特征图在0°方向下的隶属度延伸游程矩阵内，游程越集中在长游程处，说明钢筋杂质点在0°方向延伸程度越大；
为钢筋杂质隶属度特征图在90°方向下的隶属度延伸游程
矩阵内依据游程长度加权而得到的延伸程度，即钢筋杂质隶属度特征图在90°方向下的隶属度延伸游程矩阵内，游程越集中在长游程处，说明钢筋杂质点在90°方向延伸程度越大。
[0082] 根据钢筋杂质隶属度特征图在45度方向下的隶属度延伸游程矩阵和在135度方向下的隶属度延伸游程矩阵中的参数值和游程长度，得到钢筋杂质隶属度特征图中各隶属度级别对应的第二延伸偏向；根据如下公式计算钢筋杂质隶属度特征图中各隶属度级别对应的第二延伸偏向：
[0083]
[0084] 其中， 为钢筋杂质隶属度特征图中第g个隶属度级别对应的第二延伸偏向，为钢筋杂质隶属度特征图对应的最大游程长度， 为钢筋杂质隶属度特征图在45度方向下的隶属度延伸游程矩阵中第g行、第d列的元素值， 为钢筋杂质隶属度特征图在135度方向下的隶属度延伸游程矩阵中第g行、第d列的元素值； 的值越大，表明钢筋杂质隶属度特征图在45度和135度两个相互垂直方向的延伸特征差异越大。
[0085] 本实施例中通过上述过程可以得到钢筋杂质隶属度特征图中各隶属度级别对应的第一延伸偏向和钢筋杂质隶属度特征图中各隶属度级别对应的第二延伸偏向；同理，按照上述计算钢筋杂质隶属度特征图中各隶属度级别对应的第一延伸偏向和对应的第二延伸偏向的方式，可以得到木片塑料杂质隶属度特征图中各隶属度级别对应的第一延伸偏向和木片塑料杂质隶属度特征图中各隶属度级别对应的第二延伸偏向。
[0086] 由于钢筋杂质的特性为两个相互垂直方向的延伸特征差异不同，即钢筋杂质隶属度特征图中各隶属度级别对应的第一延伸偏向和对应的第二延伸偏向的值越大，表明再生骨料中钢筋杂质越多。本实施例计算钢筋杂质隶属度特征图中各像素点隶属度的均值，记为钢筋杂质隶属度特征图对应的平均隶属度；根据钢筋杂质隶属度特征图对应的平均隶属度、钢筋杂质隶属度特征图中各隶属度级别对应的第一延伸偏向和钢筋杂质隶属度特征图中各隶属度级别对应的第二延伸偏向，得到钢筋杂质隶属度特征图对应的延伸特性；根据如下公式计算钢筋杂质隶属度特征图对应的延伸特征指标：
[0087]
[0088] 其中， 为钢筋杂质隶属度特征图对应的延伸特征指标， 为钢筋杂质隶属度特征图对应的平均隶属度； 越大，表明再生骨料中存在越多的钢筋杂质。
[0089] 由于木片塑料杂质的特性为两个相互垂直方向的延伸特征差异相同，即木片塑料隶属度特征图中各隶属度级别对应的第一延伸偏向和对应的第二延伸偏向的值越小，表明再生骨料中木片塑料杂质越多。因此本实施例计算木片塑料杂质隶属度特征图中各像素点隶属度的均值，记为木片塑料杂质隶属度特征图对应的平均隶属度；根据木片塑料杂质隶属度特征图对应的平均隶属度、木片塑料杂质隶属度特征图中各隶属度级别对应的第一延伸偏向和木片塑料杂质隶属度特征图中各隶属度级别对应的第二延伸偏向，得到木片塑料杂质隶属度特征图对应的延伸特征指标；根据如下公式计算木片塑料杂质隶属度特征图对应的延伸特征指标：
[0090]
[0091] 其中， 为木片塑料杂质隶属度特征图对应的延伸特征指标， 为木片塑料杂质隶属度特征图对应的平均隶属度， 为木片塑料杂质隶属度特征图中第g个隶属度级别对应的第一延伸偏向， 为木片塑料杂质隶属度特征图中第g个隶属度级别对应的第二延伸偏向； 越大，表明再生骨料中存在越多的木片塑料杂质越多。
[0092] 步骤S004，根据所述延伸特性和所述再生骨料平铺图像上像素点数量，得到再生骨料中各类型杂质的占比；根据所述再生骨料中各类型杂质的占比，得到所述再生骨料对应的级别。
[0093] 由于木片塑料杂质隶属度特征图对应的延伸特征指标可以反映再生骨料中木片塑料杂质的占比，钢筋杂质隶属度特征图对应的延伸特征指标可以反映再生骨料中钢筋杂质的占比；因此本实施例通过对木片塑料杂质隶属度特征图和钢筋杂质隶属度特征图对应的延伸特征指标进行分析，得到再生骨料中木片塑料杂质的占比和再生骨料中钢筋杂质的占比。具体为：
[0094] 获取再生骨料平铺灰度图像上像素点总数量；将木片塑料杂质隶属度特征图对应的延伸特征指标与再生骨料平铺灰度图像上像素点总数量的比值记为再生骨料中木片塑料杂质的占比；将钢筋杂质隶属度特征图对应的延伸特征指标与再生骨料平铺灰度图像上像素点总数量的比值记为再生骨料中钢筋杂质的占比。然后获取再生骨料中木片塑料杂质的占比和再生骨料中钢筋杂质的占比的和，记为再生骨料中杂质总占比；若再生骨料中杂质总占比小于预设第二阈值，则判定该再生骨料的质量为第一级别，即最优级别；若再生骨料中钢筋杂质的占比小于预设第一阈值，但是再生骨料中木片塑料杂质的占比大于第二阈值小于第三阈值，则判定该再生骨料的质量为第二级别，即次优级别；若再生骨料中钢筋杂质的占比大于第一阈值，但是再生骨料中木片塑料杂质的占比大于第三阈值小于第四阈值，则判定该再生骨料的质量为第三级别；若再生骨料中杂质总占比大于预设第四阈值，则判定该再生骨料的质量为第四级别。
[0095] 本实施例设置预设第一阈值为3%，预设第二阈值为5%，预设第三阈值为10%，预设第四阈值为20%；作为其它的实施方式，也可以根据实际情况设置不同的再生骨料的质量分级方式。
[0096] 有益效果：本实施例依据再生骨料平铺灰度图像对应的灰度直方图和EM算法，得到再生骨料平铺灰度图像上各像素点对应的各隶属度；然后通过对各像素点对应的各隶属度进行分析，得到再生骨料平铺灰度图像对应的各目标隶属度特征图；并且隶属度特征图中一个像素点对应一个隶属度；又根据各目标隶属度特征图上各像素点对应的隶属度，得到各目标隶属度特征图对应的各隶属度延伸游程矩阵；之后根据隶属度延伸游程矩阵，得到各目标隶属度特征图对应的延伸特征指标；最后根据延伸特征指标和再生骨料平铺灰度图像上像素点总数量，得到再生骨料中各类型杂质的占比；根据得到的再生骨料中各类型杂质的占比，得到再生骨料对应的质量级别。本实施例提供的用于建筑垃圾回收的再生骨料分级评价方法能够相对准确的得到再生骨料中杂质含量，进而也能够相对准确的得到再生骨料的质量评价级别，并且该方法与需要复杂采样才能检测得到再生骨料中杂质含量的方式相比，具有快速、便捷的特点。
[0097] 以上所述实施例仅用以说明本申请的技术方案，而非对其限制；尽管参照前述实施例对本申请进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本申请各实施例技术方案的精神和范围，均应包含在本申请的保护范围之内。
