pyngl绘图命令解析

喵星观天

惊闻pyngl取代了ncl的绘图功能，打算把几年前制作的pyng命令手册发布出来。

Ngl全称：Ncl's Graphic Library
中文含义：ncl移植到了python的绘图程序


ngl绘图效果：

                               
登录/注册后可看大图

示例代码：

1. 
   
2. from __future__ import print_function
   
3. import numpy as np
   
4. import Nio, Ngl, os, sys
   
5. from wrf import getvar, latlon_coords, to_np
   
6. 
   
7. 
   
8. # Read data
   
9. filename = "wrfout_d03_2012-04-22_23_00_00"
   
10. if(not os.path.exists(filename)):
    
11.   print("You do not have the necessary '{}' file to run this example.".format(filename))
    
12.   print("You need to supply your own WRF output file")
    
13.   print("WRF output files usually have names like '{}'".format(filename))
    
14.   sys.exit()
    
15. 
    
16. # Read some WRF data
    
17. a  = Nio.open_file(filename+".nc")  # Must add ".nc" suffix for Nio.open_file
    
18. ua = getvar(a,"ua")
    
19. va = getvar(a,"va")
    
20. 
    
21. # First timestep, lowest (bottommost) level, every 5th lat/lon
    
22. nl    = 0
    
23. nt    = 0
    
24. nstep = 5     # a stride to cull some of the streamlines
    
25. u     = ua[nl,::nstep,::nstep]
    
26. v     = va[nl,::nstep,::nstep]
    
27. spd   = np.sqrt(u**2+v**2)
    
28. 
    
29. # Get the latitude and longitude points
    
30. lat, lon = latlon_coords(ua)
    
31. lat = to_np(lat)
    
32. lon = to_np(lon)
    
33. 
    
34. # Open file for graphics
    
35. wks_type = "png"
    
36. wks = Ngl.open_wks(wks_type,"wrf4")
    
37. 
    
38. res                   = Ngl.Resources()
    
39. 
    
40. res.mpDataBaseVersion = "MediumRes"         # Better map outlines
    
41. res.mpLimitMode       = "LatLon"            # Zoom in on map area of interest
    
42. res.mpMinLatF         = np.min(lat[:])-0.1
    
43. res.mpMaxLatF         = np.max(lat[:])+0.1
    
44. res.mpMinLonF         = np.min(lon[:])-0.1
    
45. res.mpMaxLonF         = np.max(lon[:])+0.1
    
46. 
    
47. res.mpFillOn                = True
    
48. res.mpLandFillColor         = "gray85"
    
49. res.mpOceanFillColor        = "transparent"
    
50. res.mpInlandWaterFillColor  = "transparent"
    
51. res.mpGridLatSpacingF       = 1
    
52. res.mpGridLonSpacingF       = 1
    
53. #res.mpGridAndLimbOn        = False
    
54. 
    
55. res.stLineThicknessF  = 3.0          # Thicker streamlines
    
56. res.stMonoLineColor   = False        # Use multiple colors for streamlines
    
57. 
    
58. res.tiMainString      = "U10/V10 streamlines color by wind speed"
    
59. 
    
60. # Necessary to overlay on map correctly.
    
61. res.vfYArray          = lat[::nstep,::nstep]
    
62. res.vfXArray          = lon[::nstep,::nstep]
    
63. 
    
64. plot = Ngl.streamline_scalar_map(wks,u,v,spd,res)
    
65. 
    
66. Ngl.end()
    

复制代码

转帖教材请征求本人同意

kc121

你这命令听起来比较强大啊！学习一下。

chiqu296

非常实用，感谢分享！

huyuhengwan

感谢分享！！

beesfly

已经下载，谢谢楼主分享。以后就要走python的路了

astiny

8说了，就差一个as-yet-to-be-named Python package了。NCAR gkd。

cqcswky

能否介绍一下pyngl是什么？

喵星观天

cqcswky 发表于 2019-2-21 19:29
能否介绍一下pyngl是什么？

ncl移植到了python的绘图程序

sreveryone

十分感谢您对ncl、python在地球科学领域的应用的传播做的贡献

yqcxm

感谢分享！！