Triple Collocation(TC)方法可用于估计同一地球物理变量的不同数据集中的随机误差方差,特别在卫星反演产品(反射率、植被指数、土壤水分等)星地校验中应用广泛。
TC方法原理参见相关论文,这里给出Matlab代码实现。
1. 生成数据集
我们可以基于定义的误差模型生成三个时间序列的数据集,
其中,Θ是真值。
%% 1-data generation
2. 归一化数据集
以任一数据集的均值和标准差为参考(例如以x),将三个数据集进行归一化。
%% 将y和z归一化到x(y和z的均值和标准差一样)
3. 求解不同数据集的误差
%% 求解误差
三个数据集的误差分别为:0.0167、0.0670、 0.0425和设置的误差0.02、0.07、0.04较为接近。
方法2
% In approach 2 we can estimate the triple collocation errors, the scaling parameter β and the signal to noise ratio directly from the covariances of the dataset. For a general overview and how approaches 1 and 2 are related please see [Gruber_2015].
%% Estimation of the error variances from the covariances of the datasets