05-GRACE工具箱提取流域的质量变化时间序列，并分析周年、半周年振幅、相位、趋势

利用冯伟老师的工具箱，我们可以提取得到以等效水高形式表示的全球质量变化格网数据。在前一篇专栏，我提供了GRACE改进的数据处理和分析的工具箱，注意：我最近修改了部分代码，具体内容见评论。

一、提取时间序列

在【04-冯伟老师GRACE工具箱改进版-可读取RL06数据】运行下面的程序后：

在变量区存在以下的一些变量：

下面的代码是提取长江流域2002-2017年163个月的时间序列，yangtze.bln在冯伟老师的数据包里面。

运行结果：

二、进行振幅相位、趋势的提取

在完成时间序列提取后，进行全球的振幅、相位和趋势的计算，具体的函数调用如下：

其中的time为时间变量：

tt为缺失的月份：

运行以下的代码，得到

运行结果【质量变化的线性趋势】

【周年振幅】

【周年相位】