标注站点图

柿子柿子柿子

最近总有人问我站点图怎么画，问多了有点烦。
正好借此机会我把我用的.ncl共享出来。
当然，这个就只当做是大家学习交流使用的啦~

1. load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
   
2. load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
   
3. load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
   
4. load "$NCARG_ROOT/lib/ncarg/nclscripts/cnmap/cnmap.ncl"
   
5. 
   
6. begin
   
7. 
   
8. wks = gsn_open_wks("eps","Station")
   
9. gsn_define_colormap( wks ,"gui_default")
   
10. 
    
11. ;=================Set Map=================;
    
12. res              = True            
    
13. res@tiMainString ="Stations Map"
    
14. res@gsnMaximize  = True
    
15. res@gsnDraw      = False
    
16. res@gsnFrame     = False
    
17. 
    
18. res@mpMinLatF               = 0.
    
19. res@mpMaxLatF               = 70.
    
20. res@mpMinLonF               = 70.
    
21. res@mpMaxLonF               = 140.
    
22. 
    
23. res@mpFillOn                = True
    
24. res@mpOutlineOn             = False                          ; Use outlines from shapefile
    
25. res@mpDataBaseVersion       = "MediumRes"
    
26. res@mpDataSetName           = "Earth..4"
    
27. res@mpAreaMaskingOn         = True
    
28. res@mpMaskAreaSpecifiers    = (/"China"/)
    
29. res@mpLandFillColor         = "white"
    
30. res@mpInlandWaterFillColor  = "white"
    
31. res@mpOceanFillColor        = "white"
    
32. res@mpOutlineBoundarySets   = "NoBoundaries"
    
33. 
    
34. map = gsn_csm_map(wks,res)
    
35. 
    
36. ;=================Draw Map=================;
    
37. 
    
38. cnres           = True
    
39. cnres@china     = True        ; draw china map or not
    
40. cnres@river     = True         ; draw changjiang&huanghe or not
    
41. cnres@province  = True   ; draw province boundary or not
    
42. cnres@nanhai    = False     ; draw nanhai or not
    
43. cnres@diqu      = True         ; draw diqujie or not
    
44. 
    
45. chinamap = add_china_map(wks,map,cnres)
    
46. 
    
47. ;=================Draw Station=================;
    
48. 
    
49. nn  = 109
    
50. lon = new((/nn/), double)
    
51. lat = new((/nn/), double)
    
52. lon = asciiread("lon.txt", (/nn/), "double")
    
53. lat = asciiread("lat.txt", (/nn/), "double")
    
54. 
    
55. mkres = True
    
56. mkres@gsMarkerIndex = 1
    
57. mkres@gsMarkerSizeF = 0.05
    
58. mkres@gsMarkerColor = 230
    
59. stnmap = gsn_add_polymarker(wks, map, lon, lat, mkres)
    
60. 
    
61. draw(map)
    
62. frame(wks)
    
63. 
    
64. end
    

复制代码

其中需要修改的地方有：经纬度范围（res@mpMinLatF、res@mpMaxLatF、res@mpMinLonF、res@mpMaxLonF），地图（res@mpMaskAreaSpecifiers）、站数（nn）、站点经纬度数据文件（lon = asciiread("lon.txt", (/nn/), "double")、lat = asciiread("lat.txt", (/nn/), "double")）、标注的图形大小及颜色（mkres@gsMarkerIndex、mkres@gsMarkerSizeF、mkres@gsMarkerColor）。
而对于cnmap的文件，我用的是在本站里面提供的那个~你的NCL没有的话就去搜一下然后根据他的安装说明弄就好了。

另外，GrADS里面标注站点的方式可以参考课本（《GrADS软件基础教程》气象出版社）里标注台风（大约P88左右）的那个例子。

stepdance

好东西！正需要

梁苏洁

学习了

FrankieLJY

谢谢楼主分享

cc4936

真的很有用，画出来了，比官网的例子容易懂

沧梧

谢谢楼主！

帅蛋fly

楼主。我的库load "$NCARG_ROOT/lib/ncarg/nclscripts/cnmap/cnmap.ncl"加载不出来是为什么？出的图没有中国地图耶

jianglegejiang

楼主 我想在上面加上各个站点的值 请问用什么函数呢？

柿子柿子柿子

jianglegejiang 发表于 2018-8-20 17:08
楼主 我想在上面加上各个站点的值 请问用什么函数呢？

gsn_add_text
不要用循环，把所有的数值做成数组

柿子柿子柿子

帅蛋fly 发表于 2018-8-13 11:07
楼主。我的库load "$NCARG_ROOT/lib/ncarg/nclscripts/cnmap/cnmap.ncl"加载不出来是为什么？出的图没有中 ...

这个是气象家园顶置贴里面的地图上数据，因为原版的中国地图有问题