GF6 WFV卫星视角影响、与Sentinel2协同及精度分析

GF6 WFV卫星视角影响、与Sentinel2协同及精度分析

李国春

摘要：Sentinel2和GF6 WFV都是共享免费的重要数据源，定量化处理精度是人们广泛关心的问题。本文用同地区同日时间相差7分钟的两个Sentinel2和GF6 WFV数据集进行处理分析比较，认为使用RSD平台经本文流程处理后，两种数据相同（近）波段表现一致，可以协同交互混合使用。注意到并提出了GF6 WFV大视角扫描对数据的影响及实现了在RSD平台的订正处理。

别以为我一本正的摘要是要写文章发表呢，没那事儿。和我原来写的东西一样，这不是空谈，是已经代码实现的RSD功能，数据和软件都是现成的，点点鼠标就可以验证。不是说好的：Talk is cheap，show me the code? 好，言归正传。

GF6 WFV除了有庞大的数据体量外，在几何位置关系上也有一些自己的特点。主要表现在宽扫描刈幅介于中分卫星和高分卫星之间，扫描数据兼有一些中分卫星数据特征，处理数据时需要进行特别关注。

一个GF6 WFV宽视角会带来下面几个问题：

一. 太阳天顶角 通常高分数据可以只使用一个平均太阳天顶角就可以，这是因为视场角比较小，星下点和边缘太阳天顶角差异带来的误差小到可以忽略。中分数据视场角很大（例如±55.4°），多提供逐像元（或者块）的太阳天顶角数据。GF6 WFV 扫描跨度达7~8个纬距（约±32.2°），因此需要逐像元计算太阳天顶角而不能使用平均太阳天顶角。可以根据时间和经纬度位置计算逐像元太阳天顶角和方位角（或者为节省时间计算部分点然后插值求得）。

二. 反射率定标系数 一些处理技术（软件）直接在算得反射率定标系数除以太阳天顶角的余弦，这样在使用线性式进行反射率定标时得到的数据直接是天顶方向的反射率，这里记为TOAup。不能使用同一个太阳天顶角时不除以这个太阳天顶角的余弦，算得的反射率时到达传感器的，记为TOAsensor。中分卫星数据以及Landsat、Sentinel2等常用共享免费数据反射波段定标后反射率均为TOAsensor。GF6 WFV同样，在RSD中直接使用反射率定标系数算得的结果是到达传感器的。如果希望得到天顶方向的反射率，不但需要太阳角度订正，还需要卫星天顶角方向的订正。

三. 卫星天顶角 图1是卫星探测角度的示意图，α是我们想求的卫星天顶角。645km的卫星高度有811km的扫描宽度，视场角这么大不能再使用天底接近0°的α角。GF6 WFV加载切分后可能会丢失像元卫星天顶角信息，所以RSD在加载时就进行了卫星天顶角修正。

图1 GF1 WFV 的卫星天顶角示意图

卫星天顶角引起的数据差异有多大呢？我们用下面给出的数据集在左侧裁剪除一块，进行卫星视角订正和不订正的比较，图2为剪切出来的两块数据，一块进行了卫星视角订正（层1），另一块保持原数据不变（层2）。

图2 GF1 WFV 卫星视角订正比较区域

星下点轨迹附近卫星视角订正前后没有差异，因为是接近垂直向下探测卫星天顶角近似为0。图2是扫描带边缘的，卫星视角订正前后差异很大，图3是卫星视角订正前后的比较。

图3 GF6 WFV 扫描带边缘1、4波段卫星视角订正前后比较

图3a是1波段蓝光的比较结果，图3b是4波段比较结果。两种均值差异均达19.8%。

从星下点的无差异逐渐到边缘的20%左右的差异，说明GF6 WFV不进行卫星视角订正可能会造成很大的处理误差。问题是这个卫星视角不像众所周知的太阳订正，还没有引起广泛重视，至少在高分数据定量处理是这样（侧摆也多是几何应用）。

四. 太阳反辉区 虽然有精心设计的准太阳同步轨道，这样大范围扫描有时还是无法避免太阳反辉区。RSD有菜单能够识别太阳反辉区并加以修正，还有菜单项计算太阳反射（入射）矢量与卫星矢量的夹角（这个夹角小到某值可认为像元落入太阳反辉区），但是这些是为中分数据特别设计的，GF6 WFV暂时还没有，稍等一段我抽出空来给GF6 WFV也加上。暂时先这样用着，中高纬问题不大，低纬地区小心绕开使用。

下面用一个例子来说明GF6 WFV 的处理过程，并且与一个同时的Seninel2数据进行比较。

一、 准备数据。

GF6_WFV_E124.0_N40.2_20190529_L1A1119885324.tar.gz

S2A_MSIL1C_20190529T023551_N0207_R089_T51TVF_20190529T042555.zip

上述GF6 WFV和哨兵2两种数据为同一天相同区域，时间相差7分钟。GF6 WFV数据区域完全覆盖哨兵2数据区域。

上述数据可以从网盘 百度网盘-链接不存在 下载。

二、 确定数据区域

打开哨兵2 L1C 10m数据，处理区域以这个数据范围为界。一个问题是哨兵2 数据是10m，GF6 WFV是16m，我们统一按照10m处理，这样便于比较。打开Sentinel2 L1C数据的MTD_MSIL1C.xml文件，选择10m数据，结果见图4。

图4  Sentinel2 L1C数据范围