凸包のライブラリ

Question

私はC＃を初めて使用しており、凸包を計算するのに苦労しています。 C＃にはこれのための何らかの種類の数学ライブラリがありますか？そうでなければ、私は自分自身を実装しなければならないと思う。

Answer 1

MIConvexHull -https://designengrlab.github.io/MIConvexHull/ - は、C＃の高性能凸包実装であり、高次元凸包もサポートしています。 LGPLライセンス。

Answer 2

ここでは、Monotone Chainアルゴリズム、a.k.a Andrew's Algorithmを使用して書いた2D凸包アルゴリズムです。それが存在すると仮定されているいくつかのものに依存している

public static IListSource<Point> ComputeConvexHull(List<Point> points, bool sortInPlace = false) 
{ 
    if (!sortInPlace) 
     points = new List<Point>(points); 
    points.Sort((a, b) => 
     a.X == b.X ? a.Y.CompareTo(b.Y) : a.X.CompareTo(b.X)); 

    // Importantly, DList provides O(1) insertion at beginning and end 
    DList<Point> hull = new DList<Point>(); 
    int L = 0, U = 0; // size of lower and upper hulls 

    // Builds a hull such that the output polygon starts at the leftmost point. 
    for (int i = points.Count - 1; i >= 0 ; i--) 
    { 
     Point p = points[i], p1; 

     // build lower hull (at end of output list) 
     while (L >= 2 && (p1 = hull.Last).Sub(hull[hull.Count-2]).Cross(p.Sub(p1)) >= 0) { 
      hull.RemoveAt(hull.Count-1); 
      L--; 
     } 
     hull.PushLast(p); 
     L++; 

     // build upper hull (at beginning of output list) 
     while (U >= 2 && (p1 = hull.First).Sub(hull[1]).Cross(p.Sub(p1)) <= 0) 
     { 
      hull.RemoveAt(0); 
      U--; 
     } 
     if (U != 0) // when U=0, share the point added above 
      hull.PushFirst(p); 
     U++; 
     Debug.Assert(U + L == hull.Count + 1); 
    } 
    hull.RemoveAt(hull.Count - 1); 
    return hull; 
} 

は、私のblog post詳細については参照してください。

Answer 3

以下は、Qwertieの回答で使用された同じJavaソースのC＃への翻字ですが、非標準クラスには依存しません。

class ConvexHull 
{ 
    public static double cross(Point O, Point A, Point B) 
    { 
     return (A.X - O.X) * (B.Y - O.Y) - (A.Y - O.Y) * (B.X - O.X); 
    } 

    public static List<Point> GetConvexHull(List<Point> points) 
    { 
     if (points == null) 
      return null; 

     if (points.Count() <= 1) 
      return points; 

     int n = points.Count(), k = 0; 
     List<Point> H = new List<Point>(new Point[2 * n]); 

     points.Sort((a, b) => 
      a.X == b.X ? a.Y.CompareTo(b.Y) : a.X.CompareTo(b.X)); 

     // Build lower hull 
     for (int i = 0; i < n; ++i) 
     { 
      while (k >= 2 && cross(H[k - 2], H[k - 1], points[i]) <= 0) 
       k--; 
      H[k++] = points[i]; 
     } 

     // Build upper hull 
     for (int i = n - 2, t = k + 1; i >= 0; i--) 
     { 
      while (k >= t && cross(H[k - 2], H[k - 1], points[i]) <= 0) 
       k--; 
      H[k++] = points[i]; 
     } 

     return H.Take(k - 1).ToList(); 
    } 
} 

Answer 4

私は提供されたすべてのコードで多くの凸包アルゴリズム/実装を比較しました。すべてはCodeProjectの記事に含まれています。

アルゴリズムを比較：

グラハムが

チャンをスキャン

• モノトーンチェーン
• MiConvexHull（ドロネー三角形分割をし、ボロノイメッシュ）
• Ouellet（鉱山）

の記事：

• 2017年10月13日 - 5月アルゴリズム/実装とテストベンチ：Fast and improved 2D Convex Hull algorithm and its implementation in O(n log h)
• 2014年5月20日 - 私自身のアルゴリズム説明：A Convex Hull Algorithm and its implementation in O(n log h)