基于AE+C#的栅格IDW插值的实现

  查阅了许多参考资料，代码五花八门，能用的确少之又少，即便百般调试之下也会出些意料之外且无能为力的BUG,其中艰难不足为外人道也。但是最终仍是实现了插值的效果，于是写下这篇博文一是给有需求的小伙伴们一点参考，二是纪念下自己在这个点投入的并不对称的汗水和收获。。。。

 在我这个系统中的IDW插值实现,主要可以分为两个部分，一是对目标点图层进行插值，二是对插值结果（栅格图层）进行栅格渲染（代码会在下面给出）。在这之前先给大家提供实现此功能所需添加的引用，每次参考哪些大神的博客时，都在在这个上面话不少功夫。

using ESRI.ArcGIS.Carto;
using ESRI.ArcGIS.Geodatabase;
using ESRI.ArcGIS.DataSourcesFile;
using Microsoft.VisualBasic;

using ESRI.ArcGIS.GeoAnalyst;

  进行插值时主要可分为以下几个步骤，一是创建目标点图层的要素图层，IDW插值的环境设置，包括像元大小、搜索半径等（如果程序在这里报错，建议将目标点图层在ARCMAP中先利用IDW插值功能运行一遍，如果能成功出图，记录下各类参数），代码如下：

#region 插值
            IInterpolationOp pInterpolationOp;
            pInterpolationOp = new RasterInterpolationOpClass();
            IFeatureClass pFeatureClass;
            IFeatureLayer pFeaturelayer;
            pFeaturelayer = this.axMapControl1.Map.get_Layer(0) as IFeatureLayer;//点图层在axMapControl控件的第一个，索引为0
            pFeatureClass = pFeaturelayer.FeatureClass; 
            IRasterRadius pRadius;
            pRadius = new RasterRadiusClass();
            object maxDistance = Type.Missing;
            pRadius.SetVariable(12, ref maxDistance);//定义搜索半径，本例为12
            IFeatureClassDescriptor pFCDescriptor;
            pFCDescriptor = new FeatureClassDescriptorClass();
            pFCDescriptor.Create(pFeatureClass, null, "value");//设置进行插值的属性字段
            object dCellSize = 113.039027413432;//设置像元值（参考ARCMAP)
            IRasterAnalysisEnvironment pEnv;
            pEnv = new RasterAnalysis();
            pEnv = pInterpolationOp as IRasterAnalysisEnvironment;
            pEnv.SetCellSize(esriRasterEnvSettingEnum.esriRasterEnvValue, ref dCellSize);
            object objectbarrier = Type.Missing;
            IGeoDataset rasDataset;
            rasDataset = pInterpolationOp.IDW((IGeoDataset)pFCDescriptor,
                2, pRadius, ref objectbarrier);
            IRaster pOutRaster;
            pOutRaster = rasDataset as IRaster;
            IRasterLayer pOutRasLayer;
            pOutRasLayer = new RasterLayerClass();
            pOutRasLayer.CreateFromRaster(pOutRaster);
            this.axMapControl1.AddLayer(pOutRasLayer, 0);// 输出结果至axMapControl1,并放置在第一层
            axMapControl1.ActiveView.Refresh();//刷新视图
            DialogResult dr = MessageBox.Show("插值成功,请选择是否进行分级渲染", "分级渲染选择",

           MessageBoxButtons.OKCancel, MessageBoxIcon.Information);//建立一个文件对话框的实例，用于接收用户指令。

  如果用户选择进行分级渲染，则进入一下代码模块：

 if (dr == DialogResult.OK)
            {
                //用户选择确认的操作（分级渲染）
                IRasterLayer pRasterLayer = axMapControl1.get_Layer(i) as IRasterLayer;//获取需要分级渲染的栅格图层，i为图层在axMapControl中的索引i
                int number = Convert.ToInt32(Interaction.InputBox("请输入栅格影像分类数量(默认为10)",
                    "字符串", "", 500, 250)); //调用VB的输入框获取用户需要的分级层数，默认为10
                if (number == 0)
                    number = 10;
                funColorForRaster_Classify(pRasterLayer, number);

                axMapControl1.ActiveView.Refresh();

             } 

    funColorForRaster_Classify(IRasterLayer pRasterLayer,int number)是一个自定义的方法，用于为输入的栅格图层按照指定的分级数量进行分级渲染，其代码如下：

            IRasterClassifyColorRampRenderer pRClassRend = new
            RasterClassifyColorRampRenderer() as IRasterClassifyColorRampRenderer;
            IRasterRenderer pRRend = pRClassRend as IRasterRenderer;
            IRaster pRaster = pRasterLayer.Raster;//栅格图层转换为栅格图
            IRasterBandCollection pRBandCol = pRaster as IRasterBandCollection;
            IRasterBand pRBand = pRBandCol.Item(0);
            if (pRBand.Histogram == null)
            {
                pRBand.ComputeStatsAndHist();
            }
            pRRend.Raster = pRaster;
            pRClassRend.ClassCount = number;//设置分级层次，默认为10
            pRRend.Update();
            //设置初始渐变色和结束渐变色
            MessageBox.Show("请选择初始颜色", "颜色选择", MessageBoxButtons.OK, MessageBoxIcon.Information);
            IRgbColor pFromColor = new RgbColor() as IRgbColor;
            selectColor(pFromColor);//自定义的一个方法，用于打开颜色选取对话框，为IRgbColor类的实例赋值
            MessageBox.Show("请选择结束颜色", "颜色选择", MessageBoxButtons.OK, MessageBoxIcon.Information);
            IRgbColor pToColor = new RgbColor() as IRgbColor;
            selectColor(pToColor);
            //进行栅格渲染并输出至axMapControl
            IAlgorithmicColorRamp colorRamp = new AlgorithmicColorRamp() as
            IAlgorithmicColorRamp;
            colorRamp.Size = 10;
            colorRamp.FromColor = pFromColor;
            colorRamp.ToColor = pToColor;
            bool createColorRamp;
            colorRamp.CreateRamp(out createColorRamp);
            IFillSymbol fillSymbol = new SimpleFillSymbol() as IFillSymbol;
            for (int i = 0; i < pRClassRend.ClassCount; i++)
            {
                fillSymbol.Color = colorRamp.get_Color(i);
                pRClassRend.set_Symbol(i, fillSymbol as ISymbol);
                pRClassRend.set_Label(i,
                pRClassRend.get_Break(i).ToString("0.00"));
            }
            pRasterLayer.Renderer = pRRend;
            IActiveView activeview = this.axMapControl1.ActiveView;
            activeview.ContentsChanged();
            //刷新窗口
            activeview.PartialRefresh(esriViewDrawPhase.esriViewGeography, null,
            null);

     selectColor(IRgbColor)的代码如下：

 //打开颜色选择对话框 并未IRGB类的三原色赋值
        private void selectColor(IRgbColor selColor)
        {
            ColorDialog colorDialog = new ColorDialog();
            colorDialog.AllowFullOpen = true;
            colorDialog.FullOpen = true;
            colorDialog.ShowHelp = true;
            colorDialog.Color = Color.Black;  //获取或设置用户所选定的颜色  
            colorDialog.ShowDialog();
            selColor.Red = colorDialog.Color.R;
            selColor.Green = colorDialog.Color.G;
            selColor.Blue = colorDialog.Color.B;
            MessageBox.Show("您选择的各颜色分量值分别为" + "\nR:" + colorDialog.Color.R.ToString()
                + "\nG:" + colorDialog.Color.G.ToString() + "\nB:" + colorDialog.Color.B.ToString());

        }

  本次分享到此结束，希望对有相关需要的小伙伴们有所帮助。