求助：风场流线图

geoeco

想画一个县域的地面风场月变化图，下载了U/V 风场两个文件：uwnd.mon.ltm.nc
vwnd.mon.ltm.nc

请问Meteoinfo里是否可以实现用这两个文件作流线图，如何做,有没有教程可以参考？

一般风玫瑰图就可以了，做流线变化图主要目的是想做的更直观一些，试了一下，没有成功，所以在坛子里求助，
先谢谢了

MeteoInfo

可以利用脚本程序来做。

1. #-----------------------------------------------------
   
2. # Author: Yaqiang Wang
   
3. # Date: 2014-11-24
   
4. # Purpose: Read U/V data from two files and plot streamline
   
5. # Note: Sample
   
6. #-----------------------------------------------------
   
7. 
   
8. import os
   
9. from org.meteoinfo.data.meteodata import MeteoDataInfo, DrawMeteoData
   
10. from org.meteoinfo.legend import LegendManage
    
11. from org.meteoinfo.shape import ShapeTypes
    
12. from java.awt import Color
    
13. 
    
14. #---- Set directories
    
15. dataDir = 'D:/Temp/nc'
    
16. 
    
17. #---- Open data
    
18. print 'Open data...'
    
19. DataU = MeteoDataInfo()
    
20. DataV = MeteoDataInfo()
    
21. ufn = os.path.join(dataDir, 'uwnd.sig995.2010.nc')
    
22. vfn = os.path.join(dataDir, 'vwnd.sig995.2010.nc')
    
23. DataU.openNetCDFData(ufn)
    
24. DataV.openNetCDFData(vfn)
    
25. 
    
26. #---- Get grid data
    
27. print 'Get grid data...'
    
28. uwnd = DataU.getGridData('uwnd')
    
29. vwnd = DataV.getGridData('vwnd')
    
30. 
    
31. #---- Draw streamline
    
32. print 'Draw streamline...'
    
33. ls = LegendManage.createSingleSymbolLegendScheme(ShapeTypes.Polyline, Color.blue, 1)
    
34. aLayer = DrawMeteoData.createStreamlineLayer(uwnd, vwnd, 4, ls, 'Streamline_UV', True)
    
35. miapp.getMapDocument().getActiveMapFrame().addLayer(aLayer)
    
36. 
    
37. print 'Finished!'

#-----------------------------------------------------<br /> # Author: Yaqiang Wang<br /> # Date: 2014-11-24<br /> # Purpose: Read U/V data from two files and plot streamline<br /> # Note: Sample<br /> #-----------------------------------------------------<br /> <br /> import os<br /> from org.meteoinfo.data.meteodata import MeteoDataInfo, DrawMeteoData<br /> from org.meteoinfo.legend import LegendManage<br /> from org.meteoinfo.shape import ShapeTypes<br /> from java.awt import Color<br /> <br /> #---- Set directories<br /> dataDir = 'D:/Temp/nc'<br /> <br /> #---- Open data<br /> print 'Open data...'<br /> DataU = MeteoDataInfo()<br /> DataV = MeteoDataInfo()<br /> ufn = os.path.join(dataDir, 'uwnd.sig995.2010.nc')<br /> vfn = os.path.join(dataDir, 'vwnd.sig995.2010.nc')<br /> DataU.openNetCDFData(ufn)<br /> DataV.openNetCDFData(vfn)<br /> <br /> #---- Get grid data<br /> print 'Get grid data...'<br /> uwnd = DataU.getGridData('uwnd')<br /> vwnd = DataV.getGridData('vwnd')<br /> <br /> #---- Draw streamline<br /> print 'Draw streamline...'<br /> ls = LegendManage.createSingleSymbolLegendScheme(ShapeTypes.Polyline, Color.blue, 1)<br /> aLayer = DrawMeteoData.createStreamlineLayer(uwnd, vwnd, 4, ls, 'Streamline_UV', True)<br /> miapp.getMapDocument().getActiveMapFrame().addLayer(aLayer)<br /> <br /> print 'Finished!'

geoeco

谢谢群主快速回复，34行的参数 4
aLayer = DrawMeteoData.createStreamlineLayer(uwnd, vwnd, 4, ls, 'Streamline_UV', True)

是指4月吗？

不好意思，不是这个专业，全是基础的问题。

MeteoInfo

geoeco 发表于 2014-11-24 11:15
谢谢群主快速回复，34行的参数 4
aLayer = DrawMeteoData.createStreamlineLayer(uwnd, vwnd, 4, ls, 'Str ...

这个是流线密度的参数

geoeco

谢谢，正在看教程六。
如何画出不同月的平均流线呢？

MeteoInfo

geoeco 发表于 2014-11-24 11:36
谢谢，正在看教程六。
如何画出不同月的平均流线呢？

1. 需要下载最新的MeteoInfo Java 1.1.6R2版本（http://www.meteothinker.com/Downloads.aspx）

2. 利用下面的脚本计算一月份的平均风场U/V分量，然后绘制流场图，其它月份可以参考自己改改。

1. #-----------------------------------------------------
   
2. # Author: Yaqiang Wang
   
3. # Date: 2014-11-24
   
4. # Purpose: Read Monthly averaged U/V data from two files and plot streamline
   
5. # Note: Sample
   
6. #-----------------------------------------------------
   
7. 
   
8. import os
   
9. from org.meteoinfo.data import GridData
   
10. from org.meteoinfo.data.meteodata import MeteoDataInfo, DrawMeteoData
    
11. from org.meteoinfo.legend import LegendManage
    
12. from org.meteoinfo.shape import ShapeTypes
    
13. from java.awt import Color
    
14. from java.util import Calendar
    
15. 
    
16. #---- Set directories
    
17. dataDir = 'D:/Temp/nc'
    
18. 
    
19. #---- Open data
    
20. print 'Open data...'
    
21. DataU = MeteoDataInfo()
    
22. DataV = MeteoDataInfo()
    
23. ufn = os.path.join(dataDir, 'uwnd.sig995.2010.nc')
    
24. vfn = os.path.join(dataDir, 'vwnd.sig995.2010.nc')
    
25. DataU.openNetCDFData(ufn)
    
26. DataV.openNetCDFData(vfn)
    
27. 
    
28. #---- Get January average U/V data
    
29. print 'Get grid data...'
    
30. cal = Calendar.getInstance()
    
31. tn = DataU.getDataInfo().getTimeNum()
    
32. u_sum = GridData()
    
33. v_sum = GridData()
    
34. n = 0
    
35. i = 0
    
36. for t in DataU.getDataInfo().getTimes():
    
37.         cal.setTime(t)
    
38.         if cal.get(Calendar.MONTH) + 1 == 1:
    
39.                 print t
    
40.                 DataU.setTimeIndex(i)
    
41.                 DataV.setTimeIndex(i)
    
42.                 u = DataU.getGridData('uwnd')
    
43.                 v = DataV.getGridData('vwnd')
    
44.                 if n == 0:
    
45.                         u_sum = u
    
46.                         v_sum = v
    
47.                 else:
    
48.                         u_sum = u_sum.add(u)
    
49.                         v_sum = v_sum.add(v)
    
50.                 n += 1
    
51. 
    
52.         i += 1
    
53. 
    
54. uwnd = u_sum.div(n)
    
55. vwnd = v_sum.div(n)
    
56. 
    
57. #---- Draw streamline
    
58. print 'Draw streamline...'
    
59. layer = DrawMeteoData.createStreamlineLayer(uwnd, vwnd, 'Streamline_January', True)
    
60. miapp.getMapDocument().getActiveMapFrame().addLayer(layer)
    
61. 
    
62. print 'Finished!'

#-----------------------------------------------------<br /> # Author: Yaqiang Wang<br /> # Date: 2014-11-24<br /> # Purpose: Read Monthly averaged U/V data from two files and plot streamline<br /> # Note: Sample<br /> #-----------------------------------------------------<br /> <br /> import os<br /> from org.meteoinfo.data import GridData<br /> from org.meteoinfo.data.meteodata import MeteoDataInfo, DrawMeteoData<br /> from org.meteoinfo.legend import LegendManage<br /> from org.meteoinfo.shape import ShapeTypes<br /> from java.awt import Color<br /> from java.util import Calendar<br /> <br /> #---- Set directories<br /> dataDir = 'D:/Temp/nc'<br /> <br /> #---- Open data<br /> print 'Open data...'<br /> DataU = MeteoDataInfo()<br /> DataV = MeteoDataInfo()<br /> ufn = os.path.join(dataDir, 'uwnd.sig995.2010.nc')<br /> vfn = os.path.join(dataDir, 'vwnd.sig995.2010.nc')<br /> DataU.openNetCDFData(ufn)<br /> DataV.openNetCDFData(vfn)<br /> <br /> #---- Get January average U/V data<br /> print 'Get grid data...'<br /> cal = Calendar.getInstance()<br /> tn = DataU.getDataInfo().getTimeNum()<br /> u_sum = GridData()<br /> v_sum = GridData()<br /> n = 0<br /> i = 0<br /> for t in DataU.getDataInfo().getTimes():<br />         cal.setTime(t)<br />         if cal.get(Calendar.MONTH) + 1 == 1:<br />                 print t<br />                 DataU.setTimeIndex(i)<br />                 DataV.setTimeIndex(i)<br />                 u = DataU.getGridData('uwnd')<br />                 v = DataV.getGridData('vwnd')<br />                 if n == 0:<br />                         u_sum = u<br />                         v_sum = v<br />                 else:<br />                         u_sum = u_sum.add(u)<br />                         v_sum = v_sum.add(v)<br />                 n += 1<br /> <br />         i += 1<br /> <br /> uwnd = u_sum.div(n)<br /> vwnd = v_sum.div(n)<br /> <br /> #---- Draw streamline<br /> print 'Draw streamline...'<br /> layer = DrawMeteoData.createStreamlineLayer(uwnd, vwnd, 'Streamline_January', True)<br /> miapp.getMapDocument().getActiveMapFrame().addLayer(layer)<br /> <br /> print 'Finished!'

geoeco

谢谢群主，已经搞好了。
还有一个和编程无关的问题，就是月风场平均数据每日有多个数据 （6小时间隔），
对于月平均值， 选哪个时段的数据最有代表性？目前我每天只取6：00 的数据做的月平均。

MeteoInfo

geoeco 发表于 2014-11-24 14:23
谢谢群主，已经搞好了。
还有一个和编程无关的问题，就是月风场平均数据每日有多个数据 （6 ...

这个你根据自己的需求定吧。

cs1204

感谢分享
！

cs1204

感谢分享
！