用Micaps数据画流场一直出错

adam_5

用MICAPS第11类数据画流场，但是一直画不出来，请大神帮忙看看吧，不知道怎么去除下载数据扣金钱的规则，不好意思
错误信息

                               
登录/注册后可看大图

1. import matplotlib.pyplot as plt
   
2. from mpl_toolkits.basemap import Basemap
   
3. from matplotlib.collections import PathCollection
   
4. from matplotlib.path import Path
   
5. from matplotlib import cm
   
6. import numpy as np
   
7. class MicapsData11:
   
8. 
   
9.    def ReadData(self,filename):
   
10.         file000=open(filename)
    
11.         lines=file000.readlines()
    
12.         foredata=lines[0].split()+lines[1].split()+lines[2].split()
    
13.         dataflag = foredata[0]
    
14.         style = foredata[1]
    
15.         title = foredata[2]
    
16.         yy = int(foredata[3])
    
17.         mm = int(foredata[4])
    
18.         dd = int(foredata[5])
    
19.         hh = int(foredata[6])
    
20.         forehh = int(foredata[7])
    
21.         level = foredata[8]
    
22.         deltalon = float(foredata[9])
    
23.         deltalat = float(foredata[10])
    
24.         beginlon = float(foredata[11])
    
25.         endlon = float(foredata[12])
    
26.         beginlat = float(foredata[13])
    
27.         endlat = float(foredata[14])
    
28.         sumlon = int(foredata[15])
    
29.         sumlat = int(foredata[16])
    
30. 
    
31.         value=[]
    
32.         for i in range(len(lines)-3):
    
33.             line= lines[i+3].split()
    
34.             value+=line
    
35. 
    
36. 
    
37.         self.U = np.array(value[:sumlat*sumlon],dtype=float).reshape(sumlat,sumlon)
    
38.         self.V = np.array(value[sumlat*sumlon:],dtype=float).reshape(sumlat,sumlon)
    
39.         self.Z = np.zeros((sumlat, sumlon))
    
40.         for i in range(sumlon):
    
41.             for j in range(sumlat):
    
42.                 self.Z[j,i]=np.sqrt(self.U[j,i]**2+self.V[j,i]**2)
    
43. 
    
44.         x = np.arange(beginlon, endlon + deltalon, deltalon)
    
45.         y = np.arange(beginlat, endlat + deltalat, deltalat)
    
46.         self.X, self.Y= np.meshgrid(x,y)
    
47. 
    
48. def Draw(self):
    
49.         fig=plt.figure()
    
50.         ax=fig.add_subplot(111)
    
51. 
    
52. 
    
53. 
    
54.         bmap = Basemap(width=6000000, height=5000000,
    
55.                        resolution='l', projection='lcc',
    
56.                        lat_1=45., lat_2=25, lat_0=35, lon_0=105., ax=ax)
    
57.         bmap.drawcoastlines()
    
58.         bmap.fillcontinents(color='ivory', lake_color='skyblue', zorder=0)
    
59. 
    
60.         shp_info = bmap.readshapefile("E:\\python\\projects\\MicapsDataDraw-master\\shapefile\\bou2_4l", 'states',
    
61.                                       drawbounds=True)
    
62.         paths = []
    
63.         for line in shp_info[4]._paths:
    
64.             paths.append(Path(line.vertices, codes=line.codes))
    
65.         PathCollection(paths, linewidths=5, facecolors=cm.jet(np.random.rand(1)), zorder=1)
    
66. 
    
67.         
    
68.         bmap.drawmeridians(range(30, 160, 10), linewidth=0.5, dashes=[4, 2], labels=[0, 0, 0, 1], color='k', zorder=1,
    
69.                            fontsize=10)
    
70.         bmap.drawparallels(range(10, 80, 10), linewidth=0.5, dashes=[4, 2], labels=[1, 1, 0, 0], color='k', zorder=1,
    
71.                            fontsize=10)
    
72.         
    
73.         bmap.drawmapboundary(fill_color='skyblue')
    
74. 
    
75.         U_rot,V_rot, X,Y,= bmap.rotate_vector(self.U,self.V,self.X,self.Y,returnxy=True)
    
76.         ax.streamplot(X,Y,U_rot,V_rot,density=[1,1],zorder=0)
    
77. 
    
78.        plt.show()
    
79. 
    
80. 
    
81. if __name__ == '__main__':
    
82. 
    
83.     micapsData11 = MicapsData11()
    
84.     micapsData11.ReadData('E:\\micaps\\qxj\\micaps\\high\\uv\\500\\16110908.000')
    
85. 
    
86.     micapsData11.Draw()
    
87. 
    
88. 
    
89. 
    
90. 
    

复制代码

adam_5

前面的问题去掉rotate_vector命令就不会再出x行不相等的问题，但又出现了新的问题，

po_po1

看提示是说数据不匹配，看看main()函数中用到的X,Y,U_rot,V_rot的shape

adam_5

po_po1 发表于 2017-5-23 15:31
看提示是说数据不匹配，看看main()函数中用到的X,Y,U_rot,V_rot的shape

你是说第二个问题吗？都看了我把shape都print出来，都是一样的。最后把行列数都改成双数就合适，单数就出问题，但是在MATPLOTLIB的实例中改成单数就能出来，真是搞不懂了