尝试EOF分析时遇到问题

dnh

根据气象家园之前发的帖子，自己尝试用Python做EOF分析，但出现

  File "E:/PythonScripts/sstEOF.py", line 40, in <module>
    solver = Eof(sst_ann, weights=wgts)
  File "D:\Program\Anaconda3\lib\site-packages\eofs\standard.py", line 136, in __init__
    raise ValueError('weight array dimensions are incompatible')

ValueError: weight array dimensions are incompatible
本人是菜鸟，不知道出现了什么问题，还望大佬指教。

1. 
   

复制代码

dnh

1. import matplotlib.pyplot as plt
   
2. import numpy as np
   
3. from eofs.standard import Eof
   
4. from eofs.examples import example_data_path
   
5. from netCDF4 import Dataset
   
6. 
   
7. filename = example_data_path('F:/School/NCEPNCAR/SST/sst.mnmean.nc')
   
8. ncin = Dataset(filename, 'r')
   
9. 
   
10. lon = ncin.variables['lon'][:]
    
11. lat = ncin.variables['lat'][:]
    
12. time = ncin.variables['time'][:]
    
13. sst = ncin.variables['sst'][1201:1800,:,:]
    
14. ncin.close()
    
15. 
    
16. lat_bounds = [-30, 30]
    
17. lon_bounds = [100, 270]
    
18. 
    
19. lat_li = np.argmin(np.abs(lat - lat_bounds[0])) #
    
20. lat_ui = np.argmin(np.abs(lat - lat_bounds[1]))
    
21. # longitude lower and upper index
    
22. lon_li = np.argmin(np.abs(lon - lon_bounds[0]))
    
23. lon_ui = np.argmin(np.abs(lon - lon_bounds[1]))
    
24. sstSubset = sst[:, lat_ui:lat_li, lon_li:lon_ui]
    
25. 
    
26. sst[sst <=-9.96921e+36 ] = np.nan
    
27. 
    
28. #计算sst年均
    
29. sst_ann = np.zeros((49, lat_li-lat_ui, lon_ui-lon_li), dtype=np.float32)
    
30. for i in range(0, 49):
    
31.     for j in range(0, 11):
    
32.         sst_ann[i,:, :] = sst_ann[i,:, : ] + sstSubset[12*i+j,:,:]
    
33.     sst_ann[i,:, :] = sst_ann[i,:, :]/12.
    
34. 
    
35. lats = np.linspace(-30, 30, 60, endpoint=False)
    
36. lons = np.linspace(100, 270, 85, endpoint=False)
    
37. coslat = np.cos(np.deg2rad(lats[:]))
    
38. wgts = np.sqrt(coslat)[..., np.newaxis]
    
39. solver = Eof(sst_ann, weights=wgts)

import matplotlib.pyplot as plt<br /> import numpy as np<br /> from eofs.standard import Eof<br /> from eofs.examples import example_data_path<br /> from netCDF4 import Dataset<br /> <br /> filename = example_data_path('F:/School/NCEPNCAR/SST/sst.mnmean.nc')<br /> ncin = Dataset(filename, 'r')<br /> <br /> lon = ncin.variables['lon'][:]<br /> lat = ncin.variables['lat'][:]<br /> time = ncin.variables['time'][:]<br /> sst = ncin.variables['sst'][1201:1800,:,:]<br /> ncin.close()<br /> <br /> lat_bounds = [-30, 30]<br /> lon_bounds = [100, 270]<br /> <br /> lat_li = np.argmin(np.abs(lat - lat_bounds[0])) #<br /> lat_ui = np.argmin(np.abs(lat - lat_bounds[1]))<br /> # longitude lower and upper index<br /> lon_li = np.argmin(np.abs(lon - lon_bounds[0]))<br /> lon_ui = np.argmin(np.abs(lon - lon_bounds[1]))<br /> sstSubset = sst[:, lat_ui:lat_li, lon_li:lon_ui]<br /> <br /> sst[sst <=-9.96921e+36 ] = np.nan<br /> <br /> #计算sst年均<br /> sst_ann = np.zeros((49, lat_li-lat_ui, lon_ui-lon_li), dtype=np.float32)<br /> for i in range(0, 49):<br />     for j in range(0, 11):<br />         sst_ann[i,:, :] = sst_ann[i,:, : ] + sstSubset[12*i+j,:,:]<br />     sst_ann[i,:, :] = sst_ann[i,:, :]/12.<br /> <br /> lats = np.linspace(-30, 30, 60, endpoint=False)<br /> lons = np.linspace(100, 270, 85, endpoint=False)<br /> coslat = np.cos(np.deg2rad(lats[:]))<br /> wgts = np.sqrt(coslat)[..., np.newaxis]<br /> solver = Eof(sst_ann, weights=wgts)

cookie-o-o

还不太会python。可能是sst_ann的维度中没有和wgts维度匹配的，你可以print(wgts.shape)和print(sst_ann.shape)，看看是不是这个方面出问题了。

dnh

cookie-o-o 发表于 2019-10-22 08:12
还不太会python。可能是sst_ann的维度中没有和wgts维度匹配的，你可以print(wgts.shape)和print(sst_ann.sh ...

已解决，谢谢