关于meteoinfo类库开发中的疑问

孤蓝et

        近期总算着手从VB6想C#转变，主要是有meteoinfo强大的气象分析类库做为后盾，使我有信心将以前做过和后续的开发向Meteoinfo类库转换。而且斑竹的Meteoinfo类库完整的demo程序，基本能给大家一个着手写基于meteoinfo类库开发的引子。不过作为C#初学者，对代码的阅读和理解毕竟有限，还是希望能通过学习，多做一些进步吧。之前在demo的基础上修改程序的过程中，也遇到不少问题，大部分是我不精通C#和meteoinfo的说明文档，但是斑竹还是给予细心的解答。所以希望我此处的提问、技术疑问能得到大家的帮助，也把自己的问题和解决方案写出来与大家分享。


1、demo程序中的default.mip中出现中文字符乱码情况。
解决方法（修改部分代码，在程序中加入<?xml version='1.0' encoding='gb2312'?>，或者在default.mip第一行加入<?xml version="1.0" encoding="gb2312"?>）：

1. 
   
2.          public void SaveProjFile(string aFile)
   
3.         {
   
4.             m_FileName = aFile;
   
5.             XmlDocument doc = new XmlDocument();
   
6.             doc.LoadXml("<?xml version='1.0' encoding='gb2312'?><MeteoInfo name='" + Path.GetFileNameWithoutExtension(aFile) + "' type='projectfile'></MeteoInfo>");
   
7.             XmlElement root = doc.DocumentElement;            
   
8.             //Add MapView content
   
9.             frmMain.G_LayerLegend.MapView.ExportProjectXML(ref doc, root, m_FileName);               
   
10.             //Save project file            
    
11.             doc.Save(aFile);
    
12.         }

复制代码

2、图层名称自定义
        demo程序中，比如绘制第一类数据的示例程序，默认情况，左侧图层列表中的名词前部写为“stationmodel_”，这部分字段假如我不需要，或者改为其他名称，能实现吗？      
          通过下述两句代码，即可实现图层名称的自定义。很简单。      
            aLayer.IsMaskout = true;
            aLayer.LayerName = "测试数据";

1. 
   
2.         private void TSMI_MICAPS1_StationModel_Click(object sender, EventArgs e)
   
3.         {
   
4.             MICAPS1DataInfo aDataInfo = new MICAPS1DataInfo();
   
5.             string aFile = Application.StartupPath + @"\Sample\09031417.000";
   
6.             aDataInfo.ReadDataInfo(aFile);
   
7.             LegendScheme aLS = LegendManage.CreateSingleSymbolLegendScheme(ShapeTypes.Point, Color.Blue, 12);
   
8.             double[,] stationModelData = new double[10, 1];
   
9.             Extent aExtent = new Extent();
   
10.             stationModelData = aDataInfo.GetStationModelData(ref aExtent);
    
11.             VectorLayer aLayer = new VectorLayer(ShapeTypes.Point);
    
12.             aLayer = DrawMeteoData.CreateStationModelLayer(stationModelData,
    
13.                     aDataInfo.UNDEF, aLS, "MICAPS1");
    
14.             aLayer.IsMaskout = true;
    
15.             aLayer.LayerName = "测试数据";
    
16.             int aLayerHandle = G_LayerLegend.AddLayer(aLayer, true);
    
17.         }

复制代码

3、地名显示【没有实现】
      目前meteoinfo支持丰富的地图文件，暂时不太熟需shp文件。所以我在着手考虑实现micaps 第17类数据的支持，来实现台站名的显示，这个在业务平台的开发中还算重要，能让预报人员明白相关台站的信息。此部分参照GIS相关知识或http://bbs.06climate.com/forum.p ... &extra=page%3D1中介绍的方法，目前本人未能实验成功。如有解决方案，立即呈上。

4、利用MICAPS第一类数据进行离散点分析。
读入数组的使用，stationmodel和stationdata的使用？以及等值线分析方法。
      今天看了许久的demo和meteoinfo class library help，忙活了一晚上也没有实现我需要的结果，为了实现打开micaps第一类数据后，自动分析温度等值线。我想参照TSMI_LLStShaded_Click，TSMI_MICAPS1_StationModel_Click这两部分的代码，看能不能合并实现我的想法，可是在数据读入就发现很没有头绪：
micaps 第一类数据的读入方法
stationModelData = aDataInfo.GetStationModelData(ref aExtent);
经纬度的数据
StationData stationData = aDataInfo.GetStationData(0);
             //Read micaps_1 data info,将micaps 第一类数据使用linlatstationdata方法读入。
            //stationModelData[0, 0] stationModelData[1, 0] 经纬度
            //stationModelData[2, 0] stationModelData[3, 0] 风向 风俗
            //stationModelData[7, 0] 温度
            //stationModelData[8, 0] 露点
            //stationModelData[9, 0] 修正海平面气压

            //GetStationData,将micaps 第一类数据使用GetStationData方法读入。
            //GetStationData(1)  风向
            //GetStationData(2)  风速
            //GetStationData(3)  气压
            //GetStationData(11) 露点/相对湿度
            //GetStationData(14) 温度
            //stationData       得到符合要求的经纬度、数据记录。这个是程序用来绘制LATLON数据。

      我怎么能从stationModelData = aDataInfo.GetStationModelData(ref aExtent);里面找到经纬度、温度数据在数组的位置，然后赋值给S = ContourDraw.FilterDiscreteData_Radius(S, gridInterp.Radius,gridInterp.GridDataParaV.dataExtent, stationData.UNDEF);用来绘图？

以下代码即为我修改后的代码。大家可以试试。

1.         private void TSMI_MICAPS1_StationModel_Click(object sender, EventArgs e)
   
2.         {
   
3.             MICAPS1DataInfo aDataInfo = new MICAPS1DataInfo();
   
4.             string aFile = Application.StartupPath + @"\Sample\09031417.000";
   
5.             aDataInfo.ReadDataInfo(aFile);
   
6. 
   
7.             LegendScheme aLS = LegendManage.CreateSingleSymbolLegendScheme(ShapeTypes.Point, Color.Blue, 12);
   
8. 
   
9.             double[,] stationModelData = new double[10, 1];
   
10.             Extent aExtent = new Extent();
    
11.             stationModelData = aDataInfo.GetStationModelData(ref aExtent);
    
12. 
    
13.             //Read micaps_1 data info,将micaps 第一类数据使用linlatstationdata方法读入。
    
14.             //stationModelData[0, 0] stationModelData[1, 0] 经纬度
    
15.             //stationModelData[2, 0] stationModelData[3, 0] 风向 风俗
    
16.             //stationModelData[7, 0] 温度
    
17.             //stationModelData[8, 0] 露点
    
18.             //stationModelData[9, 0] 修正海平面气压
    
19. 
    
20.             VectorLayer aLayer = new VectorLayer(ShapeTypes.Point);
    
21.             aLayer = DrawMeteoData.CreateStationModelLayer(stationModelData,
    
22.                     aDataInfo.UNDEF, aLS, "MICAPS1");
    
23.             aLayer.IsMaskout = true;
    
24.             aLayer.LayerName = "测试数据";
    
25.             int aLayerHandle = G_LayerLegend.AddLayer(aLayer, true);
    
26. 
    
27. 
    
28.             //GetStationData,将micaps 第一类数据使用GetStationData方法读入。
    
29.             //GetStationData(1)  风向
    
30.             //GetStationData(2)  风速
    
31.             //GetStationData(3)  气压
    
32.             //GetStationData(11) 露点/相对湿度
    
33.             //GetStationData(14) 温度
    
34.             //stationData       得到符合要求的经纬度、数据记录。这个是程序用来绘制LATLON数据。
    
35.             StationData stationData = aDataInfo.GetStationData(14);
    
36.             double[,] S = stationData.Data;
    
37. 
    
38. 
    
39.             
    
40. 
    
41.             //Interpolate
    
42.             GridDataPara aGDP = new GridDataPara();
    
43.             aGDP.dataExtent.minX = 70;
    
44.             aGDP.dataExtent.maxX = 100;
    
45.             aGDP.dataExtent.minY = 33;
    
46.             aGDP.dataExtent.maxY = 50;
    
47.             aGDP.xNum = 30;
    
48.             aGDP.yNum = 17;
    
49.             GridInterpolation gridInterp = new GridInterpolation();
    
50.             gridInterp.GridDataParaV = aGDP;
    
51. 
    
52.             gridInterp.GridInterMethodV = GridInterMethod.IDW_Radius;
    
53.             gridInterp.Radius = 2;
    
54.             gridInterp.MinPointNum = 1;
    
55. 
    
56.             double[] X = new double[1];
    
57.             double[] Y = new double[1];
    
58.             ContourDraw.CreateGridXY(gridInterp.GridDataParaV, ref X, ref Y);
    
59. 
    
60.             S = ContourDraw.FilterDiscreteData_Radius(S, gridInterp.Radius,
    
61.                 gridInterp.GridDataParaV.dataExtent, stationData.UNDEF);
    
62.             GridData gridData = ContourDraw.InterpolateDiscreteData_Radius(S,
    
63.                 X, Y, gridInterp.MinPointNum, gridInterp.Radius, stationData.UNDEF);
    
64. 
    
65. 
    
66.             //Create legend scheme  等值线
    
67.             bool hasNoData = true;
    
68.             LegendScheme aLS_2 = LegendManage.CreateLegendSchemeFromGridData(gridData, LegendType.GraduatedColor,
    
69.                 ShapeTypes.Polyline, ref hasNoData);
    
70.             
    
71. 
    
72.             //Create layer
    
73.             VectorLayer aLayer_2 = new VectorLayer(ShapeTypes.Polyline);
    
74.             aLayer_2 = DrawMeteoData.CreateContourLayer(gridData, aLS_2, "Rain");
    
75.             aLayer_2.IsMaskout = true;
    
76.             G_LayerLegend.AddLayer(aLayer_2, true);
    
77. 
    
78.             //Create legend scheme  网格点 point
    
79.             //Create legend scheme  填色 shade
    
80. 
    
81. 
    
82. 
    
83.             this.Cursor = Cursors.Default;
    
84. 
    
85.         }
    

复制代码

后续完成这部分的修改，我会把我修改的demo程序共享出来，也想把自己使用meteoinfo做二次开发的过程写下来，能为其他人做一个借鉴把。也希望更多的人能使用这个类库开发更实用的业务平台。
感谢meteoinfo斑竹无私的奉献。

MeteoInfo

赞一个先！

1、是正解，也很简单就不多说了。
2、文件名是可以用代码控制的，在aLayer.IsMaskout = true; 后加一行aLayer.LayerName = "..."即可。
3、其实就是个站点文件，可以参照我的站点文件绘制图形的帖子生成站点shape文件。

zsmnantong

感谢楼主和MeteoInfo 斑竹，很给力！

孤蓝et

MeteoInfo 发表于 2011-9-4 17:39

                               
登录/注册后可看大图

赞一个先！

1、是正解，也很简单就不多说了。

关于台站shp文件的制作，斑竹应该是推荐这个帖子吧？http://bbs.06climate.com/forum.p ... &extra=page%3D1
我看了下，通过这个方法是可以制作shp文件，但是输出的shp仅能实现台站的点的显示，我还需要显示台站的中文名？这个用这个方法实现不了吧？
好像看了以前的帖子，也的确没有找到您的这个软件有直接显示台站的名称的方法。

MeteoInfo

孤蓝et 发表于 2011-9-5 00:43

                               
登录/注册后可看大图

关于台站shp文件的制作，斑竹应该是推荐这个帖子吧？http://bbs.06climate.com/forum.php?mod=viewthread ...

当然可以显示中文名，你的站点文本文件中不是有中文站名吗。用Excel编辑成MeteoInfo的Lon/Lat Station文本格式（具体见你提到的帖子），中文站名也作为一列就ok了，在MeteoInfo中绘图时选择Station Info，这样所有字段都会在shape文件中了。

或者你把micaps第17类站点文件作为附件放上，我来做一个shape文件共享给大家。

孤蓝et

MeteoInfo 发表于 2011-9-5 08:19

                               
登录/注册后可看大图

当然可以显示中文名，你的站点文本文件中不是有中文站名吗。用Excel编辑成MeteoInfo的Lon/Lat Station文本 ...

我尝试用excel转换了一下格式，但是在meteoinfo中打开数据出错咯。呵呵。我上传一下我的数据文件。
又麻烦您了。

==================
由于站点数据涉及国家保密数据，问题解决后我删除了原始站点数据

MeteoInfo

孤蓝et 发表于 2011-9-5 10:54

                               
登录/注册后可看大图

我尝试用excel转换了一下格式，但是在meteoinfo中打开数据出错咯。呵呵。我上传一下我的数据文件。
又 ...

大家尝试自己转换一下吧，这个涉及了国家保密数据，因此不便在这里公开分享，请大家的理解

孤蓝et

      二次开发的讨论还是太少？惆怅，这C#还挺难入门。。。怪不得叫C#，横竖都是二。唉。。单纯发个牢骚。
4、读入数组的使用，stationmodel和stationdata的使用？以及等值线分析方法。
      今天看了许久的demo和meteoinfo class library help，忙活了一晚上也没有实现我需要的结果，为了实现打开micaps第一类数据后，自动分析温度等值线。我想参照TSMI_LLStShaded_Click，TSMI_MICAPS1_StationModel_Click这两部分的代码，看能不能合并实现我的想法，可是在数据读入就发现很没有头绪：
micaps 第一类数据的读入方法
stationModelData = aDataInfo.GetStationModelData(ref aExtent);
经纬度的数据
StationData stationData = aDataInfo.GetStationData(0);
      我怎么能从stationModelData = aDataInfo.GetStationModelData(ref aExtent);里面找到经纬度、温度数据在数组的位置，然后赋值给S = ContourDraw.FilterDiscreteData_Radius(S, gridInterp.Radius,gridInterp.GridDataParaV.dataExtent, stationData.UNDEF);用来绘图？
      
      莫名，说明我还没有彻底看透C#中的数组赋值问题，同样也没有真正读懂里面的程序思路。求版主的解答～_～，万分感谢啊。。。

1.         private void TSMI_LLStShaded_Click(object sender, EventArgs e)
   
2.         {
   
3.             this.Cursor = Cursors.WaitCursor;
   
4. 
   
5.             //Read data info
   
6.             LonLatStationDataInfo aDataInfo = new LonLatStationDataInfo();
   
7.             string aFile = Application.StartupPath + "\\Sample\\rain_2008072220.csv";
   
8.             aDataInfo.ReadDataInfo(aFile);
   
9. 
   
10.             //Get station data
    
11.             StationData stationData = aDataInfo.GetStationData(0);
    
12. 
    
13.             //Interpolate
    
14.             GridDataPara aGDP = new GridDataPara();
    
15.             aGDP.dataExtent.minX = 60;
    
16.             aGDP.dataExtent.maxX = 140;
    
17.             aGDP.dataExtent.minY = -20;
    
18.             aGDP.dataExtent.maxY = 60;
    
19.             aGDP.xNum = 80;
    
20.             aGDP.yNum = 80;
    
21.             GridInterpolation gridInterp = new GridInterpolation();
    
22.             gridInterp.GridDataParaV = aGDP;
    
23. 
    
24.             gridInterp.GridInterMethodV = GridInterMethod.IDW_Radius;
    
25.             gridInterp.Radius = 2;
    
26.             gridInterp.MinPointNum = 1;
    
27. 
    
28.             double[] X = new double[1];
    
29.             double[] Y = new double[1];
    
30.             ContourDraw.CreateGridXY(gridInterp.GridDataParaV, ref X, ref Y);
    
31.             double[,] S = stationData.Data;
    
32.             S = ContourDraw.FilterDiscreteData_Radius(S, gridInterp.Radius,
    
33.                 gridInterp.GridDataParaV.dataExtent, stationData.UNDEF);
    
34.             GridData gridData = ContourDraw.InterpolateDiscreteData_Radius(S,
    
35.                 X, Y, gridInterp.MinPointNum, gridInterp.Radius, stationData.UNDEF);
    
36. 
    
37.             //Create legend scheme
    
38.             bool hasNoData = true;
    
39.             LegendScheme aLS = LegendManage.CreateLegendSchemeFromGridData(gridData, LegendType.GraduatedColor,
    
40.                 ShapeTypes.Polygon, ref hasNoData);
    
41.             ((PolygonBreak)aLS.breakList[0]).DrawFill = false;
    
42. 
    
43.             //Create layer
    
44.             VectorLayer aLayer = new VectorLayer(ShapeTypes.Polygon);
    
45.             aLayer = DrawMeteoData.CreateShadedLayer(gridData, aLS, "Rain");
    
46.             aLayer.IsMaskout = true;
    
47.             G_LayerLegend.AddLayer(aLayer, true);
    
48. 
    
49.             this.Cursor = Cursors.Default;
    
50.         }
    

复制代码

MeteoInfo

孤蓝et 发表于 2011-9-11 01:47

                               
登录/注册后可看大图

二次开发的讨论还是太少？惆怅，这C#还挺难入门。。。怪不得叫C#，横竖都是二。唉。。单纯发个牢骚。 ...

如果你只是要对micaps第一类数据中的某个变量进行操作，不要用MICAPS1DataInfo的GetStationModelData方法（这个是为了站点填图用的），用GetStationData(vIdx)就可以了，参数是变量的次序（也就是第几个变量）。变量名称及次序可以在MeteoInfo桌面软件中打开一个micaps第一类数据查一下。后续的插值及生成等值线图层就和TSMI_LLStShaded_Click中的代码没什么两样了。

孤蓝et

继续研究地名显示。