在地图上绘制站点bar图

MeteoInfo

先读取所有站点的相关数据。在地图上绘制站点bar图有两种方法，一种方法是为每一个站点添加一个坐标系（Axes），坐标系的建立目前还不能直接利用地图坐标，所以需要用data2pixel函数将站点坐标转为图像坐标，然后再生成bar图所需的坐标系（需要使用units='pixels'参数），再用bar函数绘制直方图即可。需要注意的这种方式成图后图形大小交互式改变会导致坐标系位置偏差，需要重新运行脚本。

1. #Read data
   
2. fn = os.path.join(migl.get_sample_folder(), 'ASCII', 'wp.csv')
   
3. df = DataFrame.read_table(fn, delimiter=',', format='%9f', index_col=0)
   
4. lons = df['lon'].values
   
5. lats = df['lat'].values
   
6. models = df['model%'].values
   
7. ois = df['oi%'].values
   
8. enois = df['enoi%'].values
   
9. n = len(lons)
   
10. maxv = max([models.max(), ois.max(), enois.max()])
    
11. 
    
12. #Creat a map axes
    
13. #figure(figsize=[700,550], newfig=False)
    
14. proj = projinfo(proj='lcc', lon_0=105, lat_1=25, lat_2=47)
    
15. ax = axesm(outerposition=[0,0,1,1], projinfo=proj, gridline=False)
    
16. geoshow('cn_province', edgecolor='gray')
    
17. geoshow('china', edgecolor='b')
    
18. scatter(lons, lats, color='k', edge=False, size=4)
    
19. axism([78, 130, 14, 53])
    
20. 
    
21. #Draw figure so we can use data2pixel function
    
22. draw()
    
23. h = 25
    
24. w = 20
    
25. cols = ['r','b','g']
    
26. for lon,lat,model,oi,enoi in zip(lons,lats,models,ois,enois):
    
27.     x, y = ax.data2pixel(lon, lat)
    
28.     print lon, lat
    
29.     data = [model, oi, enoi]
    
30.     #Draw fixed pixels size pie charts at specific longitude/latitude
    
31.     ax1 = axes(position=[x, y, w, h], units='pixels', axis=False)   
    
32.     left = 0
    
33.     width = .1
    
34.     for i in range(3):
    
35.         bar(left, data, width, color=cols, edgecolor='gray')
    
36.         left += width  
    
37.     ylim(0, maxv)  
    
38. 
    
39. ax1 = axes(position=[150, 100, w * 3, h], units='pixels', axis=False)
    
40. left = 0
    
41. width = 1
    
42. labels = ['model (%)','oi (%)','enoi (%)']
    
43. bb = []
    
44. data = [60, 80, 100]
    
45. for i in range(3):
    
46.     b = bar(left, data, width, color=cols, edgecolor='gray', label=labels)
    
47.     bb.append(b)
    
48.     left += width  
    
49. text(left+0.5, 50, '100', xalign='left', yalign='center', fontsize=12)
    
50. ylim(0, maxv)
    
51. xlim(0, 9)
    
52. ax2 = axes(position=[130, 30, 100, 100], units='pixels', axis=False)
    
53. legend(ax1.get_legend(), loc='lower left', fontsize=12, frameon=False)
    
54. 
    
55. #Add south China Sea
    
56. axesm(position=[0.8,0.08,0.15,0.2], projinfo=proj, axison=False, gridlabel=False, frameon=True)
    
57. geoshow('cn_border', facecolor=(0,0,255))
    
58. axism([106, 123, 2, 23])

#Read data<br /> fn = os.path.join(migl.get_sample_folder(), 'ASCII', 'wp.csv')<br /> df = DataFrame.read_table(fn, delimiter=',', format='%9f', index_col=0)<br /> lons = df['lon'].values<br /> lats = df['lat'].values<br /> models = df['model%'].values<br /> ois = df['oi%'].values<br /> enois = df['enoi%'].values<br /> n = len(lons)<br /> maxv = max([models.max(), ois.max(), enois.max()])<br /> <br /> #Creat a map axes<br /> #figure(figsize=[700,550], newfig=False)<br /> proj = projinfo(proj='lcc', lon_0=105, lat_1=25, lat_2=47)<br /> ax = axesm(outerposition=[0,0,1,1], projinfo=proj, gridline=False)<br /> geoshow('cn_province', edgecolor='gray')<br /> geoshow('china', edgecolor='b')<br /> scatter(lons, lats, color='k', edge=False, size=4)<br /> axism([78, 130, 14, 53])<br /> <br /> #Draw figure so we can use data2pixel function<br /> draw()<br /> h = 25<br /> w = 20<br /> cols = ['r','b','g']<br /> for lon,lat,model,oi,enoi in zip(lons,lats,models,ois,enois):<br />     x, y = ax.data2pixel(lon, lat)<br />     print lon, lat<br />     data = [model, oi, enoi]<br />     #Draw fixed pixels size pie charts at specific longitude/latitude<br />     ax1 = axes(position=[x, y, w, h], units='pixels', axis=False)    <br />     left = 0<br />     width = .1<br />     for i in range(3):<br />         bar(left, data<i>, width, color=cols<i>, edgecolor='gray')<br />         left += width  <br />     ylim(0, maxv)  <br /> <br /> ax1 = axes(position=[150, 100, w * 3, h], units='pixels', axis=False)<br /> left = 0<br /> width = 1<br /> labels = ['model (%)','oi (%)','enoi (%)']<br /> bb = []<br /> data = [60, 80, 100]<br /> for i in range(3):<br />     b = bar(left, data<i>, width, color=cols<i>, edgecolor='gray', label=labels<i>)<br />     bb.append(b)<br />     left += width  <br /> text(left+0.5, 50, '100', xalign='left', yalign='center', fontsize=12)<br /> ylim(0, maxv) <br /> xlim(0, 9)<br /> ax2 = axes(position=[130, 30, 100, 100], units='pixels', axis=False)<br /> legend(ax1.get_legend(), loc='lower left', fontsize=12, frameon=False)<br /> <br /> #Add south China Sea<br /> axesm(position=[0.8,0.08,0.15,0.2], projinfo=proj, axison=False, gridlabel=False, frameon=True)<br /> geoshow('cn_border', facecolor=(0,0,255))<br /> axism([106, 123, 2, 23])

另一种方式是先用scatter函数生成站点图层，并将数据加入站点图层的属性表中，然后用图层的add_charts函数添加所有站点的直方图。

1. #Read data
   
2. fn = 'D:/Temp/ascii/wp.csv'
   
3. df = DataFrame.read_table(fn, delimiter=',', format='%9f', index_col=0, index_format='%s')
   
4. ids = df.index.data
   
5. lons = df['lon'].values
   
6. lats = df['lat'].values
   
7. models = df['model%'].values
   
8. ois = df['oi%'].values
   
9. enois = df['enoi%'].values
   
10. n = len(lons)
    
11. maxv = max([models.max(), ois.max(), enois.max()])
    
12. 
    
13. #Creat a map axes
    
14. figure(figsize=[1000,656], newfig=False)
    
15. proj = projinfo(proj='lcc', lon_0=105, lat_1=25, lat_2=47)
    
16. ax = axesm(outerposition=[0,0,1,1], projinfo=proj, gridline=False)
    
17. geoshow('china', edgecolor='b')
    
18. layer = scatter(lons, lats, color='k', edge=False, size=4)
    
19. layer.addfield('ID', 'string', ids)
    
20. layer.addfield('model', 'float', models)
    
21. layer.addfield('oi', 'float', ois)
    
22. layer.addfield('enoi', 'float', enois)
    
23. cols = ['r', 'g', 'b']
    
24. labels = ['model (%)','oi (%)','enoi (%)']
    
25. ls = makesymbolspec('polygon', dict(color=cols[0], edgecolor='gray', caption=labels[0]), \
    
26.     dict(color=cols[1], edgecolor='gray', caption=labels[1]),
    
27.     dict(color=cols[2], edgecolor='gray', caption=labels[2]))
    
28. layer.add_charts(['model','oi','enoi'], ls, avoidcoll=False, align='none')
    
29. layer.addlabels('ID', avoidcoll=False, yoffset=-5)
    
30. layer.movelabel('58606', x=10)
    
31. axism([78, 130, 14, 53])
    
32. 
    
33. # Add legend axes
    
34. w = 24
    
35. h = 50
    
36. ax1 = axes(position=[150, 100, w * 3, h], units='pixels', axis=False)
    
37. left = 0
    
38. width = 1
    
39. bb = []
    
40. data = [60, 80, 100]
    
41. for i in range(3):
    
42.     b = bar(left, data, width, color=cols, edgecolor='gray', label=labels)
    
43.     bb.append(b)
    
44.     left += width  
    
45. text(left+0.5, 50, '100', xalign='left', yalign='center', fontsize=12)
    
46. ylim(0, maxv)
    
47. xlim(0, 9)
    
48. ax2 = axes(position=[130, 30, 100, 100], units='pixels', axis=False)
    
49. legend(ax1.get_legend(), loc='lower left', markerscale=1.2, fontsize=12, frameon=False)
    
50. 
    
51. #Add south China Sea
    
52. axesm(position=[0.78,0.08,0.15,0.2], projinfo=proj, axison=False, gridlabel=False, frameon=True)
    
53. geoshow('cn_border', facecolor=(0,0,255))
    
54. axism([106, 123, 2, 23])

#Read data<br /> fn = 'D:/Temp/ascii/wp.csv'<br /> df = DataFrame.read_table(fn, delimiter=',', format='%9f', index_col=0, index_format='%s')<br /> ids = df.index.data<br /> lons = df['lon'].values<br /> lats = df['lat'].values<br /> models = df['model%'].values<br /> ois = df['oi%'].values<br /> enois = df['enoi%'].values<br /> n = len(lons)<br /> maxv = max([models.max(), ois.max(), enois.max()])<br /> <br /> #Creat a map axes<br /> figure(figsize=[1000,656], newfig=False)<br /> proj = projinfo(proj='lcc', lon_0=105, lat_1=25, lat_2=47)<br /> ax = axesm(outerposition=[0,0,1,1], projinfo=proj, gridline=False)<br /> geoshow('china', edgecolor='b')<br /> layer = scatter(lons, lats, color='k', edge=False, size=4)<br /> layer.addfield('ID', 'string', ids)<br /> layer.addfield('model', 'float', models)<br /> layer.addfield('oi', 'float', ois)<br /> layer.addfield('enoi', 'float', enois)<br /> cols = ['r', 'g', 'b']<br /> labels = ['model (%)','oi (%)','enoi (%)']<br /> ls = makesymbolspec('polygon', dict(color=cols[0], edgecolor='gray', caption=labels[0]), \<br />     dict(color=cols[1], edgecolor='gray', caption=labels[1]),<br />     dict(color=cols[2], edgecolor='gray', caption=labels[2]))<br /> layer.add_charts(['model','oi','enoi'], ls, avoidcoll=False, align='none')<br /> layer.addlabels('ID', avoidcoll=False, yoffset=-5)<br /> layer.movelabel('58606', x=10)<br /> axism([78, 130, 14, 53])<br /> <br /> # Add legend axes<br /> w = 24<br /> h = 50<br /> ax1 = axes(position=[150, 100, w * 3, h], units='pixels', axis=False)<br /> left = 0<br /> width = 1<br /> bb = []<br /> data = [60, 80, 100]<br /> for i in range(3):<br />     b = bar(left, data<i>, width, color=cols<i>, edgecolor='gray', label=labels<i>)<br />     bb.append(b)<br />     left += width  <br /> text(left+0.5, 50, '100', xalign='left', yalign='center', fontsize=12)<br /> ylim(0, maxv) <br /> xlim(0, 9)<br /> ax2 = axes(position=[130, 30, 100, 100], units='pixels', axis=False)<br /> legend(ax1.get_legend(), loc='lower left', markerscale=1.2, fontsize=12, frameon=False)<br /> <br /> #Add south China Sea<br /> axesm(position=[0.78,0.08,0.15,0.2], projinfo=proj, axison=False, gridlabel=False, frameon=True)<br /> geoshow('cn_border', facecolor=(0,0,255))<br /> axism([106, 123, 2, 23])

Linda.W

谢谢王老师

gaorui123

谢谢分享，值得学习

                               
登录/注册后可看大图

                               
登录/注册后可看大图

                               
登录/注册后可看大图

                               
登录/注册后可看大图

                               
登录/注册后可看大图

                               
登录/注册后可看大图

来自安卓APP客户端