基于多源高维特征的城市土地利用分类确定方法和系统与流程

1.本发明涉及数据处理技术领域，特别是涉及一种基于多源高维特征的城市土地利用分类确定方法和系统。


背景技术：

2.城市土地利用是人类改造地球环境以适宜自身生存发展的顶级表现。由于人口增长和城市化进程的加快，城市土地利用正发生着迅猛变化。高质量的精细化城市土地动态制图，对于把握城市空间分布格局、预测城市发展规模与变化趋势，实现城市土地资源优化配置与集约利用至关重要。
3.早期的城市土地利用分类往往以实地调查、人工判读为主，耗费人力物力且时效性差；而基于遥感影像的半自动分类方法仅利用了影像的光谱信息，忽视了遥感影像中的空间信息和地面辅助信息，只能对城市土地利用信息进行粗略提取和估计，分类效果较差。因此，城市土地利用制图在空间尺度、数据精度以及制图手段上仍存在不足，难以满足城市精细化管理和城市建设更新的需求。针对这一问题，宫鹏等人提出基本城市土地利用类型(essential urban land use categories，euluc)的分类体系，采用了机器学习、多源时空数据，实现了城市地块尺度的城市土地利用分类，但是基本制图单元内部土地利用混合现象极大制约了分类精度。


技术实现要素：

4.本发明的目的是提供一种基于多源高维特征的城市土地利用分类确定方法和系统，提高土地利用分类精度。
5.为实现上述目的，本发明提供了如下方案：
6.一种基于多源高维特征的城市土地利用分类确定方法，包括：
7.获取城市空间数据；所述城市空间数据包括：多时相的哨兵二号遥感影像、poi数据、osm路网数据和不透水面数据；
8.对所述城市空间数据进行预处理，得到预处理数据；所述预处理数据包括：预处理后的多时相的哨兵二号遥感影像、预处理后的poi数据、预处理后的osm路网数据和预处理后的不透水面数据；
9.提取所述预处理数据中的多时相高维特征，构成多时相高维特征图层；
10.基于所述预处理后的osm路网数据对所述不透水面进行分割得到基本制图单元；
11.基于所述基本制图单元和所述多时相高维特征图层构建样本地块；
12.获取所述样本地块的类型标签，并结合所述样本地块与其对应的所述类型标签生成原始样本集；
13.基于所述原始样本集生成扩充样本集；
14.采用随机森林模型，以所述扩充样本集为输入得到扩充样本集的多样点分类结果；
15.对所述扩充样本集的多样点分类结果进行“投票”选择，获得改进后的城市土地利用初步分类结果；对所述改进后的城市土地利用初步分类结果映射处理后得到城市土地利用分类图。
16.优选地，所述对所述城市空间数据进行预处理，得到预处理数据，具体包括：
17.对所述多时相的哨兵二号遥感影像进行第一预处理得到预处理后的多时相的哨兵二号遥感影像；所述第一预处理包括：拼接、配准和裁剪；
18.对所述poi数据进行第二预处理得到预处理后的poi数据；所述第二预处理包括：清洗、去重和重分类；
19.对osm路网数据进行第三预处理得到预处理后的osm路网数据；所述第三预处理包括重分类、道路采样和建立道路缓冲区。
20.优选地，所述提取所述预处理数据中的多时相高维特征，构成多时相高维特征图层，具体包括：
21.对所述预处理后的多时相的哨兵二号遥感影像进行特征提取，得到多时相的纹理特征和遥感指数；
22.对所述预处理后的poi数据进行核密度分析和归一化处理得到核密度值；
23.对所述纹理特征、所述遥感指数和所述核密度值进行特征融合，得到所述多时相高维特征。
24.优选地，所述基于所述原始样本集生成扩充样本集，具体包括：
25.将所述原始样本集分隔为多个样本子集；
26.以在所述基本制图单元中布设多样点的方式对每一所述样本子集进行扩充；
27.基于扩充后的样本子集生成所述扩充样本集。
28.根据本发明提供的具体实施例，本发明公开了以下技术效果：
29.本发明提供的基于多源高维特征的城市土地利用分类确定方法，在对获取的城市空间数据进行预处理后，提取预处理数据中的多时相高维特征，构成多时相高维特征图层，基于预处理后的osm路网数据对不透水面进行分割得到基本制图单元，然后，基于基本制图单元和多时相高维特征图层构建样本地块后，获取样本地块的类型标签，并结合样本地块与其对应的类型标签生成原始样本集，接着基于原始样本集生成扩充样本集，最后，采用随机森林模型以扩充样本集为输入得到的城市土地利用分类结果，对城市土地利用分类结果进行投票选择和映射处理后得到城市土地利用分类图，进而更好地解决了土地利用混杂所造成的影响，达到提高土地利用分类精度的目的。
30.对应于上述提供的基于多源高维特征的城市土地利用分类确定方法，本发明还对应提供了以下实施系统，具体为：
31.一种基于多源高维特征的城市土地利用分类确定系统，包括：
32.数据获取模块，用于获取城市空间数据；所述城市空间数据包括：多时相的哨兵二号遥感影像、poi数据、osm路网数据和不透水面数据；
33.预处理模块，用于对所述城市空间数据进行预处理，得到预处理数据；所述预处理数据包括：预处理后的多时相的哨兵二号遥感影像、预处理后的poi数据、预处理后的osm路网数据和预处理后的不透水面数据；
34.特征提取模块，用于提取所述预处理数据中的多时相高维特征，构成多时相高维
特征图层；
35.单元分割模块，用于基于所述预处理后的osm路网数据对所述不透水面进行分割得到基本制图单元；
36.单元构建模块，用于基于所述基本制图单元和所述多时相高维特征图层构建样本地块；
37.样本生成模块，用于获取所述样本地块的类型标签，并结合所述样本地块与其对应的所述类型标签生成原始样本集；
38.样本扩充模块，用于基于所述原始样本集生成扩充样本集；
39.分类结果确定模块，用于采用随机森林模型，以所述扩充样本集为输入得到扩充样本集的多样点分类结果；
40.选择
‑
映射模块，用于对所述扩充样本集的多样点分类结果进行“投票”选择，获得改进后的城市土地利用初步分类结果，还用于对所述改进后的城市土地利用初步分类结果映射处理后得到城市土地利用分类图。
41.优选地，所述预处理模块具体包括：
42.第一预处理单元，用于对所述多时相的哨兵二号遥感影像进行第一预处理得到预处理后的多时相的哨兵二号遥感影像；所述第一预处理包括：拼接、配准和裁剪；
43.第二预处理单元，用于对所述poi数据进行第二预处理得到预处理后的poi数据；所述第二预处理包括：清洗、去重和重分类；
44.第三预处理单元，用于对osm路网数据进行第三预处理得到预处理后的osm路网数据；所述第三预处理包括重分类、道路采样和建立道路缓冲区。
45.优选地，所述特征提取模块具体包括：
46.特征提取单元，用于对所述预处理后的多时相的哨兵二号遥感影像进行特征提取，得到多时相的纹理特征和遥感指数；
47.poi处理单元，用于对所述预处理后的poi数据进行核密度分析和归一化处理得到核密度值；
48.融合单元，用于对所述纹理特征、所述遥感指数和所述核密度值进行特征融合，得到所述多时相高维特征。
49.优选地，所述样本扩充模块具体包括：
50.分割单元，用于将所述原始样本集分隔为多个样本子集；
51.扩充单元，用于以在所述基本制图单元中布设多样点的方式对每一所述样本子集进行扩充；
52.扩充样本生成单元，用于基于扩充后的样本子集生成所述扩充样本集。
53.因本发明提供的基于多源高维特征的城市土地利用分类确定系统达到的技术效果与上述提供的城市土地利用分类确定方法达到的技术效果相同，故在此不再进行赘述。
附图说明
54.为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图
获得其他的附图。
55.图1为本发明提供的基于多源高维特征的城市土地利用分类确定方法的流程图；
56.图2为本发明实施例提供的基于多源高维特征的城市土地利用分类确定方法实施架构图；
57.图3为本发明实施例提供的基于多源高维特征的城市土地利用分类确定方法实施架构简图；
58.图4为本发明实施例提供的2018年郑州市归一化植被指数分布图；
59.图5为本发明实施例提供的2018年郑州市归一化建筑物指数分布图；
60.图6为本发明实施例提供的各要素核密度与归一化处理结果图；
61.图7为本发明实施例提供的基本制图单元示意图；
62.图8为本发明实施例提供的样本地块单元示意图；
63.图9为本发明实施例提供的郑州市城市土地利用一级分类图；
64.图10为本发明实施例提供的郑州市城市土地利用二级分类图；
65.图11为本发明提供的基于多源高维特征的城市土地利用分类确定系统的结构示意图。
具体实施方式
66.下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
67.本发明的目的是提供一种基于多源高维特征的城市土地利用分类确定方法和系统，提高土地利用分类精度。
68.为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。
69.如图1和2所示，本发明提供的基于多源高维特征的城市土地利用分类确定方法，包括：
70.步骤100：获取城市空间数据。城市空间数据包括：多时相的哨兵二号遥感影像、poi数据、osm路网数据和不透水面数据。
71.步骤101：对城市空间数据进行预处理，得到预处理数据。预处理数据包括：预处理后的多时相的哨兵二号遥感影像、预处理后的poi数据、预处理后的osm路网数据和预处理后的不透水面数据。该步骤具体实施过程为：
72.对多时相的哨兵二号遥感影像进行第一预处理得到预处理后的多时相的哨兵二号遥感影像。第一预处理包括：拼接、配准和裁剪。
73.对poi数据进行第二预处理得到预处理后的poi数据。第二预处理包括：清洗、去重和重分类。
74.对osm路网数据进行第三预处理得到预处理后的osm路网数据。第三预处理包括重分类、道路采样和建立道路缓冲区。
75.步骤102：提取预处理数据中的多时相高维特征，构成多时相高维特征图层。可依
据多时相高维特征生成多时相高维特征库。该步骤具体实施过程为：
76.对预处理后的多时相的哨兵二号遥感影像进行特征提取，得到多时相的纹理特征和遥感指数。遥感指数包括归一化植被指数和归一化建筑指数。
77.对预处理后的poi数据进行核密度分析和归一化处理得到核密度值。
78.对纹理特征、遥感指数和核密度值进行特征融合，得到多时相高维特征。具体的，本发明基于arcgis10.2、envi5.3软件分别计算了遥感影像的纹理特征、遥感指数与poi数据的核密度值，建立高维特征库。
79.在上述步骤102中，纹理特征能够包含物质表面形态信息及其与周围环境的关系，反映影像局部模式重复、排列的规则，是影像的重要属性之一。归一化后的植被指数ndvi能反映植被的空间分布情况，可以对识别城市内部公园绿地等类别有较大帮助；归一化后的建筑指数ndbi用来反映建筑用地的信息，对提取城市建筑用地有较大帮助。而poi数据直观的反映了人类活动的轨迹和范围，可以有效识别各个地类具体分布的空间布局。
80.步骤103：基于预处理后的osm路网数据对不透水面进行分割得到基本制图单元。
81.步骤104：基于基本制图单元和多时相高维特征图层构建样本地块。具体的，在基本制图单元内部随机布设随机点，并采用实地调查、google earth平台和数据比对等方式获取随机点所在地块内部的地物混合情况，再通过高维特征库获取特征值，构成样本地块。若地块内部土地利用类别较为单一，则作为样本地块，反之，则不用。
82.步骤105：获取样本地块的类型标签，并结合样本地块与其对应的类型标签生成原始样本集。具体的，采用实地调查、数据比对等方式获取样本地块的类型标签(即土地利用类型)；结合样本地块与其对应的样本标签，生成原始样本集；
83.步骤106：基于原始样本集生成扩充样本集。该步骤的具体实施过程为：
84.将原始样本集分隔为多个样本子集。以在基本制图单元中布设多样点的方式对每一样本子集进行扩充。
85.基于扩充后的样本子集生成扩充样本集。
86.具体的，将原始样本集随机分为不同比例的实验组，一个实验组即为一个样本子集。通过在基本制图单元内部的布设多样点来扩充样本集。这一过程的实施框架如图3所示。
87.步骤107：采用随机森林模型，以扩充样本集为输入得到扩充样本集的多样点分类结果。其中，随机森林模型是基于rstudio软件，采用r编程语言搭建形成的。具体搭建过程为：首先，引入random forest程序包；其次，加载训练数据与验证数据的csv文件，通过random forest()函数对初始模型进行训练、predict()函数预测验证数据的分类结果；最后，将预测结果写入csv文件中，形成搭建好的随机森林模型。
88.步骤108：对扩充样本集的多样点分类结果进行“投票”选择，获得改进后的城市土地利用初步分类结果。对改进后的城市土地利用初步分类结果映射处理后得到城市土地利用分类图。具体的，通过arcgis10.2软件中“连接”工具将分类结果映射至地理单元，在制图软件中输出并保存制图得到的城市土地利用分类图。
89.下面以采用本发明上述提供的基于多源高维特征的城市土地利用分类确定方法确定郑州市土地分类为例，对本发明所提供技术方案的优点进行说明。在实际应用过程中，本发明可应用于任意城市的土地分类确定。
90.郑州市位于河南中部偏北地区，位于东经112
°
42'
‑
114
°
13'，北纬34
°
16'
‑
34
°
58'。北临黄河，西依嵩山，东南为广阔的黄淮平原。是中部地区中心城市以及国家重要的枢纽城市。
91.获取的郑州市的土地数据包括：2018年哨兵二号遥感影像数据、openstreetmap数据、兴趣点(poi)数据和不透水面数据。
92.其中，2018年哨兵二号遥感影像数据基于usgs平台下载，影像的获取时间分别为2月、6月、9月、11月，云量小于5％，如表1所示：
93.表1遥感影像来源表
[0094][0095]
鉴于模型运算效率和软硬件设施限制，选取波段2、波段3、波段4、波段5、波段6、波段7、波段8和波段8a，并进行重采样，将空间分辨率统一至10m。如表2所示：
[0096]
表2 sentinel
‑
2参数表
[0097][0098][0099]
其中，openstreetmap(网址：http://www.openstreetmap.org)是一种由众多志愿者采集、加工生产的开源地理数据，用户都可以在其平台上上传和编辑数据。具有更新快、成本低、免费获取等优点，在划分城市地块、城市街区等领域得到广泛应用。由于osm平台是由众多志愿者自发采集、编辑、上传数据，数据之间存在重复区域。因此，对osm数据进行拓扑检查、修改或删除错误要素以及重分类等预处理，将道路分为高速公路、城市快速路、一级主干道、二级次要道路、三级次要道路、小路以及特殊类型道路等七种类别。其次，对不同等级道路进行采样，并建立缓冲区，生成道路网，道路采样距离见表3。
[0100]
表3不同等级道路采样距离表
[0101][0102]
其中，poi数据基于百度电子地图爬取和bigemap平台下载两种获取方式，共有72061条poi数据。该数据具有更新速度快、数据量大、分类体系较细、容易获取等优势。为满足研究需求，需对poi数据进行清洗、去重、属性检查以及重分类等处理，如表4所示：
[0103]
表4 poi重分类结果表
[0104][0105]
具体的，首先，纹理包含了物质表面形态信息及其与周围环境的关系，反映了影像局部模式重复、排列的规则，是影像的重要属性之一。其次，ndvi指数能反映植被的空间分布情况，可以对识别城市内部公园绿地等类别有较大帮助；归一化建筑指数ndbi用来反映建筑用地的信息，对提取城市建筑用地有较大帮助。最后，poi数据直观的反映了人类活动的轨迹和范围，可以有效识别各个地类具体分布的空间布局。因此，选取特征有：纹理特征、遥感指数、poi核密度和归一化处理结果。
[0106]
其中，纹理特征基于灰度共生矩阵(gray
‑
level co
‑
occurrence matrix,glcm)进行提取，包括：波段的均值、对比度、方差、协同度、相关性、相异性、熵和二阶矩。
[0107]
遥感指数选取了归一化植被指数(ndvi)和归一化建筑指数(ndbi)。归一化植被指数用来监测和反映植被覆盖状况，通过不同波段的计算而得到：
[0108]
ndvi＝(nir
‑
r)/(nir r)，r指的是红光波段的波段反射率，nir是近红外波段的波段反射率；
[0109]
归一化建筑指数反映建筑用地信息和城市化发展强度，同样基于不同波段间的计算而得到：
[0110]
ndbi＝(swir
‑
nir)/(swir nir)，nir是近红外波段的波段反射率，swir指中红外波段的波段反射率。
[0111]
通过arcgis10.2软件的“地图代数”工具分别计算遥感影像的ndvi与ndbi值，如图4和图5所示。
[0112]
poi的核密度值和归一化处理结果。通过arcgis10.2对预处理后的土地数据进行核密度分析，同时为了将不同类别的poi数据统一至相同量纲下，对poi数据进行了归一化处理。通过arcgis10.2软件的“核密度分析”与“地图代数”工具分别计算各类别poi数据的核密度值，如图6中(a)
‑
(e)所示。
[0113]
因此，融合上述提取的多时相遥感影像的纹理特征、归一化植被指数、归一化建筑指数和核密度值，得到了2个光谱特征、256个纹理特征以及5个poi辅助特征作为输入模型特征值，构建高维特征库，具体见表5：
[0114]
表5
[0115][0116]
通过预处理后的osm路网数据分割郑州市不透水面，并对分割后面积较小且破碎的地块进行整合，降低基本制图单元的破碎程度，最终获得12466个基本制图单元，如图7所示。
[0117]
采用随机选取的方式来确定样本地块：通过获得的基本制图单元内部布设随机点，采用实地调查、数据比对以及google earth平台判别地块内部混合度，筛选地物类型较为单一的样本地块，最终得到337个样本地块单元，如图8所示。
[0118]
基于rstudio软件搭建随机森林模型。首先，加载“randomforest”程序包，调用随机森林模型，读取训练集与验证集。其次，构建参数最佳组合。先确定mtry值，将ntree默认为500，mtry范围为2到20，步长为1，遍历模型。观察oob值，可以发现：当mtry为16时，袋外误差值最小。其次，当确定mtry的最佳值为16后，设置ntree范围为100
‑
2000，步长为100，可以发现：当ntree≥600时，模型的误差值波动较小，并逐渐趋于稳定状态。为了防止实验结果过拟合，本发明采用ntree值为600作为最优值。由此得到模型参数的最佳组合。最后，对数据集进行训练预测，记录各组试验结果。
[0119]
对获取的样本地块进行特征采样，并通过google earth、实地调查等确定土地利用类型，生成原始样本集；
[0120]
将生成的原始样本集分为不同比例的实验组，各比例的样本数量如表6所示。
[0121]
表6“原始样本集”实验组样本数量构成表
[0122][0123][0124]
首先将原始样本集输入搭建好的随机森林模型中，获得城市土地利用分类结果，如表7所示。
[0125]
表7“原始样本集”各组分类精度表
[0126][0127]
结果表明：当训练样本数量增多时，一级分类精度提升，二级分类精度在训练样本占比达到70％以后降低，实验结果可能与二级地类各类别验证样本数量较少、样本内部各类别分布差异较大有关。
[0128]
然后，通过对基本制图单元的多点特征采样的方式在样本地块内部布设多个样本点。经过多次试验，本发明选取3n(n＝10，为分类的类别数)作为布设的样本点个数，实现样
本集的扩充，并代入各实验组中进行训练预测，得到一、二级地物的分类精度，如表8所示。
[0129]
表8“扩充样本集”各组分类精度表
[0130][0131][0132]
结果表明：对基本制图单元内部多个样点的预测结果进行总体精度统计，发现一级分类精度与二级分类精度降低。这是由于多个样点的预测结果平均化造成精度下降，但是在其个别地类分类精度有较大提升。
[0133]
最后，将基于“扩充样本集”的多样点分类预测结果进行基于基本制图单元的“投票”机制选择，即将每一单元中被预测次数最多的土地利用类型作为该单元最终的土地利用类型。代入各实验组中，得到不同比例下一级、二级地物的分类精度，如表9所示。
[0134]
表9基于“投票”机制选择的各组分类精度表
[0135][0136]
可见，将多时相高维特征、基本制图单元的多样点扩充、基本制图单元的多样点预测“投票”机制结合使用，可以有效提高地物总体分类精度，得到了郑州市城市土地利用的精细制图结果，如图9和图10所示。
[0137]
综上，本发明上述提供的技术方案相对于现有技术具有以下优点：
[0138]
1、本发明通过结合多时相高维特征、基本制图单元的多样点扩充和多样点预测“投票”机制的方法，有效提升了城市土地利用分类精度；
[0139]
2、本发明将“投票”选择机制与多时相高维特征信息相结合，较好的降低了基本制图单元内部土地混杂的影响，为城市土地利用分类提供了新的研究思路。
[0140]
此外，对应于上述提供的基于多源高维特征的城市土地利用分类确定方法，本发明还对应提供了一种基于多源高维特征的城市土地利用分类确定系统。如图11所示，该城市土地利用分类确定系统包括：数据获取模块1、预处理模块2、特征提取模块3、单元分割模块4、单元构建模块5、样本生成模块6、样本扩充模块7、分类结果确定模块8和选择
‑
映射模块9。
[0141]
数据获取模块1用于获取城市空间数据。城市空间数据包括：多时相的哨兵二号遥感影像、poi数据、osm路网数据和不透水面数据。
[0142]
预处理模块2用于对城市空间数据进行预处理，得到预处理数据。预处理数据包括：预处理后的多时相的哨兵二号遥感影像、预处理后的poi数据、预处理后的osm路网数据和预处理后的不透水面数据。
[0143]
特征提取模块3用于提取预处理数据中的多时相高维特征，构成多时相高维特征图层。
[0144]
单元分割模块4用于基于预处理后的osm路网数据对不透水面进行分割得到基本制图单元。
[0145]
单元构建模块5用于基于基本制图单元和多时相高维特征图层构建样本地块。
[0146]
样本生成模块6用于获取样本地块的类型标签，并结合样本地块与其对应的类型标签生成原始样本集。
[0147]
样本扩充模块7用于基于原始样本集生成扩充样本集。
[0148]
分类结果确定模块8用于采用随机森林模型，以扩充样本集为输入得到扩充样本集的多样点分类结果。
[0149]
选择
‑
映射模块9用于对扩充样本集的多样点分类结果进行“投票”选择，获得改进后的城市土地利用初步分类结果，还用于对改进后的城市土地利用初步分类结果映射处理后得到城市土地利用分类图。
[0150]
其中，上述预处理模块2具体包括：第一预处理单元、第二预处理单元和第三预处理单元。
[0151]
第一预处理单元用于对多时相的哨兵二号遥感影像进行第一预处理得到预处理后的多时相的哨兵二号遥感影像。第一预处理包括：拼接、配准和裁剪。
[0152]
第二预处理单元用于对poi数据进行第二预处理得到预处理后的poi数据。第二预处理包括：清洗、去重和重分类。
[0153]
第三预处理单元用于对osm路网数据进行第三预处理得到预处理后的osm路网数据。第三预处理包括重分类、道路采样和建立道路缓冲区。
[0154]
上述特征提取模块3具体包括：特征提取单元、poi处理单元和融合单元。
[0155]
特征提取单元用于对预处理后的多时相的哨兵二号遥感影像进行特征提取，得到多时相的纹理特征和遥感指数。
[0156]
poi处理单元用于对预处理后的poi数据进行核密度分析和归一化处理得到核密度值。
[0157]
融合单元用于对纹理特征、遥感指数和核密度值进行特征融合，得到多时相高维特征。
[0158]
上述样本扩充模块7具体包括：分割单元、扩充单元和扩充样本生成单元。
[0159]
分割单元用于将原始样本集分隔为多个样本子集。
[0160]
扩充单元用于以在基本制图单元中布设多样点的方式对每一样本子集进行扩充。
[0161]
扩充样本生成单元用于基于扩充后的样本子集生成扩充样本集。
[0162]
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。
[0163]
本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不应理解为对本发明的限制。