已有元数据的grd和ctl文件，现要用grads做eof分解，但不知道怎么去做

阿朝睁眼了 · 发表于 2022-4-11 11:34:53

已经有处理好的grd数据，且楼主只会使用grads进行绘图，对eof分解了解不多，现在需要利用处理好的grd数据进行eof分解，很急!

阿朝睁眼了 · 发表于 2022-4-11 11:37:34

现在想知道用grads和grd数据能不能进行eof分解，以及去做eof分解的步骤。

雨落森林 · 发表于 2022-4-11 11:54:46

grads肯定不行啊， EOF分解，核心就是求矩阵的特征值和特征向量，至少也要用fortran吧

阿朝睁眼了 · 发表于 2022-4-11 13:37:10

雨落森林发表于 2022-4-11 11:54
grads肯定不行啊， EOF分解，核心就是求矩阵的特征值和特征向量，至少也要用fortran吧

那我接下来是再fortran中用已有的grd数据进行处理吗

阿朝睁眼了 · 发表于 2022-4-11 15:11:41

如何使用Fortran对grd文件数据进行EOF分解？

雨落森林 · 发表于 2022-4-11 21:28:59

阿朝睁眼了发表于 2022-4-11 15:11
如何使用Fortran对grd文件数据进行EOF分解？

读取数据矩阵，然后，求距平，然后，计算协方差矩阵，然后计算矩阵的特征值和特征向量，然后，得到的每一个特征向量就是一个pattern，大概这样。不过不建议用fortran（除非你真的愿意从底层一步一步实现来加深自己对EOF的印象，建议用NCL或者python或者matlab，都很方便。

		自动登录	找回密码
密码			立即注册