MeteoInfo二次开发教程（十二）

MeteoInfo

介绍一下绘制剖面图。剖面图和水平面图一样都是二维图，可以用同样的控件来实现剖面图的绘制。

Demo程序运行后的结果如下：

                               
登录/注册后可看大图

新建一个Form（frmSectionPlot），加入一个SplitContainer，左边添加一个LayersLegend（layersLegend1），右边添加一个MapLayout（mapLayout1），将layersLegend1的MapLayout属性设为mapLayout1。


打开layersLegend1的MapFrames属性，添加一个MapFrame并将其Active属性设为true。


在frmSectionPlot的构建方法中将layersLegend1.ActiveMapFrame.MapView的IsGeoMap属性设为false（这个很重要，用于区别水平平面图）。

1. using System;
   
2. using System.Collections.Generic;
   
3. using System.ComponentModel;
   
4. using System.Data;
   
5. using System.Drawing;
   
6. using System.Linq;
   
7. using System.Text;
   
8. using System.Windows.Forms;
   
9. 
   
10. namespace MeteoInfoDemo.Forms
    
11. {
    
12.     public partial class frmSectionPlot : Form
    
13.     {
    
14.         public frmSectionPlot()
    
15.         {
    
16.             InitializeComponent();
    
17. 
    
18.             this.layersLegend1.ActiveMapFrame.MapView.IsGeoMap = false;
    
19.             this.layersLegend1.ActiveMapFrame.LayoutBounds = new Rectangle(40, 36, 606, 420);
    
20.             this.mapLayout1.MouseMode = MeteoInfoC.Layout.MouseMode.Select;
    
21.         }
    
22.     }
    
23. }

复制代码

在主窗体frmMain中添加一个菜单项，在其Click事件中加入绘制剖面图的代码。这里以model.ctl数据为例，新建MeteoDataInfo对象（aDataInfo），打开model.ctl文件，将aDataInfo的DimensionSet设为PlotDimension.Level_Lat，以读取某一经度上随纬度变化的剖面图格点数据。读取数据后需要将格点数据的X/Y坐标归一化（SetXYCoords方法），图形的X/Y坐标标注需要自己设定（MapView.XGridStrs, MapView.YGridStrs）。其它和绘制水平平面图基本一样。具体代码如下：

1.        private void TSMI_LevelLat_Click(object sender, EventArgs e)
   
2.         {
   
3.             //Create a MeteoDataInfo object
   
4.             MeteoDataInfo aDataInfo = new MeteoDataInfo();
   
5. 
   
6.             //Open GrADS data file
   
7.             string aFile = Application.StartupPath + "\\Sample\\model.ctl";
   
8.             aDataInfo.OpenGrADSData(aFile);
   
9. 
   
10.             //Set dimensions
    
11.             aDataInfo.DimensionSet = PlotDimension.Level_Lat;
    
12.             aDataInfo.LonIndex = 0;
    
13.             aDataInfo.TimeIndex = 0;
    
14. 
    
15.             //Get GridData
    
16.             GridData qdata = aDataInfo.GetGridData("Q");
    
17.             SetXYCoords(qdata);
    
18. 
    
19.             //Create a legend scheme
    
20.             bool hasUndefData = false;
    
21.             LegendScheme aLS = LegendManage.CreateLegendSchemeFromGridData(qdata,
    
22.                         LegendType.GraduatedColor, ShapeTypes.Polygon, ref hasUndefData);
    
23. 
    
24.             //Create a contour layer
    
25.             VectorLayer aLayer = DrawMeteoData.CreateShadedLayer(qdata, aLS, "Shaded_PS", "PS");
    
26. 
    
27.             //Create a sction plot form object
    
28.             frmSectionPlot sPlot = new frmSectionPlot();
    
29.             sPlot.layersLegend1.ActiveMapFrame.MapView.XGridStrs = getLatLabels(aDataInfo);
    
30.             sPlot.layersLegend1.ActiveMapFrame.MapView.YGridStrs = getLevelLabels(aDataInfo);
    
31.             sPlot.layersLegend1.ActiveMapFrame.AddLayer(aLayer);
    
32.             Extent aExtent = new Extent();
    
33.             aExtent.minX = aLayer.Extent.minX - 0.1;
    
34.             aExtent.maxX = aLayer.Extent.maxX + 0.1;
    
35.             aExtent.minY = aLayer.Extent.minY - 0.1;
    
36.             aExtent.maxY = aLayer.Extent.maxY + 0.1;
    
37.             sPlot.layersLegend1.ActiveMapFrame.MapView.ZoomToExtent(aExtent);
    
38.             sPlot.layersLegend1.Refresh();
    
39. 
    
40.             //Change title of the layout
    
41.             //LayoutGraphic aTitle = sPlot.layersLegend1.MapLayout.GetTexts()[0];
    
42.             //aTitle.SetLabelText("MeteoInfo - Section Plot - Level_Lat");
    
43. 
    
44.             //Add or change the legend in layout
    
45.             LayoutLegend aLegend;
    
46.             if (sPlot.layersLegend1.MapLayout.GetLegends().Count > 0)
    
47.                 aLegend = sPlot.layersLegend1.MapLayout.GetLegends()[0];
    
48.             else
    
49.                 aLegend = sPlot.layersLegend1.MapLayout.AddLegend(650, 100);
    
50.             aLegend.LegendStyle = LegendStyleEnum.Bar_Vertical;
    
51.             aLegend.LegendLayer = aLayer;
    
52.             sPlot.layersLegend1.MapLayout.PaintGraphics();
    
53. 
    
54.             sPlot.Show(this);
    
55.         }
    
56. 
    
57.         private List getLatLabels(MeteoDataInfo aDataInfo)
    
58.         {
    
59.             List latLabels = new List();
    
60.             double[] lats = aDataInfo.GetY();
    
61.             string drawStr;
    
62.             foreach (double lat in lats)
    
63.             {
    
64.                 if (lat > 0)
    
65.                 {
    
66.                     drawStr = lat.ToString() + "N";
    
67.                 }
    
68.                 else if (lat < 0)
    
69.                 {
    
70.                     drawStr = (-lat).ToString() + "S";
    
71.                 }
    
72.                 else
    
73.                 {
    
74.                     drawStr = "EQ";
    
75.                 }
    
76.                 latLabels.Add(drawStr);
    
77.             }
    
78. 
    
79.             return latLabels;
    
80.         }
    
81. 
    
82.         private List getLevelLabels(MeteoDataInfo aDataInfo)
    
83.         {
    
84.             List levelLabels = new List();
    
85.             List levels = aDataInfo.GetLevels();
    
86.             foreach (double level in levels)
    
87.             {
    
88.                 levelLabels.Add(level.ToString());
    
89.             }
    
90. 
    
91.             return levelLabels;
    
92.         }
    
93. 
    
94.         private void SetXYCoords(GridData aGridData)
    
95.         {
    
96.             int i;           
    
97. 
    
98.             for (i = 0; i < aGridData.XNum; i++)
    
99.             {
    
100.                 aGridData.X[i] = i;
     
101.             }
     
102.             for (i = 0; i < aGridData.YNum; i++)
     
103.             {
     
104.                 aGridData.Y[i] = i;
     
105.             }
     
106.         }

复制代码

西山守望

,好好学习

西山守望

SetXYCoords(qdata); getLatLabels(aDataInfo);好像没有这两个个命令啊，再有是不是应该把截面图中的layerslegend1和maplayout1的属性设为public啊？要不实例化后调用不了这两个控件

MeteoInfo

西山守望 发表于 2012-12-8 09:01

                               
登录/注册后可看大图

SetXYCoords(qdata); getLatLabels(aDataInfo);好像没有这两个个命令啊，再有是不是应该把截面图中的layers ...

那两个函数不是类库里的，二次开发时自己加上的。

可以设为public，或者internal等，能外部访问就可以了。

yuanliang110

学习学习中

wolfwang

很好很强大的东西啊，谢谢，学习中！

???еС?

纠结中，有谁能回答我的好多好多问题呢？

wuxudan

准备用来做航线剖面图

默存

学习中，请教王老师一个问题，如果把MapView或者MapLayout的内容保存为图片呢？看到您的MeteoInfo软件有这个功能的。请指导。

默存

默存 发表于 2013-8-16 14:47

                               
登录/注册后可看大图

学习中，请教王老师一个问题，如果把MapView或者MapLayout的内容保存为图片呢？看到您的MeteoInfo软件有这个 ...

终于找到这个函数了，ExportToPicture。