一种残膜回收机防缠绕挑膜装置的制 一种秧草收获机用电力驱动行走机构

一种建材装备钣金件缺陷检测方法与流程

2023-01-15 06:00:58 来源:中国专利 TAG:


1.本发明涉及数据处理领域,具体涉及一种建材装备钣金件缺陷检测方法。


背景技术:

2.在建材装备钣金件生产及装配过程中,建材装备钣金件表面可能会出现一些缺陷,钣金件表面缺陷不仅影响产品的美观和舒适度,而且一般也会对其使用性能带来不良影响,所以生产企业对钣金件产品表面缺陷检测非常重视,以便及时发现缺陷产品,从而有效控制产品质量。因此对钣金件产品表面进行缺陷检测是质量检测中必不可少的步骤。
3.以往钣金件的缺陷检测方法为通过图像采集设备获取钣金件的表面图像,并对钣金件的表面图像进行灰度化处理得到其表面图像的灰度图像,通过比较灰度图像中像素点的灰度值与预设阈值的大小,将钣金件的表面图像中的像素点分为正常像素点和异常像素点,由异常像素点得到表面图像中的缺陷区域。但是由于钣金件的表面图像的采集环境影响,例如光线问题,会导致表面图像中像素点的灰度值存在采集误差,使得仅根据预设阈值划分的正常像素点和异常像素点存在偏差,从而不能实现对像素点的精确类别划分,导致提取出的缺陷区域不准确,难以有效的判断钣金件产品的质量。


技术实现要素:

4.为了解决上述技术问题,本发明的目的在于提供一种建材装备钣金件缺陷检测方法,所采用的技术方案具体如下:本发明提出一种建材装备钣金件缺陷检测方法,该方法包括:采集钣金件的表面图像,得到表面图像对应的灰度图像;构建标准高斯模型,将灰度图像中每个像素点的灰度值代入所述标准高斯模型中得到的结果作为所述表面图像中对应像素点的置信度;基于每个像素点的置信度将所述表面图像中的所有像素点划分为正常像素点、异常像素点和待定像素点;以任意一个待定像素点为中心获取对应待定像素点的局部区域,根据所述局部区域中每个像素点的置信度获取对应待定像素点的局部特征值,基于局部特征值将待定像素点划分为第一类待定像素点和第二类待定像素点;对每个第一类待定像素点的置信度进行优化得到优化置信度;根据第一类待定像素点的优化置信度将第一类待定像素点划分为正常像素点和异常像素点;根据每个第二类待定像素点的置信度计算所有第二类待定像素点之间的置信度分布指标;基于归一化置信度分布指标将第二类待定像素点划分正常像素点和异常像素点;将所述表面图像中的异常像素点的灰度值置为1,正常像素点的灰度值置为0,对灰度值为1的像素点进行连通域分析,得到至少一个缺陷连通域,基于缺陷连通域中每个像素点的置信度计算钣金件的表面缺陷评估值。
5.进一步地,所述基于每个像素点的置信度将所述表面图像中的所有像素点划分为正常像素点、异常像素点和待定像素点包括:
设置第一置信度阈值和第二置信度阈值,将置信度低于第一置信度阈值的像素点归为所述表面图像的异常像素点,将置信度高于第二置信度阈值的像素点归为所述表面图像的正常像素点,将置信度处于所述第一置信度阈值与所述第二置信度阈值之间的像素点归为所述表面图像的待定像素点。
6.进一步地,所述根据所述局部区域中每个像素点的置信度构建对应待定像素点的局部特征值包括:局部特征值的公式为:式中,c为待定像素点,为待定像素点c的局部区域内的像素点类别的数量,为待定像素点c的局部区域内第n种像素点类别,为待定像素点c的局部区域内第n种像素点类别的像素点数量在局部区域中所有像素点数量的占比,为待定像素点c的局部特征值,为以常数e为底数的对数函数。
7.进一步地,所述基于局部特征值将待定像素点划分为第一类待定像素点和第二类待定像素点包括:设置局部特征值阈值,将局部特征值低于等于局部特征值阈值对应的所述待定像素点归为所述表面图像的第一类待定像素点,将局部特征值高于局部特征值阈值对应的所述待定像素点归为所述表面图像的第二类待定像素点。
8.进一步地,所述对每个第一类待定像素点的置信度进行优化得到优化置信度包括:通过置信度优化模型对第一类待定像素点的置信度进行优化,所述置信度优化模型的公式为:式中,a为第一类待定像素点,为第一类待定像素点a的局部区域内的像素点总数,为局部区域内像素点i的置信度,为第一类待定像素点a和局部区域内像素点i之间的欧氏距离,为第一类待定像素点a和局部区域内像素点i之间的置信度的差值绝对值,为第一类待定像素点a的优化置信度,为常数。
9.进一步地,所述根据第一类待定像素点的优化置信度将第一类待定像素点划分为正常像素点和异常像素点包括:将每个第一类待定像素点的优化置信度分别与第一置信度阈值和第二置信度阈值作差得到对应的第一差值和第二差值,若第一差值的绝对值较小,则将所述第一类待定像素点归为异常像素点;若第二差值的绝对值较小,则将所述第一类待定像素点归为正常像素点。
10.进一步地,所述根据每个第二类待定像素点的置信度计算所有第二类待定像素点之间的置信度分布指标包括:根据每个第二类待定像素点的置信度分别获取最大置信度、最小置信度和置信度方差,计算最大置信度和最小置信度之间的差值绝对值以及最小置信度和设定的极小值之间的相加结果,将以差值绝对值为分子、相加结果为分母的比值与置信度方差相乘得到置信度分布指标。
11.进一步地,基于归一化置信度分布指标将第二类待定像素点划分正常像素点和异常像素点包括:对置信度分布指标进行归一化处理,得到归一化置信度分布指标;当所述归一化置信度分布指标低于等于预设的置信度分布指标阈值时,分别获取正常像素点、异常像素点和第二类待定像素点的类别特征向量;将第二类待定像素点的类别特征向量和正常像素点的类别特征向量之间的余弦相似度作为第一余弦相似度,将第二类待定像素点的类别特征向量和异常像素点的类别特征向量之间的余弦相似度作为第二余弦相似度;若所述第一余弦相似度较大,则将第二类待定像素点归为正常像素点;若所述第二余弦相似度较大,则将第二类待定像素点归为异常像素点;当所述归一化置信度分布指标高于预设的置信度分布指标阈值时,使用聚类算法将第二类待定像素点分为第一类像素点和第二类像素点,分别获取正常像素点、异常像素点、第一类像素点和第二类像素点的类别特征向量;将第一类像素点的类别特征向量和正常像素点的类别特征向量之间的余弦相似度作为第三余弦相似度,将第一类像素点的类别特征向量和异常像素点的类别特征向量之间的余弦相似度作为第四余弦相似度,若所述第三余弦相似度较大,则将第一类像素点归为正常像素点,若所述第四余弦相似度较大,则将第一类像素点归为异常像素点;将第二类像素点的类别特征向量和正常像素点的类别特征向量之间的余弦相似度作为第五余弦相似度,将第一类像素点的类别特征向量和异常像素点的类别特征向量之间的余弦相似度作为第六余弦相似度,若所述第五余弦相似度较大,则将第二类像素点归为正常像素点,若所述第六余弦相似度较大,则将第二类像素点归为异常像素点;其中类别特征向量的获取方法为:将当前类别中所有像素点的置信度均值作为类别特征向量的起始点,将当前类别中像素点灰度均值代入所述标准高斯模型后所得到的置信度作为类别特征向量的终点,根据所述起始点和终点构建当前类别的类别特征向量;其中当前类别包括:正常像素点构成正常类别,异常像素点构成异常类别,第二类待定像素点构成的第二类待定类别,第一类像素点构成的第一类别和第二类像素点构成的第二类别。
12.进一步地,所述基于缺陷连通域中每个像素点的置信度计算钣金件的表面缺陷评估值包括:所述钣金件的表面缺陷评估值的计算公式为:
式中,u为钣金件的表面缺陷评估值;q为钣金件的表面图像中缺陷连通域的总数;为缺陷连通域q的缺陷程度:,为设定的极小值,为缺陷连通域q内所有异常像素点的置信度均值,若缺陷连通域q内异常像素点为第一类待定像素点,则利用优化置信度计算所述置信度均值,若缺陷连通域q内异常像素点为第二类待定像素点,则利用置信度计算所述置信度均值;为常数。
13.本发明具有如下有益效果:本发明实施例采集钣金件的表面图像得到对应的灰度图像,将灰度图像中每个像素点的灰度值代入标准高斯模型获取每个像素点对应的置信度,基于置信度对表面图像中的所有像素点进行初步划分,即分为正常像素点、异常像素点和待定像素点;为了准确确认待定像素点属于正常像素点和异常像素点,通过待定像素点的局部区域内像素点的置信度获取对应待定像素点的局部特征值,通过局部特征值对待定像素点进行划分,获得第一类待定像素点和第二类待定像素点;根据第一类待定像素点的周围像素点的置信度对第一类待定像素点的置信度进行优化,以使得第一类待定像素点划分的划分结果更加严谨;根据第二类待定像素点的置信度分析其归一化置信度分布指标以根据归一化置信度分布指标对应划分标准将第二类待定像素点进行划分,根据不同的划分标准对像素点进行层层划分,缩小待分类的像素点的范围,进而基于待分类的像素点的置信度特征,使用不同的方法处理不同类别的像素点,在对所有像素点进行划分时,使得对像素点的划分更准确,存在误差的情况更小,进而根据最终划分的正常像素点和异常像素点能够准确检测出钣金件的缺陷。
附图说明
14.为了更清楚地说明本发明实施例或现有技术中的技术方案和优点,下面将对实施例或现有技术描述中所需要使用的附图作简单的介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其它附图。
15.图1为本发明一个实施例所提供的一种建材装备钣金件缺陷检测方法的步骤流程图。
具体实施方式
16.为了更进一步阐述本发明为达成预定发明目的所采取的技术手段及功效,以下结合附图及较佳实施例,对依据本发明提出的一种建材装备钣金件缺陷检测方法,其具体实施方式、结构、特征及其功效,详细说明如下。在下述说明中,不同的“一个实施例”或“另一个实施例”指的不一定是同一实施例。此外,一或多个实施例中的特定特征、结构或特点可由任何合适形式组合。
17.除非另有定义,本文所使用的所有的技术和科学术语与属于本发明的技术领域的技术人员通常理解的含义相同。
18.下面结合附图具体的说明本发明所提供的一种建材装备钣金件缺陷检测方法的具体方案。
19.请参阅图1,其示出了本发明一个实施例提供的一种建材装备钣金件缺陷检测方法的步骤流程图,该方法包括:步骤s1:采集钣金件的表面图像,得到表面图像对应的灰度图像;构建标准高斯模型,将灰度图像中每个像素点的灰度值代入标准高斯模型中得到的结果作为表面图像中对应像素点的置信度;基于每个像素点的置信度将表面图像中的所有像素点划分为正常像素点、异常像素点和待定像素点。
20.通过相机等图像采集设备对钣金件的表面进行图像采集,相机的具体部署及视角等实施者可根据实际情况自行设置,获取钣金件的表面图像。对钣金件的表面图像进行去噪以及增强操作,以避免噪声数据对缺陷检测的影响,同时提高图像对比度,加大缺陷区域与钣金件表面正常区域的区别特征,提高后续的检测精度。
21.将经过去噪和增强处理后的钣金件的表面图像灰度化处理得到对应的灰度图像。利用正常钣金件的表面图像中像素点的灰度值进行拟合操作,得到高斯模型的期望和方差,根据期望和方差构建标准高斯模型的公式为,规定标准高斯模型的函数值的取值范围为(0,1),其中,x为像素点的灰度值,为期望,为方差,为像素点的灰度值为x所对应的高斯值,函数值越大,则像素点是钣金件正常像素点的可能性越高。
22.需要说明的是,图像的灰度化处理方法,以及高斯模型的期望和方差的获取可选用现有技术中数据拟合的方法,具体方法在此不做限定,均为本领域技术人员熟知的技术手段。
23.将钣金件的灰度图像中每个像素点的灰度值分别代入标准高斯模型,获取每个像素点所对应的高斯值,代表钣金件的表面图像中像素点k的高斯值,将每个像素点的高斯值记为其置信度,根据灰度图像中每个像素点的置信度获取表面图像的置信度分布图,以初步对缺陷区域进行识别。
24.设置第一置信度阈值,第二置信度阈值;将钣金件的表面图像中置信度低于第一置信度阈值的像素点归为表面图像的异常像素点,将置信度高于第二置信度阈值的像素点归为表面图像的正常像素点,后续不需要再对异常像素点和正常像素点进行缺陷分析,在保证检测精度的同时降低系统检测量;将置信度处于第一置信度阈值与第二置信度阈值之间的像素点归为表面图像的待定像素点,在后续步骤中对待定像素点进行进一步的缺陷分析,以实现对钣金件表面缺陷区域的准确提取。
25.优选的,本方案中第一置信度阈值取经验值0.3,第二置信度阈值取经验值0.7。
26.步骤s2:以任意一个待定像素点为中心获取对应待定像素点的局部区域,根据局部区域中每个像素点的置信度获取对应待定像素点的局部特征值,基于局部特征值将待定像素点划分为第一类待定像素点和第二类待定像素点;对每个第一类待定像素点的置信度进行优化得到优化置信度;根据第一类待定像素点的优化置信度将第一类待定像素点划分
为正常像素点和异常像素点。
27.为提高对钣金件的缺陷检测精度,对于待定像素点,本发明将对其进行详细的分析:本发明以待定像素点c为例,将待定像素点c作为中心点,获取待定像素点c的局部区域内每个像素点的置信度,将置信度相同的像素点作为一个像素点类别,获取待定像素点c的局部区域内的像素点类别数量。在本发明实施例中,局部区域的尺寸为7*7,实施者可根据实际情况自行设置。
28.获取局部区域的置信度后,基于熵的思想对局部区域内的置信度分布情况进行分析,基于熵的思想,熵值越大,则局部区域内对应待定像素点的置信度分布越杂乱,本方案利用局部特征值用于表征局部区域内的置信度分布均匀程度,局部特征值越大,局部区域的置信度越不规律,其中,局部特征值的构建公式为:式中,c为待定像素点,为待定像素点c的局部区域内的像素点类别的数量,为待定像素点c的局部区域内第n种像素点类别,为待定像素点c的局部区域内第n种像素点类别的像素点数量在局部区域中所有像素点数量的占比,为待定像素点c的局部特征值,为以常数e为底数的对数函数。
29.需要说明的是,局部特征值表示待定像素点局部区域内的置信度分布情况,待定像素点c的局部区域内第n种像素点类别的像素点数量在局部区域中所有像素点数量的占比越大即越大,表示局部区域内的置信度分布越均匀,则局部特征值越小,表示待定像素点c的局部区域内的置信度越规律;反之,则认为待定像素点c的局部区域内的置信度杂乱。
30.设置局部特征值阈值,将局部特征值低于等于局部特征值阈值对应的待定像素点归为表面图像的第一类待定像素点,第一类待定像素点的局部区域分布较为规律;将局部特征值高于局部特征值阈值对应的待定像素点归为表面图像的第二类待定像素点,第二类待定像素点的局部区域分布较为杂乱。
31.优选的,本方案中局部特征值阈值取经验值0.4。
32.通过置信度优化模型对第一类待定像素点的置信度进行优化,以获取第一类待定像素点的优化置信度。其中,置信度优化模型的公式为:式中,a为第一类待定像素点,为第一类待定像素点a的局部区域内的像素点总数,为局部区域内像素点i的置信度,为第一类待定像素点a和局部区域内像素点i之间的欧氏距离,为第一类待定像素点a和局部区域内像素点i之间的置信度的差值绝
对值,为第一类待定像素点a的优化置信度,为常数。
33.需要说明的是,第一类待定像素点a和局部区域内像素点i之间的欧氏距离越大,则最终得到的第一类待定像素点a的优化置信度越小;反之,则最终得到的第一类待定像素点a的优化置信度越大。若第一类待定像素点a和局部区域内像素点i之间的置信度的差值绝对值越大,则最终得到的第一类待定像素点a的优化置信度越小;反之,则最终得到的第一类待定像素点a的优化置信度越大。若局部区域内像素点i的置信度越大,则最终得到的第一类待定像素点a的优化置信度越大;反之,则最终得到的第一类待定像素点a的优化置信度越小。并使通过置信度优化模型后得到的优化置信度处于(0,1),以达到归一化作用。
34.由于第一类待定像素点的局部区域分布较为规律,因此基于第一类待定像素点的优化置信度对第一类待定像素点进行分析,并将第一类待定像素点进行归类:以第一类待定像素点b为例,对于第一类待定像素点b的优化置信度与第一置信度阈值的第一差值绝对值和第一类待定像素点b的优化置信度与第二置信度阈值的第二差值绝对值;若第一差值绝对值较小,则将第一类待定像素点归为异常像素点;若第二差值绝对值较小,则将第一类待定像素点归为正常像素点。
35.步骤s3:根据每个第二类待定像素点的置信度计算所有第二类待定像素点之间的置信度分布指标;基于归一化置信度分布指标将第二类待定像素点划分正常像素点和异常像素点;将所述表面图像中的异常像素点的灰度值置为1,正常像素点的灰度值置为0,对灰度值为1的像素点进行连通域分析,得到至少一个缺陷连通域,基于缺陷连通域中每个像素点的置信度计算钣金件的表面缺陷评估值。
36.第二类待定像素点的局部区域分布较为杂乱,需要通过置信度分布指标对第二类待定像素点进行进一步分析,其中,置信度分布指标的公式为:式中,为所有第二类待定像素点的置信度方差,为第二类待定像素点的置信度,为第二类待定像素点中的最大置信度,第二类待定像素点中的最小置信度,为极小值,避免分母为零,实施者可自行设置,本发明将其设置为,为第二类待定像素点的置信度分布指标。
37.需要说明的是,置信度指标表示集合内像素点的相似程度,第二类待定像素点中的最大置信度与第二类待定像素点中的最小置信度的差值越小,则代表第二类待定像素点中的置信度分布越集中,置信度分布指标越小,所以第二类待定
像素点的相似度越高;其中所有第二类待定像素点的置信度方差越小,表示第二类待定像素点中的置信度离散程度越小,则可以认为当前集合内像素点的相似度较高;反之,则认为当前集合内像素点的相似度较低。
38.为便于后续分析,对置信度分布指标进行归一化处理,得到归一化置信度分布指标,保证归一化置信度分布指标处于(0,1)。
39.设置置信度分布指标阈值,若所求的第二类待定像素点的归一化置信度分布指标低于等于置信度分布指标阈值,则认为第二类待定像素点的相似度较高,将所有第二类待定像素点作为一个集合;若所求的第二类待定像素点的归一化置信度分布指标高于置信度分布指标阈值,则认为第二类待定像素点的差异较大,进而使用聚类算法对第二类待定像素点进行聚类分析,通过聚类算法将第二类待定像素点分为第一类像素点和第二类像素点。
40.优选的,本方案中置信度分布指标阈值取经验值0.5。
41.需要说明的是,本发明选取k-means聚类算法对第二类待定像素点进行聚类分析,具体方法在此不做介绍,为本领域技术人员熟知的技术手段。
42.若所求的第二类待定像素点的归一化置信度分布指标低于等于置信度分布指标阈值时,依据类别特征向量的获取方法分别获取正常像素点的类别特征向量,异常像素点的类别特征向量和第二类待定像素点的类别特征向量。计算第二类待定像素点的类别特征向量和正常像素点的类别特征向量之间的第一余弦相似度,第二类待定像素点的类别特征向量和异常像素点的类别特征向量之间的第二余弦相似度。对于第一余弦相似度和第二余弦相似度,如果第一余弦相似度较大,则将第二类待定像素点归为正常像素点;如果第二余弦相似度较大,则将第二类待定像素点归为异常像素点。
43.若所求的第二类待定像素点的归一化置信度分布指标高于置信度分布指标阈值时,首先使用k-means聚类算法将第二类待定像素点分为第一类像素点和第二类像素点,之后依据类别特征向量的获取方法分别获取正常像素点的类别特征向量,异常像素点的类别特征向量,第一类像素点的类别特征向量和第二类像素点的类别特征向量。计算第一类像素点的类别特征向量和正常像素点的类别特征向量之间的第三余弦相似度,第一类像素点的类别特征向量和异常像素点的类别特征向量之间的第四余弦相似度;对于第三余弦相似度和第四余弦相似度,如果第三余弦相似度较大,则将第一类像素点归为正常像素点;如果第四余弦相似度较大,则将第一类像素点归为异常像素点。计算第二类像素点的类别特征向量和正常像素点的类别特征向量之间的第五余弦相似度,第二类像素点的类别特征向量和异常像素点的类别特征向量之间的第六余弦相似度;对于第五余弦相似度和第六余弦相似度,如果第五余弦相似度较大,则将第二类像素点归为正常像素点;如果第六余弦相似度较大,则将第二类像素点归为异常像素点。
44.类别特征向量的获取方法为:将当前类别中所有像素点的置信度均值作为类别特
征向量的起始点,将当前类别中像素点灰度均值代入标准高斯模型后所得到的置信度作为类别特征向量的终点,根据起始点和终点构建当前类别的类别特征向量。
45.以第二类待定像素点对应的类别为例,获取其类别特征向量的过程为:根据所有第二类待定像素点的置信度计算置信度均值,将置信度均值作为类别特征向量的起始点,根据所有的第二类待定像素点的灰度值计算灰度值均值,将灰度值均值代入标准高斯模型后所得的置信度作为类别特征向量的终点,由起始点和终点构建的第二类待定像素点的类别特征向量为。
46.优选的,本方案中余弦相似度的计算方法为:以第一余弦相似度为例,第一余弦相似度是由第二类待定像素点的类别特征向量和正常像素点的类别特征向量这两个向量通过计算得到的,第一余弦相似度的计算公式为,其中为所有正常像素点的置信度均值,为将所有正常像素点的灰度值均值代入标准高斯模型后所得的置信度。
47.通过对第一类待定像素点和第二类待定像素点的分类,将表面图像中的所有像素点划分为了正常像素点和异常像素点,实现了对像素点类别的精确划分,降低了错检误检的状况。将表面图像中的异常像素点的灰度值置为1,正常像素点的灰度值置为0,可以准确得到钣金件的表面图像中灰度值为1的像素点显示的各个缺陷连通域,通过缺陷连通域进而实现对钣金件的缺陷检测分析。基于缺陷连通域中每个像素点的置信度计算钣金件的表面缺陷评估值,钣金件的表面缺陷评估值的计算公式为:式中,q为钣金件的表面图像中缺陷连通域的总数;为缺陷连通域q的缺陷程度:,为设定的极小值,本发明设置为,为缺陷连通域q内所有异常像素点的置信度均值,若缺陷连通域q内异常像素点为第一类待定像素点,则利用优化置信度计算置信度均值,若缺陷连通域q内异常像素点为第二类待定像素点,则利用置信度计算置信度均值,为常数;u为钣金件的表面缺陷评估值,若u越大,则认为钣金件的表面缺陷越严重。本发明将对表面缺陷评估值函数进行归一化处理,保证钣金件的表面缺陷评估值处于(0,1)。
48.需要说明的是,表面缺陷评估值表示钣金件的表面缺陷程度,缺陷连通域q内所有异常像素点的置信度均值越小,表示异常像素点所对应的表面图像位置有缺陷,则缺陷连通域q的缺陷程度越大,表示钣金件的表面程度越严重,使得表面缺陷评估值u越大,
则认为钣金件表面效果越不好,缺陷状况越严重;反之,则认为钣金件的表面合格。
49.设置表面缺陷评估值阈值,当所求钣金件的表面缺陷评估值高于表面缺陷评估值阈值时,则该钣金件的表面效果不佳,表面出现较为严重的缺陷问题,需要对钣金件进行再次加工,以保证钣金件的表面质量;反之,则认为该钣金件合格。
50.优选的,本方案中表面缺陷评估值阈值取经验值0.5。
51.需要说明的是:上述本发明实施例先后顺序仅仅为了描述,不代表实施例的优劣。在附图中描绘的过程不一定要求示出的特定顺序或者连续顺序才能实现期望的结果。在某些实施方式中,多任务处理和并行处理也是可以的或者可能是有利的。
52.本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。
53.以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
再多了解一些

本文用于创业者技术爱好者查询,仅供学习研究,如用于商业用途,请联系技术所有人。

发表评论 共有条评论
用户名: 密码:
验证码: 匿名发表

相关文献