C#+AE 三维分析

红太狼 2023-02-09 14:06 97阅读 0赞

因为自己对三维的模型比较感兴趣，所以就好好的学习了一下AE中的一些三维分析的功能，因此也就做一下笔记记录一下学习的过程。

### 文章目录 ###

*  一、三维数据
 *  二、TIN数据分析
 *   *  2.1属性信息
     *  2.2创建TIN数据
     *  2.3TIN表面分析
     *  2.4坡度坡向图
 *  三、小结

# 一、三维数据 #

> 简单来讲，三维数据就是在传统的二维数据(x,y)中添加了一个z值,也就变为了三维数据(x,y,z)，这个z值并不局限于高程，它也可以指很多其他的东西，如亮度、人口数量等等。所以说虽然几十年前没有我们现在常说的基于高程的三维数据，但是我们不能说之前就没有三维数据，只不过那时的三维数据的z值可能不是高程而已。

除此之外呢，在开始三维分析的学习之前，首先就要了解一下三维数据的模型，因为不同的地理信息软件产品有着不同的数据格式和数据模型，了解该产品的数据模型我们才能更好的学习怎样使用他们。

> AE的三维数据模型有以下几种：  
> 1、3D要素模型：三维点数据、三维线数据和多面体(MultiPatch)数据。而这其中特别要注意的就是多面体数据，因为多面体数据是有多个面组成的，而这些面的组合方式又有以下三种情况：  
> (1)Triangle Strip：条带式三角形，它由一系列相邻三角形连接而成，每两个相邻的三角形有一条公共边和两个公共点。

![在这里插入图片描述][watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2RheXVoYWl0YW5nMQ_size_16_color_FFFFFF_t_70]

> (2)Triangle Fan：直译是扇形三角形，但是有感觉又不是很准确。该组织方式是由一系列点定义的三角形形成的，并且所有的三角形共享一个顶点。

![在这里插入图片描述][watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2RheXVoYWl0YW5nMQ_size_16_color_FFFFFF_t_70 1]

> (3)Ring：环对象，是一种封闭的路径对象，它的起点和终点是同一个点，主要用在多边形对象中。

![在这里插入图片描述][20200519110053419.png]

> 2、表面数据：就是指具有空间连续特征的地理要素的集合，表示地球表面某部分或整体范围内的地理要素或现象。这种数据我们常用的有栅格表面、不规则三角网(TIN)和Terrain数据集等。  
> 因为栅格表面和TIN我们并不陌生，所以就不需要过多介绍了，值得一说的就是Terrain数据集：它是一种基于TIN的数据集，是多分辨率的TIN。Terrain具有一系列TIN，其中每个TIN都在特定的地图比例尺范围内使用。当地图范围较大时，使用粗颗粒TIN；而当放大地图并将实现集中于特定的地图范围是，则使用细颗粒的TIN，这样可以更为详细的描述地物的特征。

# 二、TIN数据分析 #

因为我们最常使用的可能就是TIN数据了，所以我就实现了一下TIN数据的一些相关分析。

## 2.1属性信息 ##

代码如下：

TINAttributeInfo .cs

public partial class TINAttributeInfo : DevExpress.XtraEditors.XtraForm
         { 
            public TINAttributeInfo()
            { 
                InitializeComponent();
            }
    
            public List<string> layerPaths {  get; set; }
    
            private void TINAttributeInfo_Load(object sender, EventArgs e)
            { 
                comboBoxEdit_Layer.Properties.Items.Clear();
                if (layerPaths == null)
                { 
                    return;
                }
    
                for (int i = 0; i < layerPaths.Count; i++)
                { 
                    comboBoxEdit_Layer.Properties.Items.Add(layerPaths[i]);
                }
            }
    
            public string X {  get; set; }
            public string Y {  get; set; }
            public string Elevation {  get; set; }
            public string Slope {  get; set; }
            public string Aspect {  get; set; }
            public int layerIndex {  get; set; }
    
            public void ShowInfo()
            { 
                textBox_elevation.Text = todouble(Elevation);
                textBox_slope.Text = todouble(Slope);
                textBox_aspect.Text = todouble(Aspect);
                textBox_y.Text = todouble(Y);
                textBox_x.Text = todouble(X);
                
            }
    
            public string todouble(string s)
            { 
                double f = Convert.ToDouble(s);
                return f.ToString("0.01");
            }
    
            private void comboBoxEdit_Layer_SelectedIndexChanged(object sender, EventArgs e)
            { 
                layerIndex = comboBoxEdit_Layer.SelectedIndex;
            }
        }

mianForm.cs

//三维属性信息
            private void barButtonItem_SurfaceAnalysis_ItemClick(object sender, DevExpress.XtraBars.ItemClickEventArgs e)
            { 
                TINAttr = new TINAttributeInfo();
                List<string> layerPaths = new List<string>();  //建立图层路径集合
                for (int i = 0; i < axMapControl_MapView.Map.LayerCount; i++)
                { 
                    layerPaths.Add(axMapControl_MapView.Map.Layer[i].Name);
                }
    
                TINAttr.layerPaths = layerPaths;
                TINAttr.Show();     //可以支持多线程操作
                axMapControl_MapView.MousePointer = ESRI.ArcGIS.Controls.esriControlsMousePointer.esriPointerCrosshair;
            }

实现效果：  
![在这里插入图片描述][watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2RheXVoYWl0YW5nMQ_size_16_color_FFFFFF_t_70 2]  
![在这里插入图片描述][watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2RheXVoYWl0YW5nMQ_size_16_color_FFFFFF_t_70 3]

## 2.2创建TIN数据 ##

CreateTIN.cs

public partial class CreateTIN : DevExpress.XtraEditors.XtraForm
     { 
         public CreateTIN()
         { 
             InitializeComponent();
         }
    
         public AxSceneControl axSceneControl {  get; set; }
    
         private IFeatureLayer pFeatureLayer;
         private IFeatureClass pFeatureClass;
    
         //获取并加载矢量数据集
         private void button_InputPath_Click(object sender, EventArgs e)
         { 
             OpenFileDialog openFileDialog = new OpenFileDialog();
             openFileDialog.Title = "打开Shp点图层文件";
             openFileDialog.Filter = "Shape文件(*.shp)|*.shp";
             openFileDialog.ShowDialog();
             textBox_InputPath.Text = openFileDialog.FileName;
    
             string path;
             string filname;
    
             FileInfo fileinfo = new FileInfo(textBox_InputPath.Text);
             path = fileinfo.DirectoryName;
             filname = System.IO.Path.GetFileName(textBox_InputPath.Text);
             try
             { 
                 IWorkspaceFactory pWSF = new ShapefileWorkspaceFactory();
                 IWorkspace pWS = pWSF.OpenFromFile(path, 0);
                 IFeatureWorkspace pFWS = pWS as IFeatureWorkspace;
                 pFeatureClass = pFWS.OpenFeatureClass(filname);
                 pFeatureLayer = new FeatureLayerClass();
                 pFeatureLayer.FeatureClass = pFeatureClass;
                 pFeatureLayer.Name = pFeatureClass.AliasName;
                 ILayer layer = pFeatureLayer as ILayer;
                 axSceneControl.Scene.AddLayer(layer, false);
             }
             catch
             { 
                 MessageBox.Show("Error");
             }
    
         }
    
         //保存文件
         private void simpleButton_OutputPath_Click(object sender, EventArgs e)
         { 
             SaveFileDialog saveTIN = new SaveFileDialog();
             saveTIN.Title = "选择TIN文件保存路径";
             saveTIN.Filter = "TIN(*.tin)|*.tin|All files (*.*)|*.*";
             saveTIN.RestoreDirectory = true;
    
             if (saveTIN.ShowDialog() == DialogResult.OK)
             { 
                 textBox_OutputPath.Text= saveTIN.FileName;
             }
         }
    
         //矢量数据转换并显示
         private void simpleButton_FeatureTransferAndShow_Click(object sender, EventArgs e)
         { 
             try
             { 
                 string path = textBox_OutputPath.Text.ToString();
                 if (path == "" || path == null)
                 { 
                     MessageBox.Show("结果保存路径不能为空！！");
                     return;
                 }
                 CreateTinFromFeature(pFeatureClass, path);
                 this.Close();
             }
             catch {  }
         }
    
         //矢量数据创建TIN
         public void CreateTinFromFeature(IFeatureClass featureclass, string savepath)
         { 
             IGeoDataset pGeodataset = featureclass as IGeoDataset;
             IEnvelope extent = pGeodataset.Extent;
             extent.SpatialReference = pGeodataset.SpatialReference;
    
             IFields fields = featureclass.Fields;
             IField heightfield = fields.get_Field(4);
    
             ITinEdit tinedit = new TinClass();
             tinedit.InitNew(extent);
    
             try
             { 
                 tinedit.AddFromFeatureClass(featureclass, null, heightfield, null, esriTinSurfaceType.esriTinMassPoint);
             }
             catch
             { 
                 MessageBox.Show("创建TIN失败");
             }
    
             tinedit.SaveAs(savepath);
             tinedit.StopEditing(false);
    
             AddTinData(savepath);
         }
    
         //添加并显示
         public void AddTinData(string path)
         { 
             IWorkspaceFactory pWSF = new TinWorkspaceFactory();
             ITinWorkspace TinWS;
             ITin tin;
             FileInfo fileinfo = new FileInfo(path);
             if (pWSF.IsWorkspace(fileinfo.DirectoryName))
             { 
                 TinWS = pWSF.OpenFromFile(fileinfo.DirectoryName, 0) as ITinWorkspace;
                 tin = TinWS.OpenTin(fileinfo.Name);
                 ITinLayer tinlayer = new TinLayerClass();
                 tinlayer.Dataset = tin;
                 tinlayer.Name = fileinfo.Name;      //给图层添上名字
                 ILayer layer = tinlayer as ILayer;
                 axSceneControl.Scene.AddLayer(layer, true);
                 axSceneControl.SceneGraph.RefreshViewers();
             }
         }
    
         /// <summary>
         /// 栅格数据转TIN
         /// </summary>
         /// <param name="sender"></param>
         /// <param name="e"></param>
    
         //栅格输出路径
         private void simpleButton_RasterOutputPath_Click(object sender, EventArgs e)
         { 
             SaveFileDialog saveTIN = new SaveFileDialog();
             saveTIN.Title = "选择TIN文件保存路径";
             saveTIN.Filter = "TIN(*.tin)|*.tin|All files (*.*)|*.*";
             saveTIN.RestoreDirectory = true;
    
             if (saveTIN.ShowDialog() == DialogResult.OK)
             { 
                 textBox_RasterOutputPath.Text = saveTIN.FileName;
             }
         }
    
         //获取栅格输入数据集
         private void simpleButton_RasterInputPath_Click(object sender, EventArgs e)
         { 
             OpenFileDialog openFileDialog = new OpenFileDialog();
             openFileDialog.Title = "打开.tif文件";
             openFileDialog.Filter = ".tif文件(*.tif)|*.tif|PNG(*.png)|*.png";
             openFileDialog.ShowDialog();
             textBox_RasterDataInput.Text = openFileDialog.FileName;
    
             ILayer layer = GetDEMLayer(textBox_RasterDataInput.Text);
             ISceneGraph pSceneGraph = this.axSceneControl.SceneGraph;
             pSceneGraph.Scene.AddLayer(layer);
             pSceneGraph.RefreshViewers();
         }
    
         //获取DEM图层
         private ILayer GetDEMLayer(string path)
         { 
             ILayer layer;
             IRasterLayer rasterLayer = null;
             IWorkspaceFactory rasterWorkspaceFactory = new RasterWorkspaceFactoryClass();
             IRasterWorkspace rasterWorkspace;
             FileInfo file = new FileInfo(path);
             if (rasterWorkspaceFactory.IsWorkspace(file.DirectoryName))
             { 
                 rasterWorkspace = rasterWorkspaceFactory.OpenFromFile(file.DirectoryName, 0) as IRasterWorkspace;
                 IRasterDataset rasterDataset = rasterWorkspace.OpenRasterDataset(file.Name);
                 rasterLayer = new RasterLayerClass();
                 rasterLayer.CreateFromDataset(rasterDataset);
    
             }
             layer = rasterLayer as ILayer;
             return layer;
         }
    
         //由栅格数据创建TIN
         public void CreateTinFromRaster(IRaster raster, string savepath)
         { 
             double m_zTolerance = 0.001;
             IGeoDataset pGeoData = raster as IGeoDataset;
             IEnvelope pExtent = pGeoData.Extent;
             IRasterBandCollection pRasBC = raster as IRasterBandCollection;
             IRasterBand pRasBand = pRasBC.Item(0);
             IRawPixels pRawPixels = pRasBand as IRawPixels;
             IRasterProps pProps = pRawPixels as IRasterProps;
    
             int iWid = pProps.Width;
             int iHei = pProps.Height;
    
             double w = iWid / 1000.0f;
             double h = iHei / 1000.0f;
    
             IPnt pBlockSize = new DblPntClass();
             bool IterationFlag;
    
             if (w < 1 && h < 1) //横纵都小于1000个像素
             { 
                 pBlockSize.X = iWid;
                 pBlockSize.Y = iHei;
                 IterationFlag = false;
             }
             else
             { 
                 pBlockSize.X = 1001.0f;
                 pBlockSize.Y = 1001.0f;
                 IterationFlag = true;
             }
    
             double cellsize = 0.0f;     //栅格大小
             IPnt pPnt1 = pProps.MeanCellSize(); //栅格平均大小
             cellsize = pPnt1.X;
    
             ITinEdit pTinEdit = new TinClass() as ITinEdit;
             pTinEdit.InitNew(pExtent);
    
             ISpatialReference pSpatial = pGeoData.SpatialReference;
             pExtent.SpatialReference = pSpatial;
    
             IPnt pOrigin = new DblPntClass();
             IPnt pPixelBlockOrigin = new DblPntClass();
    
             //栅格左上角像素中心坐标
             double bX = pBlockSize.X;
             double bY = pBlockSize.Y;
    
             pBlockSize.SetCoords(bX, bY);
             IPixelBlock pPixelBlock = pRawPixels.CreatePixelBlock(pBlockSize);
    
             object nodata = pProps.NoDataValue;     //无值标记
             ITinAdvanced2 pTinNodeCount = pTinEdit as ITinAdvanced2;
             int nodeCount = pTinNodeCount.NodeCount;
    
             object vtMissing = Type.Missing;
    
             object vPixels = null;      //格子
             if (IterationFlag)   //当为一个处理单元格子时
             { 
                 pPixelBlockOrigin.SetCoords(0.0f, 0.0f);
                 pRawPixels.Read(pPixelBlockOrigin, pPixelBlock);
    
                 vPixels = pPixelBlock.get_SafeArray(0);
                 double xMin = pExtent.XMin;
                 double yMax = pExtent.YMax;
                 pOrigin.X = xMin + cellsize / 2;
                 pOrigin.Y = yMax - cellsize / 2;
                 bX = pOrigin.X;
                 bY = pOrigin.Y;
    
                 pTinEdit.AddFromPixelBlock(bX, bY, cellsize, cellsize, nodata, vPixels, m_zTolerance, ref vtMissing, out vtMissing);
             }
             else  //当有多个处理单元格时，依次循环处理每个单元格
             { 
                 int i = 0, j = 0, count = 0;
                 int FirstGoNodeCount = 0;
                 while (nodeCount != FirstGoNodeCount)
                 { 
                     count++;
                     nodeCount = pTinNodeCount.NodeCount;
                     //依次循环处理
                     for (i = 0; i < h + 1; i++)
                     { 
                         for (j = 0; j < w + 1; j++)
                         { 
                             double bX1, bY1, xMin1, yMax1;
                             bX1 = pBlockSize.X;
                             bY1 = pBlockSize.Y;
    
                             pPixelBlockOrigin.SetCoords(j * bX1, i * bY1);
                             pRawPixels.Read(pPixelBlockOrigin, pPixelBlock);
                             vPixels = pPixelBlock.get_SafeArray(0);
    
                             xMin1 = pExtent.XMin;
                             yMax1 = pExtent.YMax;
    
                             bX1 = pBlockSize.X;
                             bY1 = pBlockSize.Y;
    
                             pOrigin.X = xMin1 + j * bX1 * cellsize + cellsize / 2.0f;
                             pOrigin.Y = yMax1 + i * bY1 * cellsize - cellsize / 2.0f;
    
                             bX1 = pOrigin.X;
                             bY1 = pOrigin.Y;
    
                             pTinEdit.AddFromPixelBlock(bX1, bY1, cellsize, cellsize, nodata, vPixels, m_zTolerance, ref vtMissing, out vtMissing);
    
                             FirstGoNodeCount = pTinNodeCount.NodeCount;
                         }
                     }
                 }
             }
    
             //保存TIN文件
             pTinEdit.SaveAs(savepath);
             pTinEdit.StopEditing(true);
             MessageBox.Show("转换完成", "消息");
             AddTinData(savepath);
         }
    
         //栅格转换为TIN
         private void simpleButton_RasterTransAndShow_Click(object sender, EventArgs e)
         { 
             try
             { 
                 string path = textBox_RasterOutputPath.Text.ToString();
                 if (path == "" || path == null)
                 { 
                     MessageBox.Show("结果保存路径不能为空！！");
                     return;
                 }
                 IScene scene = axSceneControl.Scene;
                 IRaster raster = null;
                 for (int i = 0; i < scene.LayerCount; i++)
                 { 
                     ILayer lyr = scene.get_Layer(i);
                     if (lyr is IRasterLayer)
                     { 
                         IRasterLayer rLayer = lyr as IRasterLayer;
                         raster = rLayer.Raster;
                         break;
                     }
                 }
    
                 CreateTinFromRaster(raster, path);
                 this.Close();
             }
             catch {  }
         }
     }

mainForm.cs

//创建TIN
    	private void barButtonItem_CreateTIN_ItemClick(object sender, ItemClickEventArgs e)
    	{ 
    	    CreateTIN tin = new CreateTIN();
    	    tin.axSceneControl = axSceneControl;
    	    tin.Show();
    	}

实现效果：  
![在这里插入图片描述][watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2RheXVoYWl0YW5nMQ_size_16_color_FFFFFF_t_70 4]  
![在这里插入图片描述][watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2RheXVoYWl0YW5nMQ_size_16_color_FFFFFF_t_70 5]  
![在这里插入图片描述][watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2RheXVoYWl0YW5nMQ_size_16_color_FFFFFF_t_70 6]

## 2.3TIN表面分析 ##

代码分析：  
TINSurfaceAnalysis.cs:

public partial class TINSurfaceAnalysis : DevExpress.XtraEditors.XtraForm
     { 
         public TINSurfaceAnalysis()
         { 
             InitializeComponent();
         }
    
         public AxSceneControl axSceneControl {  get; set; }
         private IScene scene;
    
         private void TINSurfaceAnalysis_Load(object sender, EventArgs e)
         { 
             comboBoxEdit_Layer.Properties.Items.Clear();
             scene= axSceneControl.Scene;
             if (scene==null)
             { 
                 return;
             }
    
             for (int i = 0; i < scene.LayerCount; i++)
             { 
                 comboBoxEdit_Layer.Properties.Items.Add(scene.Layer[i].Name);
             }
    
         }
         //全局变量
         private ILayer layer;
         private ITinLayer tinlayer;
         private ITin tin;
         private ISurface surface;
         private IPoint point;
         private double volume;
         private double zfactor;
         private bool reftype = true;
    
         private void simpleButton_Volumn_Click(object sender, EventArgs e)
         { 
             zfactor = Convert.ToDouble(textEdit_Height.Text);
             if (reftype)
             { 
                 volume = surface.GetVolume(zfactor, esriPlaneReferenceType.esriPlaneReferenceAbove);
                 memoEdit_ResultShow.Text += "在面高度" + zfactor.ToString() + "之上的体积为:\r\n" + volume.ToString() + "\r\n";
             }
             else
             { 
                 volume = surface.GetVolume(zfactor, esriPlaneReferenceType.esriPlaneReferenceBelow);
                 memoEdit_ResultShow.Text += "在面高度" + zfactor.ToString() + "之下的体积为:\r\n" + volume.ToString() + "\r\n";
             }
         }
    
         private void comboBoxEdit_Layer_SelectedIndexChanged(object sender, EventArgs e)
         { 
             //获取TIN图层，并得到TIN表面
             layer = axSceneControl.Scene.get_Layer(comboBoxEdit_Layer.SelectedIndex);       //根据选择的索引来确定图层
             tinlayer = layer as ITinLayer;
             tin = tinlayer.Dataset;
             surface = tin as ISurface;
         }
    
         //选择基础面上下
         private void radioGroup_SelectBasicface_SelectedIndexChanged(object sender, EventArgs e)
         { 
             if (radioGroup_SelectBasicface.SelectedIndex==0)
             { 
                 reftype = true;
             }
             else
             { 
                 reftype = false;
             }
         }
    
         //投影面积
         private void simpleButton_ProjectArea_Click(object sender, EventArgs e)
         { 
             zfactor = Convert.ToDouble(textEdit_Height.Text);
             if (reftype)
             { 
                 volume = surface.GetProjectedArea(zfactor, esriPlaneReferenceType.esriPlaneReferenceAbove);
                 memoEdit_ResultShow.Text += "在面高度" + zfactor.ToString() + "之上的投影面积积为:\r\n" + volume.ToString() + "\r\n";
             }
             else
             { 
                 volume = surface.GetProjectedArea(zfactor, esriPlaneReferenceType.esriPlaneReferenceBelow);
                 memoEdit_ResultShow.Text += "在面高度" + zfactor.ToString() + "之下的投影面积为:\r\n" + volume.ToString() + "\r\n";
             }
         }
    
         //表面积
         private void simpleButton_SurfaceArea_Click(object sender, EventArgs e)
         { 
             if (reftype)
             { 
                 volume = surface.GetSurfaceArea(zfactor, esriPlaneReferenceType.esriPlaneReferenceAbove);
                 memoEdit_ResultShow.Text += "在面高度" + zfactor.ToString() + "之上的表面积积为:\r\n" + volume.ToString() + "\r\n";
             }
             else
             { 
                 volume = surface.GetSurfaceArea(zfactor, esriPlaneReferenceType.esriPlaneReferenceBelow);
                 memoEdit_ResultShow.Text += "在面高度" + zfactor.ToString() + "之下的表面积为:\r\n" + volume.ToString() + "\r\n";
             }
         }
    
         //清楚结果
         private void simpleButton_Clear_Click(object sender, EventArgs e)
         { 
             memoEdit_ResultShow.Text = "";
         }
    
         //关闭
         private void simpleButton_Close_Click(object sender, EventArgs e)
         { 
             this.Close();
         }
     }

mainForm.cs:

//TIN表面分析
    	private void barButtonItem_SurfaceAnalysis_ItemClick(object sender, ItemClickEventArgs e)
    	{ 
    	    TINSurfaceAnalysis TINsurf = new TINSurfaceAnalysis();
    	    TINsurf.axSceneControl = axSceneControl;
    	    TINsurf.Show();
    	}

实现效果：  
![在这里插入图片描述][watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2RheXVoYWl0YW5nMQ_size_16_color_FFFFFF_t_70 7]  
![在这里插入图片描述][watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2RheXVoYWl0YW5nMQ_size_16_color_FFFFFF_t_70 8]

## 2.4坡度坡向图 ##

代码如下：

public partial class SlopeAndAspect : DevExpress.XtraEditors.XtraForm
    	{ 
    	    public SlopeAndAspect()
    	    { 
    	        InitializeComponent();
    	    }
    	
    	    public AxSceneControl axSceneControl {  get; set; }
    	    private IScene scene;
    	
    	    //窗体生成事件
    	    private void SlopeAndAspect_Load(object sender, EventArgs e)
    	    { 
    	        comboBoxEdit_Layer.Properties.Items.Clear();
    	        scene = axSceneControl.Scene;
    	
    	        if (scene==null)
    	        { 
    	            return;
    	        }
    	
    	        for (int i = 0; i < scene.LayerCount; i++)
    	        { 
    	            comboBoxEdit_Layer.Properties.Items.Add(scene.Layer[i].Name);
    	        }
    	
    	    }
    	
    	    //坡度分析
    	    private void simpleButton_SlopeAnalysis_Click(object sender, EventArgs e)
    	    { 
    	        try
    	        { 
    	            ILayer layer = axSceneControl.Scene.Layer[comboBoxEdit_Layer.SelectedIndex];
    	            string filepath = getLayerPath(layer);
    	            if (filepath == null || textBox_Output.Text==null)
    	            { 
    	                return;
    	            }
    	            string outputPath = textBox_Output.Text;
    	            ExecuteSlope(filepath, outputPath);
    	            ShowResult(outputPath);
    	            MessageBox.Show("生成成功！");
    	
    	        }
    	        catch (Exception ex)
    	        { 
    	            MessageBox.Show(ex.Message);
    	        }
    	    }
    	
    	    //坡度分析
    	    private void ExecuteSlope(string input, string output)
    	    { 
    	        Geoprocessor gp = new Geoprocessor();
    	        gp.OverwriteOutput = true;
    	        ESRI.ArcGIS.Analyst3DTools.SurfaceSlope Slope = new ESRI.ArcGIS.Analyst3DTools.SurfaceSlope(input, output);
    	        IGeoProcessorResult result = gp.Execute(Slope, null) as IGeoProcessorResult;
    	        if (result.Status != esriJobStatus.esriJobSucceeded)
    	        { 
    	            MessageBox.Show("Failed to Aspect layer");
    	        }
    	    }
    	
    	    //保存路径
    	    private void simpleButton_SaveResult_Click(object sender, EventArgs e)
    	    { 
    	        SaveFileDialog saveDlg = new SaveFileDialog();
    	        saveDlg.CheckPathExists = true;
    	        saveDlg.Filter = "Shapefile (*.shp)|*.shp";
    	        saveDlg.OverwritePrompt = true;
    	        saveDlg.Title = "Output Layer";
    	        saveDlg.RestoreDirectory = true;
    	
    	        DialogResult dr = saveDlg.ShowDialog();
    	        if (dr == DialogResult.OK)
    	        { 
    	            textBox_Output.Text = saveDlg.FileName;
    	        }
    	    }
    	    //坡向分析
    	    private void ExecuteAspect(string input, string output)
    	    { 
    	        Geoprocessor gp = new Geoprocessor();
    	        gp.OverwriteOutput = true;
    	        ESRI.ArcGIS.Analyst3DTools.SurfaceAspect Aspect = new ESRI.ArcGIS.Analyst3DTools.SurfaceAspect(input, output);
    	        IGeoProcessorResult result = gp.Execute(Aspect, null) as IGeoProcessorResult;
    	        if (result.Status != esriJobStatus.esriJobSucceeded)
    	        { 
    	            MessageBox.Show( "Failed to Aspect layer");
    	        }
    	    }
    	
    	    private void simpleButton_AspectAnalysis_Click(object sender, EventArgs e)
    	    { 
    	        try
    	        { 
    	            ILayer layer = axSceneControl.Scene.Layer[comboBoxEdit_Layer.SelectedIndex];
    	            string filepath = getLayerPath(layer);
    	            if (filepath == null || textBox_Output.Text == null)
    	            { 
    	                return;
    	            }
    	            string outputPath = textBox_Output.Text;
    	            ExecuteAspect(filepath, outputPath);
    	            ShowResult(outputPath);
    	            MessageBox.Show("生成成功！");
    	
    	        }
    	        catch (Exception ex)
    	        { 
    	            MessageBox.Show(ex.Message);
    	        }
    	    }
    	
    	    //显示结果
    	    private void ShowResult(string filepath)
    	    { 
    	        FileInfo fileinfo = new FileInfo(filepath);
    	        string path = fileinfo.DirectoryName;
    	        string filename = System.IO.Path.GetFileName(filepath);
    	        IWorkspaceFactory pWSF = new ESRI.ArcGIS.DataSourcesFile.ShapefileWorkspaceFactoryClass();
    	        IWorkspace pWS = pWSF.OpenFromFile(path, 0);
    	        IFeatureWorkspace pFWS = pWS as IFeatureWorkspace;
    	        IFeatureClass pFC = pFWS.OpenFeatureClass(filename);
    	        IFeatureLayer pFLayer = new FeatureLayerClass();
    	        pFLayer.FeatureClass = pFC;
    	        pFLayer.Name = pFC.AliasName;
    	        ILayer pLayer = pFLayer as ILayer;
    	        axSceneControl.Scene.AddLayer(pLayer, true);
    	        axSceneControl.SceneGraph.RefreshViewers();
    	    }
    	
    	    private static string getLayerPath(ILayer pLayer)
    	    { 
    	        IDatasetName pDatasetName = (pLayer as IDataLayer2).DataSourceName as IDatasetName;
    	        IWorkspaceName pWorkspaceName = pDatasetName.WorkspaceName;
    	        return pWorkspaceName.PathName +"\\"+ pLayer.Name;
    	    }
    	
    	    private void simpleButton_Close_Click(object sender, EventArgs e)
    	    { 
    	        this.Close();
    	    }
    	}

mainForm.cs:

//TIN表面坡度坡向图
    	private void barButtonItem_SlopeAndAspect_ItemClick(object sender, ItemClickEventArgs e)
    	{ 
    	    SlopeAndAspect sa = new SlopeAndAspect();
    	    sa.axSceneControl = axSceneControl;
    	    sa.Show();
    	}

实现效果：  
![在这里插入图片描述][watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2RheXVoYWl0YW5nMQ_size_16_color_FFFFFF_t_70 9]  
![在这里插入图片描述][20200524092039184.png]  
![在这里插入图片描述][watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2RheXVoYWl0YW5nMQ_size_16_color_FFFFFF_t_70 10]

# 三、小结 #

整个过程走下来，自我感觉AE的封装程度其实已经很高了，需要真正学习的部分其实主要就是要了解三维数据在AE中的组织方式和熟练使用更为基础的接口对象\*~\*。

> 参考书籍：《ArcgisEngine地信信息系统开发教程》

[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2RheXVoYWl0YW5nMQ_size_16_color_FFFFFF_t_70]: https://img-blog.csdnimg.cn/20200519104655341.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2RheXVoYWl0YW5nMQ==,size_16,color_FFFFFF,t_70
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2RheXVoYWl0YW5nMQ_size_16_color_FFFFFF_t_70 1]: https://img-blog.csdnimg.cn/20200519105818843.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2RheXVoYWl0YW5nMQ==,size_16,color_FFFFFF,t_70
[20200519110053419.png]: https://img-blog.csdnimg.cn/20200519110053419.png
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2RheXVoYWl0YW5nMQ_size_16_color_FFFFFF_t_70 2]: https://img-blog.csdnimg.cn/20200524082341969.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2RheXVoYWl0YW5nMQ==,size_16,color_FFFFFF,t_70
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2RheXVoYWl0YW5nMQ_size_16_color_FFFFFF_t_70 3]: https://img-blog.csdnimg.cn/20200524082331424.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2RheXVoYWl0YW5nMQ==,size_16,color_FFFFFF,t_70
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2RheXVoYWl0YW5nMQ_size_16_color_FFFFFF_t_70 4]: https://img-blog.csdnimg.cn/20200524090009489.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2RheXVoYWl0YW5nMQ==,size_16,color_FFFFFF,t_70
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2RheXVoYWl0YW5nMQ_size_16_color_FFFFFF_t_70 5]: https://img-blog.csdnimg.cn/20200524090026715.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2RheXVoYWl0YW5nMQ==,size_16,color_FFFFFF,t_70
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2RheXVoYWl0YW5nMQ_size_16_color_FFFFFF_t_70 6]: https://img-blog.csdnimg.cn/20200524090054821.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2RheXVoYWl0YW5nMQ==,size_16,color_FFFFFF,t_70
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2RheXVoYWl0YW5nMQ_size_16_color_FFFFFF_t_70 7]: https://img-blog.csdnimg.cn/20200524090355349.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2RheXVoYWl0YW5nMQ==,size_16,color_FFFFFF,t_70
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2RheXVoYWl0YW5nMQ_size_16_color_FFFFFF_t_70 8]: https://img-blog.csdnimg.cn/20200524090309717.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2RheXVoYWl0YW5nMQ==,size_16,color_FFFFFF,t_70
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2RheXVoYWl0YW5nMQ_size_16_color_FFFFFF_t_70 9]: https://img-blog.csdnimg.cn/20200524092006975.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2RheXVoYWl0YW5nMQ==,size_16,color_FFFFFF,t_70
[20200524092039184.png]: https://img-blog.csdnimg.cn/20200524092039184.png
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2RheXVoYWl0YW5nMQ_size_16_color_FFFFFF_t_70 10]: https://img-blog.csdnimg.cn/20200524092130169.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2RheXVoYWl0YW5nMQ==,size_16,color_FFFFFF,t_70