REOF 小结

wuyejincao

在写这篇帖子之前，我已经被REOF折磨了148个小时，所以别怪我附件要卖点贡献了！

---

我一共实验了7个不同的REOF程序，最后实验表明，我们论坛里面的REOF程序最全面，最好用。下面我简单说一下实验结论：
1.李建平老师的那个reof肯定是有问题，前面eof结果是对的，后面的reof部分还是有待改进。
2.吴洪宝老师的那个程序总觉得乱七八糟的，特别是那个write（6，*），由于当初run一下窗口就闪没了，故舍弃~
3.论坛里的这个REOF是一个叫海燕的选手分享的，也有些乱，我整理了下，把程序改的更简单易懂一点。

这是7个不同版本的REOF程序： 7个版本的REOF.rar (35.76 KB, 下载次数: 639, 售价: 2 贡献)

2012-10-17 21:33 上传

点击文件名下载附件

---

接下来是我run通的一个REOF程序，里面包含了我的测试数据文件等等整个过程所有需要的文件，方便大家测试。 test.rar (722.61 KB, 下载次数: 990, 售价: 2 贡献)

2012-10-17 21:51 上传

点击文件名下载附件
简单介绍下流程：
1.首先用spring.f90对我的资料进行了预处理，这一步是为了剔除缺测（实际上我是用插值的方法），没有缺测的可以跳过这一步。
2.运行reof_luntan.f。
3.REOF的载荷向量、经纬度存进了tzxl.txt这个文件中。
4.接下来就是清风大侠教导我们的那样，利用站点数据画图了。我的测试数据是226个站（有一个站有些问题，经纬度是0,0，没有关系，别管它），取10个载荷向量。一时之间搜不到清风的帖子，所以我手动传一下清风大侠的这篇神作： GrADS站点文件作图详细解决方案.pdf (1.35 MB, 下载次数: 1085)

2012-10-17 21:45 上传

点击文件名下载附件

---

FAQ：
1.Q:很多人run这个程序的时候，会出现——math的错误。
A:个人觉得这是因为数据中有缺测，所以你先把数据进行预处理吧，我的预处理方法是对undef进行插值（反正最后画图时还是要插值的！）。
2.Q:旋转经验正交函数个数p怎么取？（这里借用@wode_q_x的回答！）
A:1、由经验正交函数的累积方差贡献来确定。一般可取累积方差贡献率达85%为标准来确定旋转特征向量的个数p。方差贡献百分率根据具体问题适当增减。
2、通过特征值对数曲线变化来确定旋转特征向量的个数。如果某个特征值之后的直线的斜率明显变小，即以该点特征值的个数作为旋转特征向量的个数p。
3、用North等（1982）特征值误差范围来确定旋转特征向量的个数p。

---

感谢这段时间来帮助我做REOF的人，感谢论坛给予我的宝贵信息，特别感谢@mofangbao@言深深@麻花先生
求各种表扬~@管理员@版主@超级版主@实习版主   

PS：后来才发现，很多程序是可以run通的，当时没有run通是因为没把缺测处理好。但lijianping老师的那个REOF是肯定有问题，至今没有run通。

yy456

楼主测试数据里面用的数据是春季的什么数据？spring226文件里面原始数据时怎么摆放的？

dr-liu

感谢楼主分享，原来想用matlab做reof分区，但是资源很少。fortran一点不会，慢慢学吧

youjianlaodao

我实验了一下，真是不错的程序哦，但是怎么没有看到旋转后的特征值和解释方差啊？
还有
Sum of variance of squared element of each colum of rotated loding martix
at each circulation
    0.189   0.371   0.399   0.403   0.406   0.408

     Sum of squared element of each of the first NP rotated loading vector
      3.39932     2.99185     1.95800     1.86290     1.58330     1.32448     1.37368     1.65991     1.56700     2.64141

这两个是什么a ?

wuyejincao

兰溪之水 发表于 2012-10-17 23:08

                               
登录/注册后可看大图

write（6，*）和write(*,*)效果一样，1，2，5，6有默认的通道，话不要说得太满。。。

兰溪说的对，我忘了6是默认的屏幕输出了，我又把吴洪宝老师的程序试验了下，结果跟论坛那个REOF结果一致，我贴出err这个文件，以示结果，感谢@兰溪之水。顺便把改过的吴老师的程序再贴一下，不卖贡献了： reof_wu.f (12.47 KB, 下载次数: 624)

2012-10-18 14:35 上传

点击文件名下载附件
EOF EXPAND RESULTS:


  Eigenvalue(lambad)
      43.1394     32.6709     23.8412     18.6576     17.6527     15.2149     13.1467     11.0495      9.7548      8.9113


The percentage of lambad
      19.0882     14.4562     10.5492      8.2556      7.8109      6.7323      5.8171      4.8892      4.3163      3.9430


The accumulated percentage of lambad
      19.0882     33.5444     44.0936     52.3492     60.1601     66.8924     72.7095     77.5986     81.9149     85.8580

REOF EXPAND RESULTS:

     Sum of squared element of each of the first NP rotated loading vector
      33.7379     31.2132     21.2260     25.3954     16.4844     14.4596     17.5117     12.5635     10.0760     11.3714
REOF PERCENTAGE
      14.9283     13.8112      9.3921     11.2369      7.2940      6.3980      7.7485      5.5591      4.4584      5.0316
REOF accumulated percentage
      14.9283     28.7395     38.1315     49.3684     56.6624     63.0604     70.8090     76.3680     80.8265     85.8580


The correlation coefficients martix between unrotated  PC and rotated PC
   0.78  0.42  0.05 -0.44 -0.09 -0.02  0.10 -0.02 -0.03 -0.04
  -0.15  0.77  0.13  0.45 -0.09  0.01 -0.26 -0.23 -0.03  0.19
   0.03 -0.27  0.81 -0.15 -0.26  0.17 -0.32 -0.11  0.03  0.21
   0.58 -0.31 -0.05  0.63  0.13 -0.19 -0.32  0.07  0.08  0.01
   0.04  0.04  0.42  0.12  0.72  0.02  0.42 -0.31  0.06 -0.15
   0.15 -0.08 -0.09  0.32 -0.34  0.73  0.42 -0.13  0.11 -0.01
  -0.02  0.02  0.27  0.24 -0.40 -0.54  0.58  0.25  0.07  0.10
  -0.04  0.17  0.26  0.13  0.05  0.23 -0.09  0.66 -0.26 -0.57
  -0.06  0.12  0.02 -0.07  0.00  0.01 -0.12  0.15  0.95 -0.20
   0.03  0.08  0.00 -0.03  0.33  0.22  0.08  0.55  0.06  0.73

wuyejincao

做REOF的人确实不多，接下来我把分区做好了再跟大家分享利用REOF做分区图的经验~

兰溪之水

write（6，*）和write(*,*)效果一样，1，2，5，6有默认的通道，话不要说得太满。。。

river

不错，没有白被折磨

wode_q_x

恭喜恭喜！

topmad

很全面，测试了各种程序，连测试数据也带上了，很有心，值得学习

mofangbao

总结了以后思路应该更清晰了，很不错的总结~
另外，你说的有些不行的程序有一部分可能到你手上的时候被前面的师兄师姐修改过了，所以。。。

做个霸气的木头

历经坎坷，楼主终于出结果了，恭喜恭喜，还这么细致的奉献出来了，辛苦辛苦~

蔚蓝情怀

实验做的很卖力