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
{}