1

私は CGAL に慣れていませんが、.locate() 関数を使用して三角形内の点を見つけると、いくつかの問題があります。まず、3 点 (-1,0) (1,0) (0,1) で三角形を作成します。次に、いくつかの点の相対位置を三角形で確認します。点は次のとおりです。

p0(0.8, 0.1) ---- locate in the triangle；
p1(0.95, 0) ---- locate on the edge (-1,0) ---- (1,0)；
p2(0.8, 0.2) ---- locate on the edge (1,0) ---- (0,1)；
p3(0.6, 0.4) ---- locate on the edge (1,0) ---- (0,1)；
p4(0.7,0.3) ---- locate on the edge (1,0) ---- (0,1)；

しかし、検索結果は次のとおりです。

p0 ---- FACE, right!
p1 ---- EDGE, right!
p2 ---- OUTSIDE_CONVEX_HULL, wrong!
p3 ---- EDGE, right!
p4 ---- FACE, wrong!

点 p2、p3、p4 が同じエッジ (1,0) ---- (0,1) に位置している間

問題の理由は何ですか?

コード：

typedef CGAL::Exact_predicates_inexact_constructions_kernel K;
typedef CGAL::Delaunay_triangulation_2<K> Delaunay_triangulation;

void main()
{
    //generate the triangles
    std::vector<Point2> pts;
    pts.push_back(Point2(1,0));
    pts.push_back(Point2(-1,0));
    pts.push_back(Point2(0,1));
    Delaunay_triangulation dt;
    dt.insert(pts.begin(), pts.end());

    //locate point in the triangle
    Point2 pt0 = Point2(0.8, 0.1);//the point located in the triangle.(on the FACE ?)
    Point2 pt1 = Point2(0.95, 0);//the point located on the segment ( (1,0) ---- (-1,0) ) (on the EDGE ?)
    Point2 pt2 = Point2(0.8, 0.2);//the point locate on the segment ( (1,0) ---- (0,1) ) (on the EDGE ?)
    Point2 pt3 = Point2(0.6, 0.4);//the point locate on the segment ( (1,0) ---- (0,1) ) (on the EDGE ?)
    Point2 pt4 = Point2(0.7,0.3);//the point locate on the segment ( (1,0) --- (0,1) ) (on the EDGE ?)
    Face_handle fh0, fh1, fh2, fh3, fh4;
    Delaunay_triangulation::Locate_type lt0, lt1, lt2, lt3, lt4;
    int li;
    fh0 = dt.locate(pt0, lt0, li);//for pt0, lt0 is equal to FACE, right!
    fh1 = dt.locate(pt1, lt1, li);//for pt1, lt1 is equal to EDGE, right!
    fh2 = dt.locate(pt2, lt2, li);//for pt2, lt2 is equal to OUTSIDE_CONVEX_HULL, I think it's wrong!
    fh3 = dt.locate(pt3, lt3, li);//for pt3, lt3 is equal to EDGE, right!
    fh4 = dt.locate(pt4, lt4, li);//for pt4, lt3 is equal to FACE, I thing it's wrong!
}
4

1 に答える 1

2

浮動小数点数は、書き込めるすべての 10 進数を表すことはできません。丸め誤差があります。そのため、ポイントの座標はあなたが思うものではありません。次のプログラムを実行します。

#include <iostream>
#include <CGAL/Point_2.h>
#include <CGAL/Simple_cartesian.h>

typedef CGAL::Simple_cartesian<double> Kernel;
typedef Kernel::Point_2 Point;

int main()
{
  std::cout.precision(17);
  Point p0(0.8, 0.1);
  Point p1(0.95, 0);
  Point p2(0.8, 0.2);
  Point p3(0.6, 0.4);
  Point p4(0.7,0.3);
  std::cout << p0 << std::endl;
  std::cout << p1 << std::endl;
  std::cout << p2 << std::endl;
  std::cout << p3 << std::endl;
  std::cout << p4 << std::endl;
  return 0;
}

ポイントの座標を高精度で表示します。出力は次のとおりです。

0.80000000000000004 0.10000000000000001
0.94999999999999996 0
0.80000000000000004 0.20000000000000001
0.59999999999999998 0.40000000000000002
0.69999999999999996 0.29999999999999999

これらの座標は、 によって返されるすべての場所を説明していますlocate()。

于 2013-10-22T09:49:34.233 に答える