海南岛高分一号卫星遥感深加工数据集

江威，何国金，程博，龙腾飞，冷宛春，尹然宇

1.中国科学院遥感与数字地球研究所，北京 100094

2.中国科学院大学，北京 100049

3.海南省地球观测重点实验室，海南三亚 572029

4.三亚中科遥感研究所，海南三亚 572029

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

数据库（集）名称 海南岛高分一号卫星遥感深加工数据集数据作者 江威、何国金、程博、龙腾飞、冷宛春、尹然宇数据通信作者 何国金（hegj@radi.ac.cn）数据时间范围 2015-2016年地理区域 海南全岛（地理范围包括北纬18°7′-20°7′，东经108°30′-111°10′）空间分辨率 2 m数据量 105.2 GB数据格式 *.tif ( GeoTIFF, 8bit float )数据服务系统网址ftp://124.16.184.141/GF1_hainan/ftp://10.7.2.2/GF1_hainan/http://www.sciencedb.cn/dataSet/handle/679基金项目海南省重大科技计划项目（ZDKJ2016021）；国家重点研发计划全球变化及应对专项课题项目（2016YFA0600302）；中国科学院A类战略性先导科技专项课题（XDA190090300）；国家发展和改革委员会促进大数据发展重大工程（2016-999999-65-01-000696-01）。

数据库（集）组成本数据集主要由以下五个部分组成：第一部分是海南岛高分一号卫星遥感镶嵌产品，格式为GeoTIFF，第二部分是海南岛高分一号卫星遥感镶嵌产品缩略图，第三部分是海南岛高分一号卫星遥感镶嵌分幅产品，第四部分是海南岛高分一号卫星遥感镶嵌分幅产品索引图，第五部分是海南岛高分一号深加工原始影像列表。此外，提供了一份样例数据，包含正射校正、影像融合和影像镶嵌数据。

引 言

海南岛是21世纪海上丝绸之路战略的重要支点，也是我国面积最大的经济特区，具有特殊的区位优势和战略地位[1]。新时期海南自由贸易试验区的建设过程中，及时掌握资源环境动态信息对于海南岛资源环境和社会经济可持续发展具有重要意义。

卫星遥感技术对于开展大尺度资源环境调查具有独特优势，这类调查必须依靠高质量的卫星遥感基础数据[2]。何国金等人首次提出了即得即用（Ready To Use，RTU）卫星遥感产品理念[3-4]，该理念是将卫星遥感数据加工成用户能够直接使用的级别，免去用户在使用过程中繁琐的预处理过程，目前已经发布了基于Landsat系列数据研制的中国区域RTU系列产品，包括正射校正产品、地表反射率产品、地表温度产品、镶嵌产品和标准分幅产品。随着国产卫星的发射，江威等人进一步利用国产高分卫星研制了全国16 m分辨率的卫星遥感RTU产品，进一步丰富了原有RTU数据类型[5]。尽管这些产品包含了海南岛范围，但数据的空间分辨率较粗，需要研制覆盖海南岛更高分辨率的卫星遥感深加工产品，以期满足海南岛资源环境精细调查的需求。

本数据集选用2015-2016年的国产高分一号2 m分辨率卫星遥感数据，是较早针对海南岛研制的高分卫星遥感深加工产品。本数据集特点是云量少、定位精度高和色彩均衡，能够直接下载使用。

1 数据采集和处理方法

1.1 数据源

高分一号卫星于2013年4月发射，目前在轨5年多，2018年3年又成功发射了高分一号系统星的02、03、04星，目前已经组网形成星座。高分一号搭载了2台2 m分辨率全色8 m分辨率多光谱相机，4台16 m分辨率多光谱相机，其中2 m分辨率卫星参数如表1[5]。

本数据采用的高分一号2 m分辨率卫星数据来源于海南省高分数据与应用中心，下载的数据产品为 1A级，尚未进行正射校正，数据文件中包含了缩略图、有理函数系数（Rational Polynomial Coefficients，RPC）文件、原始数据以及影像元数据等文件，全色影像和多光谱影像分为2个TIFF格式存储[6]。选择2015-2016年高质量的高分一号2 m分辨率卫星数据，影像的云量在10%以下，优选夏季影像，下载后检查影像，确保没有缺失、条带和异常像元等情况。根据以上要求，共选取覆盖海岛67景高分一号2 m分辨率卫星影像。

除了影像数据外，还收集了海南岛30 m 数字高程模型（DEM）数据和Landsat 8 全色15 m分辨率参考影像，如图1所示，DEM高程数据来源于地理空间数据云[7]，Landsat 8影像来源于中国遥感卫星地面站[8]。

表1 GF-1卫星2 m分辨率传感器主要载荷技术指标

图1 辅助数据

1.2 数据处理

针对多源卫星遥感数据，需要对卫星遥感数据开展数据预处理，主要包含影像配准、平差、几何校正、影像融合和影像镶嵌，流程如图2。以下主要对影像配准、控制点采集、平差、正射校正、影像融合和匀色镶嵌等关键步骤进行说明。

图2 高分一号影像数据处理流程

1.2.1 影像配准

在影像成像过程中，多种误差会造成全色和多光谱位置存在错位现象，需要首先对影像进行配准，以高分辨率全色影像为参考，将多光谱数据进行配准，配准后影像偏差要求不超过1个像元。

1.2.2 控制点采集

利用高分一号影像自带的RPC系数校正影像，影像定位精度较低，需要借助参考影像采集高精度控制点，优化影像自带的RPC系数，从而实现高精度正射校正[9]。传统的地面控制点都是人工选择，主要是道路交叉点等不易变动的位置，本次控制点选取采用频域相位相关匹配法，该方法匹配精度高，而且速度较快，能够识别特征不明显的控制点[5]。地面控制点（GCP）采集搜索范围为100像素×100像素，每景影像控制点数阈值设置为100，通过网格控制，使控制点在空间上分布均匀，在特定缺少GCP区域，进行手动增加。

1.2.3 平差

高分一号影像之间重叠较大，需要对影像采集连接点进行平差，以保证影像具有较高的接边精度[5]。通过对影像进行连接点采集，进一步优化每景影像的RPC系数，连接点采集的算法也为频率域相位相关，连接点搜索半径设置为50像素×50像素。

将各时期土地利用类型面积与社会经济相关指标进行多元回归分析，找出影响山东省土地利用变化的主要因素。通过SPSS软件，对这些因子进行降维分析，构建多元逐步回归分析模型，分析各地类面积变化与社会经济指标的相关性[3]。

1.2.4 正射校正

卫星在采集过程受诸多成像因素影像，影像会发生一定的几何畸变，正射校正的目的是消除这种畸变，使得采集的影像能够准确代表地面实际位置。有理函数模型是目前高分影像校正通用的方法，模型将像元坐标表示以相应地面点空间坐标为自变量的有理多项式的比值[9]，通过引入大量的参数，能够实现较高的定位精度，该模型表达式如下[10]：

式中（Tn，Sn）和（Xn，Yn，Zn）分别为像素坐标（T，S）和地面坐标（X，Y，Z）经平移和缩放后的正则化坐标，进一步地，该方程形式可以展开为：

ai,bi,ci,di(i=0,1…19)即为RPC参数，其中b0和d0的值为1。

由于高分一号影像自带的RPC系数精度不够高，通过前面步骤采集的控制点和连接点可对RPC系数进行优化，并利用DEM高程数据进一步提高复杂地形区域影像定位精度。针对单景影像，首先利用有优化的理函数模型对全色影像进行校正，然后将多光谱影像与已经纠正好的全色影像进行配准，这样能够同时保证全色和多光谱影像定位精度。

1.2.5 影像融合

高分辨率遥感数据的全色图像具有较高的空间分辨率，而多光谱图像具有丰富的光谱特征，对全色和多光谱影像进行融合可获取高空间分辨率的多光谱影像。目前遥感商业软件提供了多种影像融合方法，如Brovey、Gram-Schmidt、HSV、PCA和Pansharp。江威等人从目视效果和定量指数的角度评价了这5种方法用于高分二号影像融合的效果，结果表明Pansharp 法在图像信息、细节以及光谱都具有较好的保持效果[11]。本次产品是采用Pansharp方法进行影像融合。

由于影像成像条件不一致，融合后的影像会有较大的色差，需要对影像进行匀色，使其整体上色调能够保持一致。本次匀色镶嵌采用的是捆绑匀色方法，该方法首先在每个影像及其重叠影像间使用“束捆绑”来计算每个影像的均值和西格玛全局调整，通过该调整使影像之间的色调更加均衡；其次，系统会自动采集一系列特征点，在影像对之间进行小范围的局部调整。该方法的优势是不考虑影像叠加次序，对影像是全局到局部精细化调整。另外，由于部分影像存在少许云量覆盖，可以通过人工编辑镶嵌线编辑，用无云影像替代有云影像，最后生成云量较少的影像产品[5]。

2 数据样本描述

本数据集为海南岛高分一号深加工产品，为全岛镶嵌影像产品，能够直观展示海南全岛资源环境空间分布。影像空间分辨率为2 m，GeoTIFF格式，格式为阿尔伯斯圆锥等面积投影，坐标系为WGS1984，海南全岛镶嵌图如图3所示。

3 数据质量控制和评估

3.1 正射校正精度

影像采集的控制点和连接点精度如表2。全区域共采集控制点4079个，平均每景控制点61个，控制点的均方根误差（RMS）为1.35个像元，每影像影像采集的控制点个数及其均方根误差如图4。此外，全区共采集了连接点数量为27685个（图5），连接点水平均方根误差为0.18个像元，分别选取上下和左右邻接影像展示接边精度（图6），发现影像的接边精度很高，道路和河流线状地物接边精度较高。

图3 海南岛高分一号影像镶嵌图

表2 影像采集的控制点和连接点精度

图4 单景影像控制点数量和控制点均方根误差

图5 影像连接点分布图

图6 影像接边图

3.2 影像融合质量

影像融合后会得到分辨率较高的多光谱影像，图7显示了融合前后的影像效果，从图中可以看出，融合后的影像相比融合前的多光谱影像，对于地物细节信息有较大的提升，同时对于地物的色彩也保持得较好，说明本文采用的融合方法能够取得较高的影像融合质量。

图7 影像融合前后效果对比

3.3 影像的匀色效果

影像的匀色镶嵌效果是对于不同色度的影像进行颜色调整，从而在全区范围优化，使得镶嵌后的影像在接边处颜色过度自然。图8显示了镶嵌前后影像的效果，可以看出镶嵌后的影像在影像接边过度上有较好的处理效果，影像的色调基本均匀。

图8 影像匀色镶嵌前后效果对比

4 数据价值

高质量的遥感卫星数据产品对于海南岛资源环境动态信息监测极为重要。本研究收集了国产高分一号2 m分辨率卫星遥感数据，采用卫星遥感数据深加工流程技术，研制覆盖海南全岛的卫星遥感深加工产品。本产品云量少，定位精度高且色彩均衡，使用者可直接进行遥感应用。本数据集包含全岛镶嵌影像、镶嵌产品缩略图、镶嵌影像分幅产品和镶嵌分幅产品索引图等，可方便不同用户下载使用。

海南省既是国家生态文明试验区，也将建设成为自由贸易试验区。本数据集可作为海南岛基础数据，将在岛内重大工程选址、城市发展规划、农田保护和海岸带监测等领域发挥重要作用。

5 数据使用方法和建议

海南岛高分一号2m分辨率卫星遥感深加工数据集保存为GeoTIFF格式，包含4个波段，能够在ENVI、PCI、ArcGIS和QGIS等遥感和GIS主流软件中读取和操作。为便于下载使用，对镶嵌影像进行分块处理，并提供了影像的分块边界（图9），可以快速查询获取局部数据。与此同时，在样例数据中提供了海南省全岛60 m分辨率镶嵌产品，便于快速查看全岛数据。所采用的单景高分一号数据列表命名为海南岛高分一号深加工数据列表，格式为XLS。使用过程中，建议采用影像3（R）4（G）1（B）波段组合，能够直观反映植被和城市等信息。

图9 海南岛镶嵌影像分块示意图