基于C#的ArcEngine二次开发26：等高线点线矛盾检查

最新推荐文章于 2023-05-04 09:28:29 发布

小薛引路

最新推荐文章于 2023-05-04 09:28:29 发布

阅读量1.7k

点赞数

分类专栏：基于C#的ArcEngine二次开发 C#编程学习

本文链接：https://blog.csdn.net/m1m2m3mmm/article/details/105134908

版权

C#编程学习同时被 2 个专栏收录

62 篇文章 109 订阅

订阅专栏

基于C#的ArcEngine二次开发

61 篇文章 358 订阅

订阅专栏

1 理论基础与数据说明

1.1 要解决的问题及思路

对于点线矛盾的检查，精确的方法是建立不规则三角网；但是三角网的建立比较复杂，其质量的好坏直接决定检查结果的好坏

1.2 等高线的特性及数据成果介绍

基于上述分析，得出如下结论：

主要等高线的编码是以71开头的，因此在进行等高线点线矛盾检查时，要首先提取出TERL层中以71为开头的要素；
高程点层的GB码统一为720100

我们的数据检查要在这两个图层之间开展

2 由FeatureClass生成TIN

2.1 由点要素生成TIN代码

引文来源：ArcEngine由点生成TIN

根据高程点生成TIN的过程：

（1）根据高程点文件（Excel）生成点shapefile

（2）使用ITinEdit的AddFromFeatureClass方法生成TIN

2.1.1 数据说明

2.1.2 代码实现

   #region 创建TIN
        private void button_ok_Click(object sender , EventArgs e)
        {
            try
            {
                strTinName=textBox_TINName.Text;
                outFolder=textBox_outpath.Text;
                if (strTinName=="")
                {
                    MessageBox.Show("请输入TIN名称！");
                } 
                else if (outFolder=="")
                {
                    MessageBox.Show("请选择TIN保存路径");
                }
                else
                {
                    //生成点shape
                    Excel2Shape excel2shape = new Excel2Shape(textBox_point.Text);
                    excel2shape.CreateShapeFromExcel();
 
                    //获取shapefile
                    string path = excel2shape.path;
                    string name = excel2shape.name;
 
                    IWorkspaceFactory pWSFac = new ShapefileWorkspaceFactoryClass();
                    IFeatureWorkspace pFeatureWS = pWSFac.OpenFromFile(path , 0) as IFeatureWorkspace;
                    IFeatureClass pFeatureClass = pFeatureWS.OpenFeatureClass(name);
                    IField pField = pFeatureClass.Fields.get_Field(pFeatureClass.FindField("Z"));
                    //创建TIN
                    ITin pTin = Create_TIN(pFeatureClass , pField , outFolder);
                    ITinLayer pTinlayer = new TinLayerClass();
                    pTinlayer.Dataset = pTin;
                    pTinlayer.Name = strTinName;
 
                    pSceneControl.SceneGraph.Scene.AddLayer(pTinlayer as ILayer , true);
 
                    this.Close();
                }
              
            }
            catch (System.Exception ex)
            {
                MessageBox.Show(ex.Message + "创建TIN失败！");
            }
 
        }
        /// <summary>
        /// 创建TIN
        /// </summary>
        /// <param name="pFeatureClass">点要素类</param>
        /// <param name="pField">Z字段</param>
        /// <param name="pPath">路径</param>
        public ITin Create_TIN(IFeatureClass pFeatureClass , IField pField , string pPath)
        {
            IGeoDataset pGeoDataset = pFeatureClass as IGeoDataset;
            ITinEdit pTinEdit = new TinClass();
            pTinEdit.InitNew(pGeoDataset.Extent);
            object pObj = Type.Missing;
 
            pTinEdit.AddFromFeatureClass(pFeatureClass , null , pField , null ,
            esriTinSurfaceType.esriTinMassPoint , ref pObj);
            if (System.IO.File.Exists(pPath))
            {
                _3DCampus.helper.FolderHelper.DeleteFolder(pPath);
            }
            pTinEdit.SaveAs(pPath , ref pObj);
            pTinEdit.Refresh();
            return pTinEdit as ITin;
 
        }
        #endregion

2.2 TIN转Raster

未必能有，先放这里以供参考

public IRasterDataset tin2raster(string tempBathyTIN, string geoPath, string gridName)
        {
            string tinFolder = System.IO.Path.GetDirectoryName(tempBathyTIN);
            string tinName = System.IO.Path.GetFileName(tempBathyTIN);
            IRasterDataset rasterDataset = new RasterDatasetClass();
            try
            {
                string rasterPath = System.IO.Path.GetDirectoryName(geoPath);
                IWorkspaceFactory TinWF = new TinWorkspaceFactory();
                ITinWorkspace TinWK = TinWF.OpenFromFile(tinFolder, 0) as ITinWorkspace;
                ITinAdvanced2 tinAd = TinWK.OpenTin(tinName) as ITinAdvanced2;
                IEnvelope extent = tinAd.Extent;
                IPoint origin = extent.LowerLeft;
                origin.X = origin.X - (5 * 0.5);
                origin.Y = origin.Y - (5 * 0.5);
                int nCol = (int)Math.Round(extent.Width / 5) + 1;
                int nRow = (int)Math.Round(extent.Height / 5) + 1;

                ISpatialReference2 spatialRef = (ISpatialReference2)extent.SpatialReference;

                IWorkspaceFactory rasterWF = new RasterWorkspaceFactoryClass();
                IRasterWorkspace2 workSpace = (IRasterWorkspace2)rasterWF.OpenFromFile(rasterPath, 0);

                rasterDataset = workSpace.CreateRasterDataset(gridName, "IMAGINE Image", origin, nCol, nRow, 5, 5, 1, ESRI.ArcGIS.Geodatabase.rstPixelType.PT_FLOAT, spatialRef, true);


                IRasterBandCollection bandColl = (IRasterBandCollection)rasterDataset;
                IRasterBand rasterBand = bandColl.Item(0);
                IRawPixels rawPixels = (IRawPixels)rasterBand;
                IPnt blockSize = new DblPntClass();
                blockSize.X = nCol;
                blockSize.Y = nRow;
                IPixelBlock3 pixelBlock = (IPixelBlock3)rawPixels.CreatePixelBlock(blockSize);
                ITinSurface tinSurface = (ITinSurface)tinAd;
                IRasterProps rasterProps = (IRasterProps)rawPixels;
                object nodataFloat;
                //long nodataInt;
                object val = pixelBlock.get_PixelDataByRef(0);          
                double cellsize = 5;
                origin.X = origin.X + (5 * 0.5);
                origin.Y = origin.Y + (5 * nRow) - (5 * 0.5);

                nodataFloat = Convert.ToSingle(rasterProps.NoDataValue.ToString());
                tinSurface.QueryPixelBlock(origin.X, origin.Y, cellsize, cellsize, esriRasterizationType.esriElevationAsRaster, nodataFloat, val);
                IPnt offset = new DblPntClass();
                offset.X = 0;
                offset.Y = 0;
                rawPixels.Write(offset, pixelBlock as IPixelBlock);


            }
            catch (Exception ex)
            {
              

            }
       
            return rasterDataset;

        }

2.3 读取TIN的代码

百度抄的，没有验证

/// <summary>
/// 打开 TIN 文件
/// </summary>
/// <param name="aFolderName">TIN文件夹 名</param>
/// <param name="axMapControl">axMapControl（可已将其换为 Scene或 Globe 控件）</param>
/// <returns></returns>
public void OpenFile_TINFolder(string aFolderName, AxMapControl axMapControl)
{
    string pathToWorkspace = System.IO.Path.GetDirectoryName(aFolderName);
    string tinName = System.IO.Path.GetFileName(aFolderName);
    IWorkspaceFactory pWorkspaceFactory = new TinWorkspaceFactoryClass();
    IWorkspace pWorkspace = pWorkspaceFactory.OpenFromFile(pathToWorkspace, 0);
    ITinWorkspace pTinWorkspace = pWorkspace as ITinWorkspace;
    ITinLayer pTinLayer = new TinLayerClass();
    if (pTinWorkspace.get_IsTin(tinName))
    {
        ITin pTin = pTinWorkspace.OpenTin(tinName);
        pTinLayer.Dataset = pTin;
        pTinLayer.ClearRenderers();
        pTinLayer.Name = tinName;
    }
    axMapControl.AddLayer(pTinLayer);
    axMapControl.ActiveView.Refresh();
}

自己写的点线矛盾程序查看如下图公众号文章[ArcEngine等高线与高程点点线矛盾检查]

X 参考资料

ArcGIS三维入门(2-14)ArcGIS等高线转换DEM的方法