格点插值到站点的问题

zhksoft

1.             string m4file = @"T:\17010208.000";
   
2.             MICAPS4DataInfo m4data = new MICAPS4DataInfo();
   
3.             m4data.ReadDataInfo(m4file);
   
4.             richTextBox1.AppendText(m4data.Description.ToString());
   
5.             
   
6.             string sta_file = Path.Combine(AppDomain.CurrentDomain.BaseDirectory, @"altsta.csv");
   
7.             string out_file=Path.Combine(AppDomain.CurrentDomain.BaseDirectory, @"tmp.csv");
   
8.            
   
9.             GridData gd=new GridData();
   
10.             gd.Data =m4data.GridData;
    
11.             gd.ToStation(sta_file,out_file);
    
12.             //((GridData)m4data.GridData).ToStation(sta_file, out_file);

复制代码

@MeteoInfo
王老师，我想用格点插到站点，可是似乎MICAPS4DataInfo 和GridData之间传递数据出了问题，用类型强制转换也不行，请问有没有快捷的办法把MICAPS4DataInfo里的GridData 传到GridData变量里？  还有一个问题就是第2类和第11类的UV要插到站点上成为站点的风向和风速，是先把U和V分别插到站点然后通过站点的U和V用向量的办法计算出风向和风速，还是有其它的快捷的方法呢？

MeteoInfo

参考这个帖子：格点插值为站点数据批量处理
http://bbs.06climate.com/forum.p ... 191&fromuid=106
(出处: 气象家园)

zhksoft

MeteoInfo 发表于 2017-1-3 18:05
参考这个帖子：格点插值为站点数据批量处理
http://bbs.06climate.com/forum.php?mod=viewthread&tid=3519 ...

好的，谢谢。

孤蓝et

MeteoInfo 发表于 2017-1-3 18:05
参考这个帖子：格点插值为站点数据批量处理
http://bbs.06climate.com/forum.php?mod=viewthread&tid=3519 ...

王老师，我也遇到类似问题，
读取micaps第四类数据，想获取距离站点最近的格点数据，这个我用meteoInfo Lab可以实现，那在C#类库可以做吗？我需要用c#类库实现这个功能，盼指教。

孤蓝et

1.         /// <summary>
   
2.         /// Interpolate grid data to a station point
   
3.         /// </summary>
   
4.         /// <param name="x">x coordinate of the station</param>
   
5.         /// <param name="y">y coordinate of the station</param>
   
6.         /// <returns>interpolated value</returns>
   
7.         public double ToStation(double x, double y)
   
8.         {
   
9.             double iValue = MissingValue;
   
10.             if (x < X[0] || x > X[XNum - 1] || y < Y[0] || y > Y[YNum - 1])
    
11.                 return iValue;
    
12. 
    
13.             //Get x/y index
    
14.             int xIdx = 0, yIdx = 0;
    
15.             xIdx = (int)((x - X[0]) / XDelt);
    
16.             yIdx = (int)((y - Y[0]) / YDelt);

复制代码

在类库里面找到上面一段代码，这里面的xIdx   yIdx是不是结合网格点数据及经纬网格间隔，直接找到距离站点经纬度最近的序号？那我在这里，增加一个函数，不去做插值，直接输出value，就是距离站点最近的格点数据了吧？

孤蓝et

老师的dll的确很有用。