ArcGIS二次开发基础教程(10)：三维分析

//DEM数据的坡度分析 将分析结果添加到地图上
//首先获取DEM数据，方法有很多例如从个人地理数据库获取，也可直接获取文件数据，此处采用第二种方法
IWorkspaceFactory workspaceFactory = new ShapefileWorkspaceFactoryClass();
//从文件处打开工作空间
IRasterWorkspace rasterWorkspace = workspaceFactory.OpenFromFile(System.IO.Path.GetDirectoryName(path)) as IRasterWorkspace;
//打开数据集
IRasterDataset rasterDatset = rasterWorkspace.OpenRasterDataset(System.IO.Path.GetFileName(path));
//转换为坡度分析需要的地理数据集
IGeoDataset rasterGeoDataset = rasterDataset as IGeoDataset;
//表面分析接口 请查阅帮助文档
ISurfaceOp surfaceOp = new RasterSurfaceOpClass();
//坡度分析方法
IGeoDataset resultDataset = surfaceOp.Slope(rasterGeoDataset,esriGeoAnalysisSlopeEnum.esriGeoAnalysisSlopeDegrees,Type.Missing);
//将生成的数据集转换为栅格图层加到地图中
IRasterLayer rasterLayer = new RasterLayerClass();
//当栅格数据为单个数据栅格时
rasterLayer.CreateFromRaster(resultDataset as IRaster);
//当栅格数据为栅格数据集时
//rasterLayer.CreateFromDataset((IRasterDataset)reaultDataset);
rasterLayer.Name = "Slope";
axMapControl1.AddLayer(rasterLayer as ILayer);
axMapControl1.Refresh();
axTOCControl1.Update();

//通视分析多需要一个要素数据集表示观察点
//首先获取DEM数据，方法有很多例如从个人地理数据库获取，也可直接获取文件数据，此处采用第二种方法
IWorkspaceFactory workspaceFactory = new ShapefileWorkspaceFactoryClass();
//从文件处打开工作空间  path1是DEM数据的路径
IRasterWorkspace rasterWorkspace = workspaceFactory.OpenFromFile(System.IO.Path.GetDirectoryName(path1)) as IRasterWorkspace;
//打开数据集
IRasterDataset rasterDatset = rasterWorkspace.OpenRasterDataset(System.IO.Path.GetFileName(path1));
//转换为坡度分析需要的地理数据集
IGeoDataset rasterGeoDataset = rasterDataset as IGeoDataset;
//获取要素数据集  path2是观察点要素数据集
IFeatrueWorkspace featureWorkspace = workspaceFactory.OpenFromFile(System.IO.Path.GetDirectoryName(path2)) as IFeatureWorkspace;
IFeatrueDataset featureDataset = featureWorkspace.OpenFeatureDataset(System.IO.Path.GetFileName(path2));
IGeoDataset featureGeoDataset = featureDataset as IGeoDataset;
//表面分析接口 请查阅帮助文档
ISurfaceOp surfaceOp = new RasterSurfaceOpClass();
//通视分析方法
IGeoDataset resultDataset = surfaceOp.Visibility(rasterGeoDataset,featureGeoDataset,esriGeoAnalysisVisibilityEnum.esriGeoAnalysisVisibilityFrequency,Type.Missing);
//将生成的数据集转换为栅格图层加到地图中
IRasterLayer rasterLayer = new RasterLayerClass();
rasterLayer.CreateFromRaster(resultDataset as IRaster);
rasterLayer.Name = "Visibility";
axMapControl1.AddLayer(rasterLayer as ILayer);
axMapControl1.Refresh();
axTOCControl1.Update();