摘要&关键词
摘要：中巴经济走廊是“一带一路”战略大局的重要组成部分，对“一带一路”战略实施发挥着重大的示范和推动作用。本数据集由高分一号、高分二号国产卫星影像制作而成，空间范围为23°54′N–39°12′N、71°24′E–76°48′E，中巴公路两侧约60 km，中国喀什地区至巴基斯坦伊斯兰堡段，时间范围为2013–2017年，全色与多光谱2 m融合影像，TIFF格式。数据集正射校正精度最大误差X方向为0.35个像元，Y方向为0.4个像元；均方根误差X方向为0.42个像元，Y方向为0.38个像元。从图像融合效果看，高分影像采用PanSharpening方法融合效果较好。本数据集可以丰富本区域空间基础数据资源，能够应用于中巴经济走廊基础设施建设规划、自然灾害预警与生态安全评价等领域。
关键词：中巴经济走廊；正射影像数据集；数据预处理；正射矫正；数据融合

Abstract & Keywords
Abstract:?The China–Pakistan Economic Corridor is an important part of the Belt and Road. Its implementation plays an important role in the demonstration and promotion of the Belt and Road Initiative. The production of Digital Orthophoto Map(DOM) dataset is made up of high-resolution GF-1 and GF-2 satellite images. This dataset has a spatial coverage of 23°54'N – 39°12'N, 71°24'E – 76°48'E, which is within 60km on both sides of KKH (Kashgar to Islamabad), and a temporal span from 2013 to 2017. It contains panchromatic and multispectral fusion images in TIFF format, with a 2-m spatial resolution. The maximum error of orthorectification correction is 0.35 pixel in the X direction, and 0.4 pixel in the Y direction, while that of root mean square is 0.42 pixel in the X direction, and 0.38 pixel in the Y direction. PanSharpening method is used to fuse these high-resolution images. This dataset can enrich the spatial basic data resources of the region and can be used as basic data for infrastructure planning, natural disaster early-warning and ecological safety evaluation along the China–Pakistan Economic Corridor.
Keywords:?China–Pakistan Economic Corridor;?DOM dataset;?data preprocessing;?orthorectification;?data fusion

数据库（集）基本信息简介

数据库（集）名称	中巴经济走廊（喀什至伊斯兰堡段）高分正射影像数据集
数据作者	韩立钦、张耀南、田德宇、康建芳
数据通信作者	张耀南（yaonan@lzb.ac.cn）
数据时间范围	2013–2017年
地理区域	23°54′N–39°12′N、71°24′E–76°48′E
空间分辨率	2 m
数据量	674.8 GB
数据格式	TIFF
数据服务系统网址	http://www.crensed.ac.cn/portal/metadata/f960952b-00f3-42a9-9bd8-07e75c22f347
基金项目	国家科技基础条件平台（Y719H71006），中国科学院信息化专项（XXH13506），甘肃省高校科技转化项目（2017D-27）
数据集组成	数据集共6种数据文件，主体由两部分组成，一是正射影像样例数据，二是DEM数据，另外两个说明文件分别是中巴经济走廊正射影像数据集列表、数据文件说明。正射影像命名示例为GF1_PMS1_E75.2_N38.6_20150916_L1A0001042400-MSS1_ORTHO_PSH.tif，DEM数据命名为China-Pakistan Economic Corridor DEM.tif，数据集列表命名为China-Pakistan Corridor DOM list.xls，数据文件说明命名为data file description.docx。

Dataset Profile

Title	A high-resolution DOM dataset of the China-Pakistan Economic Corridor (Kashgar to Islamabad)
Data corresponding author	Zhang Yaonan (yaonan@lzb.ac.cn)
Data authors	Han Liqin, Zhang Yaonan, Tian Deyu，Kang Jianfang
Time range	From 2013 to 2017
Geographical scope	23°54′ N – 39°12′ N, 71°24′ E– 76°48′ E
Spatial resolution	2 m
Data volume	674.8 GB
Data format	TIFF
Data service system	<http://www.crensed.ac.cn/portal/metadata/f960952b-00f3-42a9-9bd8-07e75c22f347>
Sources of funding	Data Sharing Fundamental Program for the Construction of the National Science and Technology Infrastructure Platform (Y719H71006); The 13th Five-Year Informatization Plan of Chinese Academy of Sciences (XXH13506); Science and Technology Transformation Project of Gansu High Education (2017D-27).
Dataset composition	The dataset consists of two parts of data, one for sample DOM and the other for DEM .It comprises six data files in total. The data files are recorded as: GF1_PMS1_E75.2_N38.6_20150916_L1A0001042400-MSS1_ORTHO_PSH.tif and China-Pakistan Economic Corridor DEM.tif and China-Pakistan Economic Corridor DOM lists.xls and data file description.docx.

引 言
中巴经济走廊是“一带一路”战略大局的重要组成部分，北起中国喀什、南至巴基斯坦瓜达尔港，主要沿中巴公路展开，包括公路、铁路、油气和光缆通道在内的一条经济大动脉^[1]。走廊穿越喜马拉雅山、喀喇昆仑山、兴都库什山，地势北高南低（海拔460–4750 m），峡谷深切、冰川广布、雪峰林立^[2]。气候垂直分带明显，植被受气候、海拔、地形影响，不均匀发育显著，特殊的自然地理条件使研究区内冰崩、雪崩、冻融、滑坡、泥石流、岩崩、洪水、冰湖溃决等各类地质灾害频发，给基础设施规划建设带来了极大的挑战^{[3,4,5 ]}。
近年来，随着SPOT等一系列卫星遥感数据的推出，高分影像在区域制图^[6]、灾害监测^[7]和资源调查^[8-9]等领域迅速普及。以高分对地观测系统宏观、快速、动态、高空间、高时间、高光谱的优势开展中巴经济走廊地质灾害专题制图、灾情监测，同时融合气象、水文、地质、冰川、冻土、土壤类型等数据，建立孕灾认知、演变分析的机理模型，将对该区域基础设施建设提供有力的数据支持和保障^[10-11]。
本数据集全部采用国产高分专项数据资源，也是国内大规模利用国产高分数据制备中巴经济走廊基础影像数据的一次尝试，实现了中巴经济走廊（喀什至伊斯兰堡段）2 m分辨率，中巴公路两侧约60 km的数字正射影像图（DOM，Digital Orthophoto Map），其中中国境内415公里（G314国道），巴基斯坦境内809公里（N35国道）。本数据集丰富了此区域空间数据资源，同时拓展了高分数据的应用领域。

1 ? 数据采集和处理方法
1.1 ? 数据来源与预处理
数据源主要有卫星影像和数字高程模型，如表1。中巴经济走廊（喀什至伊斯兰堡段）正射影像数据全部来源于我国高分辨率对地观测系统高分一号和高分二号卫星，高分一号为2013–2017年成像的全色2 m、多光谱8 m数据，高分二号为2015–2017年成像的全色1 m、多光谱4 m数据。在成像条件较好的中巴经济走廊中国段及巴基斯坦北部地区，影像云量控制在5%以下，巴基斯坦中部地区云量控制在10%以下，并进行相应的去云处理。其中高分一号数据占97.2%，高分二号为影像空白区补充数据。
表1 ? 研究采用数据列表

序号	名称	时间	来源	类型
1	中巴经济走廊（喀什至伊斯兰堡段）高分一号影像	2013–2017年	高分辨率对地观测系统甘肃数据与应用中心	栅格
2	中巴经济走廊（喀什至伊斯兰堡段）高分二号影像	2015–2017年	高分辨率对地观测系统甘肃数据与应用中心	栅格
3	ASTER GDEM数据集	2009年	美国航天局与日本产经省	栅格

本数据集通过ENVI和ArcGIS软件工具，利用数字高程模型对高分一号、高分二号影像数据进行正射校正、全色与多光谱融合、匀色、镶嵌等处理。由于高分一号和高分二号卫星影像采集的时间分辨率限制，研究区影像资源横跨5个年份，考虑到不同成像时间的季节差异，对高分一号所有影像数据进行统一匀色。对于高分一号未覆盖区域，利用高分二号影像资源作为补充，独立进行匀色处理。最终按照一定范围裁切生成数字正射影像分块数据集，供用户分块检索使用（图1）。其中，正射校正、数据融合的精度与效果是数据生产的关键。




图1 ? 正射影像生产流程
ASTER GDEM^[12]数据是NASA根据新一代对地观测卫星Terra星载热辐射和反射计（即ASTER）3N和3B波段垂直向下观测成像，历时近10年制作完成。数据空间分辨率为1″×1″（约30 m×30 m）,全球范围内置信度为95%时，垂直精度约为20 m。本研究区域多为高寒山区，地面高精度控制数据几乎为空白，DEM精度很难保证。
赵国松^[13]等在华中地区3°×3°的研究区分析发现，ASTER GDEM产品高程值比STRM DEM平均低5.42 m。刘勇^[14]指出山区的STRM DEM精度要比平原地区低很多，但是总体精度可以接受。如图2所示，在对公格尔–慕士塔格地区的ASTER GDEM、STRM DEM和ASTER L1T单景DEM统计分析发现，ASTER L1T单景DEM有明显异常，ASTER GDEM和STRM DEM无明显异常，但是STRM DEM最大值比慕士塔格峰主峰海拔多出200 m以上，而ASTER GDEM与实地测量值基本相符，因此，本数据集利用ASTER GDEM作为正射校正的高程数据源。




图2 ? 高程数据统计分析结果

1.2 ? 数据处理方法
1.2.1 ? RPC正射校正模型
卫星影像正射校正模型分为严格的几何纠正模型和近似纠正模型。本数据集无法精确地获知高分一号和高分二号的传感器参数与轨道参数等，不能建立共线方程，而多项式近似纠正模型适用于地面平坦地区。因此，面对地形较为复杂的中巴经济走廊地区，选择RPC有理函数模型法进行影像正射校正。
有理函数模型^[15]是将影像中的像点在像方坐标系中的坐标（r,c）和其对应的地面点，经过空间转换投影后在物方坐标系中的空间坐标（X,Y,Z），用比值多项式关联起来。为了增强参数求解的稳定性，将地面坐标和影像坐标正则化到?1.0～1.0之间。
定义如下：
（1）
在式（1）中，（X_n,Y_n,Z_n ）为正则化的地面点坐标，（r_n,c_n ）为正则化的像点坐标。多项式P₁,P₂,P₃,P₄ 如式（2）：
（2）
多项式系数0≤n₁ ≤3；0≤n₂ ≤3；0≤n₃ ≤3，且n₁ +n₂ +n₃ =3，每个多项式P₁,P₂,P₃,P₄ 都是包含了20个系数的三次多项式，如式（3）:
P₁=a₀+a₁X+a₂Y+a₃Z+a₄XY+a₅XZ+a₆YZ+a₇X²+a₈Y²+a₉Z²+a₁₀XYZ+a₁₁X³+a₁₂X²Y+
a₁₃X²Z+a₁₄XY²+a₁₅XZ²+a₁₆Y³+a₁₇Y²Z+a₁₈YZ²+a₁₉Z³;
P₂=b₀+b₁X+b₂Y+b₃Z+b₄XY+b₅XZ+b₆YZ+b₇X²+b₈Y²+b₉Z²+b₁₀XYZ+b₁₁X³+b₁₂X²Y+
b₁₃X²Z+b₁₄XY²+b₁₅XZ²+b₁₆Y³+b₁₇Y²Z+b₁₈YZ²+b₁₉Z³; （3）
P₃=c₀+c₁X+c₂Y+c₃Z+c₄XY+c₅XZ+c₆YZ+c₇X²+c₈Y²+c₉Z²+c₁₀XYZ+c₁₁X³+c₁₂X²Y+
c₁₃X²Z+c₁₄XY²+c₁₅XZ²+c₁₆Y³+c₁₇Y²Z+c₁₈YZ²+c₁₉Z³;
P₄=d₀+d₁X+d₂Y+d₃Z+d₄XY+d₅XZ+d₆YZ+d₇X²+d₈Y²+d₉Z²+d₁₀XYZ+d₁₁X³+d₁₂X²Y+
d₁₃X²Z+d₁₄XY²+d₁₅XZ²+d₁₆Y³+d₁₇Y²Z+d₁₈YZ²+d₁₉Z³;
且满足：
（4）
式（4）中的X₀,Y₀,Z₀,X_s,Y_s,Z_s 为地面坐标进行位移的比例系数；r₀,c₀,r_s,c_s 为影像坐标进行位移的比例系数。式（3）中的a,b,c,d是多项式系数，且a₀,b₀,c₀,d₀ 通常取值为1。

1.2.2 ? 影像融合方法
遥感影像融合^[16-17]是对多个传感器影像进行多级别、多层次的处理过程，以达到提高影像分辨率、增强目标特征、提高分类精度的目的，是提升遥感影像应用能力的重要手段。本数据集采用最邻近扩散全色锐化法（Nearest Neighbor Diffusion pan-sharpening，PanSharpening）完成高分一号和高分二号卫星影像全色与多光谱数据融合。与比值变换法（Brovey Transform，Brovey）、主成份变换法（Principal Component Analysis，PCA）和正交变换法（Gram-Schmidt，GS）等方法相比较，PanSharpening方法能够较好地保持多光谱影像的色彩、纹理和光谱特征，通过反映影像质量的均值、标准差、平均梯度、光谱质量的偏差指数、相关系数和影像信息量的交叉熵几个方面进行影像质量定量评价。
（1）均值(μ)
影像的全部像元灰度的算术平均值，反映了影像中地物的平均反射率。其中，F(i,j)为融合影像F在像素点(i,j)处的灰度值，M和N为影像F的大小。均值越高，则影像整体亮度越高。
μ =\(\frac{1}{M×N}\sum _{i=1}^{M}\sum _{j=1}^{N}F\left(i,j\right)\) （5）
（2）标准差(std)
标准差由均值间接得到，表示影像灰度像素值与平均值的离散程度。其中，F(i,j)为融合影像F在像素点(i,j)处的灰度值，M和N为影像F的大小，\(\mu \mathrm{为}\mathrm{灰}\mathrm{度}\mathrm{平}\mathrm{均}\mathrm{值}\)，标准差越大，灰度离散程度越大，影像反差越大，视觉效果越好。
\(std=\sqrt{\frac{1}{M\times N}\sum_{i=1}^{M}\sum_{j=1}^{N}\left ( F\left ( i,j \right )-\mu \right )}\) （6）
（3）平均梯度(g)
平均梯度反映了影像的平均灰度变化率，即清晰度。F(i,j)为融合影像F在像素点(i,j)处的灰度值，M和N为影像F的大小，其值表示融合影像中的微小细节反差和纹理变化特征，在融合影像中，平均梯度越大，影像清晰度越高。
g=\(\frac{1}{M×N}\sum _{i=1}^{M}\sum _{j=1}^{N}\sqrt{\frac{\left(\left(\frac{\partial F\left(i,j\right)}{\partial i}\right)2+\left(\frac{\partial F\left(i,j\right)}{\partial j}\right)2\right)}{2}}\) （7）
（4）偏差指数(dc)
光谱扭曲度直接反映了融合影像对原光谱影像的失真程度。其值表示融合影像与原多光谱影像像元灰度值的差异和匹配程度。其中，F(i,j)为融合影像F在像素点(i,j)处的灰度值，M和N为影像F的大小，A(i,j)表示原多光谱影像在像素点(i,j)处的灰度值，偏差指数越大，影像失真越强烈。
dc=\(\frac{1}{M×N}\sum _{i=1}^{M}\sum _{j=1}^{N}\frac{|F\left(i,j\right)-\mathrm{A}\left(i,j\right)|}{\mathrm{A}\left(i,j\right)}\) （8）
（5）相关系数(cc)
相关系数反映了融合影像与源影像之间光谱特征的相关程度，以及融合影像光谱信息的保持能力。其中，F(i,j)为融合影像F在像素点(i,j)处的灰度值，M和N为影像F的大小，A(i,j)表示原多光谱影像在像素点(i,j)处的灰度值，\(\mu F\)与\(\mu A\)分别表示融合影像与源影像的灰度平均值。相关系数越大，融合影像从源影像中获取的信息越多，融合效果越好。
cc=\(\frac{\sum _{i=1}^{M}\sum _{j=1}^{N}\left(F\left(i,j\right)-\mu F\right)\left(A\left(i,j\right)-\mu A\right)}{\sum _{i=1}^{M}\sum _{j=1}^{N}\left(F\left(i,j\right)-\mu F\right)2\left(A\left(i,j\right)-\mu A\right)2}\) （9）
（6）影像信息量( )
影像的熵值反映了影像信息的丰富程度。交叉熵（ce）用来衡量A、F两幅影像灰度分布的差异。对于一幅单一影像，各像素的灰度值是相互独立的，则影像灰度分布为P=｛P₀,P₁,…P_i,…P_n ｝，P_i 表示影像像素灰度值为i的概率，即灰度值为i的像素与影像总像素之比，l为影像总的灰度级数,,表示两幅影像像元灰度为i的概率。交叉熵越小，融合影像灰度分布与源影像的差异越小，即融合影像包含源影像的信息量越多，融合效果越好。
（10）



2 ? 数据样本描述
本数据集分为中巴经济走廊（喀什至伊斯兰堡段）正射影像分幅数据、正射影像分幅索引数据和ASTER DEM数据。正射影像为2 m分辨率，部分1 m分辨率，TIFF格式，分幅索引数据采用SHP格式保存，所有数据的坐标系为WGS1984。局部数据结果如图3。




图3 ? 中巴经济走廊（喀什至伊斯兰堡段）局部正射影像数据集

3 ? 数据质量控制和评估
3.1 ? 正射校正精度评估
本数据集采用基于RPC参数构建有理函数模型的基础上，在正射校正过程中，利用优化后的ASTER GDEM进行第一次纠正，并对不正常变形区域，勾画边界进行掩膜处理，逐一将掩膜区域进行二维DEM可视化编辑，进行高程平滑，最后将修正后的DEM再导入正射校正模块，以此保证正射校正精度。
以研究区盖孜河谷区域样例影像（成像时间2016年）为例，从影像像元误差进行正射校正精度分析。在目标影像和参考影像上均匀选择25个检查点，其中参考影像来源于天地图第16级基础影像，空间分辨率为1.88 m。以参考影像为基准，目标影像检查点对应的最大误差X方向为0.35个像元，Y方向为0.4个像元；均方根误差X方向为0.42个像元，Y方向为0.38个像元。
从地面实测点与目标影像地物目标进行点位绝对定位精度误差分析，实测检查点来源于2017年中巴经济走廊科学考察中利用差分GPS实测的控制点，控制点分布如图4所示。由于研究区地处帕米尔高原与喀喇昆仑山交汇的无人区，实测控制点布设在中巴公路两侧标志明显的路基上，样例区段共7个，以其绝对定位误差与影像中对应标志点进行解算，统计影像标志点的平均平面中误差，结果表明，单点定位平面中误差最大为3.36 m，平均平面中误差为1.18 m。




图4 ? 地面检查点分布

3.2 ? 影像融合评价
从目视角度上定性评价融合效果，Brovey方法颜色失真较大；GS方法清晰度较差；PanSharpening方法和PCA方法整体效果较好，但PanSharpening方法纹理较细腻，色彩、光谱等信息保持较好，高分影像采用PanSharpening方法融合效果较好，如图5所示。

（a）


（b）


（c）


（d）


图5 ? 四种融合方法效果（a）PanSharpening融合；（b）Brovey融合；（c）GS融合；（d）PCA融合

基于融合效果定量评价指标的计算公式，主要针对蓝、绿、红和近红外4个波段，分别对4种融合方法进行定量统计，结果如表2所示。
表2 ? 融合效果评价

融合方法	波段	μ	std	g	dc	cc
Brovey	Blue	0.0561	0.0382	2.4011	0.6231	0.3715	10.6124
	Green	0.0440	0.0300	1.9015	0.9760	0.6630	4.0227
	Red	0.0333	0.0226	4.0919	0.5948	0.4017	5.4216
	NIR	0.0747	0.0589	2.2800	0.3190	0.7173	5.9013
GS	Blue	0.07136	0.0503	1.8762	0.6176	0.4282	2.1786
	Green	0.0940	0.0656	1.803	0.4761	0.3132	3.7913
	Red	0.1203	0.0848	5.8418	0.8074	0.3128	4.5201
	NIR	0.1863	0.1309	4.0825	0.7014	0.5710	5.9514
PCA	Blue	0.0714	0.0509	4.1137	0.1691	0.5179	2.1947
	Green	0.0940	0.0662	3.9371	0.2285	0.7060	4.7489
	Red	0.1204	0.0858	6.4918	0.2271	0.5710	4.4702
	NIR	0.1864	0.1287	3.9272	0.3384	0.6257	3.7513
PanSharpenig	Blue	0.0725	0.0515	4.2794	0.1541	0.6248	3.3985
	Green	0.0955	0.0669	3.8371	0.2080	0.7136	2.9453
	Red	0.1222	0.0867	8.2592	0.3536	0.5981	4.3971
	NIR	0.1892	0.1297	4.1293	0.3227	0.6929	2.5296

Brovey变化的结果是均值、标准差及平均梯度值都较低，影像质量效果较差，影像细节不清晰。虽然其影像信息量和相关系数与其他方法相当，但其偏差指数过高，光谱失真现象较严重。
GS变化的结果是均值、标准差与PCA和PanSharpening方法相当，但平均梯度较低，混合像元较多，地物分辨能力较差，其偏差指数值过高，光谱失真现象较严重。
PCA变化的结果是均值、标准差和评价梯度值较高，影像质量效果较好，虽比PanSharpening方法稍低，但比Brovey和GS方法优势明显。另外，其偏差指数较低，光谱失真较小，信息量保持不丰富。
PanSharpening变化的结果是均值、标准差和平均梯度与PCA方法能力相当，但比Brovey和GS方面优势明显，影像质量较好。光谱偏差指数值整体较小，光谱信息保持失真较小。相关系数值整体偏大，融合信息保持能力强。影像信息熵比其他方法均小，融合效果细腻，信息量丰富。
因此，从实验结果看，PanSharpening方法较适合与本区域高分影像融合处理。

3.3 ? 镶嵌接边评价
其中影像镶嵌采用平均值作为栅格单元值，把多景影像按镶嵌线进行无缝拼接，主要包括计算轮廓线、颜色校正、影像羽化/调和、接边等处理步骤，镶嵌接边精度如表3，影像接边处色彩过渡自然，地物合理接边，无重影现象，纹理清晰，效果如图4所示。
表3 ? 影像镶嵌接边精度评价

序号	类型	内容
1	平均误差/像元	1.24
2	标准偏差/像元	1.13
3	接边处是否清晰	清晰
4	纹理是否清晰	清晰
5	色差差异情况	过渡自然

4 ? 数据价值
中巴经济走廊主要穿越多年冻土、冰川、积雪区，气候环境恶劣，人烟稀少，基础影像数据匮乏，国产米级分辨率影像数据几乎为空白，开展本区域冰川、积雪、冻土等特殊环境要素调查迫切需要高分辨率DOM数据。
本文基于国产高分一号、高分二号遥感数据，利用遥感图像处理软件开展中巴走廊高分DOM数据生产。针对正射纠正中的复杂地形影响，采用RPC有理函数模型进行纠正，针对全色与多光谱影像融合，采用PanSharpening工具，并与Brovey、GS、PCA等方法进行对比分析，最终实现了正射影像半个像元纠正精度，获得纹理清晰、光谱损失较小的正射数据集。
中巴经济走廊正射影像数据是走廊内工程选址、设计和规划的重要基础数据，也可作为全球气候变化背景下，走廊资源调查的本底数据。在正射影像数据基础上，结合气候、水文、生态等数据综合分析，对区域内工程长期安全运营和生态环境可持续发展有着重要意义。

5 ? 数据使用方法和建议
中巴经济走廊（喀什至伊斯兰堡段）正射影像数据集保存为栅格TIF格式。ArcGIS、QGIS、ENVI、ERDAS等常用的GIS与遥感软件可支持该数据的读取和操作。

致 谢
感谢高分辨率对地观测系统甘肃数据与应用中心提供的高分影像数据，感谢美国航空航天局与日本经济产业省提供的ASTER GDEM V2高程数据，感谢美国DigitalGlobe商业卫星与IKONOS商业卫星公司提供参考影像数据。


[1]
习近平《在纳扎尔巴耶夫大学的演讲》[EB/OL]. (2013-09-06) [2018-08-01]. http://news.xinhuanet.com/world/2013-09/08/c_117273079.htm.

+?CSCD?·?Baidu Scholar

[2]
朱颖彦, 杨志全, 廖丽萍. 中巴喀喇昆仑公路冰川地貌地质灾害[J]. 灾害学, 2014, 2014(3): 81-90.

+?CSCD?·?Baidu Scholar

[3]
方成杰, 钱德玲, 徐士彬, 等. 基于可拓学和熵权的中巴公路泥石流易发性评价[J]. 自然灾害学报, 2016, 25(6): 18-26.

+?CSCD?·?Baidu Scholar

[4]
郭万钦, 刘时银, 许君利, 等. 木孜塔格西北坡鱼鳞川冰川跃动遥感监测[J]. 冰川冻土, 2012, 34(4): 765-774.

+?CSCD?·?Baidu Scholar

[5]
GUL C, KANG S C, GHAURI B, et al. Using Landsat images to monitor changes in the snow-covered area of selected glaciers in northern Pakistan[J]. Journal of Mountain Science, 2017, 14(10): 1-15.

+?CSCD?·?Baidu Scholar

[6]
MANTOVANI F, GRACIA F J, DE COSMO P D, et al. A new approach to landslide Geomorphological mapping using the Open Source software in the Olvera area(Cadiz,Spain)[J]. Landslides, 2010, 7: 69-74.

+?CSCD?·?Baidu Scholar

[7]
LOMBARDO L, BACHOFER F, CAMA M, et al. Exploiting Maximum Entropy method and ASTER data for assessing debris flow and debris slide susceptibility for the Giampilieri catchment (north‐eASTERn Sicily, Italy)[J]. Earth Surface Processes & Landforms, 2016, 41(12): 1776-1789.

+?CSCD?·?Baidu Scholar

[8]
MEER F D V D, WERFF H M A V D, RUITENBEEK F J A V, et al. Multi- and hyperspectral geologic remote sensing: A review[J]. International Journal of Applied Earth Observation & Geoinformation, 2012, 14(1):112-128.

+?CSCD?·?Baidu Scholar

[9]
MONTGOMERY D R,Dietrich W.E.Landscape dissection and drainage area-slope thresholds[J]. Process Models & Theoretical Geomorphology, 1994: 221–246.

+?CSCD?·?Baidu Scholar

[10]
RICHTER R. Correction of atmospheric and topographic effects for high spatial resolution satellite imagery[J]. International Journal of Remote Sensing, 1997, 18(5): 1099-1111.

+?CSCD?·?Baidu Scholar

[11]
ROCKWELL, B.W. Automated mapping of mineral groups and green vegetation from Landsat Thematic Mapper imagery with an example from the San Juan Mountains, Colorado: U.S[EB/OL]. (2013-09-06)[2018-08-01]. Geological Survey Scientific Investigations Map 3252, 25 p. pamphlet, 1 map sheet, scale 1:325,000, https://pubs.usgs.gov/sim/3252.

+?CSCD?·?Baidu Scholar

[12]
TACKIKAWA, T., HATO, M., KAKU, M., & IWASAKI, A.?Characteristics of ASTER GDEM version 2[J/OL]. International Geoscience and Remote Sensing Symposium (IGARSS),2011, 3657-3660. DOI:?10.1109/IGARSS.2011.6050017.

+?CSCD?·?Baidu Scholar

[13]
赵国松, 杜耘, 凌峰, 等. ASTER GDEM与SRTM3高程差异影响因素分析[J]. 测绘科学, 2012, 37(4): 167-170.

+?CSCD?·?Baidu Scholar

[14]
刘勇. 中国西北部高山高原地区SRTM3数据质量评价[C].中国地理信息系统协会GIS理论与方法专业委员会2007年学术研讨会暨地理元胞自动机和应用研讨会, 广州, 2007.

+?CSCD?·?Baidu Scholar

[15]
魏春. 基于改进的有理函数模型的卫星影像正射校正方案及实现[D]. 北京: 中国科学院大学, 2015.

+?CSCD?·?Baidu Scholar

[16]
AIAZZI B, ALPARONE L, BARONTI S, et al. Twenty-five years of Pansharpening: a critical review and new developments[J]. Signal and Image Processing for Remote Sensing, 2012: 533-548.

+?CSCD?·?Baidu Scholar

[17]
VIVONE G, ALPARONE L, CHANUSSOT JET et al. A critical comparison among pansharpening algorithms[J]. IEEE Transactions on Geoscience and Remote Sensing, 2015, 53(5): 2565–2586.

+?CSCD?·?Baidu Scholar


数据引用格式
韩立钦, 张耀南, 田德宇,等. 中巴经济走廊（喀什至伊斯兰堡段）高分正射影像数据集[DB/OL].国家特殊环境、特殊功能观测研究台站共享平台, 2018. (2018-03-21). DOI: 10.12072/casnw.048.2018.db.

稿件与作者信息

论文引用格式
韩立钦, 张耀南, 田德宇, 等. 中巴经济走廊（喀什至伊斯兰堡段）高分正射影像数据集[J/OL]. 中国科学数据, 2019, 4(3). (2019-07-11). DOI: 10.11922/csdata.2018.0052.zh.
韩立钦hanliqin

主要承担工作：高分影像处理，地形校正。
（1983—），男，河南范县人，博士研究生，研究方向为定量遥感、遥感地质灾害监测。

张耀南zhangyaonan

主要承担工作：数据处理流程设计。
yaonan@lzb.ac.cn
（1966—），男，甘肃省天水市人，博士，研究员，研究方向为地学大数据。

田德宇tiandeyu

主要承担工作：数据批处理实现。
（1993—），男，内蒙古自治区四子王旗人，硕士研究生，研究方向为寒旱区遥感应用。

康建芳kangjianfang

主要承担工作：数据质量控制与管理。
（1981—），女，甘肃秦安人，硕士，工程师，研究方向为寒旱区大数据应用。