GIS C#+AE9.1疊加求交代碼

 

try
   {
    //分析層
    ILayer pLayer=this.axMapControl2.get_Layer(0);
    IFeatureLayer pInputFeatLayer=pLayer as IFeatureLayer;  
    ITable pInputTable=pLayer as ITable;
    IFeatureClass pInputFeatClass=pInputFeatLayer.FeatureClass;

    //疊加表
    pLayer=this.axMapControl2.get_Layer(1);
    ITable pOverlayTable=pLayer as ITable;
    //疊加分析表
    IFeatureClassName pFeatClassName=new FeatureClassNameClass();
    pFeatClassName.FeatureType=esriFeatureType.esriFTSimple;
    pFeatClassName.ShapeFieldName="shape";
    pFeatClassName.ShapeType=pInputFeatClass.ShapeType;

    //工作空間名稱
    IWorkspaceName pNewWSName=new WorkspaceNameClass();
    pNewWSName.WorkspaceFactoryProgID ="esriDataSourcesFile.ShapefileWorkspaceFactory";
    pNewWSName.PathName = @"E:/temp";

    //資料集名稱
               IDatasetName pDatasetName;
               pDatasetName = pFeatClassName.FeatureDatasetName;
               pDatasetName.Name = "Intersect_result";
               pDatasetName.WorkspaceName = pNewWSName ;
               double  _tol=0;
                IBasicGeoprocessor pBGP;
                pBGP =new BasicGeoprocessorClass();
                IFeatureClass  pOutputFeatClass ;
             pOutputFeatClass = pBGP.Intersect(pInputTable, false, pOverlayTable, false, _tol, pFeatClassName);
                IFeatureLayer      pOutputFeatLayer;
                pOutputFeatLayer = new  FeatureLayerClass();
    pOutputFeatLayer.FeatureClass = pOutputFeatClass;
                pOutputFeatLayer.Name = pOutputFeatClass.AliasName;
                axMapControl2.AddLayer(pOutputFeatLayer,1);
   }
   catch

{}