中巴经济走廊专题 I 区论文(评审中) 版本 ZH2
下载
2016年中巴经济走廊冻土分布数据集
A dataset of permafrost distribution along the China-Pakistan Economic Corridor in 2016
 >>
: 2018 - 08 - 15
: 2018 - 10 - 15
: 2018 - 10 - 15
545 6 0
摘要&关键词
摘要:中巴经济走廊穿越帕米尔高原和喀喇昆仑山系,在海拔4000 m以上区域广泛发育和分布着多年冻土冻融作用导致的多种地质灾害。研究中巴经济走廊冻土分布和制图是解决其实际工程问题的基础,也对水资源利用、生态安全和边防建设有重要意义。研究区大致空间范围大致为23°47′N–40°55′N,60°20′E–80°16′E,包括中国新疆喀什地区、克孜勒苏柯尔克孜自治州以及巴基斯坦全境。本文收集了2016年MODIS地表温度数据,2009年中国帕米尔高原冰川编目数据,2003–2004年巴基斯坦冰川编目数据和2008年世界土壤数据库(HWSDv1.2),基于多年冻土顶板温度模型(TTOP)得到中巴经济走廊多年冻土分布数据(GeoTIFF格式,空间分辨率1 km)。采用统计学中的确定系数对数据制备方法进行分析和评估,并通过现有文献资料对数据结果质量进行验证。本数据集可作为中巴经济走廊多年冻土变化的本底调查资料,为该区域工程建设中冻胀融沉研究提供基础数据支撑,也可与气候、水文等数据综合分析,揭示中巴经济走廊喀喇昆仑山区水–土–气–生之间定量联系,深化对该区域气候变化背景下生态环境和可持续发展的科学认识。
关键词:中巴经济走廊;多年冻土分布;TTOP模型;MODIS地表温度
Abstract & Keywords
Abstract:  The China-Pakistan Economic Corridor, an iconic project of the Belt and road Initiative, matters greatly to the economical and cultural exchanges of China-Pakistan region. Passing through the Pamirs and Karakoram Mountain System, the Corridor has been frequented by various types of geological disasters at altitudes above 4,000 meters caused by freeze-thaw cycles in permafrost. Studies on permafrost distribution and mapping constitute the basis for solving practical engineering problems along the Corridor, and they are of great importance to water resource utilization, ecological security and border defence. The study has a spatial coverage of approximately 23°47′ N – 40°55′ N, 60°20′ E – 80°16′ E, including Kashgar in Xinjiang, Kizilsu Kirghiz Autonomous Prefecture and the whole Pakistan. The data of permafrost distribution along the Corridor (format: Tiff; spatial resolution: 1 km) were acquired on the basis of TTOP Model, including MODIS surface temperature data in 2016, glacial cataloging data for the Pamirs of China in 2009, glacier cataloging for Pakistan in 2003 – 2004 and World Soil Database for 2008 (HWSD v1.2). The data were then analyzed through statistical methods and validated against existed literature. This dataset can serve as fundamental survey material for studying permafrost changes along the Corridor, providing basic data support for researches on frost heaving and thawing during the construction of the region. Besides, the dataset could be analyzed together with climate, hydrological, and other data to reveal the quantitative relation among the hydrology-soil-atmosphere-ecology of the Corridor, thereby deepening scientific understanding of the region’s ecological environment and sustainable development against the background of climate changes.
Keywords: China-Pakistan Economic Corridor; permafrost distribution; TTOP Model; MODIS land surface temperature
数据库(集)基本信息简介
数据库(集)名称2016年中巴经济走廊冻土分布数据集
数据作者艾鸣浩、张耀南、康建芳、冯克庭、田德宇
数据通信作者张耀南(yaonan@lzb.ac.cn)
数据时间范围2016年1月1日至2016年12月27日,时间间隔为8天
地理区域23°47′ N–40°55′ N, 60°20′ E–80°16′ E
空间分辨率1 km
数据量8.42 MB
数据格式TIF、SHP
数据服务系统网址http://www.crensed.ac.cn/portal/metadata/95e66f86-f868-429f-b62e-a526df42c566
基金项目国家科技基础条件平台“特殊环境特殊功能观测研究台站共享服务平台”(Y719H71006),中国科学院信息化专项“寒旱区环境演变研究‘科技领域云’的建设与应用”(XXH13506)
数据集组成1. 中巴经济走廊多年冻土分布栅格图;2. 中巴经济走廊多年冻土分布矢量图。
Dataset Profile
TitleA dataset of permafrost distribution along the China-Pakistan Economic Corridor in 2016
Data corresponding authorZhang Yaonan (yaonan@lzb.ac.cn)
Data authorsAi Minghao, Zhang Yaonan, Kang Jianfang, Feng keting, Tian Deyu
Time rangeFrom January 1, 2016 to December 27, 2016, at an interval of 8 days
Geographical scope23°47′ N–40°55′ N, 60°20′ E–80°16′ E
Spatial resolution1 km
Data volume8.42 MB
Data formatTIF, SHP
Data service system<http://www.crensed.ac.cn/portal/metadata/95e66f86-f868-429f-b62e-a526df42c566>
Sources of fundingNational R&D Infrastructure and Facility Development Program of China ”National Special Environment and Function of Observation and Research Station Shared Service Platform” (Y719H71006); Information Program of the Chinese Academy of Sciences “Construction and Application of ‘Technology Cloud’ in Cold and Arid Regions Environment Evolution” (XXH13506).
Dataset composition1.China-Pakistan Economic Corridor Permafrost Distribution Raster;
2.China-Pakistan Economic Corridor Permafrost Distribution Vetor.
引 言
冻土是北半球冰冻圈分布最广泛的因子,占陆地面积的56%,其中多年冻土约占24%[1]。多年冻土具有热稳定性差、水热活动强烈、厚层地下冰和高含冰量冻土所占比重大等特点,在全球气候变暖背景下,多年冻土对工程活动的热扰动极为敏感。在多年冻土区工程建设中,施工期不可避免地会改变包括地表覆盖、地形因子及浅层地表能量平衡等局地因素,加上开通运营后动载荷对冻融循环水热环境的改变,直接导致冻土温度升高和上限抬升,进而增加了冻胀、融沉和冻拔等冻土灾害风险,对多年冻土区工程施工和安全运营构成严重威胁,也对区域水文地质条件和生态环境演化产生巨大影响[2,3,4 ]
中巴经济走廊作为“一带一路”倡议的重要组成部分,是中国与中亚、西亚、南亚地区的交通要道,对中国和巴基斯坦以及周边国际的经济文化交流有着重要意义。走廊穿越喀喇昆仑山系,在海拔4300–4750 m以上高寒山区,冻胀草丘,热融塘湖,石环、石链等冻土地貌广泛分布。正确的评价中巴走廊多年冻土分布问题是规划解决中巴经济走廊公路、铁路、通信和输油管线“四位一体”实际工程问题的基础,也对走廊内水资源利用、生态环境保护具有重要意义,尤其对边防建设有重要的战略意义。
传统的冻土调查通常采用野外工程勘查进行,但中巴走廊沿线地域广阔,实测资料较少,现有的多年冻土考察主要在中巴公路沿线,冻土资料基本一片空白。由于冻土与气候系统,特别是温度条件关系密切,国内外依据冻土与气温相互作用特性发展了一系列冻土分布模型。Simth等在分析气候-冻土关系基础上,考虑局部地形和土壤条件提出一种联系地温和气温的,半经验半物理的多年冻土顶板温度模型(TTOP)[5]。TTOP模型目前已应用于加拿大全国尺度和区域尺度如Machensize河流域、布拉多高原和昂加瓦湾等的冻土分布模拟制图上,取得了良好的效果[6,7,8 ]
随着遥感技术的快速发展,通过遥感技术可以大面积多时相地获得与多年冻土分布有关的地面信息进行空间分析,结合现有的冻土空间分布模型,可以实现宏观、动态、快速的实现大尺度上的冻土分布模拟。本研究以中巴经济走廊为研究区,包括中国新疆的喀什地区、克孜勒苏柯尔克孜自治州以及巴基斯坦全境(图1)。利用MODIS地表温度产品作为数据源,代替实测数据输入到TTOP模型中,制备中巴走廊多年冻土分布数据集。


图1   研究区范围地理位置示意图
1   数据采集和处理方法
1.1   数据来源与预处理
中巴经济走廊多年冻土分布数据制备基于4种数据:2016年MODIS地表温度产品MOD11A2 LST数据,2009年中国帕米尔高原第二次冰川编目数据集(v1.0),2003–2004年巴基斯坦冰川编目数据集,2008年世界土壤数据库(Harmonized World Soil Database version 1.2,即HWSDv1.2),如表1所示。
地表温度作为地表能量平衡中的主要参数以及气候系统的主要因子,能够表征地气间能量和水分交换的程度,是影响冻土发育、分布以及演化的关键因素之一,也是冻土建模的上边界条件。MOD11A2采用劈窗法从MOIDS数据中反演得到地表温度,为8天合成L3级产品,空间分辨率1 km,包括白天和夜间数据,来源于美国国家航空航天局(NASA)陆地过程分布式数据档案中心(LPDAAC)。数据以HDF-EOS格式存储,并采用正弦曲线投影坐标系。覆盖整个中巴走廊需要6景MOD11A2数据。本研究获取2016年全年MOD11A2 L3数据共5544景。首先使用MRT工具对数据进行镶嵌和重采样,将所有数据坐标系转为WGS-84坐标系,然后利用中巴走廊矢量边界和开源栅格空间数据转换库GDAL中的warp工具进行数据裁剪,最后得到462景研究区LST数据。
利用地表温度计算冻土分布会将部分冰川错分类成冻土,需要根据研究区内冰川分布数据将其去除。中国帕米尔高原第二次冰川编目数据集和巴基斯坦冰川编目均采用Landsat TM/ETM+数据经过校正后通过自动提取和专家干预修订获得,数据来源于寒区旱区科学数据中心(WestDC)[9-10]
TTOP模型需要不同类型土壤的冻融热导系数作为输入参数。土壤数据来源于2008年维也纳国际应用系统分析研究所(IIASA)和联合国粮农组织(FAO)所构建的HWSDv1.2。HWSDv1.2包括不同种类土壤的全球分布情况,土壤分类标准为FAO-90,数据分辨率1 km,数据投影WGS84。
将冰川编目数据和土壤类型数据使用中巴经济走廊矢量边界裁剪,形成研究区内的冰川分布数据和土壤类型分布数据。
表1   研究采用数据列表
序号名称时间来源类型
1MOD11A2 LST2016年LPDAAC栅格
2中国帕米尔高原第二次冰川编目数据集 (v1.0)2009年WestDC矢量
3巴基斯坦冰川编目数据集2003–2004年WestDC矢量
4世界土壤数据库(v1.2)2008年FAO栅格
1.2   数据处理步骤
1.2.1   TTOP模型
TTOP模型其实质是采用活动层顶板温度来估计多年冻土的热状况。模型假定土壤基质匀质,在年尺度上,地–气热交换界面处于热平衡状态。TTOP模型定义如下:
\(\mathrm{T}\mathrm{T}\mathrm{O}\mathrm{P}=\mathrm{ }\frac{\frac{{K}_{t}}{{K}_{f}}({n}_{t}\bullet DDT-{n}_{f}\bullet DDF)}{P}\) (1)
式中,\(\mathrm{T}\mathrm{T}\mathrm{O}\mathrm{P}\)为多年冻土判别指标,\(\mathrm{T}\mathrm{T}\mathrm{O}\mathrm{P}<0\)表示存在多年冻土,\(\mathrm{T}\mathrm{T}\mathrm{O}\mathrm{P}>0\)表示季节性冻土或无冻土;P为周期,取365天;\({K}_{t}\)\({K}_{f}\)分别为土壤融化和冻结时的导热系数;\(\mathrm{D}\mathrm{D}\mathrm{T}\)为融化指数,指一年中连续高于0℃的持续时间和其数值乘积的总和,以\(\mathrm{℃}\bullet day\)表示;\(\mathrm{D}\mathrm{D}\mathrm{F}\)为冻结指数,指一年中低于0℃的持续时间和其数值乘积的总和,以\(\mathrm{℃}\bullet day\)表示。\({n}_{t}\)\({n}_{f}\)分别为融化和冻结因子,是由于地表植被和雪盖会影响地气热交换而引入融化和冻结因子修正模型结果。中巴走廊多年冻土发育地区处于高寒山区,也是极端干旱区,植被稀疏,本研究忽略雪盖和植被影响,将\({n}_{t}\)\({n}_{f}\)都赋值为1。
1.2.2   地表温度数据日尺度重建
\(\mathrm{T}\mathrm{T}\mathrm{O}\mathrm{P}\) 模型输入数据\(\mathrm{D}\mathrm{D}\mathrm{T}\)\(\mathrm{D}\mathrm{D}\mathrm{F}\)均为日尺度指数,因此需要将MODIS LST的8天合成地表温度产品重建到日尺度。地表温度通常表现出明显的周期性,李述训等采用余弦函数对青藏高原时间序列地表温度进行模拟[11],建立了地表温度与时间的函数关系,表达式如下:
\({T}_{s}\left(t\right)=\stackrel{-}{T}+Acos(\omega t+\phi )\) (2)
其中,\({T}_{s}\)为地表温度,\(t\)为时间,\(\stackrel{-}{T}\)为时间序列内的平均温度;\(A\)为时间序列内的温度振幅;\(\omega =2\pi /T\mathrm{中},T\)为时间序列周期,\(\phi \)为相位角。
上述模拟方法虽能够反映地表温度周期变化情况,但参数计算不够准确。以平均温度和温度振幅作为余弦函数系数使函数拟合后整体偏差较大。本研究以\(\stackrel{-}{T}\)\(A\)\(\phi \)为待定系数,利用最小二乘法对全年46景数据每一景数据对应像元进行数值拟合,得到拟合函数,再离散化到日尺度,完成地表温度数据的重建。
1.2.3   导热系数
土壤融化和冻结时的导热系数是\(\mathrm{T}\mathrm{T}\mathrm{O}\mathrm{P}\)模型中与土壤性质有关的参数,其值取决于土壤颗粒级配、干容重、矿物成分和含水量等因素。由于缺少中巴走廊土壤性质的详细参数,难以计算每一网格上的土壤导热系数。王之夏等在研究青藏高原多年冻土分布时基于青藏高原第四纪地图,将土壤按地质成因分类,然后针对每一类采用Johansen方法计算其冻融状态下的导热系数[12,13,14 ]。由于此方法基于土壤本身性质,地区差异性相对较小,因此本研究认为可用在中巴经济走廊。HWSD数据包含土壤分类信息,可获取栅格上土壤分类名称,但采用的FAO土壤分类体系是以土壤诊断层和诊断特性为基础,与王之夏等以土壤的工程地质成因的分类方式并不完全相同,经过对比筛选后得到基于HWSD数据的冻融导热系数,见表2。
表2   不同类型土壤的冻融导热系数表
序号土壤类型FAO分类缩写\({K}_{t }(W/m\bullet ℃)\)\({K}_{f} (W/m\bullet ℃)\)
1薄层土LP1.342.32
2冲积土FL1.552.17
3有机土HS0.521.70
4冰碛物GG1.652.78
5其他CL、RK等1.502.20
1.2.4   多年冻土分布模拟流程
首先重建地表温度数据,逐像元计算\(\mathrm{D}\mathrm{D}\mathrm{T}\)\(\mathrm{D}\mathrm{D}\mathrm{F}\);然后制作研究区土壤融化导热系数图和土壤冻结导热系数图;最后将以上数据输入\(\mathrm{T}\mathrm{T}\mathrm{O}\mathrm{P}\)模型中,剔除冰川覆盖区域,通过二值分类得到中巴走廊多年冻土分布图。数据处理流程如图2所示:


图2   地表温度数据时间序列拟合示意图
2   数据样本描述
中巴经济走廊多年冻土分布数据集分为两个文件夹存储,分别为冻土分布栅格数据和冻土分布矢量数据。前者文件夹包括冻土分布栅格图CPEC_PermafrostDistribution2016,空间分辨率1 km,保存格式为tif。后者文件夹内为多年冻土分布矢量图CPEC_PermafrostDistribution2016,采用SHP格式保存。所有数据地理坐标系为WGS1984。数据结果展示如图3。


图3   中巴经济走廊多年冻土分布图
3   数据质量控制和评估
3.1   数据重建质量评估
以余弦函数和最小二乘法拟合的MODIS地表温度数据可用确定系数来评价。确定系数表示一个随机变量与多个随机变量关系的数字特征,用来反映模型对样本数据的近似程度,常以R表示。R值域为\((0,\mathrm{ }1)\),越接近1表明拟合的优度越高,其表达式如下:
\(R=\frac{\sum_{i=1}^{n}(\hat{y_{l}}-\bar{y})^{2} }{\sum_{i=1}^{n}(y_{i}-\bar{y})^{2}}\) (3)
其中,\(\stackrel{-}{y}\)为样本数据均值,\(\hat{y_{l}}\)为模型预测值,\({y}_{i}\)为样本值。
中巴经济走廊边界区域内共有1061933像素点,按时间序列拟合后R最大值0.99,最小值0.43。R值大于0.9的像素点占总数的69.71%,处于0.8和0.9之间的占总数27.64%,小于0.8的占总数2.65%,频率直方图如图4所示。因此,拟合重建的方式在中巴走廊冻土分布模拟中是适用的。


图4   拟合R值频率分布直方图
3.2   多年冻土分布模拟结果评估
由于缺少冻土分布的实地勘测资料,因此以文献所提中巴公路沿线冻土数据作为验证。张学进等人提出中巴多年冻土分布于海拔4500 m以上的红其拉普路段K807+000至K811+343.165路段内,连续分布长度4.3 km[15]。朱颖彦等人提到中巴走廊多年冻土分布在红其拉普附近,海拔4300 m以上路段[16]。两篇文献提到多年冻土位置接近,可互为验证。本研究以遥感数据和模型模拟的方式得到多年冻土分布图,在中巴公路沿线所得结果如图5所示,与以上两篇文献吻合。


图5   红其拉普冻土分布图
4   数据价值
传统冻土图绘制主要以野外实测数据为基础,如钻孔地温数据,常规气象数据等,虽然结果精度较高,但其必须建立在时间序列足够长、观测数据足够可靠的基础上。中巴经济走廊多年冻土分布区域自然环境严苛,地缘政治环境复杂,在此区域进行冻土的调查和原位观测具有较多限制。
本文基于遥感数据MODIS地表温度数据产品,利用\(\mathrm{T}\mathrm{T}\mathrm{O}\mathrm{P}\)模型模拟中巴走廊冻土分布情况。针对数据与模型时间尺度不匹配的情况,采用拟合方式对数据进行时间序列重建。使用世界土壤数据集数据估算不同种类土壤冻融导热系数作为\(\mathrm{T}\mathrm{T}\mathrm{O}\mathrm{P}\)模型参数,最后获得走廊多年冻土分布图,并采用决定系数对数据重建过程中的拟合优度进行评价。野外实测是调查冻土分布的最佳手段,本文搜集了文献中关于多年冻土分布位置的描述,通过和模拟结果进行比较可知,基于MODIS地表温度数据和\(\mathrm{T}\mathrm{T}\mathrm{O}\mathrm{P}\)模型的多年冻土分布模拟结果与文献中的野外实测结果具有较高的一致性。
中巴经济走廊多年冻土分布数据是走廊内工程选址、设计和规划的重要基础数据,也可作为全球气候变化背景下,走廊多年冻土变化的本底数据。在多年冻土分布数据基础上,结合气候、水文、生态等数据综合分析,对区域内工程长期安全运营和生态环境可持续发展有着重要意义。
5   数据使用方法和建议
2016中巴经济走廊冻土分布数据集保存为栅格TIF格式和矢量SHP格式。 ArcGIS, QGIS, ENVI, ERDAS等常用的GIS与遥感软件可支持该数据的读取和操作。
致 谢
感谢LPDAAC提供的MODIS数据,感谢FAO提供的世界土壤数据库HWSD v1.2,感谢寒区旱区科学数据中心提供的中国帕米尔高原第二次冰川编目数据和巴基斯坦冰川编目数据。
[1]
ZHANG T, BARRY R G, KONWLES K, et al. Statistics and characteristics of permafrost and ground-ice distribution in the Northern Hemisphere[J]. Polar Geogr, 1999, 23: 132-154.
[2]
SHENG Y, ZHANG J M, LIU Y Z, et al. Thermal regime in the embankment of Qinghai-Tibet Highway in permafrost regions[J]. Cold Regions Sci Technol, 2002, 35: 25-44.
[3]
WU Q B, CHENG G D, MA W, et al. Technical approaches on permafrost thermal stability for Qinghai-Tibet Railway[J]. Geomech Geoeng: Int J, 2006, 1: 119-127.
[4]
MA W, MU Y, LI G, et al. Responses of embankment thermal regime to engineering activities and climate change along the Qinghai-Tibet Railway(in Chinese)[J]. Scientia Sinica Terrae, 2013, 43:478-489.
[5]
SMITH M W, RISEBOROUGH D W. Permafrost monitoring and detection of climate change[J]. Permafrost and Periglacial Process, 1996, 7: 301-309.
[6]
HENRY K, AND M W. Smith 2001. A Model-Based Map of Ground Temperatures for the Permafrost Regions of Canada[J]. Permafrost and Periglacial Processes, 2001, 12: 389-398.
[7]
WRIGHT JF, DUCHESNE C, COTE M M. Regional-scale permafrost mapping using the TTOP ground temperature model[C]. In Proceedings of the Eighth International Conference on permafrost. Zurich, Switzerland, Lisse, 2003.
[8]
WAY R G, LEWKOWICZ A G. Modelling the spatial distribution of permafrost in Labrador–Ungava using the temperature at the top of permafrost[J/OL]. Canadian Journal of Earth Sciences, 2016, 53(10): 1010-1028. https://doi.org/10.1139/cjes-2016-0034.
[9]
刘时银, 郭万钦. 中国帕米尔高原第二次冰川编目数据集(V1.0)[J/OL]. 寒区旱区科学数据中心, 2014. DOI:10.3972/glacier.004.2013.db.
[10]
MOOL P K, BAJRACHARYA S R, SHRESTHA B, et al. Inventory of Glaciers and Glacial Lakes and the Identification of Potential Glacial Lake Outburst Floods (GLOFs) Affected by Global Warming in the Mountains of Himalayan Region(Indus Basin Pakistan Himalaya) [M]. Kathmandu, Nepal: ICIMOD, 2005.
[11]
李述训, 程国栋. 冻融土中的水热运输问题[M]. 兰州: 兰州大学出版社, 1995.
[12]
王之夏. 多年冻土分布模型在青藏高原的应用研究[D]. 兰州: 兰州大学, 2010.
[13]
JOHANSEN O. Thermal Conductivity of Soil[D]. Trondheim: University of Trondheim, Norway, 1975.
[14]
LI X, CHENG G D, WU Q B et al. Modeling Chinese cryospheric change by using GIS technology[J]. Cold Regions Science and Technology, 2003, 36(1-3): 1-9.
[15]
张学进. 中巴公路沿线地质灾害分布特征及防治[J]. 重庆工商大学学报(自然科学版). 2013, 30(2): 45-50.
[16]
朱颖彦, 杨志全, 廖丽萍. 中巴喀喇昆仑公路冰川地貌地质灾害[J]. 灾害学, 2014, 2014(3): 81-90.
数据引用格式
艾鸣浩, 张耀南, 康建芳, 等. 2016年中巴经济走廊冻土分布数据集[DB/OL]. 国家特殊环境、特殊功能测研究台站共享平台, 2018. (2018-06-13). DOI: 10.12072/casnw.049.2018.db.
稿件与作者信息
论文引用格式
艾鸣浩, 张耀南, 康建芳, 等. 2016年中巴经济走廊冻土分布数据集[J/OL]. 中国科学数据, 2018. (2018-10-12). DOI: 10.11922/csdata.2018.0053.zh.
艾鸣浩
Ai Minghao
主要承担工作:时间序列地表温度重建, 模型流程设计与实现。
aimh@lzb.ac.cn
(1986—),男,山东省济宁市人,硕士,工程师,研究方向为寒旱区大数据应用。
张耀南
Zhang Yaonan
主要承担工作:数据处理流程设计。
(1966—),男,甘肃省天水市人,博士,研究员,研究方向为地学大数据。
康建芳
Kang Janfang
主要承担工作:数据搜集与分析。
(1980—),女,甘肃省天水市人,硕士,工程师,研究方向为遥感数据制备方法。
冯克庭
Feng Keting
主要承担工作:数据前处理方案设计。
(1980—),男,宁夏省回族自治区中卫市人,博士研究生,工程师,研究方向为遥感干旱。
田德宇
Tian Deyu
主要承担工作:数据批处理实现。
(1993—),男,内蒙古自治区四子王旗人,硕士研究生,学生,研究方向为寒旱区遥感应用。
出版历史
I区发布时间:2018年10月15日 ( 版本ZH2
参考文献列表中查看
中国科学数据
csdata