meteoinfo地面填图问题

♂雨已~○

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)

2013-1-25 20:17 上传

点击文件名下载附件

MeteoInfo

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

                               
登录/注册后可看大图

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

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

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

♂雨已~○

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

                               
登录/注册后可看大图

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

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

♂雨已~○

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

                               
登录/注册后可看大图

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

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

MeteoInfo

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

                               
登录/注册后可看大图

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

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

♂雨已~○

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

                               
登录/注册后可看大图

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

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

2013-7-30 18:54 上传

点击文件名下载附件

MeteoInfo

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

                               
登录/注册后可看大图

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

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

♂雨已~○

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

                               
登录/注册后可看大图

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

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

♂雨已~○

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)

2013-8-4 17:51 上传

点击文件名下载附件

MeteoInfo

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

                               
登录/注册后可看大图

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

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