巴基斯坦土地覆盖数据产品季节性精度评价
1.
2.
3.
Seasonal Accuracy Assessments of Three Land Cover Datasets in Pakistan
1.
2.
3.
通讯作者:
收稿日期: 2019-03-25 修回日期: 2020-05-21 网络出版日期: 2020-07-03
基金资助: |
|
Received: 2019-03-25 Revised: 2020-05-21 Online: 2020-07-03
作者简介 About authors
郭紫燕(1995-),女,江苏连云港人,博士研究生,主要从事土地利用/土地覆盖遥感监测研究E⁃mail:
关键词:
Keywords:
本文引用格式
郭紫燕, 杨康, 刘畅, 程亮, 李满春.
Guo Ziyan, Yang Kang, Liu Chang, Cheng Liang, Li Manchun.
1 引 言
全球高精度土地覆盖数据产品对指导资源开发,预警生态环境破坏,应对气候变化以及监测地理国情均具有重要意义 [1,2,3,4,5]。目前,随着数据处理技术的发展,已有越来越多的全球高精度土地覆盖数据产品可供免费下载 [6],并用于各方面的研究 [7,8,9,10,11,12,13,14]。这些产品间的全球或区域精度比较作为产品选择与使用的前提,得到越来越多的重视 [15,16,17]。然而,目前这些产品均为年际产品,仅采用了单时相遥感影像用于制图,难以反映年内土地覆盖类型的动态变化,对描述不同季节土地覆盖情况的适应性有待商榷。近年来已有研究开始关注土地覆盖的季节性差异 [18,19],但是土地覆盖年际产品在较长时间内仍将被广泛使用。因此,研究分析评价了几套全球高精度土地覆盖数据年际产品的季节性精度。
选取巴基斯坦作为研究区,评价3套全球高精度土地覆盖数据产品在巴基斯坦夏冬两季的精度差异。巴基斯坦总面积796 095 km2,南北气候差异显著,印度河自北向南贯穿全境[20],境内土地覆盖类型齐全,含耕地、森林、草地、灌木、水体、人造地表、冰川积雪和裸地等多个类型,土地覆盖分类精度评价更为全面客观。同时,巴基斯坦部分土地覆盖类型呈现显著的季节性转换特征,例如受耕作周期影响的耕地和受气候条件影响的冰川等,这对研究季节性土地覆盖十分有益。此外,“中巴经济走廊”是“一带一路”旗舰产品 [21],分析巴基斯坦土地覆盖类型及其动态变化,对于监测巴基斯坦地理国情、掌握“中巴经济走廊”态势具有重要意义,是高精度土地覆盖数据产品应用的重要领域。
2 数据与方法
2.1 研究数据
研究选择GlobCover-2009(300 m)、FROM-GLC-2010(30 m)和GlobeLand30-2010(30 m)3套全球高精度土地覆盖数据产品进行精度评价,其分类体系见表 1。
表 1 3套土地覆盖数据产品分类体系及统一分类体系
Tab.1
产品 | GlobCover-2009 | FROM-GLC-2010 | GlobeLand30-2010 | 统一分类 |
---|---|---|---|---|
类型 | 水田 雨林农田 镶嵌型农田(50%~70%)/植被(20%~50%) | 农田 | 耕地 | 耕地 |
镶嵌型植被(草原/灌木林/森林) (50%~70%)/农田(20%~50%) 开阔阔叶常绿阔叶林或半落叶林(>15%)(>5 m) 阔叶落叶林(>40%)(>5 m) 封闭型针刺常绿林(>40%)(>5 m) 开阔混合阔叶林和针刺林(>15%)(>5 m) 镶嵌型森林或灌木丛(50%~70%)/草地(20%~50%) 镶嵌型草原(50%~70%)/森林或灌木林(20%~50%) (阔叶或针刺,常绿或落叶)灌木丛(>15%)(<5 m) 开放型草本植被(草原,热带稀树草或地衣/苔藓)(>15%) 稀疏植被(<15%) 定期淹没或淹水的土壤封闭开放草原或木质植被 新鲜,咸水或盐水(>15%) | 森林 草地 灌木 | 森林 草地 灌木地 苔原 | 植被 | |
人造表面及相关区域(城市面积>50%) | 不透水面 | 人造地表 | 人造地表 | |
裸地 | 裸地 | 裸地 | 裸地 | |
水体 封闭阔叶林或灌木丛永久性淹水 - 盐水或微咸水(>40%) | 水体 | 水体 湿地 | 水体 | |
永久冰雪 | 积雪冰川 | 冰川和永久积雪 | 冰川积雪 |
GlobCover-2009(Global Land Cover Map for 2009)数据产品由欧洲空间局(European Space Agency,ESA)、联合国环境规划署(United Nations Environment Programme,UNEP)等组织共同参与完成。GlobCover-2009产品以2009年ENVISAT卫星的MERIS(MEdium Resolution Imaging Spectrometer)影像作为数据源,空间分辨率为300 m,使用了LCCS (Land Cover Classification System)分类体系,采用分层分区的分类方法进行制图,共包含22种土地覆盖类型,全球数据总精度约为73% [22]。
FROM-GLC(Finer Resolution Observation and Monitoring of Global Land Cover)-2010数据产品由清华大学生产,以Landsat TM /ETM + 影像为数据源,使用最大似然分类器(Maximum Likelihood Classifier,MLC)、随机森林(Random Forest,RF)、支持向量机(Support Vector Machine,SVM)等分类器,通过监督分类的方法制图,是首套30 m空间分辨率的全球土地覆盖数据产品,包含8种土地覆盖类型,全球数据总精度约为71.5% [23]。
2.2 研究方法
2.2.1 统一土地覆盖类型
本研究评价的3套高精度土地覆盖数据产品分类体系和类型存在差异,不能直接进行精度比较,因此研究首先统一不同数据产品的土地覆盖类型,归并对应的土地覆盖类型[26,27,28,29],得到耕地、裸地、植被、人造地表、水体和冰川积雪6个主要土地覆盖类型(表1)。GlobCover-2009数据产品中植被类型详细,其中“镶嵌型植被”等类型无法具体划分森林、草地或灌木,因此,在巴基斯坦植被面积较小,植被内部发生季节间转换的可能性小的前提下,研究将各产品中的森林、草地、灌木和苔原均合并为“植被”类型。此外,FROM-GLC-2010数据产品将湿地划分为水体,为了保持一致,本文将GlobCover-2009和GlobeLand30-2010中的永久性淹水和湿地划分为水体。
2.2.2 土地覆盖分类精度评价
研究共选取2009~2011年30 m Landsat-5、Landsat-7多光谱遥感影像共122景,其中夏季、冬季影像各61景,均覆盖巴基斯坦全境。以样本点为参考数据对土地覆盖数据产品进行精度评价是目前最常用的方法 [30]。研究选择了997个样本点(1 000个随机点中去掉3个空值点),利用夏冬两季的遥感影像分别对验证样本点进行目视解译,确定夏冬两季的土地覆盖类型属性,形成两套独立的验证点数据集。此外,为了排除验证样本点个数对精度评价结果的影响,本研究对样本点进行9次分层抽样,以100为间隔,抽取点数为100~900共9组验证样本点,对3套产品、2个季节重复进行精度评价。
研究建立了3套土地覆盖数据与夏冬两季验证点间的混淆矩阵,计算了3套数据在2个季节的总精度(Overall Accuracy,OA,产品与真值类型相同的点数/样本点数)与Kappa系数,以反映3套数据在不同季节用来描述巴基斯坦土地覆盖真实状况的适应程度。除了总精度,还分别计算了各种土地覆盖类型在2个季节的生产者精度(Produce Accuracy,PA,产品与真值均为某类型的点数/真值中为某类型的点数)及用户精度(User Accuracy,UA,产品与真值均为某类型的点数/产品中为某类型的点数),以分析各类型间的混淆情况[31]。以上4个精度评价指标中,OA和Kappa系数描述待评价数据的总体情况,UA和PA反映各种土地覆盖类型的精度及混淆状况。由UA和PA的概念可知,UA更适合比较同一数据产品不同季节间的差异,而PA更适合比较同一季节不同数据产品间的差异。
3 研究结果与讨论
3.1 巴基斯坦土地覆盖特征
研究经过数据预处理和统一分类,得到3套巴基斯坦土地覆盖数据产品,获得了巴基斯坦土地覆盖基本情况(图 1)。虽然3套产品对巴基斯坦土地覆盖的描述存在差别,但整体上各地类的比例较为一致,这与戴昭鑫等在南美地区得出的结论具有相似性 [32]。巴基斯坦土地覆盖以裸地为主,大约占总面积的40%,显著高于全球水平(GlobCover-2009、FROM-GLC-2010和GlobeLand30-2010依次为50%、30%、44%,全球裸地覆盖率约为16.5%),主要分布在巴基斯坦西南部。裸地、植被(约25%)和耕地(约28%)之和约占巴基斯坦总面积的90%,空间上自西向东的优势土地覆盖类型依次为裸地、植被、耕地。水体大部分集中在印度河和沿海湿地,低于全球水平(约1% vs. 3. 6%),人造地表面积较小,低于1%,低于全球水平(GlobCover-2009、FROM-GLC-2010和GlobeLand30-2010依次为0.32%、0.34%、0.71%,全球人造地表覆盖率约为0.66%),集中分布在东部伊斯兰堡和南部沿海卡拉奇附近的城市群。
图1
研究根据2010年夏冬两季遥感影像,对巴基斯坦境内1 000个验证样本点进行目视解译,这些样本点在夏冬两季的属性在反映了巴基斯坦土地覆盖的季节差异,因此研究构建冬季—夏季转移矩阵(表2)。在巴基斯坦夏冬两季,土地覆盖类型空间变化主要有耕地—裸地、植被—裸地、水体—裸地和冰川积雪—裸地几种类型。由此可知,巴基斯坦土地覆盖的季节差异主要变现为裸地与其他土地覆盖类型的转变。
表2 巴基斯坦土地覆盖样本点冬季—夏季转移矩阵
Table 2
夏冬 | 耕地 | 裸地 | 植被 | 水体 | 人造地表 | 冰川积雪 | 夏季总计 |
---|---|---|---|---|---|---|---|
耕地 | 91 | 26 | 0 | 0 | 0 | 0 | 117 |
裸地 | 3 | 437 | 9 | 0 | 0 | 13 | 462 |
植被 | 2 | 95 | 160 | 2 | 0 | 7 | 266 |
水体 | 0 | 5 | 0 | 34 | 1 | 0 | 40 |
人造地表 | 0 | 0 | 0 | 0 | 52 | 0 | 52 |
冰川积雪 | 0 | 2 | 2 | 0 | 0 | 56 | 60 |
冬季总计 | 96 | 565 | 171 | 36 | 53 | 76 | 997 |
3.2 土地覆盖数据产品分类精度评价
研究生成了3套土地覆盖数据产品GlobCover-2009、FROM-GLC-2010和GlobeLand30-2010在夏冬两季的混淆矩阵,进而评价了3套产品各自的夏冬两季土地覆盖分类精度(表 3)。
表3 巴基斯坦3套土地覆盖数据产品2季精度分析结果
Table 3
精度指标 | 季节 | GlobCover -2009 | FROM-GLC-2010 | GlobeLand 30-2010 |
---|---|---|---|---|
总精度 | 夏 | 59.1 | 61.2 | 65.6 |
冬 | 59.5 | 59.0 | 63.9 | |
Kappa系数 | 夏 | 0.41 | 0.44 | 0.52 |
冬 | 0.38 | 0.40 | 0.48 |
3套产品夏季与冬季的总精度差别不大,GlobCover-2009的总精度冬季略高于夏季,FROM-GLC-2010和GlobeLand30-2010的总精度夏季略高于冬季(表3),因此GlobCover-2009所反映的土地覆盖特征更符合巴基斯坦冬季的真实情况,FROM-GLC-2010和GlobeLand30-2010更符合夏季的真实情况。此外,3套产品中,空间分辨率为300 m的GlobCover-2009的总精度最低(约59.3%),空间分辨率为30 m的FROM-GLC-2010总精度(约60.1%)仅略高于GlobeCover2009, GlobeLand30-2010总精度最高(约64.8%)。因此,巴基斯坦土地覆盖数据产品的精度随着空间分辨率的提高有所上升(约1%~6%),但并不显著。
尽管3套产品在夏冬两季的总精度差别不大,但各个土地覆盖类型在不同数据产品及不同季节的精度差异明显。研究对比了3套产品夏冬两季的生产者精度和用户精度(图2)。
图2
图2
3套土地覆盖数据产品夏冬两季用户精度和生产者精度
Fig.2
User Accuracy (UA)and Produceser Accuracy (PA) for three land cover datasets in summer and winter
在耕地分类中,3套产品的生产者精度GlobeLand30-2010分类效果最优,其次为GlobCover-2009, FROM-GLC-2010再次之;从用户精度来看,3套产品夏季均优冬季。在裸地的分类中,3套产品的生产者精度GlobCover-2009分类效果最优,其次为GlobeLand30-2010, FROM-GLC-2010再次之;从用户精度来看,3套产品裸地分类冬季好于夏季。在植被分类中,3套产品的生产者精度FROM-GLC-2010最优,GlobeLand30-2010优于GlobCover-2009;从用户精度的角度看,夏季明显高于冬季,季节性差异较为显著。在人造地表的分类中,GlobeLand30-2010分类效果最优,其次为FROM-GLC-2010,GlobCover-2009再次之。人造地表的用户精度季节波动很小。水体的生产精度GlobeLand30-2010最高,其次为FROM-GLC-2010,GlobCover-2009再次之;从用户精度来看,3套产品夏季均好于冬季。在冰川积雪的分类中,生产者精度FROM-GLC-2010最高,其次为GlobeLand30-2010,GlobCover-2009再次之;从用户精度的角度看,3套产品冬季优于夏季。
综上,GlobCover-2009能较好地刻画裸地的覆盖情况,FROM-GLC-2010能较好地刻画植被和冰川积雪的覆盖情况,GlobeLand30-2010对耕地、人造地表和水体有更好的表达效果。3套土地覆盖数据产品中各类型对季节的适用性较为一致,耕地、裸地和冰川积雪的分布更符合冬季的真实情况,而植被、水体更符合夏季的真实情况,人造地表分布无明显的季节差异(表 4)。
表 4 巴基斯坦不同土地覆盖类型使用产品和季节
Table 4
土地覆盖 类型 | 最优土地覆盖数据产品 | 更适用季节 | |||
---|---|---|---|---|---|
GlobCover -2009 | FROM-GLC-2010 | GlobeLand 30-2010 | 夏 | 冬 | |
耕地 | √ | √ | |||
裸地 | √ | √ | |||
植被 | √ | √ | |||
人造地表 | √ | — | — | ||
水体 | √ | √ | |||
冰川积雪 | √ | √ |
研究将3套土地覆盖数据产品在巴基斯坦的精度评价结果(夏冬两季平均精度)与全球精度进行比较。GlobCover-2009 、FROM-GLC-2010和GlobeLand30-2010在巴基斯坦的总精度均低于其全球精度(73% vs. 59.3%[22],71.5 % vs. 60.1%[23],83.5 % vs. 64.8%[24])。而从不同土地覆盖类型来看,GlobCover-2009、FROM-GLC-2010和GlobeLand30-2010人造地表在巴基斯坦的用户精度均高于全球精度(96.4% vs. 63.6%,95.0 % vs. 30.8%,95.6 % vs. 86.9%);从水体在巴基斯坦的用户精度来看,FROM-GLC-2010低于全球(75.7 % vs. 80.6%),GlobCover-2009和GlobeLand30-2010高于全球(82.6 % vs. 74.1%、95.4% vs. 92.1%);从裸地在巴基斯坦的用户精度来看,GlobCover-2009低于全球(74.0% vs. 88.1%),FROM-GLC-2010和GlobeLand30-2010高于全球(73.2 % vs. 62.7%、82.7 % vs. 77.3%);从冰川积雪在巴基斯坦的用户精度来看,FROM-GLC-2010低于全球(83.6% vs. 95.6%),GlobCover-2009和GlobeLand30-2010高于全球(92.3% vs. 74.1%、76.2% vs. 75.8%)。耕地和植被在统一分类时进行了类型归并,因此无法与产品类型的全球精度进行比较。3套全球土地覆盖数据产品在巴基斯坦的精度评价结果中,总精度均有所下降(>10%),各类型精度上也都与全球精度有所偏差(人造地表精度均有所升高,其他类别因产品而异)。巴基斯坦土地覆盖的特殊性是造成这些偏差的原因之一,例如巴基斯坦的裸地覆盖率显著高于全球水平(约40% vs. 14.6%),植被覆盖率显著低于全球水平(约25 vs. 53.2%)[23]等,正确分类不同土地覆盖类型在产品生产中的难易程度不同,因此各种土地覆盖类型的比例会影响土地覆盖数据产品的精度。巴基斯坦土地覆盖类型在不同季节变化显著,从用户精度看,植被夏季比冬季精度高约16%,冰川积雪冬季比夏季高约11%,在这种季节差异下,年际土地覆盖产品生产时使用单一时相影像带来的误差更加明显。此外,相对于全球精度验证,本研究在巴基斯坦区域选择了更加密集的验证样本点(约1个/800 km2,FROM-GLC-2010的全球验证样本共38 664个,约1个/2 000 km2)[33],这也给精度评价结果带来一定影响。
3.3 验证样本点个数对季节性精度评价结果影响
土地覆盖数据产品精度评价结果受样本点数量的影响(图3)。
图3
图3
不同分层采样点数下3套土地覆盖数据产品在夏冬两季的总精度
Fig.3
Total accuracy of three land cover datasets using different stratified sampling points in summer and winter
当样本点个数小于700时,土地覆盖数据产品总精度会随验证样本点数量的变化产生较为剧烈的波动,3套产品夏冬两季分类精度波动的平均值为3.4%,其中精度差异最大是冬季的GlobeLand30-2010(约8%)。当样本点个数大于700时,各产品总精度逐渐趋于稳定,3套产品夏冬两季分类精度波动的平均值为0.7%。由此推测,巴基斯坦土地覆盖分类精度评价合理的验证样本点个数为700个,因此研究选择的1 000个样本点足以反映各种土地覆盖数据产品在巴基斯坦的精度特征。在今后验证土地覆盖数据产品分类精度时,每1 000 km2应至少选择1个样本点。
4 结 论
从土地覆盖制图的角度,GlobeLand30-2010、FROM-GLC-2010和GlobCover-2009 3套土地覆盖数据产品在巴基斯坦地区夏冬两季总分类精度差异不大,GlobeLand30-2010(63.9% vs. 65.6%)与FROM-GLC-2010(59.0% vs. 61.2%)的总精度夏季略高于冬季,GlobCover-2009的总精度冬季略高于夏季(59.5% vs. 59.1%)。此外,土地覆盖数据产品在描述不同的土地覆盖类型时精度有所差别,除人造地表外,各种土地覆盖类型的用户精度均有明显差别,尤其是植被和裸地夏冬两季的分类精度差别很大。因此,在巴基斯坦的土地覆盖制图中,应当充分考虑土地覆盖的季节性特征。
从土地覆盖产品使用的角度,3套产品中GlobeLand30-2010对耕地、人造地表和水体有更好的分类效果,FROM-GLC-2010对植被和冰川积雪的分类更为准确,GlobCover-2009对裸地的分类更为准确;3套产品对耕地、裸地和冰川积雪分类更符合冬季的真实情况,对植被、水体分类更符合夏季的真实情况,人造地表分类无明显的季节差异。总的来说,在使用2010年巴基斯坦土地覆盖数据产品时,若需要全分类的产品可选择总精度最高的GlobeLand30-2010;若需要某一土地覆盖类型的产品可根据生产者精度评价结果进行选择。
目前,全球高精度土地覆盖数据产品均为年际产品,而本研究的精度评价结果显示年内土地覆盖类型存在季节性差异。因此,在未来的土地覆盖制图和产品使用中,应当充分考虑土地覆盖类型的季节性变化及其对产品精度的影响。在制图中,可以针对不同季节分别进行土地覆盖制图,在同一年生产两套或多套数据产品,或者在制图时设置反映季节变化的土地覆盖类型,例如常绿植被和季节型植被等。在产品使用中需要对备选产品进行季节性精度评价,再根据实际情况选用合适的土地覆盖数据产品。此外,在精度评价时,为排除验证样本点数量对结果的影响,应多次抽样以确保使用足够数量的样本点,例如在巴基斯坦,验证样本点密度应达到1个/1 000 km2以上。
参考文献
Land Cover Change and Its Driving Forces in the Republic of Korea Since the 1990s
[J].
1990~2015年韩国土地覆被变化及其驱动因素
[J].
Spatial-temporal Distribution of Cropland in China based on Geomorphologic Regionalization During 1990~2015
[J].
基于地貌分区的1990-2015年中国耕地时空特征变化分析
[J].
Comparing Data Sources in Six Global Data Sets of Historical Land Cover Historical Land Cover
[J].
六个全球历史土地覆盖数据集数据来源的对比分析
[J].
Impacts of Lucc on Soil Properties in the Riparian Zones of Desert Oasis with Remote Sensing Data: A Case Study of the Middle Heihe River Basin, China
[J].
The Status Quo and Expectation of Global and Local Land Cover and Land Use RS Research
[J].
全球及区域尺度土地覆盖土地利用遥感研究的现状和展望
[J].
An Overview of 21 Global and 43 Regional Land Cover Mapping Products
[J].
Different Characteristics of Land Cover Changes in Source Regions of the Yangtze River and the Yellow River in the Past 30 Years
[J].
近30年来长江源区与黄河源区土地覆被及其变化对比分析
[J].
Effects of Land-Use/Cover Change on Hydrological Processes Using a Gis/Rs-based Integrated Hydrological Model: Case Study of the East River, China
[J].
Linkage Analysis of Land Use/Cover Patterns and Hydro-chemical Characteristics in Different Seasons in Ebinur Lake Watershed, China
[J].
Accuracy Assessment of Wetland Categories from the Glob Cover 2009 Data Over China
[J].
全球土地覆盖Clob Cover 2009数据中的中国区域湿地数据精度评价
[J].
Streamflow Responses to Climate Change and Lucc in A Semi-arid Watershed of Chinese Loess Plateau
[J].
Analysis and Evaluation of The Ecological Sensitivity in the Middle Reaches of the Syr Darya River based on GIS
[J].
基于gis的锡尔河中游生态敏感性分析与评价
[J].
Accuracy Evaluation of the Seven Land Cover Data in Qiangtang Plateau
[J].
七套土地覆被数据在羌塘高原的精度评价
[J].
Accuracy Evaluation of the Four Remote Sensing based Land Cover Products over China
[J].
四种常用的全球1km土地覆盖数据中国区域的精度评价
[J].
Accuracy Assessment of Four Global Land Cover Datasets in China
[J].
四类全球土地覆盖数据在中国区域的精度评价
[J].
The First All-season Sample Set for Mapping Global Land Cover with Landsat8 Data
[J].
Exploring the Impact of Seasonality on Urban Land Cover Mapping Using Multi-season Sentinel-1a and GF-1 WFV Images in a Subtropical Monsoon-Climate Region
[J].
Flood Hazard Assessment in Pakistan at Grid Scale
[J].
基于格网的洪水灾害危险性评价分析——以巴基斯坦为例
[J].
Indian Ocean Passages and Their Geopolitical Impacts on China
[J].
印度洋通道及其对中国地缘环境影响
[J].
Finer Resolution Observation and Monitoring of Global Land Cover: First Mapping Results with Landsat Tm and ETM+ Data
[J].
30 Meter Global Land Cover Data Product-Globe Land 30
[J].
全球30 m地表覆盖遥感数据产品——Globe Land 30
[J].
Global Land Cover Mapping at 30m Resolution: A Pok-based Operational Approach
[J].
Accuracy Assessment and Comparative Analysis of GlobeLand30 Dataset in Henan Province
[J].
河南省globeland30数据精度评价及对比分析
[J].
Precision Validation of Multi-sources Land Cover Products Derived from Remote Sensing
[J].
国家尺度异源土地覆被遥感产品精度评价
[J].
Comparison and Assessment of Large-scale Land Cover Datasets in China and Adjacent Regions
[J].
大尺度土地覆盖数据集在中国及周边区域的精度评价
[J].
Precision Analysis and Validation of Multi-sources Landcover Products Derived From Remote Sensing in China
[J].
中国区域多源土地覆被遥感产品精度分析与验证
[J].
Sampling Designs for Accuracy Assessment of Land Cover
[J].
Regional Accuracy Assessments of the First Global Land Cover Dataset at 30-meter Resolution: A Case Study of Henan Province
[J].
首套全球30 m分辨率土地覆被产品区域尺度精度评价——以河南省为例
[J].
Agreement Analysis of Multi-source Land Cover Products Derived from Remote Sensing in South America
[J].戴昭鑫
多源卫星遥感土地覆被产品在南美洲的一致性分析
[J].
/
〈 |
|
〉 |
