样条函数插值

using System;
using System.Collections.Generic;
using System.ComponentModel;
using System.Data;
using System.Drawing;
using System.Text;
using System.Windows.Forms;
using ESRI.ArcGIS.Carto;
using ESRI.ArcGIS.Geodatabase;
using ESRI.ArcGIS.DataSourcesRaster;
using ESRI.ArcGIS.Geometry;
using ESRI.ArcGIS.SpatialAnalyst;
using ESRI.ArcGIS.GeoAnalyst;




namespace 样条函数法插值
{
    public partial class Form1 : Form
    {
        public Form1()
        {
            InitializeComponent();
        }


        private IRasterAnalysisEnvironment rasterEnv;//分析环境
        private IInterpolationOp2 interOp;//空间插值对象
        private IFeatureClass feaClass;
        private IFeatureClassDescriptor feaDes;
        private IGeoDataset inGeodataset;//输入栅格
        private IGeoDataset outGeodataset;//输出栅格


        #region 环境变量


        private double cellSize = 500;//输出像元大小
        private object Missing = Type.Missing;
        private object cellSizeObj;
        private object extentProObj;//处理范围
        private IRasterRadius radius;
        #endregion
        #region Spline变量


        private esriGeoAnalysisSplineEnum splineEnum;
        #endregion


        private void btnSpline_Click(object sender, EventArgs e)
        {
            try
            {
                interOp = rasterEnv as IInterpolationOp2;
                outGeodataset = interOp.Spline(inGeodataset, splineEnum, ref Missing, ref Missing);
                ShowResult(outGeodataset, "Spline");
                MessageBox.Show("完成");
            }
            catch { }
        }
        private void Form1_Load(object sender, EventArgs e)
        {
            rasterEnv = new RasterInterpolationOpClass();
        }




        private void cmbFields_SelectedIndexChanged(object sender, EventArgs e)
        {
            try
            {
                feaDes = new FeatureClassDescriptorClass();
                feaDes.Create(feaClass, null, cmbFields.SelectedItem.ToString());
                inGeodataset = feaDes as IGeoDataset;
            }
            catch
            {
                MessageBox.Show("请输入栅格图层！");
            };
        }


        private void txtCellSize_TextChanged(object sender, EventArgs e)
        {
            try
            {
                cellSize = Convert.ToDouble(txtCellSize.Text);
                cellSizeObj = cellSize;
                rasterEnv.SetCellSize(esriRasterEnvSettingEnum.esriRasterEnvValue, ref cellSizeObj);
            }
            catch { }
        }


     


    


        private void cmbSplineType_SelectedIndexChanged(object sender, EventArgs e)
        {


        
            int index = cmbSplineType.SelectedIndex;
            switch (index)
            {
                case 0:// 规则样条函数方法 
                    splineEnum = esriGeoAnalysisSplineEnum.esriGeoAnalysisRegularizedSpline;
                    break;
                case 1://TENSION 张力样条函数方法
                    splineEnum = esriGeoAnalysisSplineEnum.esriGeoAnalysisTensionSpline;
                    break;
            }
        }




        private void cmbLayers_MouseClick(object sender, MouseEventArgs e)
        {
            ComboBox c = sender as ComboBox;
            c.Items.Clear();


            IMap map = axMapControl1.Map;
            if (map != null)
            {
                for (int i = 0; i < map.LayerCount; i++)
                {
                    c.Items.Add(map.get_Layer(i).Name);
                }
            }
        }


        //通过图层名得到图层
        private ILayer getLayerFromName(string layerName)
        {
            ILayer layer;
            IMap map = axMapControl1.Map;
            for (int i = 0; i < map.LayerCount; i++)
            {
                layer = map.get_Layer(i);
                if (layerName == layer.Name)
                    return layer;
            }
            return null;
        }




        private void cmbLayers_SelectedIndexChanged(object sender, EventArgs e)
        {
            try
            {
                cmbFields.Items.Clear();


                ILayer layer = getLayerFromName(cmbLayers.SelectedItem.ToString());
                IFeatureLayer fl = layer as IFeatureLayer;
                feaClass = fl.FeatureClass;
                for (int i = 0; i < feaClass.Fields.FieldCount; i++)
                {
                    cmbFields.Items.Add(feaClass.Fields.get_Field(i).Name);
                }
                extentProObj = layer;
                //设置空间处理范围
                rasterEnv.SetExtent(esriRasterEnvSettingEnum.esriRasterEnvValue, ref extentProObj, Missing);//设置空间处理范围
            }
            catch
            {
                MessageBox.Show("请输入栅格图层！");
            }
        }
        //显示分析结果
        private void ShowResult(IGeoDataset geoDataset, string interType)
        {
            IRasterLayer rasterLayer = new RasterLayerClass();
            IRaster raster = new Raster();
            raster = (IRaster)geoDataset;
            rasterLayer.CreateFromRaster(raster);
            rasterLayer.Name = interType;


            axMapControl1.AddLayer((ILayer)rasterLayer, 1);
            axMapControl1.ActiveView.Refresh();
        }


        private void axMapControl1_OnMouseDown(object sender, ESRI.ArcGIS.Controls.IMapControlEvents2_OnMouseDownEvent e)
        {


        }
    }
}