AE调用GP工具的方法（转）

第一，首先要明确自己需要调用arctoolbox里面的什么工具，实现什么样的功能。

第二，按照需求看看在arctoolbox工具中是怎么实现功能的，然后确定需要的数据源。

第三，编写command或tool工具，编写自己要的功能工具。

第四，以自己编写的一个简单的克里金插值为例来说明一下。（自己曾经没有用调用GP的方式实现这个功能写了很长很长的代码，但是试过了调用GP工具的方式，只用了那么几行简单的代码就OK了，还是很方便的，呵呵）。

1）首先创建一个工程。

2）创建一个Icommand类。

3）定义工具的名称等信息

        public Kriging()
        {
            base.m_category = "GPTools";
            base.m_caption = "克吕金插值"; 
            base.m_message = "克吕金插值";
            base.m_toolTip = "克吕金插值";
            base.m_name = "Kriging";  
            try
            {
                string bitmapResourceName = GetType().Name + ".bmp";
                base.m_bitmap = new Bitmap(GetType(), bitmapResourceName);
            }
            catch (Exception ex)
            {
                System.Diagnostics.Trace.WriteLine(ex.Message, "Invalid Bitmap");
            }
        }

4）在单击的方法中填写代码

      public override void OnClick()
        {
            Geoprocessor GP = new Geoprocessor();
            ESRI.ArcGIS.Analyst3DTools.Kriging pKriging = new Kriging();
            pKriging.in_point_features = @"E:\海平面气压.shp";
            pKriging.z_field = "pslv";
            pKriging.out_surface_raster = @"E:\\Krigingnew";
            pKriging.semiVariogram_props = "Spherical";
            GP.Execute(pKriging, null);
        }

当然，可以编写个界面的形式，可以自动选择in_point_features 与out_surface_raster 的路径。这样更人性化一点，这里只是简单的写写。

5）运行程序即可，在arcmap中调用试了一下，还是挺成功的。

本来克里金插值插出来grid数据是一个规则边界的（根据插值的范围来定的Extent），上面的图是插值后经过Extraction后的效果，看起来好看多了。

Extraction工具也可以调用GP工具，代码如下：

         Geoprocessor GP = new Geoprocessor();
         ESRI.ArcGIS.SpatialAnalystTools.ExtractByMask pExtractByMask = new pExtractByMask();
         pExtractByMask.in_raster = @"C:\map\grid\dqjsl";
         pExtractByMask.in_mask_data = @"C:\map\XZJQ_ShiJi_R";
         pExtractByMask.out_raster = @"C:\map\clipresult";
         CommonClass.RunTool(GP, pExtractByMask, null);

当然还可以设置其他的细节的东西，比如各用到的参数等。