MeteoInfo二次开发教程（七）

MeteoInfo

讲讲站点数据生成的各种图层，主要包括Station Point（站点）、Station Model（站点填图）、Weather Symbol（天气现象符号）、 Station Vector（站点风场矢量）、Station Shaded（站点数据插值成格点数据再生成等值线填色图层）。

用DrawMeteoData类的一些静态方法来创建各种图层：
Station_Point:  DrawMeteoData.CreateSTPointLayer(StationData stationData, LegendScheme aLS, string layerName, string        fieldName)
Station_Model:  DrawMeteoData.CreateStationModelLayer(double[,] stationModelData, double undefine, LegendScheme aLS,    string layerName, bool isSurface)
Weather_Symbol:  DrawMeteoData.CreateWeatherSymbolLayer(StationData weatherData, LegendScheme aLS, string LName)
Station_Vector:  DrawMeteoData.CreateSTVectorLayer(StationData uData, StationData vData, StationData stData, LegendScheme aLS, bool ifColr, string LName, bool isUV)

就不多说了，大家看代码吧：

1.         private void TSMI_StationPoint_Click(object sender, EventArgs e)
   
2.         {
   
3.             //Create a MeteoDataInfo object
   
4.             MeteoDataInfo aDataInfo = new MeteoDataInfo();
   
5. 
   
6.             //Open SYNOP data file
   
7.             string aFile = Application.StartupPath + "\\Sample\\12010615.syn";
   
8.             string stFile = Application.StartupPath + "\\Sample\\SYNOP_Stations.csv";
   
9.             aDataInfo.OpenSYNOPData(aFile, stFile);
   
10. 
    
11.             //Get StationData
    
12.             StationData visData = aDataInfo.GetStationData("Visibility");
    
13. 
    
14.             //Create a legend scheme
    
15.             bool hasUndefData = false;
    
16.             LegendScheme aLS = LegendManage.CreateLegendSchemeFromStationData(visData, LegendType.GraduatedColor,
    
17.                 ShapeTypes.Point, ref hasUndefData);
    
18.             for (int i = 0; i < aLS.BreakNum; i++)
    
19.                 ((PointBreak)aLS.breakList[i]).Size = 8;
    
20. 
    
21.             //Create a contour layer
    
22.             VectorLayer aLayer = DrawMeteoData.CreateSTPointLayer(visData, aLS, "StationPoint_Vis", "Vis");
    
23. 
    
24.             //Add layer
    
25.             layersLegend1.ActiveMapFrame.AddLayer(aLayer);
    
26.             layersLegend1.Refresh();
    
27. 
    
28.             //Change title of the layout
    
29.             LayoutGraphic aTitle = mapLayout1.GetTexts()[0];
    
30.             aTitle.SetLabelText("MeteoInfo Class Library Demo - Station Point Layer");
    
31. 
    
32.             //Add or change the legend in layout
    
33.             LayoutLegend aLegend;
    
34.             if (mapLayout1.GetLegends().Count > 0)
    
35.                 aLegend = mapLayout1.GetLegends()[0];
    
36.             else
    
37.                 aLegend = mapLayout1.AddLegend(650, 100);
    
38.             aLegend.LegendStyle = LegendStyleEnum.Bar_Vertical;
    
39.             aLegend.LegendLayer = aLayer;
    
40.             if (tabControl1.SelectedIndex == 1)
    
41.                 mapLayout1.PaintGraphics();
    
42.         }
    
43.         
    
44.         private void TSMI_StationModel_Click(object sender, EventArgs e)
    
45.         {
    
46.             //Create a MeteoDataInfo object
    
47.             MeteoDataInfo aDataInfo = new MeteoDataInfo();
    
48. 
    
49.             //Open SYNOP data file
    
50.             string aFile = Application.StartupPath + "\\Sample\\12010615.syn";
    
51.             string stFile = Application.StartupPath + "\\Sample\\SYNOP_Stations.csv";
    
52.             aDataInfo.OpenSYNOPData(aFile, stFile);
    
53. 
    
54.             //Create a legend scheme
    
55.             LegendScheme aLS = LegendManage.CreateSingleSymbolLegendScheme(ShapeTypes.Point, Color.Blue, 12);
    
56. 
    
57.             //Get station model data
    
58.             double[,] stationModelData = new double[10, 1];
    
59.             Extent aExtent = new Extent();
    
60.             stationModelData = aDataInfo.GetStationModelData(ref aExtent);
    
61. 
    
62.             VectorLayer aLayer = new VectorLayer(ShapeTypes.Point);
    
63.             aLayer = DrawMeteoData.CreateStationModelLayer(stationModelData,
    
64.                     aDataInfo.UNDEF, aLS, "StationModel", true);            
    
65. 
    
66.             //Add layer
    
67.             layersLegend1.ActiveMapFrame.AddLayer(aLayer);
    
68.             layersLegend1.Refresh();
    
69. 
    
70.             //Change title of the layout
    
71.             LayoutGraphic aTitle = mapLayout1.GetTexts()[0];
    
72.             aTitle.SetLabelText("MeteoInfo Class Library Demo - Station Model Layer");
    
73.         }
    
74. 
    
75.         private void TSMI_WeatherSymbol_Click(object sender, EventArgs e)
    
76.         {
    
77.             //Create a MeteoDataInfo object
    
78.             MeteoDataInfo aDataInfo = new MeteoDataInfo();
    
79. 
    
80.             //Open SYNOP data file
    
81.             string aFile = Application.StartupPath + "\\Sample\\12010615.syn";
    
82.             string stFile = Application.StartupPath + "\\Sample\\SYNOP_Stations.csv";
    
83.             aDataInfo.OpenSYNOPData(aFile, stFile);
    
84. 
    
85.             //Get StationData
    
86.             StationData wData = aDataInfo.GetStationData("Weather");
    
87. 
    
88.             //Create a legend scheme
    
89.             LegendScheme aLS = LegendManage.CreateSingleSymbolLegendScheme(ShapeTypes.Point, Color.Blue, 15);
    
90. 
    
91.             //Create a contour layer
    
92.             VectorLayer aLayer = DrawMeteoData.CreateWeatherSymbolLayer(wData, aLS, "Weather");
    
93. 
    
94.             //Add layer
    
95.             layersLegend1.ActiveMapFrame.AddLayer(aLayer);
    
96.             layersLegend1.Refresh();
    
97. 
    
98.             //Change title of the layout
    
99.             LayoutGraphic aTitle = mapLayout1.GetTexts()[0];
    
100.             aTitle.SetLabelText("MeteoInfo Class Library Demo - Weather Symbol Layer");
     
101. 
     
102.             //Add or change the legend in layout
     
103.             LayoutLegend aLegend;
     
104.             if (mapLayout1.GetLegends().Count > 0)
     
105.                 aLegend = mapLayout1.GetLegends()[0];
     
106.             else
     
107.                 aLegend = mapLayout1.AddLegend(650, 100);
     
108.             aLegend.LegendStyle = LegendStyleEnum.Bar_Vertical;
     
109.             aLegend.LegendLayer = aLayer;
     
110.             if (tabControl1.SelectedIndex == 1)
     
111.                 mapLayout1.PaintGraphics();
     
112.         }
     
113. 
     
114.         private void TSMI_StationVector_Click(object sender, EventArgs e)
     
115.         {
     
116.             //Create a MeteoDataInfo object
     
117.             MeteoDataInfo aDataInfo = new MeteoDataInfo();
     
118. 
     
119.             //Open SYNOP data file
     
120.             string aFile = Application.StartupPath + "\\Sample\\12010615.syn";
     
121.             string stFile = Application.StartupPath + "\\Sample\\SYNOP_Stations.csv";
     
122.             aDataInfo.OpenSYNOPData(aFile, stFile);
     
123. 
     
124.             //Get StationData
     
125.             StationData wdData = aDataInfo.GetStationData("WindDirection");
     
126.             StationData wsData = aDataInfo.GetStationData("WindSpeed");
     
127. 
     
128.             //Create a legend scheme
     
129.             LegendScheme aLS = LegendManage.CreateSingleSymbolLegendScheme(ShapeTypes.Point, Color.Blue, 15);
     
130. 
     
131.             //Create a contour layer
     
132.             VectorLayer aLayer = DrawMeteoData.CreateSTVectorLayer(wdData, wsData, wdData, aLS, false, "StationVector", false);
     
133. 
     
134.             //Add layer
     
135.             layersLegend1.ActiveMapFrame.AddLayer(aLayer);
     
136.             layersLegend1.Refresh();
     
137. 
     
138.             //Change title of the layout
     
139.             LayoutGraphic aTitle = mapLayout1.GetTexts()[0];
     
140.             aTitle.SetLabelText("MeteoInfo Class Library Demo - Station Wind Vector Layer");
     
141.         }
     
142. 
     
143.         private void TSMI_StationShaded_Click(object sender, EventArgs e)
     
144.         {
     
145.             this.Cursor = Cursors.WaitCursor;
     
146. 
     
147.             //Read data info
     
148.             MeteoDataInfo aDataInfo = new MeteoDataInfo();
     
149.             string aFile = Application.StartupPath + "\\Sample\\rain_2008072220.csv";
     
150.             aDataInfo.OpenLonLatData(aFile);
     
151. 
     
152.             //Get station data
     
153.             StationData stationData = aDataInfo.GetStationData("Precipitation6h");
     
154. 
     
155.             //Interpolate
     
156.             GridDataPara aGDP = new GridDataPara();
     
157.             aGDP.dataExtent.minX = 60;
     
158.             aGDP.dataExtent.maxX = 140;
     
159.             aGDP.dataExtent.minY = -20;
     
160.             aGDP.dataExtent.maxY = 60;
     
161.             aGDP.xNum = 80;
     
162.             aGDP.yNum = 80;
     
163.             GridInterpolation gridInterp = new GridInterpolation();
     
164.             gridInterp.GridDataParaV = aGDP;
     
165. 
     
166.             gridInterp.GridInterMethodV = GridInterMethod.IDW_Radius;
     
167.             gridInterp.Radius = 2;
     
168.             gridInterp.MinPointNum = 1;
     
169. 
     
170.             double[] X = new double[1];
     
171.             double[] Y = new double[1];
     
172.             ContourDraw.CreateGridXY(gridInterp.GridDataParaV, ref X, ref Y);
     
173.             double[,] S = stationData.Data;
     
174.             S = ContourDraw.FilterDiscreteData_Radius(S, gridInterp.Radius,
     
175.                 gridInterp.GridDataParaV.dataExtent, stationData.UNDEF);
     
176.             GridData gridData = ContourDraw.InterpolateDiscreteData_Radius(S,
     
177.                 X, Y, gridInterp.MinPointNum, gridInterp.Radius, stationData.UNDEF);
     
178. 
     
179.             //Create legend scheme
     
180.             bool hasNoData = true;
     
181.             LegendScheme aLS = LegendManage.CreateLegendSchemeFromGridData(gridData, LegendType.GraduatedColor,
     
182.                 ShapeTypes.Polygon, ref hasNoData);
     
183.             ((PolygonBreak)aLS.breakList[0]).DrawFill = false;
     
184. 
     
185.             //Create layer
     
186.             VectorLayer aLayer = new VectorLayer(ShapeTypes.Polygon);
     
187.             aLayer = DrawMeteoData.CreateShadedLayer(gridData, aLS, "Rain", "Rain");
     
188.             aLayer.IsMaskout = true;
     
189. 
     
190.             //Add layer
     
191.             layersLegend1.ActiveMapFrame.AddLayer(aLayer);
     
192.             layersLegend1.ActiveMapFrame.MoveLayer(aLayer.Handle, 0);
     
193.             layersLegend1.Refresh();
     
194. 
     
195.             //Change title of the layout
     
196.             LayoutGraphic aTitle = mapLayout1.GetTexts()[0];
     
197.             aTitle.SetLabelText("MeteoInfo Class Library Demo - Station Shaded Layer");
     
198. 
     
199.             //Add or change the legend in layout
     
200.             LayoutLegend aLegend;
     
201.             if (mapLayout1.GetLegends().Count > 0)
     
202.                 aLegend = mapLayout1.GetLegends()[0];
     
203.             else
     
204.                 aLegend = mapLayout1.AddLegend(650, 100);
     
205.             aLegend.LegendStyle = LegendStyleEnum.Bar_Vertical;
     
206.             aLegend.LegendLayer = aLayer;
     
207.             if (tabControl1.SelectedIndex == 1)
     
208.                 mapLayout1.PaintGraphics();
     
209. 
     
210.             this.Cursor = Cursors.Default;
     
211.         }

复制代码

运行结果：

MeteoInfo

最近这两讲可能有些难，仔细讲的话要花很多时间准备材料，粗略地讲又会有很多东西没有提及。Demo程序已经到了一定程度了，因此将Demo程序的源代码打包上传到MeteoInfo网页，有兴趣的朋友可以去下载源代码自己琢磨。

http://www.meteothinker.com/Downloads.html 中的 MeteoInfo class library demo program。

sysun

王老师又出新教程了，关注、学习中。

fatcat916

王老师辛苦了~~

MeteoInfo

zhuvieri 发表于 2012-5-30 19:21

                               
登录/注册后可看大图

王老师 我想知道 Gis中 一个多边形，比方说 浙江这个多边形，他的中心点坐标怎么用语句写出来？好比label可 ...

在MeteoInfo中每个shape都有一个最小外接矩形（Extent），可以简单的用Extent的中心作为多边形的中心。但是对于凹多边形，这样算出来的中心有可能在多边形之外。更好的方法是计算多边形的几何中心，不过MeteoInfo目前还没有这种算法。如果你的多边形凹得不厉害，可以用我说的简化方法。或者你的多边形是固定的，比如是浙江，你可以自己手动找出它的一个大概中心，标注应该要求没那么高了。

MeteoInfo

zhuvieri 发表于 2012-5-31 08:50

                               
登录/注册后可看大图

Extent取中心点有没有语句的？还是要自己写代码算。嘻嘻。

矩形算中心点还需要函数呀，自己写吧

黄河沙

好好学习，天天向上

MeteoInfo

zhuvieri 发表于 2012-6-6 08:56

                               
登录/注册后可看大图

王老师，那个MeteoInfoC里头的一些Gis函数，又没有一个函数参考大全？或者是我可以直接下其他gis里头的函数 ...

更新了MeteoInfo类库的帮助文档，可以在此网页下载：http://www.meteothinker.com/Documentation.html

帮助文档是自动生成的，文件名：MeteoInfoC.chm

西山守望

好好学学，

geoeco

//Create legend scheme:
            bool hasNoData = true;
            LegendScheme aLS=LegendManage.CreateLegendSchemeFromGridData(gridData, LegendType.GraduatedColor,ShapeTypes.Polygon,ref hasNoData);
            ((PolygonBreak) aLS.breakList[0]).DrawFill=false;


legend 出问题了，会是哪儿的错呢？

Error      'MeteoInfoC.Legend.LegendScheme' does not contain a definition for 'breakList' and no extension method 'breakList' accepting a first argument of type 'MeteoInfoC.Legend.LegendScheme' could be found (are you missing a using directive or an assembly reference?)