利用meteoinfo实现Surfer二次开发效果，绘制等值线图，并白化（maskout）边框外内容

MeteoInfo · 发表于 2013-1-1 21:32:59

记忆碎片发表于 2013-1-1 01:12

登录/注册后可看大图

经过验证，国内网站上流传“城市经纬度”的数据本身质量有问题
我在http://www.diva-gis.org/gdata 上，选 ...

建议先在网上搜搜何为shape文件

抠数据有很多方法，手动选择、按属性选择、按位置选择等。

记忆碎片 · 发表于 2013-1-5 16:20:23

最后确认下来还是数据源的问题
找到了正确的数据，预期效果很快就出来了
上图，交作业
浙江叠加图.jpg

记忆碎片 · 发表于 2013-1-24 11:04:50

遇到一个新的问题
之前调试完好的内容，现在把核心的MeteoInfoC.dll从1.0.5升级到1.0.8，仅对部分属性名称做了对应调整。没有进行功能性修改
发现报错，而且报错很莫名其妙,说wContour缺失，经过查看，该文件存在的。
不知道如何解决好？
报错.jpg

MeteoInfo · 发表于 2013-1-24 11:44:36

记忆碎片发表于 2013-1-24 11:04

登录/注册后可看大图

遇到一个新的问题
之前调试完好的内容，现在把核心的MeteoInfoC.dll从1.0.5升级到1.0.8，仅对部分属性名称 ...

wContour.dll也需要用最新的文件。

记忆碎片 · 发表于 2013-1-25 12:06:11

MeteoInfo 发表于 2013-1-24 11:44

登录/注册后可看大图

wContour.dll也需要用最新的文件。

请问与1.0.8对应的MeteoInfo.dll对应的应该是1.6.0.0的wContour?
还对？

MeteoInfo · 发表于 2013-1-25 14:15:16

记忆碎片发表于 2013-1-25 12:06

登录/注册后可看大图

请问与1.0.8对应的MeteoInfo.dll对应的应该是1.6.0.0的wContour?
还对？

是的

记忆碎片 · 发表于 2013-2-3 00:18:18

通过如下代码得到上图效果图（源文件：

Ng值文件.txt (4.59 KB, 下载次数: 1)

2013-2-3 00:14 上传

点击文件名下载附件），可是效果很奇怪，照理应该是个方形或者圆形，至少不至于不连续吧（理想预期效果图本地贴顶楼有）？请老师诊断下问题是否处在插值参数设置不合理？还是问题处在什么地方？

this.Cursor = Cursors.WaitCursor;
//DrawStationShaded();
string sourceFile = AppDomain.CurrentDomain.BaseDirectory + @"\Sample\Ng值文件.txt";
DrawCountourLine.DrawStationShadedNew(layersLegend1, mapLayout1, sourceFile);
this.Cursor = Cursors.Default;

/// <summary>
/// 需要原文件格式“stationID,Longitude,Latitude,Value”
/// </summary>
/// <param name="_targetLayersLegend"></param>
/// <param name="_targetMapLayout"></param>
/// <param name="fileSource"></param>
public static void DrawStationShadedNew(LayersLegend _targetLayersLegend,MapLayout _targetMapLayout,string fileSource)
{
//Read data info
MeteoDataInfo aDataInfo = new MeteoDataInfo();
string aFile = fileSource;
aDataInfo.OpenLonLatData(aFile);
//Get station data
StationData stationData = aDataInfo.GetStationData("Value");
//Interpolate
GridDataSetting aGDP = new GridDataSetting();
aGDP.DataExtent.minX = 60;
aGDP.DataExtent.maxX = 140;
aGDP.DataExtent.minY = -20;
aGDP.DataExtent.maxY = 60;
aGDP.XNum = 80;
aGDP.YNum = 80;
InterpolationSetting gridInterp = new InterpolationSetting();
gridInterp.GridDataSet = aGDP;
gridInterp.InterpolationMethod = InterpolationMethods.IDW_Neighbors;
gridInterp.Radius = 2;
gridInterp.MinPointNum = 1;
double[] X = new double[1];
double[] Y = new double[1];
ContourDraw.CreateGridXY(gridInterp.GridDataSet, ref X, ref Y);
double[,] S = stationData.Data;
S = ContourDraw.FilterDiscreteData_Radius(S, gridInterp.Radius,
gridInterp.GridDataSet.DataExtent, stationData.UNDEF);
GridData gridData = ContourDraw.InterpolateDiscreteData_Radius(S,
X, Y, gridInterp.MinPointNum, gridInterp.Radius, stationData.UNDEF);
//Create legend scheme
bool hasNoData = true;
LegendScheme aLS = LegendManage.CreateLegendSchemeFromGridData(gridData, LegendType.GraduatedColor,
ShapeTypes.Polygon, ref hasNoData);
((PolygonBreak)aLS.LegendBreaks[0]).DrawFill = false;
//Create layer
VectorLayer aLayer = new VectorLayer(ShapeTypes.Polygon);
aLayer = DrawMeteoData.CreateShadedLayer(gridData, aLS, "Ng", "Ng");
//aLayer.IsMaskout = true;
//Add layer
_targetLayersLegend.ActiveMapFrame.AddLayer(aLayer);
_targetLayersLegend.ActiveMapFrame.MoveLayer(aLayer.Handle, 0);
_targetLayersLegend.Refresh();
//Change title of the layout
LayoutGraphic aTitle = _targetMapLayout.GetTexts()[0];
aTitle.SetLabelText("Station Ng Layer");
//Add or change the legend in layout
LayoutLegend aLegend;
if (_targetMapLayout.GetLegends().Count > 0)
aLegend = _targetMapLayout.GetLegends()[0];
else
aLegend = _targetMapLayout.AddLegend(650, 100);
aLegend.LegendStyle = LegendStyleEnum.Bar_Vertical;
aLegend.LegendLayer = aLayer;
//if (tabControl1.SelectedIndex == 1)
_targetMapLayout.PaintGraphics();
}

/// <summary> 
        /// 需要原文件格式“stationID,Longitude,Latitude,Value” 
        /// </summary> 
        /// <param name="_targetLayersLegend"></param> 
        /// <param name="_targetMapLayout"></param> 
        /// <param name="fileSource"></param> 
        public static void DrawStationShadedNew(LayersLegend _targetLayersLegend,MapLayout _targetMapLayout,string fileSource) 
        { 
            //Read data info 
            MeteoDataInfo aDataInfo = new MeteoDataInfo(); 
            string aFile = fileSource; 
            aDataInfo.OpenLonLatData(aFile); 
 
            //Get station data 
            StationData stationData = aDataInfo.GetStationData("Value"); 
 
            //Interpolate 
            GridDataSetting aGDP = new GridDataSetting(); 
            aGDP.DataExtent.minX = 60; 
            aGDP.DataExtent.maxX = 140; 
            aGDP.DataExtent.minY = -20; 
            aGDP.DataExtent.maxY = 60; 
            aGDP.XNum = 80; 
            aGDP.YNum = 80; 
            InterpolationSetting gridInterp = new InterpolationSetting(); 
            gridInterp.GridDataSet = aGDP; 
 
            gridInterp.InterpolationMethod = InterpolationMethods.IDW_Neighbors; 
            gridInterp.Radius = 2; 
            gridInterp.MinPointNum = 1; 
 
            double[] X = new double[1]; 
            double[] Y = new double[1]; 
            ContourDraw.CreateGridXY(gridInterp.GridDataSet, ref X, ref Y); 
            double[,] S = stationData.Data; 
            S = ContourDraw.FilterDiscreteData_Radius(S, gridInterp.Radius, 
                gridInterp.GridDataSet.DataExtent, stationData.UNDEF); 
            GridData gridData = ContourDraw.InterpolateDiscreteData_Radius(S, 
                X, Y, gridInterp.MinPointNum, gridInterp.Radius, stationData.UNDEF); 
 
            //Create legend scheme 
            bool hasNoData = true; 
            LegendScheme aLS = LegendManage.CreateLegendSchemeFromGridData(gridData, LegendType.GraduatedColor, 
                ShapeTypes.Polygon, ref hasNoData); 
            ((PolygonBreak)aLS.LegendBreaks[0]).DrawFill = false; 
 
            //Create layer 
            VectorLayer aLayer = new VectorLayer(ShapeTypes.Polygon); 
            aLayer = DrawMeteoData.CreateShadedLayer(gridData, aLS, "Ng", "Ng"); 
            //aLayer.IsMaskout = true; 
 
            //Add layer 
            _targetLayersLegend.ActiveMapFrame.AddLayer(aLayer); 
            _targetLayersLegend.ActiveMapFrame.MoveLayer(aLayer.Handle, 0); 
            _targetLayersLegend.Refresh(); 
 
            //Change title of the layout 
            LayoutGraphic aTitle = _targetMapLayout.GetTexts()[0]; 
            aTitle.SetLabelText("Station Ng Layer"); 
 
            //Add or change the legend in layout 
            LayoutLegend aLegend; 
            if (_targetMapLayout.GetLegends().Count > 0) 
                aLegend = _targetMapLayout.GetLegends()[0]; 
            else 
                aLegend = _targetMapLayout.AddLegend(650, 100); 
            aLegend.LegendStyle = LegendStyleEnum.Bar_Vertical; 
            aLegend.LegendLayer = aLayer; 
            //if (tabControl1.SelectedIndex == 1) 
                _targetMapLayout.PaintGraphics(); 
        }

MeteoInfo · 发表于 2013-2-3 10:10:14

记忆碎片发表于 2013-2-3 00:18

登录/注册后可看大图

通过如下代码得到上图效果图（源文件：），可是效果很奇怪，照理应该是个方形或者圆形，至少不至于不连续 ...

看了一下，你的数据其实是格点数据，建议你把它改成格点数据格式，比如GrADS的二进制，或者简单的Surfer文本格点格式。

你目前用站点格式处理，在绘制等值线的时候需要先插值成格点数据，插值都是一些数学假设，且参数设置不同结果会完全不同，显然你的插值参数设置是有问题的。

记忆碎片 · 发表于 2013-2-3 14:28:47

本帖最后由记忆碎片于 2013-2-3 14:31 编辑

MeteoInfo 发表于 2013-2-3 10:10

登录/注册后可看大图

看了一下，你的数据其实是格点数据，建议你把它改成格点数据格式，比如GrADS的二进制，或者简单的Surfer文 ...

果然是问题出在插值参数设置上，按照上述的代码，将差值部分代码修改为下示代码，插值后效果好多了就。
代码：

/// <summary>
/// 需要原文件格式“stationID,Longitude,Latitude,Value”
/// </summary>
/// <param name="_targetLayersLegend"></param>
/// <param name="_targetMapLayout"></param>
/// <param name="fileSource"></param>
public static void DrawStationShadedNew(LayersLegend _targetLayersLegend,MapLayout _targetMapLayout,string fileSource)
{
//Read data info
MeteoDataInfo aDataInfo = new MeteoDataInfo();
string aFile = fileSource;
aDataInfo.OpenLonLatData(aFile);
//Get station data
StationData stationData = aDataInfo.GetStationData("Value");
//Interpolate
GridDataSetting aGDP = new GridDataSetting();
aGDP.DataExtent.minX = 121;//x最小值
aGDP.DataExtent.maxX = 122;//x最大值
aGDP.DataExtent.minY = 28;//y最小值
aGDP.DataExtent.maxY = 29;//y最大值
aGDP.XNum = 81;//列数
aGDP.YNum = 81;//行数
InterpolationSetting gridInterp = new InterpolationSetting();
gridInterp.GridDataSet = aGDP;
gridInterp.InterpolationMethod = InterpolationMethods.IDW_Neighbors;//差值方法
gridInterp.Radius = 3;//差值半径
gridInterp.MinPointNum = 1;
double[] X = new double[1];
double[] Y = new double[1];
ContourDraw.CreateGridXY(gridInterp.GridDataSet, ref X, ref Y);
double[,] S = stationData.Data;
S = ContourDraw.FilterDiscreteData_Radius(S, gridInterp.Radius,
gridInterp.GridDataSet.DataExtent, stationData.UNDEF);
GridData gridData = ContourDraw.InterpolateDiscreteData_Radius(S,
X, Y, gridInterp.MinPointNum, gridInterp.Radius, stationData.UNDEF);
//Create legend scheme
bool hasNoData = true;
LegendScheme aLS = LegendManage.CreateLegendSchemeFromGridData(gridData, LegendType.GraduatedColor,
ShapeTypes.Polygon, ref hasNoData);
((PolygonBreak)aLS.LegendBreaks[0]).DrawFill = false;
//Create layer
VectorLayer aLayer = new VectorLayer(ShapeTypes.Polygon);
aLayer = DrawMeteoData.CreateShadedLayer(gridData, aLS, "Ng", "Ng");
aLayer.IsMaskout = true;
//Add layer
_targetLayersLegend.ActiveMapFrame.AddLayer(aLayer);
_targetLayersLegend.ActiveMapFrame.MoveLayer(aLayer.Handle, 0);
_targetLayersLegend.Refresh();
//Change title of the layout
LayoutGraphic aTitle = _targetMapLayout.GetTexts()[0];
aTitle.SetLabelText("Station Ng Layer");
//Add or change the legend in layout
LayoutLegend aLegend;
if (_targetMapLayout.GetLegends().Count > 0)
aLegend = _targetMapLayout.GetLegends()[0];
else
aLegend = _targetMapLayout.AddLegend(650, 100);
aLegend.LegendStyle = LegendStyleEnum.Bar_Vertical;
aLegend.LegendLayer = aLayer;
//if (tabControl1.SelectedIndex == 1)
_targetMapLayout.PaintGraphics();
}

/// <summary> 
        /// 需要原文件格式“stationID,Longitude,Latitude,Value” 
        /// </summary> 
        /// <param name="_targetLayersLegend"></param> 
        /// <param name="_targetMapLayout"></param> 
        /// <param name="fileSource"></param> 
        public static void DrawStationShadedNew(LayersLegend _targetLayersLegend,MapLayout _targetMapLayout,string fileSource) 
        { 
            //Read data info 
            MeteoDataInfo aDataInfo = new MeteoDataInfo(); 
            string aFile = fileSource; 
            aDataInfo.OpenLonLatData(aFile); 
 
            //Get station data 
            StationData stationData = aDataInfo.GetStationData("Value"); 
 
            //Interpolate 
            GridDataSetting aGDP = new GridDataSetting(); 
            aGDP.DataExtent.minX = 121;//x最小值 
            aGDP.DataExtent.maxX = 122;//x最大值 
            aGDP.DataExtent.minY = 28;//y最小值 
            aGDP.DataExtent.maxY = 29;//y最大值 
            aGDP.XNum = 81;//列数 
            aGDP.YNum = 81;//行数 
            InterpolationSetting gridInterp = new InterpolationSetting(); 
            gridInterp.GridDataSet = aGDP; 
 
            gridInterp.InterpolationMethod = InterpolationMethods.IDW_Neighbors;//差值方法 
            gridInterp.Radius = 3;//差值半径 
            gridInterp.MinPointNum = 1; 
 
            double[] X = new double[1]; 
            double[] Y = new double[1]; 
            ContourDraw.CreateGridXY(gridInterp.GridDataSet, ref X, ref Y); 
            double[,] S = stationData.Data; 
            S = ContourDraw.FilterDiscreteData_Radius(S, gridInterp.Radius, 
                gridInterp.GridDataSet.DataExtent, stationData.UNDEF); 
            GridData gridData = ContourDraw.InterpolateDiscreteData_Radius(S, 
                X, Y, gridInterp.MinPointNum, gridInterp.Radius, stationData.UNDEF); 
 
            //Create legend scheme 
            bool hasNoData = true; 
            LegendScheme aLS = LegendManage.CreateLegendSchemeFromGridData(gridData, LegendType.GraduatedColor, 
                ShapeTypes.Polygon, ref hasNoData); 
            ((PolygonBreak)aLS.LegendBreaks[0]).DrawFill = false; 
 
            //Create layer 
            VectorLayer aLayer = new VectorLayer(ShapeTypes.Polygon); 
            aLayer = DrawMeteoData.CreateShadedLayer(gridData, aLS, "Ng", "Ng"); 
            aLayer.IsMaskout = true; 
 
            //Add layer 
            _targetLayersLegend.ActiveMapFrame.AddLayer(aLayer); 
            _targetLayersLegend.ActiveMapFrame.MoveLayer(aLayer.Handle, 0); 
            _targetLayersLegend.Refresh(); 
 
            //Change title of the layout 
            LayoutGraphic aTitle = _targetMapLayout.GetTexts()[0]; 
            aTitle.SetLabelText("Station Ng Layer"); 
 
            //Add or change the legend in layout 
            LayoutLegend aLegend; 
            if (_targetMapLayout.GetLegends().Count > 0) 
                aLegend = _targetMapLayout.GetLegends()[0]; 
            else 
                aLegend = _targetMapLayout.AddLegend(650, 100); 
            aLegend.LegendStyle = LegendStyleEnum.Bar_Vertical; 
            aLegend.LegendLayer = aLayer; 
            //if (tabControl1.SelectedIndex == 1) 
                _targetMapLayout.PaintGraphics(); 
        }

插值图层后，再Maskout，然后效果图如下：

插值效果图

差值细节图

接下来的待完成工作：
1 自动化找到数据中的行数、列数、经纬度最大、最小值。

等把这部分功能完成后。抽空把内容整合起来分享给大家。现在虽然方法都有了，不过太零散呵呵再次感谢王老师一路的指导。

qq200602116 · 发表于 2014-5-21 09:05:17

记忆碎片发表于 2013-1-5 16:20
最后确认下来还是数据源的问题
找到了正确的数据，预期效果很快就出来了
上图，交作业

请问您的数据在那找到的，谢谢

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