本文将介绍 Whittaker Smoothing 时间序列数据平滑算法。Whittaker Smoothing 对于一个给定的时空图像集合,它可以生成平滑的时间序列图像。
具体而言,这些函数实现了以下过程:
- extractBits函数将从一个 QA 图像中提取一些特定的比特位,返回一个新的单波段图像,其中每个像素的值代表了提取出的比特位。
- getDifferenceMatrix函数生成一个给定阶数的差分矩阵,可以用于Whittaker Smoothing。
- unpack函数将数组图像展开为单个图像并分配波段名。
- inverseLogRatio函数计算回归结果的逆对数比率,以将其转换回百分比单位。
- whittakerSmoothing函数使用 Whittaker Smoothing 对时间序列图像集合进行平滑处理,并返回一个新的时间序列图像集合。在这个函数中,如果 isCompositional 为真,则将数据视为组成数据并对其进行预处理。然后计算给定阶数的差分矩阵并进行矩阵运算,生成一个平滑系数图像并解决平滑系数图像和原始数组图像之间的线性系统,最后重组回平滑图像集合,并在必要时将其还原为组成数据。
结果如下图所示,