循环画风玫瑰图报错求助，windows可正常使用 linux系统下报错

flyed

1. import time
   
2. def windrose2polar(a):
   
3.     r = 360 - a + 90
   
4.     r[r>360] = r - 360
   
5.     return r
   
6. 
   
7. def plotWindRole(fn, imageName, polType):
   
8.     #Read data (wind speed, weed direction, pm10)
   
9.     table = readtable(fn,delimiter=",",format='%3f')
   
10.     ws = table['WS']
    
11.     wd = table['WD']
    
12.     pm10 = table['VAL']
    
13.     N = len(ws)
    
14. 
    
15.     vmaxArray = []
    
16.     for row in pm10:
    
17.         vmaxArray.append(row)
    
18.     vmaxArray.sort()
    
19.     vmax = vmaxArray[int(len(vmaxArray)*0.99)]
    
20.    
    
21.     #Convert from windrose coordinate to polar coordinate
    
22.     rwd = windrose2polar(wd)
    
23.    
    
24.     #Degree to radians
    
25.     rwd = radians(rwd)
    
26.    
    
27.     #Calculate frequency of each wind direction bin
    
28.     wdbins = linspace(0.0, pi * 2, 9)
    
29.     dwdbins = degrees(wdbins)
    
30.     dwdbins = windrose2polar(dwdbins)
    
31.     rwdbins = radians(dwdbins)
    
32.     wdN = len(wdbins) - 1
    
33.     theta = ones(wdN + 1)
    
34.     for j in range(wdN):
    
35.         theta[j] = rwdbins[j]
    
36.     theta[wdN] = theta[0]
    
37.     wd = wd + 360./wdN/2
    
38.     wd[wd>360] = wd - 360
    
39.     wdhist = histogram(radians(wd), wdbins)[0].astype('float')
    
40.     wdhist = wdhist / N
    
41.     nwdhist = wdhist.aslist()
    
42.     nwdhist.append(nwdhist[0])
    
43.     nwdhist = array(nwdhist)
    
44.    
    
45.     #Polar coordinate to Cartesian coordinate
    
46.     rwdc, wsc = pol2cart(rwd, ws)
    
47.    
    
48.     #Get convexhull (minimum outer polygon of the wind points)
    
49.     poly = topo.convexhull(rwdc, wsc)
    
50.    
    
51.     #Get grid data
    
52.     dd = 0.5
    
53.     x = linspace(rwdc.min() - dd, rwdc.max() + dd, 5)
    
54.     y = linspace(wsc.min() - dd, wsc.max() + dd, 5)
    
55.     data = griddata((rwdc, wsc), pm10, xi=(x, y), method='idw', pointnum=5, convexhull=False)[0]
    
56.     #---------------------------------------
    
57.     #Plot figure
    
58.     pos = [0.13, 0.1, 0.775, 0.775]
    
59.    
    
60.     #Cartesian axes
    
61.     ax = axes(position=pos, aspect='equal')
    
62.     xaxis(visible=False)
    
63.     yaxis(location='right', visible=False)
    
64.     yaxis(location='left', shift = 50)
    
65.     ylabel('Wind speed (m/s)')
    
66.     levs = arange(0, vmax, vmax/7).astype(int)
    
67.     cg = contourf(x, y, data, levs, cmap='BlAqGrYeOrRe', visible=False)
    
68.     cg = cg.clip([poly])
    
69.     ax.add_graphic(cg)
    
70.     #ss = scatter(rwdc, wsc, s=6, fill=False, edgecolor='b', edgesize=1)
    
71.     #patch(poly)
    
72.     # xshift:离图的距离 shrink：colorBar的高度
    
73.     colorbar(cg, shrink=0.8, xshift=30, label=r'$\mu g/m^3'
    
74. , labelloc='bottom')
    
75.     maxv = ws.max()+1
    
76.     xlim(-maxv, maxv)
    
77.     ylim(-maxv, maxv)
    
78.     ticks = ax.get_yticks()
    
79.     ax.set_yticklabels([abs(yy) for yy in ticks])
    
80.    
    
81.     #Polar axes
    
82.     axp = axes(position=pos, polar=True)
    
83.     plot(theta, nwdhist, color='k', linewidth=2)
    
84.     axp.set_rmax(1)
    
85.     axp.set_rlabel_position(25.)
    
86.     axp.set_rtick_locations([0.2,0.4,0.6,0.8])
    
87.     #axp.set_rticks(['2','4','6','8','10 m/s'])
    
88.     axp.set_rticks(['20%','40%','60%','80%'])
    
89.     axp.set_xtick_font(size=14)
    
90.     axp.set_xtick_locations(arange(0,360,22.5))
    
91.     axp.set_xticks(['E','NEE','NE','NNE','N','NNW','NW','NWW','W','SWW',
    
92.         'SW','SSW','S','SSE','SE','SEE'])
    
93.     title(r'$Windrose \ with \ '+polType+r' \ concentrations, fontsize=18)
    
94.     savefig(imageName,640, 524)
    
95. clf()
    
96.     time.sleep(2)
    
97. 
    
98. baseDir = r"/mnt/d/pythonWorkspace/try/windrose/"
    
99. targetDir = r"/mnt/f/tmp/windRoseImages/"
    
100. 
     
101. for file in (os.listdir(baseDir)):
     
102.     if file.endswith(".txt"):
     
103.         imgName = file.replace(".txt", ".png")
     
104.         polType = imgName[imgName.rindex('-')+1:imgName.index('.')]
     
105.         if polType == 'pm25':
     
106.             polType = 'PM_{25}'
     
107.         if polType == 'o3':
     
108.             polType = 'O_{3}'
     
109.         if polType == 'pm10':
     
110.             polType = 'PM_{10}'
     
111.         polType = polType.upper()
     
112.         print(polType, baseDir, targetDir, imgName, file)
     
113.         plotWindRole(baseDir + file, targetDir + imgName, polType)
     
114.         print imgName
     
115.    
     

复制代码

如程序需求需要循环画风玫瑰图当在 clf()是否发生错误，此程序在windows下可正常运行，但在Linux中不行  linux中的报错信息为：AttributeError: 'NoneType' object has no attribute 'getFigureDock'

MeteoInfo

你是在MeteoInfolab界面里运行脚本，还是在命令行中运行的？Linux系统你是远程登录进去的吗？是否启动了图形环境？

flyed

MeteoInfo 发表于 2021-4-1 10:50
你是在MeteoInfolab界面里运行脚本，还是在命令行中运行的？Linux系统你是远程登录进去的吗？是否启动了图 ...

linux是远程进去的，没有启动图形界面

MeteoInfo

flyed 发表于 2021-4-1 10:51
linux是远程进去的，没有启动图形界面

可以尝试加上 -b 参数：  milab.sh -b test.py