查询给定区域内曲面平均高程

2019年5月23日QQ群友聊天，

谈到如何查询曲面某范围内的平均高程，

对于整个曲面的平均高程，

在曲面特性内可以直接查看，

对于给定范围内的平均高程怎没有现成的命令可以实现，

虽然可以创建剪裁曲面，

创建新曲面后进行查看，

但对于多个范围的话，

重复操作也就很烦人了。

所以就又一次体现出学习一点点二次开发知识的必要性了！

曲面有现成的API可以查询给定范围的界内体积，

体积除以底面积，

就应该是平均高程！

就这么简单，

大致估计一下代码的行数，

应该不超过100行。

刚好电脑在手边，

那就试着写一下吧，

说不准这个功能那天自己也能用上。

代码如下：

using Autodesk.AutoCAD.DatabaseServices;
using Autodesk.AutoCAD.EditorInput;
using Autodesk.AutoCAD.ApplicationServices;
using Autodesk.AutoCAD.Geometry;
using Autodesk.Civil.DatabaseServices;

namespace 翻模工具
{
    /// <summary>
    /// 2019年5月23日
    /// 计算选定范围内的平均曲面高程
    /// </summary>
    class MeanHeight
    {
        public void CalMeanHeigt()
        {
            Document doc = Autodesk.AutoCAD.ApplicationServices.Application.DocumentManager.CurrentDocument;
            Database db = doc.Database;
            TypedValue[] tv = new TypedValue[]
               {
                new TypedValue((int)DxfCode.Start,"LWPOLYLINE")
               };
            SelectionFilter sf = new SelectionFilter(tv);
            PromptSelectionResult psr = doc.Editor.GetSelection(sf);
            if (psr.Status != PromptStatus.OK) return;
            SelectionSet ss = psr.Value;
            PromptEntityOptions peo = new PromptEntityOptions("拾取曲面");
            peo.SetRejectMessage("
请选择三角网曲面");
            peo.AddAllowedClass(typeof(TinSurface), true);
            PromptEntityResult per = doc.Editor.GetEntity(peo);
            if (per.Status != PromptStatus.OK) return;
            using (Transaction tr = doc.TransactionManager.StartTransaction())
            {
                BlockTable bt = db.BlockTableId.GetObject(OpenMode.ForRead) as BlockTable;
                BlockTableRecord btr = bt[BlockTableRecord.ModelSpace].GetObject(OpenMode.ForWrite) as BlockTableRecord;
                TinSurface sur = per.ObjectId.GetObject(OpenMode.ForRead) as TinSurface;
                foreach (ObjectId id in ss.GetObjectIds())
                {
                    Polyline pl = id.GetObject(OpenMode.ForRead) as Polyline;
                    Point3dCollection pts = new Point3dCollection();
                    for (int i = 0; i < pl.NumberOfVertices; i++)
                    {
                        pts.Add(pl.GetPoint3dAt(i));
                    }
                    if (pts[0].DistanceTo(pts[pts.Count - 1]) != 0)
                    {
                        pts.Add(pts[0]);
                    }
                    double h = sur.GetBoundedVolumes(pts).Net / pl.Area;
                    Point3d pt1 = pl.GeometricExtents.MaxPoint;
                    Point3d pt2 = pl.GeometricExtents.MinPoint;
                    Point3d pt = pt2 + (pt1 - pt2)/2;
                    DBText t = new DBText();
                    t.TextString = h.ToString("0.000");
                    t.Position = pt;
                    t.HorizontalMode = TextHorizontalMode.TextCenter;
                    t.VerticalMode = TextVerticalMode.TextVerticalMid;
                    t.AlignmentPoint = pt;
                    t.Height = 2.5;
                    t.Color = Autodesk.AutoCAD.Colors.Color.FromColorIndex(Autodesk.AutoCAD.Colors.ColorMethod.ByColor, 1);
                    btr.AppendEntity(t);
                    tr.AddNewlyCreatedDBObject(t, true);
                }
                tr.Commit();
            }
        }
    }
}

当然代码还不完善，

如果要有通用性，

还要进一步完善。

测试结果如下：