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

应用于矿山环境恢复治理的数据处理方法及系统与流程

2022-05-06 10:49:11 来源:中国专利 TAG:


1.本发明涉及图像数据处理技术,具体涉及应用于矿山环境恢复治理的数据处理方法及系统。


背景技术:

2.矿山地质环境是生态环境的重要组成部分,矿产资源开发引发、产生和加剧了矿山地质环境问题,其类型、表现形式、严重程度等,与开发的矿产资源种类、开发方式、区域地质环境条件、开采规模等因素密切相关。
3.随着信息技术的发展,各种遥感数据已经开始大规模的应用在地质环境的监测、修复和治理当中。但是由于矿山地质环境的复杂性,在进行矿山地质恢复的监测过程中,计算机进行相关边缘识别时,噪声很大,季节天气等因素对遥感图像的影响也很大,所以无法大规模的推广。


技术实现要素:

4.为了至少克服现有技术中的上述不足,本技术的目的在于提供应用于矿山环境恢复治理的数据处理方法及系统。
5.第一方面,本技术实施例提供了应用于矿山环境恢复治理的数据处理方法,包括:获取目标区域矿山恢复前的遥感图像作为原始图像;从所述原始图像中识别矿区绘制矿区的闭合边界作为第一边界,并从所述原始图像中识别矿区外的至少两条山脊线作为识别标识;获取所述目标区域开始矿山恢复的遥感图像作为待识别图像;通过所述识别标识校准所述待识别图像,并通过所述第一边界识别所述待识别图像中矿区的闭合边界作为第二边界;通过所述第一边界和第二边界计算矿山恢复率。
6.现有技术中,目前进行矿山恢复情况数据分析的主要手段还是基于遥感图像进行相关分析,一般对遥感图像进行边缘识别是基于边缘两侧像素的差异进行的,例如灰度差异、色彩空间中差异、色彩通道差异等;但是在进行矿山恢复时,现有的步骤是恢复基础地貌再搭建植被网格,最后再分批次种植草本植被和木本植被,在不同的步骤情况下天气和季节等因素对遥感图像的影响是不同的。例如,在种植草本植被而没有种植木本植被的情况下,秋季的颜色差异会和春季完全不同。为了提高精度,每次遥感图像的识别都需要大量的人工干预进行降噪,这对于经常需要进行恢复程度监测的图像识别是不利的。
7.本技术实施例实施时,主要的思路是对原始图像进行较为精确的识别,再以原始图像识别的结果为约束条件进行后续图像的识别来提高精度,减少后续图像识别的人工干预。具体的,在本技术实施例中,需要先做矿区的识别,矿区一般为一个闭合区域,需要通过机器识别进行第一边界的识别,由于在矿区尚未进行恢复时,可以对遥感图像的天气和季节进行选择,从而可以获得较好的识别效果,并且也可以进行一定的人工干预降噪,获得较
为精确的第一边界。
8.而为了保证后续可以通过第一边界对后续图像进行识别,本技术实施例中还识别出了识别标识进行图像之间的校准。其中,山脊线在各种天气和季节条件下,都会比较明显;比如在有日照的时候,山脊线两侧的灰度会有明显差异;或者在光照不足时,通过山脊线两侧不同植被的色差进行识别。所以本技术实施例将至少两条山脊线作为识别标识,来提高识别的准确率。应当理解的是,这里的山脊线识别出来一般应当为多段线段的组合。
9.获取待识别图像时,该待识别图像应当是已经开始进行矿山恢复作业后的图像,以用来评估矿山恢复情况。通过识别标识将待识别图像校准至与原始图像相同或相近的角度和比例,便于进行具体的矿区边界的识别。
10.在本技术实施例中,以第一边界为约束进行第二边界的识别,通过第一边界作为参照,可以有效的进行第二边界的识别,从而过滤掉大量的噪声,提高识别的精度,再通过识别出来的数据进行矿山恢复率计算。本技术以原始图像中矿区的边界作为约束条件对第二边界进行识别,可以有效的提高第二边界识别的准确性,减少人工干预,降低了人力成本,有利于矿山恢复情况的多阶段评估。
11.在一种可能的实现方式中,通过所述识别标识校准所述待识别图像包括:将所述识别标识映射到所述待识别图像上;从所述待识别图像中识别对应所述识别标识的山脊线作为第一山脊线,并获取所述第一山脊线朝向所述矿区方向的端点作为第一端点;获取所述识别标识朝向所述矿区方向的端点作为第二端点,并将调整所述待识别图像至所有的所述第一端点和对应的所述第二端点对齐;从所述第一端点开始沿着所述第一山脊线获取至少两个山脊线的拐点作为第一拐点,并从所述第二端点开始沿着所述识别标识的山脊线获取至少两个山脊线的拐点作为第二拐点;调整所述待识别图像至所有的所述第一端点和对应的所述第二端点对齐,且所述第一拐点和对应的第二拐点对齐。
12.在一种可能的实现方式中,通过所述第一边界识别所述待识别图像中矿区的闭合边界作为第二边界包括:将所述第一边界映射到校准后的所述待识别图像上;从所述第一边界中选取扫描起点,并以所述扫描起点开始向所述第一边界围成的区域内部预设深度扫描边界数据作为待选边界数据;所述边界数据为两侧像素差异超过预设的边缘;当任意所述待选边界数据出现不连续时,抛弃该待选边界数据并继续扫描,直至存在至少一条完整连续的闭合边界。
13.当所有的所述待选边界数据均不连续时,调整所述预设深度和/或所述像素差异并再次从所述扫描起点扫描。
14.在一种可能的实现方式中,通过所述第一边界识别所述待识别图像中矿区的闭合边界作为第二边界还包括:当扫描出来待选边界数据中,完整连续的闭合边界为一个时,将所述闭合边界作为所述第二边界;
当扫描出来待选边界数据中,完整连续的闭合边界数量大于一个时,将最外圈的所述闭合边界作为所述第二边界。
15.在一种可能的实现方式中,通过所述第一边界和第二边界计算矿山恢复率包括:对所述第一边界内的区域进行积分获取所述第一边界内区域面积作为第一面积数据,并对所述第二边界内的区域进行积分获取所述第二边界内区域面积作为第二面积数据;根据所述第一面积数据和所述第二面积数据计算所述矿山恢复率。
16.第二方面,本技术实施例提供了应用于矿山环境恢复治理的数据处理系统,包括:基准获取单元,被配置为获取目标区域矿山恢复前的遥感图像作为原始图像;基准生成单元,被配置为从所述原始图像中识别矿区绘制矿区的闭合边界作为第一边界,并从所述原始图像中识别矿区外的至少两条山脊线作为识别标识;目标获取单元,被配置为获取所述目标区域开始矿山恢复的遥感图像作为待识别图像;目标识别单元,被配置为通过所述识别标识校准所述待识别图像,并通过所述第一边界识别所述待识别图像中矿区的闭合边界作为第二边界;计算单元,被配置为通过所述第一边界和第二边界计算矿山恢复率。
17.在一种可能的实现方式中,所述目标识别单元还被配置为:将所述识别标识映射到所述待识别图像上;从所述待识别图像中识别对应所述识别标识的山脊线作为第一山脊线,并获取所述第一山脊线朝向所述矿区方向的端点作为第一端点;获取所述识别标识朝向所述矿区方向的端点作为第二端点,并将调整所述待识别图像至所有的所述第一端点和对应的所述第二端点对齐;从所述第一端点开始沿着所述第一山脊线获取至少两个山脊线的拐点作为第一拐点,并从所述第二端点开始沿着所述识别标识的山脊线获取至少两个山脊线的拐点作为第二拐点;调整所述待识别图像至所有的所述第一端点和对应的所述第二端点对齐,且所述第一拐点和对应的第二拐点对齐。
18.在一种可能的实现方式中,所述目标识别单元还被配置为:将所述第一边界映射到校准后的所述待识别图像上;从所述第一边界中选取扫描起点,并以所述扫描起点开始向所述第一边界围成的区域内部预设深度扫描边界数据作为待选边界数据;所述边界数据为两侧像素差异超过预设的边缘;当任意所述待选边界数据出现不连续时,抛弃该待选边界数据并继续扫描,直至存在至少一条完整连续的闭合边界。
19.当所有的所述待选边界数据均不连续时,调整所述预设深度和/或所述像素差异并再次从所述扫描起点扫描。
20.在一种可能的实现方式中,所述目标识别单元还被配置为:当扫描出来待选边界数据中,完整连续的闭合边界为一个时,将所述闭合边界作为所述第二边界;
当扫描出来待选边界数据中,完整连续的闭合边界数量大于一个时,将最外圈的所述闭合边界作为所述第二边界。
21.在一种可能的实现方式中,所述计算还被配置为:对所述第一边界内的区域进行积分获取所述第一边界内区域面积作为第一面积数据,并对所述第二边界内的区域进行积分获取所述第二边界内区域面积作为第二面积数据;根据所述第一面积数据和所述第二面积数据计算所述矿山恢复率。
22.本发明与现有技术相比,具有如下的优点和有益效果:本发明应用于矿山环境恢复治理的数据处理方法及系统,以原始图像中矿区的边界作为约束条件对第二边界进行识别,可以有效的提高第二边界识别的准确性,减少人工干预,降低了人力成本,有利于矿山恢复情况的多阶段评估。
附图说明
23.此处所说明的附图用来提供对本发明实施例的进一步理解,构成本技术的一部分,并不构成对本发明实施例的限定。在附图中:图1为本技术实施例方法步骤示意图;图2为本技术实施例系统架构示意图;图3为本技术实施例原始图像示意图;图4为本技术实施例第一边界和识别标识示意图。
具体实施方式
24.为使本技术实施例的目的、技术方案和优点更加清楚,下面将结合本技术实施例中的附图,对本技术实施例中的技术方案进行清楚、完整地描述,应当理解,本技术中附图 仅起到说明和描述的目的,并不用于限定本技术的保护范围。另外,应当理解,示意性的附图并未按实物比例绘制。本技术中使用的流程图示出了根据本技术实施例的一些实施例实现的操作。应该理解,流程图的操作可以不按顺序实现,没有逻辑的上下文关系的步骤可以反转顺序或者同时实施。此外,本领域技术人员在本技术内容的指引下,可以向流程图添加一个或多个其它操作,也可以从流程图中移除一个或多个操作。
25.另外,所描述的实施例仅仅是本技术一部分实施例,而不是全部的实施例。通常在此处附图中描述和示出的本技术实施例的组件可以以各种不同的配置来布置和设计。因此,以下对在附图中提供的本技术的实施例的详细描述并非旨在限制要求保护的本技术的范围,而是仅仅表示本技术的选定实施例。基于本技术的实施例,本领域技术人员在没有做出创造性劳动的前提下所获得的所有其它实施例,都属于本技术保护的范围。
26.请结合参阅图1,为本发明实施例所提供的应用于矿山环境恢复治理的数据处理方法的流程示意图,所述应用于矿山环境恢复治理的数据处理方法可以应用于图2中的应用于矿山环境恢复治理的数据处理系统,进一步地,所述应用于矿山环境恢复治理的数据处理方法具体可以包括以下步骤s1-步骤s5所描述的内容。
27.s1:获取目标区域矿山恢复前的遥感图像作为原始图像;s2:从所述原始图像中识别矿区绘制矿区的闭合边界作为第一边界,并从所述原
始图像中识别矿区外的至少两条山脊线作为识别标识;s3:获取所述目标区域开始矿山恢复的遥感图像作为待识别图像;s4:通过所述识别标识校准所述待识别图像,并通过所述第一边界识别所述待识别图像中矿区的闭合边界作为第二边界;s5:通过所述第一边界和第二边界计算矿山恢复率。
28.现有技术中,目前进行矿山恢复情况数据分析的主要手段还是基于遥感图像进行相关分析,一般对遥感图像进行边缘识别是基于边缘两侧像素的差异进行的,例如灰度差异、色彩空间中差异、色彩通道差异等;但是在进行矿山恢复时,现有的步骤是恢复基础地貌再搭建植被网格,最后再分批次种植草本植被和木本植被,在不同的步骤情况下天气和季节等因素对遥感图像的影响是不同的。例如,在种植草本植被而没有种植木本植被的情况下,秋季的颜色差异会和春季完全不同。为了提高精度,每次遥感图像的识别都需要大量的人工干预进行降噪,这对于经常需要进行恢复程度监测的图像识别是不利的。
29.本技术实施例实施时,主要的思路是对原始图像进行较为精确的识别,再以原始图像识别的结果为约束条件进行后续图像的识别来提高精度,减少后续图像识别的人工干预。具体的,在本技术实施例中,需要先做矿区的识别,矿区一般为一个闭合区域,需要通过机器识别进行第一边界的识别,由于在矿区尚未进行恢复时,可以对遥感图像的天气和季节进行选择,从而可以获得较好的识别效果,并且也可以进行一定的人工干预降噪,获得较为精确的第一边界。
30.而为了保证后续可以通过第一边界对后续图像进行识别,本技术实施例中还识别出了识别标识进行图像之间的校准。其中,山脊线在各种天气和季节条件下,都会比较明显;比如在有日照的时候,山脊线两侧的灰度会有明显差异;或者在光照不足时,通过山脊线两侧不同植被的色差进行识别。所以本技术实施例将至少两条山脊线作为识别标识,来提高识别的准确率。应当理解的是,这里的山脊线识别出来一般应当为多段线段的组合。
31.获取待识别图像时,该待识别图像应当是已经开始进行矿山恢复作业后的图像,以用来评估矿山恢复情况。通过识别标识将待识别图像校准至与原始图像相同或相近的角度和比例,便于进行具体的矿区边界的识别。
32.在本技术实施例中,以第一边界为约束进行第二边界的识别,通过第一边界作为参照,可以有效的进行第二边界的识别,从而过滤掉大量的噪声,提高识别的精度,再通过识别出来的数据进行矿山恢复率计算。本技术以原始图像中矿区的边界作为约束条件对第二边界进行识别,可以有效的提高第二边界识别的准确性,减少人工干预,降低了人力成本,有利于矿山恢复情况的多阶段评估。
33.在一种可能的实现方式中,通过所述识别标识校准所述待识别图像包括:将所述识别标识映射到所述待识别图像上;从所述待识别图像中识别对应所述识别标识的山脊线作为第一山脊线,并获取所述第一山脊线朝向所述矿区方向的端点作为第一端点;获取所述识别标识朝向所述矿区方向的端点作为第二端点,并将调整所述待识别图像至所有的所述第一端点和对应的所述第二端点对齐;从所述第一端点开始沿着所述第一山脊线获取至少两个山脊线的拐点作为第一拐点,并从所述第二端点开始沿着所述识别标识的山脊线获取至少两个山脊线的拐点作为
第二拐点;调整所述待识别图像至所有的所述第一端点和对应的所述第二端点对齐,且所述第一拐点和对应的第二拐点对齐。
34.本技术实施例实施时,需要先识别具体对应的山脊线,一般来说遥感图像是卫星拍摄的图像,所以相对应的方式可以通过相对位置进行获取,这是现有技术,本技术在此不多做复述;在获取第一山脊线后,以第一端点和第二端点先执行对齐,在这个步骤中,一般的调整策略是平移和轻微的旋转即可,然后再通过拐点进行进一步的图像调整,在这个步骤中调整策略一般是轻微的缩放,调整方式一般视遥感图像源拍摄情况而定;如果采用的是航片,一般会做更多的调整。在进行调整后,待识别图像就可以第一边界相互对应的,有利于进行下一步的第二边界的识别。
35.在一种可能的实现方式中,通过所述第一边界识别所述待识别图像中矿区的闭合边界作为第二边界包括:将所述第一边界映射到校准后的所述待识别图像上;从所述第一边界中选取扫描起点,并以所述扫描起点开始向所述第一边界围成的区域内部预设深度扫描边界数据作为待选边界数据;所述边界数据为两侧像素差异超过预设的边缘;当任意所述待选边界数据出现不连续时,抛弃该待选边界数据并继续扫描,直至存在至少一条完整连续的闭合边界。
36.当所有的所述待选边界数据均不连续时,调整所述预设深度和/或所述像素差异并再次从所述扫描起点扫描。
37.本技术实施例实施时,发明人在实践中发现,一般来说在开始进行矿区修复以后,矿区的面积是会发生缩小的;利用这个特性在第一边界为基础,向第一边界围成的区域内部一定范围内进行扫描,就可以寻找到准确的第二边界,并且由于是从第一边界开始向内扫描的,符合一般的矿区修复逻辑,有效的屏蔽了矿区内部的可能存在的其他边界,如矿区道路等数据,降低了扫描运算量。
38.同时在本技术实施例中,以第一边界向内部扫描时,不可避免的也会出现一些噪声,但是这些噪声往往是不能形成闭合区间的,所以在本技术实施例中,只要边界数据出现了中断就进行抛弃,从而屏蔽噪声,最终得到闭合的第二边界。
39.在一种可能的实现方式中,通过所述第一边界识别所述待识别图像中矿区的闭合边界作为第二边界还包括:当扫描出来待选边界数据中,完整连续的闭合边界为一个时,将所述闭合边界作为所述第二边界;当扫描出来待选边界数据中,完整连续的闭合边界数量大于一个时,将最外圈的所述闭合边界作为所述第二边界。
40.本技术实施例实施时,在进行第二边界扫描时,如果只找到一个闭合的待选边界数据,那么这个数据就是第二边界;如果找到了多个闭合边界,那么就将最外侧的作为第二边界。其主要依据是发明人发现,在遥感图像中,在一些季节性明显的植被中,处于不同恢复阶段的植被之间的分界线也会被识别成第二边界,但是本技术实施例只关心已经完成恢复的部分,所以选取了最外侧的闭合边界,而如果没有扫描到相应的第二边界,则可以采用
增加扫描深度和调整边界识别参数的方式进行,即使是进行调整,相比现有技术也可以大量的降低扫描量。
41.在一种可能的实现方式中,通过所述第一边界和第二边界计算矿山恢复率包括:对所述第一边界内的区域进行积分获取所述第一边界内区域面积作为第一面积数据,并对所述第二边界内的区域进行积分获取所述第二边界内区域面积作为第二面积数据;根据所述第一面积数据和所述第二面积数据计算所述矿山恢复率。
42.本技术实施例实施时,矿山恢复率是通过已恢复的矿区面积与矿区总面积的比值计算的,其可以通过第一面积数据和第二面积数据进行计算。
43.基于上述同样的发明构思,请结合参阅图2,还提供了应用于矿山环境恢复治理的数据处理系统的功能模块框图,关于所述应用于矿山环境恢复治理的数据处理系统的详细描述如下。
44.基准获取单元,被配置为获取目标区域矿山恢复前的遥感图像作为原始图像;基准生成单元,被配置为从所述原始图像中识别矿区绘制矿区的闭合边界作为第一边界,并从所述原始图像中识别矿区外的至少两条山脊线作为识别标识;目标获取单元,被配置为获取所述目标区域开始矿山恢复的遥感图像作为待识别图像;目标识别单元,被配置为通过所述识别标识校准所述待识别图像,并通过所述第一边界识别所述待识别图像中矿区的闭合边界作为第二边界;计算单元,被配置为通过所述第一边界和第二边界计算矿山恢复率。
45.在一种可能的实现方式中,所述目标识别单元还被配置为:将所述识别标识映射到所述待识别图像上;从所述待识别图像中识别对应所述识别标识的山脊线作为第一山脊线,并获取所述第一山脊线朝向所述矿区方向的端点作为第一端点;获取所述识别标识朝向所述矿区方向的端点作为第二端点,并将调整所述待识别图像至所有的所述第一端点和对应的所述第二端点对齐;从所述第一端点开始沿着所述第一山脊线获取至少两个山脊线的拐点作为第一拐点,并从所述第二端点开始沿着所述识别标识的山脊线获取至少两个山脊线的拐点作为第二拐点;调整所述待识别图像至所有的所述第一端点和对应的所述第二端点对齐,且所述第一拐点和对应的第二拐点对齐。
46.在一种可能的实现方式中,所述目标识别单元还被配置为:将所述第一边界映射到校准后的所述待识别图像上;从所述第一边界中选取扫描起点,并以所述扫描起点开始向所述第一边界围成的区域内部预设深度扫描边界数据作为待选边界数据;所述边界数据为两侧像素差异超过预设的边缘;当任意所述待选边界数据出现不连续时,抛弃该待选边界数据并继续扫描,直至存在至少一条完整连续的闭合边界。
47.当所有的所述待选边界数据均不连续时,调整所述预设深度和/或所述像素差异
并再次从所述扫描起点扫描。
48.在一种可能的实现方式中,所述目标识别单元还被配置为:当扫描出来待选边界数据中,完整连续的闭合边界为一个时,将所述闭合边界作为所述第二边界;当扫描出来待选边界数据中,完整连续的闭合边界数量大于一个时,将最外圈的所述闭合边界作为所述第二边界。
49.在一种可能的实现方式中,所述计算还被配置为:对所述第一边界内的区域进行积分获取所述第一边界内区域面积作为第一面积数据,并对所述第二边界内的区域进行积分获取所述第二边界内区域面积作为第二面积数据;根据所述第一面积数据和所述第二面积数据计算所述矿山恢复率。
50.在一种可能的实现方式中,请参阅图3和图4,示出了第一边界和识别标识的获取示例。从图3中可以看出,一张遥感图像中,矿区的数量可能不止一个,而从图4中可以看出,对应的第一边界的数量可能不止一个,在该示例中,第一边界的数量为三个,即由虚线所围成的区域,而实线则是用于校准的识别标识,采用了山脊线识别的方式进行,从图3中可以看出,山脊线在遥感图像中有着非常明显的特点,所以识别的精度会非常高。将从图4中获取的第一边界和识别标识用于后续矿区恢复的识别,可以有效的提高矿区恢复程度的评估。
51.本领域普通技术人员可以意识到,结合本文中所公开的实施例描述的各示例的单元及算法步骤,能够以电子硬件、计算机软件或者二者的结合来实现,为了清楚地说明硬件和软件的可互换性,在上述说明中已经按照功能一般性地描述了各示例的组成及步骤。这些功能究竟以硬件还是软件方式来执行,取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能,但是这种实现不应认为超出本发明的范围。
52.在本技术所提供的几个实施例中,应该理解到,所揭露的装置和方法,可以通过其它的方式实现。例如,以上所描述的装置实施例仅仅是示意性的,例如,所述单元的划分,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式,例如多个单元或组件可以结合或者可以集成到另一个系统,或一些特征可以忽略,或不执行。另外,所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口、装置或单元的间接耦合或通信连接,也可以是电的,机械的或其它的形式连接。
53.所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显然本领域普通技术人员可以意识到,结合本文中所公开的实施例描述的各示例的单元及算法步骤,能够以电子硬件、计算机软件或者二者的结合来实现,为了清楚地说明硬件和软件的可互换性,在上述说明中已经按照功能一般性地描述了各示例的组成及步骤。这些功能究竟以硬件还是软件方式来执行,取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能,但是这种实现不应认为超出本发明的范围。
54.另外,在本发明各个实施例中的各功能单元可以集成在一个处理单元中,也可以是各个单元单独物理存在,也可以是两个或两个以上单元集成在一个单元中。上述集成的
单元既可以采用硬件的形式实现,也可以采用软件功能单元的形式实现。
55.所述集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分,或者该技术方案的全部或部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网格设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括:u盘、移动硬盘、只读存储器(rom,read-onlymemory)、随机存取存储器(ram,random access memory)、磁碟或者光盘等各种可以存储程序代码的介质。
56.以上所述的具体实施方式,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施方式而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
再多了解一些

本文用于企业家、创业者技术爱好者查询,结果仅供参考。

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

相关文献