2015, Vol. 40文章信息
- 许剑辉, 舒红
- XU Jianhui, SHU Hong
- 基于Triple-Collocation的地面观测与卫星遥感数据融合的雪深反演
- The Triple-Collocation-based Fusion of In-situ and Satellite Remote Sensing Data for Snow Depth Retrieval
- 武汉大学学报·信息科学版, 2015, 40(4): 469-473
- Geomatics and Information Science of Wuhan University, 2015, 40(4): 469-473
- http://dx.doi.org/10.13203/j.whugis20130727
-
文章历史
- 收稿日期：2013-12-02
积雪是最动态和活跃的水文学变量之一[1]，也是地球上最重要的水资源补给方式之一。雪深是积雪监测的一个重要因素，准确的雪深估计对气候、雪灾的监测非常重要。一般情况下，被动微波遥感反演的雪深具有空间分辨率低、反演精度低等问题。
针对上述问题，国内外学者提出了多种数据定量反演雪深[2, 3]，尤其是在利用多源数据融合的方法反演积雪数据的研究方面，已经取得了重要的进展[4, 5]。赵亮等结合地面观测数据和遥感数据，采用时空距离权重法计算动态反演系数，最后通过反演系数进行雪深反演[1]。Foster等发展了一种新的ANSA积雪算法融合MODIS雪盖数据、AMSR-E和QuikSCAT数据改善积雪产品[6]。Liu等根据格点与地面观测站点距离，通过计算每个站点的最优权重的方式对ANSA遥感反演的格点雪深数据进行偏差校正[7]。这种基于最优估计理论的偏差校正方法适用于地面观测数据比较密集的区域，并不适合于地面站点稀疏的区域。
本文以北疆为研究区域，北疆气象站点雪深、 AMSR-E和SSM/I遥感反演雪深为研究对象，利用一种不确定性估计的方法——Triple-Collocation方法分别估计三个雪深数据集的误差方差，在最小二乘原理基础上实现星-地多源雪深数据的融合分析。 1 基于Triple-Collocation的融合方法 1.1 Triple-Collocation方法
Triple-Collocation(TC)方法由Stoffelen提出，用于校正基于Scattermeter驱动的海洋风速，估计其误差[8]。随后，该方法作为一个不确定估计方法被用来评价土壤湿度[9, 10]、叶面积指数[11]等遥感数据的误差和不确定性。TC方法的实现需要三个相互独立的数据集。为了避免误差估计过程中的数值问题，一般认为每个独立数据集的样本数大于100[9]。假设三个相互独立的雪深(基于站点观测的普通克里金(ordinary kriging,OK)估计雪深、AMSR-E和SSM/I雪深)和假设的雪深真值间存在如下线性关系：
式中，sdptrue为真实的雪深；sdpx、sdpy和sdpz分别为OK(01,OK)估计雪深、AMSR-E和SSM/I遥感反演雪深；εx、εy和εz分别为OK、AMSR-E和SSM/I雪深估值的误差，均值为0。通过估算式(1)的均方根误差来表达每个雪深数据的误差方差。
方程两边除以βi (i=x,y,z)，sdpi*=(sdpi-αi)/βi，εi*=εi/βi并移项，整理得到：
通过式(2)中三个等式两两相减消除未知的真值sdptrue，然后进行交叉相乘并取均值，最后得到每个估计的雪深数据误差的方差。由于残差εx*、εy*和εz* 相互独立，有<εx*εy*>=<εx*εz*>=<εy*εz*>=0。
式中，<·>表示均值计算。为了估计出不同雪深数据误差的方差，需要计算出式(1)中相应的回归系数αi和βi(i =x,y,z)。由于难以精确获取真实的雪深数据sdptrue，任意选择其中一个数据集sdpOK作为参考数据，并设置αx=0，βx=1。考虑到式(1)的对称性，εx*2、εy*2和εz*2的计算与参考数据的选择无关。式(1)调整为：
对式(4)中的等式两边分别计算其均值(即计算区域雪深的空间均值)，然后将式(4)减去其均值，并且令sdp′=sdp-<sdp>，将得到的三个等式两两相乘并取均值，最终计算出相应的回归系数：
最后将计算得到的系数代入式(3)，即可得到不同雪深数据集的误差方差：
将普通克里金估计雪深、AMSR-E与SSM/I反演雪深产品进行融合，得到一个数据质量更好的雪深反演数据[12]：
式中，s pm为融合的雪深估计值；sdpx、sdpy和sdpz分别为OK估计的雪深，AMSR-E反演雪深和SSM/I反演雪深；ωx、ωy和ωz分别为OK、AMSR-E和SSM/I雪深的相对权重。为了保证融合雪深估计值是无偏的，需要满足条件ωx+ωy+ωz=1。最后根据最小二乘原理求出式(6)相应的权重系数：
以北疆为研究区域，区域空间位置为42°12′N~49°12′N，79°48′E~92°36′E，包括两大山脉(天山和阿尔泰山)和一大盆地(准格尔盆地)。北疆地区地形复杂，地貌多样，海拔最高达到6 300 m，最低仅为170 m。
研究所选用的遥感雪深为中国西部环境与生态科学数据中心提供的逐日的积雪厚度分布数据。用于反演该雪深数据集的原始数据来自美国国家雪冰数据中心处理的两个不同传感器AMSR-E和SSM/I的逐日亮温。通过对不同传感器的亮温进行交叉定标，提高亮温在时间上的一致性，然后利用文献[13]在Chang算法基础上针对中国地区进行修正的算法反演得到经纬度投影的遥感雪深。实验进行多源数据融合前，采用最邻近插值法对空间分辨率为25 km的雪深进行插值，得到0.01°空间分辨率的雪深。
实验采用同一天的站点雪深(北疆48个站点雪深日值数据集)与遥感反演雪深进行融合分析。站点雪深是对积雪过程累积量的一个测量值，和遥感瞬时观测值在天尺度下具有可比性。由于北疆地区站点稀少，为了获取较高精度的OK插值结果，选择10%的空间均匀分布的站点(5个站点：布尔津、塔城、奇台、伊宁和乌鲁木齐)的雪深数据作为验证数据，最后90%的站点(43个站点)的雪深进行普通OK插值分析，得到0.01°空间分辨率的雪深。整个北疆地区获取的0.01°空间分辨率的OK、AMSR-E和SSM/I雪深分别有450、259个像元。采用2007年1月1日(001)、1月15日(015)和1月30日(030)OK估计雪深、AMSR-E和SSM/I雪深进行融合分析。 3 结果及讨论 3.1 雪深数据间的相关系数及权重值
由于站点观测雪深、AMSR-E和SSM/I雪深的来源是相互独立的，满足TC方法需要三个数据集相互独立。本文主要采用了001、015和030共3 d OK估计雪深、AMSR-E和SSM/I雪深进行多源雪深数据的融合分析。首先根据公式sdp′=sdp-
在最小二乘原理的框架下，结合基于TC估计的误差方差，利用式(7)分别计算三个雪深数据的权重。表 2显示了不同时刻三个雪深数据所占的权重，其中AMSR-E雪深占的权重最大。OK估计雪深的权重比较大，权重最小的是SSM/I雪深。从表 1和表 2可以看出，和其他两个雪深数据的相关系数最大的雪深数据，其占的权重也最大。
图 1为星-地多源雪深的融合结果。一般情况下，被动微波不能探测到5 cm以下的浅层积雪，由于AMSR-E和SSM/I雪深反演时建立的反演模型存在误差，使得反演模型计算得到的积雪厚度出现0或者小于5 cm的情况，如AMSR-E遥感反演雪深数据在北疆中部出现一小块雪深为0的区域，AMSR-E和SSM/I反演结果都有小于5 cm的样本。本实验采用的车涛[13]反演的AMSR-E和SSM/I雪深产品保留了雪深为0和小于5 cm的情况。OK估计的雪深整体上都比较平滑，因为地统计插值过程中没有考虑其他因素如高程、气温等的影响，导致OK估计值在有些情况并不能很好地反映区域情况。特别是在北疆中部沙漠地区，几乎没有地面雪深观测站点，地统计方法会降低该区域雪深的精度，而这种情况下卫星遥感数据能够较好地反演该区域的雪深情况。AMSR-E和SSM/I遥感反演的雪深空间分辨率低，通过插值方法和星-地多源雪深数据的融合方法能够获取高空间分辨率的、高精度的雪深数据，如图 1所示。融合的雪深较好地保持了原始遥感数据的整体空间格局，同时通过融入高精度地面观测的雪深数据信息而有效地提高雪深数据的测量精度，如图 1中TC所示。从图 1中还发现，融合的雪深比AMSR-E和SSM/I雪深整体上细节不明显，这是由于OK估计雪深整体上比较平滑引起的。
从图 1还可以看出，2007年1月份属于积雪累积期，最大积雪深度主要集中在阿尔泰山、塔城、伊犁和乌鲁木齐地区。
3.3 融合结果的验证
图 2(a)~2(c)分别为43个站点对应的AMSR-E、SSM/I和基于TC融合的雪深的散点图。地面观测站点雪深与对应AMSR-E雪深的相关系数仅为0.307(001)、0.424(015)和0.314(030)，与SSM/I雪深的相关系数仅为0.312、0.377和0.298，而与融合后的TC雪深数据的相关系数分别提高到0.616、0.780和0.753。
图 3(a)~3(c)为5个站点对应的反演雪深 与站点雪深的散点图。融合后的雪深与站点雪深的相关系数分别为0.307(图 3(a))、0.568(图 3(b))和0.428(图 3(c))，与AMSR-E、SSM/I和站点雪深的相关系数相比，相关系数提高最小的为1月15日，分别提高了67.98%和69.09%。另外，与AMSR-E和SSM/I雪深数据相比较，基于TC的雪深与站点雪深的均方根误差分别减少了10.42%和27.88%(001)、19.72%和30.09% (015)、24.24%和37.63%(030)。
针对地面观测站点雪深、AMSR-E微波遥感反演雪深和SSM/I微波遥感反演雪深数据，提出了一种基于误差分布加权计算的多源数据融合方法，能有效地提高雪深的反演精度，融合的雪深有效地反映了积雪时空分布的变化。
然而，本文只对同一种类型的不同来源的数据进行空间上的融合分析研究，研究还存在一些不足，在后续的研究中，将针对不同类型的积雪数据(如雪盖、雪深)和更多来源的遥感数据，以雪深反演为目标，进行多源积雪数据的时空融合分析研究，获取高精度的雪深时空分布图。
致谢：感谢中国气象局乌鲁木齐沙漠气象研究所提供站点雪深数据，中国西部环境与生态科学数据中心(http://westdc.westgis.ac.cn)提供AMSR-E和SSM/I微波遥感反演雪深数据。
ROK-AMSR-E ROK-SSM/I RAMSR-E-SSM/I
001 0.160* 0.127* 0.713*
015 0.245* 0.204* 0.598*
030 0.179* 0.133* 0.576* 注：*p<0.000 1。 ωOK ωAMSR-E ωSSM/I
001 0.223 0.625 0.152
015 0.327 0.418 0.255
030 0.349 0.467 0.184 
图 1 不同雪深数据的融合结果Fig. 1 Merged Map of TC from Different Snow Depth Products 
图 2 AMSR-E、SSM/I和TC雪深数据的比较(43个站点)Fig. 2 Comparison of AMSR-E,SSM/I and TC-based Snow Depth Products(43 Stations) 
图 3 AMSR-E、SSM/I和TC雪深数据的比较(5个站点)Fig. 3 Comparison of AMSR-E,SSM/I and TC-based Snow Depth Products(Five Stations)
| [1] | Zhao Liang, Zhu Yuxiang, Yang Hong, et al. A Dynamic Approach to Retrieving Snow Depth Based on the Technology of Integrating Satellite Remote Sensing and in Situ Data[J]. Acta Meteorologica Sinica, 2013, 71(4):769-782 (赵亮, 朱玉祥, 杨弘, 等. 一种基于卫星遥感与地面测站数据融合技术的雪深动态反演方法[J]. 气象学报, 2013, 71(4):769-782) |
| [2] | Liu Hai, Chen Xiaoling, Song Zhen, et al. Study of Characteristic Parametric Selection and Model Construction for Snow Depth Retrieval from MODIS Image[J].Geomatics and Information Science of Wuhan University, 2011, 36(1):113-116 (刘海, 陈晓玲, 宋珍, 等. MODIS影像雪深遥感反演特征参数选择与模型研究[J]. 武汉大学学报·信息科学版, 2011, 36(1):113-116) |
| [3] | Liu Yan, Ruan Huihua, Zhang Pu, et al. Kriging Interpolation of Snow Depth at the North of Tianshan Mountains Assisted by MODIS Data[J]. Geomatics and Information Science of Wuhan University, 2012, 37(4):403-405(刘艳, 阮惠华, 张璞, 等. 利用MODIS数据研究天山北麓Kriging雪深插值[J]. 武汉大学学报·信息科学版, 2012, 37(4):403-405) |
| [4] | Foster J L, Hall D K, Eylander J, et al. Blended Visible, Passive Microwave and Scatterometer Global Snow Products[C]. The 64th Eastern Snow Conference, New Foundland, Canada, 2007 |
| [5] | Hall D K, Foster J L, Kumar S, et al. Improving the Accuracy of the AFWA-NASA (ANSA) Blended Snow-cover Product over the Lower Great Lakes Region[J]. Hydrology and Earth System Sciences Discussions, 2012, 9(1):1 141-1 161 |
| [6] | Foster J L, Hall D K, Eylander J B, et al. A Blended Global Snow Product Using Visible, Passive Microwave and Scatterometer Satellite Data[J]. International Journal of Remote Sensing, 2011, 32(5):1 371-1 395 |
| [7] | Liu Y, Peters-Lidard C D, Kumar S, et al. Assimilating Satellite-based Snow Depth and Snow Cover Products for Improving Snow Predictions in Alaska[J]. Advances in Water Resources, 2013, 54:208-227 |
| [8] | Stoffelen A. Toward the True Near-surface Wind Speed:Error Modeling and Calibration Using Triple-Collocation[J]. Journal of Geophysical Research:Oceans (1978-2012), 1998, 103(C4):7 755-7 766 |
| [9] | Scipal K, Holmes T, De Jeu R, et al. A Possible Solution for the Problem of Estimating the Error Structure of Global Soil Moisture Data Sets[J]. Geophysical Research Letters, 2008, 35(24):L24403 |
| [10] | Leroux D J, Kerr Y H, Richaume P, et al. SpatialDistribution and Possible Sources of SMOS Errors at the Global Scale[J]. Remote Sensing of Environment, 2013, 133:240-250 |
| [11] | Fang H, Wei S, Jiang C, et al. Theoretical Uncertainty Analysis of Global MODIS, CYCLOPES, and GLOBCARBON LAI Products Using a Triple-Collocation Method[J]. Remote Sensing of Environment, 2012, 124:610-621 |
| [12] | Yilmaz M T, Crow W T, Anderson M C, et al. An Objective Methodology for Merging Satellite-and Model-based Soil Moisture Products[J]. Water Resources Research, 2012, 48(11):W11502, doi:10.1029/2011WR011682 |
| [13] | Che Tao. Study on Passive Microwave Remote Sensing of Snow and Snow Data Assimilation Method[D]. Lanzhou:Cold and Arid Regions Environmental and Engineering Research Institute, CAS, 2006 (车涛. 积雪被动微波遥感反演与积雪数据同化方法研究[D].兰州:中国科学院寒区旱区环境与工程研究所, 2006) |