利用馮偉老師的工具箱，我們可以提取得到以等效水高形式表示的全球質(zhì)量變化格網(wǎng)數(shù)據(jù)。在前一篇專欄，我提供了GRACE改進的數(shù)據(jù)處理和分析的工具箱，注意：我最近修改了部分代碼，具體內(nèi)容見評論。

一、提取時間序列

在【04-馮偉老師GRACE工具箱改進版-可讀取RL06數(shù)據(jù)】運行下面的程序后：

在變量區(qū)存在以下的一些變量：

下面的代碼是提取長江流域2002-2017年163個月的時間序列，yangtze.bln在馮偉老師的數(shù)據(jù)包里面。

運行結(jié)果：

二、進行振幅相位、趨勢的提取

在完成時間序列提取后，進行全球的振幅、相位和趨勢的計算，具體的函數(shù)調(diào)用如下：

其中的time為時間變量：

tt為缺失的月份：

運行以下的代碼，得到

運行結(jié)果【質(zhì)量變化的線性趨勢】

【周年振幅】

【周年相位】