matlab 归一化_Triple Collocation方法Matlab实现

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