MeteoInfo脚本示例：格点数据提取、扩展、合并

MeteoInfo

这个脚本程序示例了格点数据提取（Extract，从一个格点数据中提取更小区域的格点数据）、扩展（Extand，将格点数据扩大到更大的空间范围，扩大 部分用缺测值）、合并（Union，两个格点数据合并为一个格点数据，合并后的格点数据以第二个格点数据为准，但如果第二个格点数据中某些格点为缺测值则以一个格点数据为准，此操作要求两个格点数据的空间范围和格点数一致）。 需要更新为最新文件（见置顶帖子）。

提取数据：GridData.Extract(int sXIdx, int sYIdx, int xNum, int yNum) ，前两个参数是x、y方向起始位置，后两个参数是提取的x、y方向的格点数

扩展数据：GridData.Extend(int leftNum, int rightNum, int topNum, int bottomNum) ，四个参数分别是左右上下扩展的格点数。

以model.ctl数据为例，获取全球气压格点数据，用Extract提取左右两个子格点数据，将左右两个格点数据分布扩展（Extand）到全球范围，再将扩展后的两个格点数据合并（Union）成全球数据。

提取的左边格点数据：


提取的右边格点数据：


合并后的全球数据：


脚本代码如下：

1. #--------------------------------------------------------        
   
2. # Author: Yaqiang Wang                                          
   
3. # Date: 2012-10-26                                                
   
4. # Purpose: Extract, extent and union grid data  
   
5. # Note: Sample                                                   
   
6. #-----------------------------------------------------------
   
7. import clr
   
8. from System.Collections.Generic import *
   
9. clr.AddReference("MeteoInfoC.dll")
   
10. from MeteoInfoC import *
    
11. from MeteoInfoC.Data import *
    
12. from MeteoInfoC.Data.MeteoData import *
    
13. 
    
14. #---- Set directories
    
15. BaseDir = "C:\\Program Files (x86)\\MeteoInfo\"
    
16. DataDir = BaseDir + "Sample\\GrADS\"
    
17. 
    
18. #---- Extract grid data
    
19. mydata = MeteoDataInfo()
    
20. mydata.OpenGrADSData(DataDir + "model.ctl")
    
21. tnum = mydata.GetTimeNumber()
    
22. datagrid = mydata.GetGridData("PS")
    
23. leftgrid = datagrid.Extract(0, 0, 36, 46)
    
24. rightgrid = datagrid.Extract(35, 0, 37, 46)        
    
25. mipy.MeteoDataset.SetDrawType("shaded")
    
26. leftlayer = mipy.MeteoDataset.Display(leftgrid)
    
27. leftlayer.LayerName = "Left_Grid"
    
28. rightlayer = mipy.MeteoDataset.Display(rightgrid)
    
29. rightlayer.LayerName = "Right_Grid"
    
30. mipy.MapDocument.ActiveMapFrame.UpdateLayerNode(leftlayer)
    
31. mipy.MapDocument.ActiveMapFrame.UpdateLayerNode(rightlayer)
    
32. mipy.MapDocument.ActiveMapFrame.MoveLayer(leftlayer, 0)
    
33. mipy.MapDocument.ActiveMapFrame.MoveLayer(rightlayer, 0)
    
34. 
    
35. #----Extend and Union grid data
    
36. leftgrid.Extend(0, 36, 0, 0)
    
37. rightgrid.Extend(35, 0, 0, 0)
    
38. leftgrid.Union(rightgrid)
    
39. unionlayer = mipy.MeteoDataset.Display(leftgrid)
    
40. unionlayer.LayerName = "Union_Grid"
    
41. mipy.MapDocument.ActiveMapFrame.UpdateLayerNode(unionlayer)
    
42. mipy.MapDocument.ActiveMapFrame.MoveLayer(unionlayer, 0)
    
43. 
    
44. print 'Finished!'

复制代码

尽头的尽头

抢个沙发，老师辛苦了

wbriver

辛苦啦王老师~
前段时间读micaps4类数据，由于格点数较多，用Shaded绘等值线图速度太慢，您建议使用Raster，速度是快了，但是图像不是很清晰，有了提取格点方法，应该能够解决Shaded绘等值线图速度太慢的问题吧？但是我简单试了下：
MICAPS4DataInfo aDataInfo = new MICAPS4DataInfo();
aDataInfo.ReadDataInfo(openFileDialog.FileName);
GridData rainTemp = aDataInfo.GetGridData();
GridData rain = rainTemp.Extract(0, 0, 400, 400);//提取格点
调试结果显示rain的值为Null，这样使用对吗？

MeteoInfo

wbriver 发表于 2012-10-28 13:39

                               
登录/注册后可看大图

辛苦啦王老师~
前段时间读micaps4类数据，由于格点数较多，用Shaded绘等值线图速度太慢，您建议使 ...

使用应该没问题。

wbriver

MeteoInfo 发表于 2012-10-28 18:44

                               
登录/注册后可看大图

使用应该没问题。

看到您的帖子一激动，没仔细看，原来二次开发中的Extract和脚本文件稍有不一样，GridData Extract(double sX, double eX, double sY, double eY)，其中参数sX、SY是x、y方向起始位置，参数eX、eY是x、y方向结束位置。。。现在搞定了，不再畏惧Shaded绘等值线图了！
还有个小问题，在使用maskout剪裁边界外数据时，若使用Raster时，剪裁可以，但是使用Shaded时，剪裁失效？请问王老师这是为啥呢？

MeteoInfo

wbriver 发表于 2012-10-28 21:26

                               
登录/注册后可看大图

看到您的帖子一激动，没仔细看，原来二次开发中的Extract和脚本文件稍有不一样，GridData Extract(double ...

脚本用的最新的库文件，Extract方法有几种重载方式，你用的库文件可能不是最新的。

MapLayer有个IsMaskout属性，设为true即可。

♂雨已~○

王老师，这个有没有对站点数据的提取呢？比如指定经纬度范围，然后进行站点数据提取？

MeteoInfo

♂雨已~○ 发表于 2012-10-29 19:12

                               
登录/注册后可看大图

王老师，这个有没有对站点数据的提取呢？比如指定经纬度范围，然后进行站点数据提取？

有。StationData的Extract(double sX, double eX, double sY, double eY)方法。

♂雨已~○

MeteoInfo 发表于 2012-10-29 22:35

                               
登录/注册后可看大图

有。StationData的Extract(double sX, double eX, double sY, double eY)方法。

这个参数里面的sX ,eX 都是经纬度哇？还是就是按照点来做的啊？希望都是经纬度的呢，有这个功能函数没？

MeteoInfo

♂雨已~○ 发表于 2012-10-29 22:57

                               
登录/注册后可看大图

这个参数里面的sX ,eX 都是经纬度哇？还是就是按照点来做的啊？希望都是经纬度的呢，有这个功能函数没？

是经纬度，你试试就知道了。