大家好,又见面了,我是你们的朋友全栈君。
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();
发布者:全栈程序员-用户IM,转载请注明出处:https://javaforall.cn/163455.html原文链接:https://javaforall.cn
【正版授权,激活自己账号】: Jetbrains全家桶Ide使用,1年售后保障,每天仅需1毛
【官方授权 正版激活】: 官方授权 正版激活 支持Jetbrains家族下所有IDE 使用个人JB账号...