meteoinfo地面填图问题

♂雨已~○ · 发表于 2013-1-25 20:21:09

MeteoInfo 发表于 2013-1-25 19:09

登录/注册后可看大图

数据有经纬度，可以用判断点是否在某区域内的函数来筛选，参考此帖：MeteoInfo脚本示例：创建点图层并判断 ...

谢谢王老师，问题已经解决！实现代码如下：
  LayerNode ID;
         ID = layersLegend1.ActiveMapFrame.GetLayerNodeByName("西南地区.shp");
         VectorLayer vLayer=(VectorLayer)ID.MapLayer;
         PointD aPoint = new PointD();
         int m = 0;
         for(int i=0;i<vLayer.ShapeList.Count;i++)
         {
            PolygonShape aPGS=(PolygonShape)vLayer.ShapeList;
               for (int j = 0; j < stationModelData1.GetLength(1);j++)
               {
                  aPoint.X=stationModelData1[0, j];
                  aPoint.Y=stationModelData1[1, j];
                     if (GeoComputation.PointInPolygon(aPGS, aPoint))
                     {
                        m++;
                     }
               }
         }
         double[,] stationModelData = new double[10, m];
         m = 0;
         for (int i = 0; i < vLayer.ShapeList.Count; i++)
         {
            PolygonShape aPGS = (PolygonShape)vLayer.ShapeList;
            for (int j = 0; j < stationModelData1.GetLength(1); j++)
            {
                  aPoint.X = stationModelData1[0, j];
                  aPoint.Y = stationModelData1[1, j];
                  if (GeoComputation.PointInPolygon(aPGS, aPoint))
                  {
                     for(int n=0;n<10;n++)
                  {
                     stationModelData[n,m] = stationModelData1[n,j];
                  }
                     m++;
                  }
            }
         }

但是我发现这么效率有点低呢，速度很慢，王老师有没有更好的办法呢？我用了一个比较笨的办法，自己是执行了两次，第一次是为了计算那个数组维数，以便申明提取后的数组呢~~~小子C#刚入门呢~~~还有一个问题王老师，就是我用Meteoinfo站点数据图位置对应不对啊，发现整体数据比底图左偏呢，我用的是micaps第一类的自动站数据画的~~~下面是我用的数据，王老师帮忙看看吧~~~

12081410.AWS (4.04 MB, 下载次数: 1)

MeteoInfo · 发表于 2013-1-25 23:26:40

♂雨已~○ 发表于 2013-1-25 20:21

登录/注册后可看大图

谢谢王老师，问题已经解决！实现代码如下：
LayerNode ID;
ID = layersLegend1.ActiveMap ...

1. 你的程序有两次大的循环，其实一次就可以了，具体你自己琢磨吧。另外需要判断的PolygonShape越多越费时，最好是一个。

2. 是软件的一个bug，谢谢你指出来，已经改了（见置顶帖子的最新文件）。

♂雨已~○ · 发表于 2013-1-26 10:14:24

MeteoInfo 发表于 2013-1-25 23:26

登录/注册后可看大图

1. 你的程序有两次大的循环，其实一次就可以了，具体你自己琢磨吧。另外需要判断的PolygonShape越多越费时 ...

恩~~~好的，王老师我再继续研究研究，嘿嘿~~~

♂雨已~○ · 发表于 2013-1-26 10:33:52

MeteoInfo 发表于 2013-1-25 23:26

登录/注册后可看大图

1. 你的程序有两次大的循环，其实一次就可以了，具体你自己琢磨吧。另外需要判断的PolygonShape越多越费时 ...

王老，还请你讲最新文件传到临时服务器或者你的官方网站上吧，现在在这个上面无法下载（自动回复：请不要使用迅雷等下载工具，使用IE右键另存为就可以）了~~~要现金，而那个临时服务器上还是1月3号的版本~~~麻烦王老师咯

MeteoInfo · 发表于 2013-1-26 17:50:45

♂雨已~○ 发表于 2013-1-26 10:33

登录/注册后可看大图

王老，还请你讲最新文件传到临时服务器或者你的官方网站上吧，现在在这个上面无法下载（自动回复：请不要 ...

放在http://www.meteothinker.com/Downloads.aspx上了，"MeteoInfo Update Files"链接。

♂雨已~○ · 发表于 2013-7-30 18:54:29

MeteoInfo 发表于 2013-1-26 17:50

登录/注册后可看大图

放在http://www.meteothinker.com/Downloads.aspx上了，"MeteoInfo Update Files"链接。

王老师，最近我用meteoinfo打开一种5分钟更新的那种地面自动观测数据一直出错，你调试一下bug哇，数据传上······

1307290620.000 (87.86 KB, 下载次数: 2)

MeteoInfo · 发表于 2013-7-30 23:11:30

♂雨已~○ 发表于 2013-7-30 18:54

登录/注册后可看大图

王老师，最近我用meteoinfo打开一种5分钟更新的那种地面自动观测数据一直出错，你调试一下bug哇，数据传上 ...

已经更新了，见置顶帖子。

♂雨已~○ · 发表于 2013-7-30 23:40:25

MeteoInfo 发表于 2013-7-30 23:11

登录/注册后可看大图

已经更新了，见置顶帖子。

好的，王老师，辛苦了啊嘿嘿

♂雨已~○ · 发表于 2013-8-4 17:53:04

MeteoInfo 发表于 2013-7-30 23:11

登录/注册后可看大图

已经更新了，见置顶帖子。

王老师，meteoinfo能不能开发支持一下NC数据的极坐标格式的显示，数据存储的模块呢？现在有一些雷达或者其他数据都有这种极坐标的NC数据格式，我这有一个nc雷达的极坐标格式的数据。nc格式的变量说明有：
dimensions:
      Azimuth = 865 ;
      Gate = 480 ;
variables:
      float Azimuth(Azimuth) ;
            Azimuth:Units = "Degrees" ;
      float Beamwidth(Azimuth) ;
            Beamwidth:Units = "Degrees" ;
      float GateWidth(Azimuth) ;
            GateWidth:Units = "Meters" ;
      float Filtered_Intensity(Horizontal)(Azimuth, Gate) ;
            Filtered_Intensity(Horizontal):Units = "dBz" ;
// global attributes:
            :TypeName = "Filtered_Intensity(Horizontal)" ;
            :DataType = "RadialSet" ;
            :Latitude = 30.5788993835449 ;
            :Longitude = 103.962112426758 ;
            :Height = 120 ;
            :Time = -1 ;
            :FractionalTime = 0. ;
            :attributes = "NyquistVelocity vcp radarName ColorMap" ;
            :NyquistVelocity-unit = "MetersPerSecond" ;
            :NyquistVelocity-value = "31.62" ;
            :vcp-unit = "dimensionless" ;
            :vcp-value = "jinjin" ;
            :radarName-unit = "dimensionless" ;
            :radarName-value = "9477CHE" ;
            :ColorMap-unit = "dimensionless" ;
            :ColorMap-value = "Reflectivity" ;
            :Elevation = 1.5 ;
            :ElevationUnits = "Degrees" ;
            :MissingData = -99900. ;
            :RangeFolded = -99901. ;
            :RadarParameters = "PRF PulseWidth MaximumRange" ;
            :PRF-unit = "Hertz" ;
            :PRF-value = 1180 ;
            :PulseWidth-unit = "MicroSeconds" ;
            :PulseWidth-value = 0.4 ;
            :MaximumRange-unit = "KiloMeters" ;
            :MaximumRange-value = 120. ;
            :ConversionPlugin = "Vol2NetCDF Conversion for EDGE 5.0" ;
}

nc极坐标雷达数据：

9477CHE--01-ZH.nc (1.6 MB, 下载次数: 1)

MeteoInfo · 发表于 2013-8-4 19:05:46

♂雨已~○ 发表于 2013-8-4 17:53

登录/注册后可看大图

王老师，meteoinfo能不能开发支持一下NC数据的极坐标格式的显示，数据存储的模块呢？现在有一些雷达或者其 ...

挺有意思的问题，不过做数据格式的详细描述是必不可少的。

		自动登录	找回密码
密码			立即注册