利用Cressman插值算法 生成surfer能用的.grd文件

changsha2012

求教版主，我想利用带经纬度的站点文件，通过Cressman插值算法 生成.grd文件。

现在我对 radlist这个参数不明白，另外如何导出一个插值完的.grd文件呢。

changsha2012

已经解决，有空补上源码。

mofangbao

changsha2012 发表于 2012-2-4 18:04

                               
登录/注册后可看大图

已经解决，有空补上源码。

期待，不要忘记哈

♂雨已~○

我也期待，先在这谢过啦哈

changsha2012

/// <summary>
        /// 由csv/txt监测站点的经纬度数据插值生成Grid文件
        /// </summary>
        /// <param name="csvFile">csv文件(或者txt)文件，要求2，3，4列依次为经度，纬度，插值数据</param>
        private string CsvToGrid(string csvFile)
        {
            LonLatStationDataInfo aDataInfo = new LonLatStationDataInfo();
            aDataInfo.ReadDataInfo(csvFile);
            //Get station data
            StationData stationData = aDataInfo.GetStationData(0);
            //Interpolate
            //范围
            Extent dataExtent = new Extent(106, 115, 24, 33);
            GridDataPara aGDP = new GridDataPara();
            aGDP.dataExtent = dataExtent;
            //格点数量
            aGDP.xNum = 300;
            aGDP.yNum = 300;
            //迭代半径
            List<double> radlist = new List<double>();
            radlist.Add(30);
            GridInterpolation gridInterp = new GridInterpolation();
            gridInterp.GridDataParaV = aGDP;
            gridInterp.GridInterMethodV = GridInterMethod.Cressman;
            gridInterp.MinPointNum = 1;
            gridInterp.RadList = radlist;
            gridInterp.UnDefData = 0.001;

            double[] X = new double[1];
            double[] Y = new double[1];
            ContourDraw.CreateGridXY(gridInterp.GridDataParaV, ref X, ref Y);
            double[,] S = stationData.Data;

            GridData gridData = new GridData();

            //Cressman
            gridData = ContourDraw.InterpolateDiscreteData_Cressman(S, X, Y, gridInterp.UnDefData, gridInterp.RadList);
            
            //导出.grd文件
            string outputGridFile = System.IO.Path.ChangeExtension(csvFile, ".grd");
            gridData.SaveAsSurferASCIIFile(outputGridFile);
            return outputGridFile;
        }

MeteoInfo

changsha2012 发表于 2012-2-9 12:24

                               
登录/注册后可看大图

///
        /// 由csv/txt监测站点的经纬度数据插值生成Grid文件
        ///

相当的绕，看来一些库函数做的还不够简洁

changsha2012

呵呵,感觉没有Surfer那么便捷~

daohua201314

谢谢楼主分享

llsshh1985

UnDefData 这个值一般设成什么啊？我想让降雨量为0的点也发挥影响。之前我把UnDefData这个设为0，好像0的站点没有发挥作用。用的IDW_Radius插值

Ocean_Baby

surfer 插值有没有IDW方法啊？