分享SEOF程序

wangji_zl

这几天正在眼机会这个，雪中送炭

jjpsl

谢谢LZ分享~~~

denghaijun

感谢楼主分享这么好的程序，非常感谢，不过原程序：test-SEOF.m中的“weight=sqrt(cos(lat(b)*pi/180))”有误，应该为：weight=sqrt(cos(lat*pi/180))。不知我的理解是否正确，麻烦楼主检查。

Aires

denghaijun 发表于 2015-6-15 22:45
感谢楼主分享这么好的程序，非常感谢，不过原程序：test-SEOF.m中的“weight=sqrt(cos(lat(b)*pi/180))”有 ...

I think the code is no problem because I have selected the SST domain, so the weighted should be performed only the latitudes of SST domain.

denghaijun

Aires 发表于 2015-6-16 19:05
I think the code is no problem because I have selected the SST domain, so the weighted should be p ...

楼主，程序中提示了“weighted for each latitude of sst, weight = sqrt(cos(lat))”，

1. ......
   
2. a=lon;
   
3. b=lat;
   
4. lon_num=length(a);
   
5. lat_num=length(b);
   
6. sst=xlsread('D:precip.xlsx');
   
7. sst=double(sst);
   
8. sst(sst<-1000)=NaN;
   
9. 
   
10. % weighted for each latitude of sst, weight = sqrt(cos(lat))
    
11. weight=sqrt(cos(lat(b)*pi/180));
    
12. ......

复制代码

其中，weight=sqrt(cos(lat(b)*pi/180));运行时出错了，错误提示：
Index exceeds matrix dimensions.
Error in test_SEOF (line 24)
weight=sqrt(cos(lat(b)*pi/180));
再根据后面程序：
for i=1:lat_num
    seof1(:,i,:)=seof1(:,i,:)/weight(i);
    seof2(:,i,:)=seof2(:,i,:)/weight(i);
end
修改成：
weight=sqrt(cos(lat*pi/180));
可以运行了，
麻烦楼主再次检查下。谢谢。

Aires

denghaijun 发表于 2015-6-16 22:25
楼主，程序中提示了“weighted for each latitude of sst, weight = sqrt(cos(lat))”，其中，weight=sqr ...

Please note that in your program, a and b refer to the lon and lat. But in my test program, a and b refer to the index of lon and lat.

denghaijun

Aires 发表于 2015-6-16 22:39
Please note that in your program, a and b refer to the lon and lat. But in my test program, a and  ...

明白了，我把程序改错了，谢谢。

幸福YOUKNOW

谢谢分享~~就是不知道哪里可以找到SEOF更具体的方法。。。

Aires

幸福YOUKNOW 发表于 2015-6-24 12:49
谢谢分享~~就是不知道哪里可以找到SEOF更具体的方法。。。

你可以参考Bin Wang (2005)年的那篇文献

swyh1105

感谢楼主分享。