Python之Anaconda编译器做EOF分析借鉴他人后小结

LemonTree007

最近在做EOF分析，原始脚本是参考了——python做EOF并绘制前3个特征向量和时间序列http://bbs.06climate.com/forum.p ... 1&fromuid=69428(出处: 气象家园)
谢谢这位大佬的分享，学到了很多！
但我要处理的数据是四维(time,level,longitude,latitude)的风场，涉及到固定一个level的问题，而参考脚本是三维的，本是个小问题，但对于python的初学者(就是我)，还不是很清楚。这里我把修改后的脚本附在这里，并加了一些必要的说明。

1. import netCDF4
   
2. import matplotlib.pyplot as plt
   
3. import numpy as np
   
4. from mpl_toolkits.basemap import Basemap
   
5. from eofs.standard import Eof
   
6. 
   
7. # read uwind
   
8. nc_obj = netCDF4.Dataset('D:\\wind_u_v\\ERA\\v.mon.mean.1958.2016.nc', 'r')
   
9. print(nc_obj.variables.keys())
   
10. lon = nc_obj.variables['longitude'][:]
    
11. lat = nc_obj.variables['latitude'][:]
    
12. time = nc_obj.variables['time'][:]
    
13. level = nc_obj.variables['level'][:]
    
14. v = nc_obj.variables['v'][:, 2, :, :] # 重点就是这里啦，我要研究850hPa，而level = 850hPa是位于’level‘的第三个位置，所以index是2
    
15. nc_obj.close()
    
16. 
    
17. # 截取研究区域和时段u风数据
    
18. lat_bounds = [-5, 30]
    
19. lon_bounds = [100, 130]
    
20. # latitude lower and upper index

复制代码

但其实现在mv-eof分析更有意义，而且地图处理的话也时兴用cartopy，所以下一步准备向这方面靠拢，做出来了再给大家分享。
这里附加一下mv-eof的官方说明网址。https://ajdawson.github.io/eofs/userguide/multivariate.html

Ｏrz

博主是什么问题，改了脚本后，不能出图吗

LemonTree007

Ｏrz 发表于 2019-3-31 15:38
博主是什么问题，改了脚本后，不能出图吗

修改之前不能出图，改了之后可以

Ｏrz

LemonTree007 发表于 2019-3-31 15:59
修改之前不能出图，改了之后可以

能把修改后的代码放出来吗，我也是刚学习的，想参考一下，谢谢

LemonTree007

Ｏrz 发表于 2019-3-31 16:13
能把修改后的代码放出来吗，我也是刚学习的，想参考一下，谢谢

我粘贴的就是修改之后的，不然我放个原来的有啥意思？？？

chongzika

我想问一下，利用站点数据如何做EOF 分析呢

LemonTree007

chongzika 发表于 2019-4-6 16:13
我想问一下，利用站点数据如何做EOF 分析呢

没有做过，家园应该有