Revit中实现GH Scan凸包算法

YaruCode

已于 2022-08-20 09:53:17 修改

阅读量200

点赞数

分类专栏： Revit二次开发文章标签： Revit 凸包计算几何算法 Revit二次开发

于 2022-08-20 09:46:20 首次发布

本文链接：https://blog.csdn.net/yasenrk/article/details/126436183

版权

Revit二次开发专栏收录该内容

27 篇文章 7 订阅

订阅专栏

该博客详细介绍了计算几何中求解二维平面上点集的凸包问题的算法，提供了一个C#实现的ConvexHull函数。首先检查点的数量，如果少于3个直接返回。接着找到最低点并按角度排序，使用单调链方法构建凸包。此外，还定义了辅助函数用于判断三点是否构成左转。整个过程高效地构建了凸包并返回结果。

摘要由CSDN通过智能技术生成

             
public static IEnumerable<XYZ> ConvexHull(IEnumerable<XYZ> points)
        {
            if (points.Count() < 3)
            {
                return points;
            }
            var inputPoints = new List<XYZ>(points);
            var llt = FindLowestThenLeftPoint(inputPoints);
            var res = new List<XYZ>();

            inputPoints = inputPoints.OrderBy(x => (x - llt).Normalize().AngleTo(XYZ.BasisX)).ToList();

            var stack = inputPoints.ToArray();

            var lowPtr = 1;
            var upPtr = 2;

            while (upPtr != stack.Length)
            {
                var p1 = stack[lowPtr];
                var p2 = stack[lowPtr - 1];

                var candidate = stack[upPtr];

                var left = ToLeft(p2, p1, candidate);
                if (left > 0)
                {
                    stack[++lowPtr] = stack[upPtr++];
                }
                else
                {
                    lowPtr--;
                }
            }
            return stack.ToList().Take(lowPtr + 1);
        }

      
private static XYZ FindLowestThenLeftPoint(IEnumerable<XYZ> points)
        {
            var res = points.FirstOrDefault();

            var enumerator = points.GetEnumerator();
            while (enumerator.MoveNext())
            {
                var temp = enumerator.Current as XYZ;
                if (temp.Y < res.Y)
                {
                    res = temp;
                }
                else if (temp.Y == res.Y)
                {
                    res = temp.X < res.X ? temp : res;
                }
            }
            return res;
        }



public static int ToLeft(XYZ p1, XYZ p2, XYZ p3, double eps = 1e-4)
        {
            var area = p1.X * p2.Y - p1.Y * p2.X
                    + p2.X * p3.Y - p2.Y * p3.X
                    + p3.X * p1.Y - p3.Y * p1.X;
            if (Math.Abs(area) < eps)
            {
                return 0;
            }

            return area > 0 ? 1 : -1;
        }