Arcengine+C#实现矢量数据的裁剪

   
最近工作需要完成图像的裁剪，因为有很多需要裁剪的图像，用arcgis做的话工作量太大，而且后期工作涉及到自动化处理，所以写了段代码用engine实现图像的裁剪功能。首先要获得你所需要被裁剪的图层的IFeatureLayer和裁剪图层的IFeatureLayer，然后调用IBasicGeoprocessor.Clip方法进行裁剪，我把该方法单独封装了一下，代码如下：

 public IFeatureClass Clip(IFeatureClass pInputFeatureClass, IFeatureClass pClipFeatureClass, string pfilePath, string pfileName)

        {

            IFeatureClassName pOutput = new FeatureClassNameClass();

            pOutput.ShapeFieldName = pInputFeatureClass.ShapeFieldName;

            pOutput.ShapeType = pInputFeatureClass.ShapeType;

            pOutput.FeatureType = esriFeatureType.esriFTSimple;

            //获取shapefile数据工作空间

            IWorkspaceName pWsN = new WorkspaceNameClass();

            pWsN.WorkspaceFactoryProgID = "esriDataSourcesFile.ShapefileWorkspaceFactory";

            pWsN.PathName = pfilePath;

            //通过IDatasetName设置输出结果相关参数

            IDatasetName pDatasetName = pOutput as IDatasetName;

            pDatasetName.Name = pfileName;

            pDatasetName.WorkspaceName = pWsN;

            //初始化IBasicGeoprocessor对象，调用clip方法

            IBasicGeoprocessor pBasicGeo = new BasicGeoprocessorClass();

            pBasicGeo.SpatialReference = axMapControl1.SpatialReference;

            //进行裁剪运算

            IFeatureClass pFeatureClass = pBasicGeo.Clip(pInputFeatureClass as ITable, true , pClipFeatureClass as ITable, true , 0.01, pOutput);

            return pFeatureClass; 

        }

        然后我添加了一个button和三个mapcontrol控件，在axmapcontrol1和axmapcontrol2的属性里直接添加了图层，代码如下：

 private void button1_Click(object sender, EventArgs e)

        {

            ILayer pInputLayer = axMapControl1.get_Layer(0);//获取在axmapcontrol1中加载的图层，是需要被裁剪的图层

            ILayer pClipLayer = axMapControl2.get_Layer(0);//获取在axmapcontrol2中加载的图层，是裁剪图层

            if (pInputLayer != null && pClipLayer != null)

            {

                IFeatureLayer pInputFeatureLayer = pInputLayer as IFeatureLayer;

                IFeatureLayer pClipFeatureLayer = pClipLayer as IFeatureLayer;

                IFeatureClass outFeatureClass = Clip(pInputFeatureLayer.FeatureClass, pClipFeatureLayer.FeatureClass, @"C:\Users\Administrator\Desktop\", "裁剪分析结果");

                if (outFeatureClass != null)

                {

                    IFeatureLayer pFeatureLayer = new FeatureLayerClass();

                    pFeatureLayer.FeatureClass = outFeatureClass;

                    pFeatureLayer.Name = outFeatureClass.AliasName;

                    this.axMapControl3.AddLayer(pFeatureLayer);

                }

            }