- html - 出于某种原因,IE8 对我的 Sass 文件中继承的 html5 CSS 不友好?
- JMeter 在响应断言中使用 span 标签的问题
- html - 在 :hover and :active? 上具有不同效果的 CSS 动画
- html - 相对于居中的 html 内容固定的 CSS 重复背景?
我正在使用 CGAL 计算两个凸多边形相交的面积。在对 this 的接受答案中发布了执行此操作的简短演示代码。问题。但是,当我修改该代码以使用我感兴趣的多边形时,CGAL 从 CGAL::intersection() 例程的深处抛出运行时异常。
这是一个简短的示例代码,它是从上面链接的 SO 问题中复制粘贴的,除了它使用我自己的多边形并打印一些关于每个多边形的诊断信息以表明它们是凸面的并使用 CCW 绕组订单。
#include <iostream>
#include <CGAL/Simple_cartesian.h>
#include <CGAL/Polygon_2.h>
#include <CGAL/Polygon_with_holes_2.h>
#include <CGAL/Boolean_set_operations_2.h>
#include <CGAL/Polygon_2_algorithms.h>
typedef CGAL::Simple_cartesian<double> K;
typedef K::Point_2 Point;
typedef CGAL::Polygon_2<K> Polygon_2;
typedef CGAL::Polygon_with_holes_2<K> Polygon_with_holes_2;
using std::cout; using std::endl;
int main(){
bool IsSimple, IsConvex, IsClockwise;
double Area;
// Define the polygons. P is long and skinny, Q is 20 equally
// spaced points around a circle of radius 0.025. The top-right
// corner of P protudes into Q.
int np=7;
int nq=20;
Point points[]={
Point( 2.67905949256324049657e-02,-2.36307305336680428809e-02),
Point( 2.19804738309115066386e-02, 1.25954939609798088895e-02),
Point( 3.24142504626934169210e-03, 8.27109808759729503436e-03),
Point(-7.07297020897392769712e-02,-1.31780902623212625713e-01),
Point(-1.72945875771868884385e+00,-4.33625989924519217311e+00),
Point(-5.64696462817379796206e+00,-2.91346747524371210147e+01),
Point(-5.66016394049875870564e+00,-4.95259849820236013329e+01)};
Point points2[]={
Point( 2.50000000000000013878e-02, 0.00000000000000000000e+00),
Point( 2.37764129073788424429e-02, 7.72542485937368662852e-03),
Point( 2.02254248593736890571e-02, 1.46946313073118266929e-02),
Point( 1.46946313073118284276e-02, 2.02254248593736890571e-02),
Point( 7.72542485937368836324e-03, 2.37764129073788424429e-02),
Point( 1.53080849893419158600e-18, 2.49999999999999979183e-02),
Point(-7.72542485937368489379e-03, 2.37764129073788424429e-02),
Point(-1.46946313073118232234e-02, 2.02254248593736890571e-02),
Point(-2.02254248593736890571e-02, 1.46946313073118284276e-02),
Point(-2.37764129073788424429e-02, 7.72542485937368923060e-03),
Point(-2.50000000000000013878e-02, 3.06161699786838317201e-18),
Point(-2.37764129073788424429e-02,-7.72542485937368402643e-03),
Point(-2.02254248593736925266e-02,-1.46946313073118232234e-02),
Point(-1.46946313073118318970e-02,-2.02254248593736890571e-02),
Point(-7.72542485937369096533e-03,-2.37764129073788389734e-02),
Point(-4.59242549680257437283e-18,-2.50000000000000013878e-02),
Point( 7.72542485937368229171e-03,-2.37764129073788424429e-02),
Point( 1.46946313073118214887e-02,-2.02254248593736925266e-02),
Point( 2.02254248593736855877e-02,-1.46946313073118318970e-02),
Point( 2.37764129073788389734e-02,-7.72542485937369183269e-03)};
Polygon_2 polyP(points, points+np);
Polygon_2 polyQ(points2, points2+nq);
// Print some properties of polygon P. It is simple, convex, and CCW
// oriented
CGAL::set_pretty_mode(cout);
cout << "created the polygon P:" << endl;
cout << polyP << endl;
cout << endl;
IsSimple = polyP.is_simple();
IsConvex = polyP.is_convex();
IsClockwise = (polyP.orientation() == CGAL::CLOCKWISE);
Area = polyP.area();
cout << "polygon P is";
if (!IsSimple) cout << " not";
cout << " simple." << endl;
cout << "polygon P is";
if (!IsConvex) cout << " not";
cout << " convex." << endl;
cout << "polygon P is";
if (!IsClockwise) cout << " not";
cout << " clockwise oriented." << endl;
cout << "the area of polygon P is " << Area << endl;
cout << endl;
// Print some properties of polygon Q. It is simple, convex, and CCW
// oriented
cout << "created the polygon Q:" << endl;
cout << polyQ << endl;
cout << endl;
IsSimple = polyQ.is_simple();
IsConvex = polyQ.is_convex();
IsClockwise = (polyQ.orientation() == CGAL::CLOCKWISE);
Area = polyQ.area();
cout << "polygon Q is";
if (!IsSimple) cout << " not";
cout << " simple." << endl;
cout << "polygon Q is";
if (!IsConvex) cout << " not";
cout << " convex." << endl;
cout << "polygon Q is";
if (!IsClockwise) cout << " not";
cout << " clockwise oriented." << endl;
cout << "the area of polygon Q is " << Area << endl;
cout << endl;
// Compute the intersection
std::list<Polygon_with_holes_2> polyI;
CGAL::intersection(polyP, polyQ, std::back_inserter(polyI));
double totalArea = 0;
typedef std::list<Polygon_with_holes_2>::iterator LIT;
for(LIT lit = polyI.begin(); lit!=polyI.end(); lit++){
totalArea+=lit->outer_boundary().area();
}
cout << "TotalArea::" << totalArea;
return 0;
}
这是测试程序的输出。
created the polygon P:
Polygon_2(
PointC2(0.0267906, -0.0236307)
PointC2(0.0219805, 0.0125955)
PointC2(0.00324143, 0.0082711)
PointC2(-0.0707297, -0.131781)
PointC2(-1.72946, -4.33626)
PointC2(-5.64696, -29.1347)
PointC2(-5.66016, -49.526)
)
polygon P is simple.
polygon P is convex.
polygon P is not clockwise oriented.
the area of polygon P is 71.1027
created the polygon Q:
Polygon_2(
PointC2(0.025, 0)
PointC2(0.0237764, 0.00772542)
PointC2(0.0202254, 0.0146946)
PointC2(0.0146946, 0.0202254)
PointC2(0.00772542, 0.0237764)
PointC2(1.53081e-18, 0.025)
PointC2(-0.00772542, 0.0237764)
PointC2(-0.0146946, 0.0202254)
PointC2(-0.0202254, 0.0146946)
PointC2(-0.0237764, 0.00772542)
PointC2(-0.025, 3.06162e-18)
PointC2(-0.0237764, -0.00772542)
PointC2(-0.0202254, -0.0146946)
PointC2(-0.0146946, -0.0202254)
PointC2(-0.00772542, -0.0237764)
PointC2(-4.59243e-18, -0.025)
PointC2(0.00772542, -0.0237764)
PointC2(0.0146946, -0.0202254)
PointC2(0.0202254, -0.0146946)
PointC2(0.0237764, -0.00772542)
)
polygon Q is simple.
polygon Q is convex.
polygon Q is not clockwise oriented.
the area of polygon Q is 0.00193136
terminate called after throwing an instance of 'CGAL::Precondition_exception'
what(): CGAL ERROR: precondition violation!
Expr: (m_traits.compare_y_at_x_2_object()(p, cv) == EQUAL) && compare_xy(cv.left(), p) == SMALLER && compare_xy(cv.right(), p) == LARGER
File: /usr/local/include/CGAL/Arr_segment_traits_2.h
Line: 706
Aborted (core dumped)
我不明白是什么导致了错误。任何帮助将不胜感激。
最佳答案
我可以重现此错误(在带有 clang++ 的 MacOS 上使用 CGAL 4.9)。据我了解,这种类型的未捕获异常不应该发生,换句话说,您发现了 CGAL 中的错误。因此,请按照错误消息中的说明提交错误报告–– 您没有发布的部分(或者可能因为版本不同而没有发布?):
Explanation: Refer to the bug-reporting instructions at http://www.cgal.org/bug_report.html
libc++abi.dylib: terminating with uncaught exception of type CGAL::Precondition_exception: CGAL ERROR: precondition violation!
Expr: (m_traits.compare_y_at_x_2_object()(p, cv) == EQUAL) && compare_xy(cv.left(), p) == SMALLER && compare_xy(cv.right(), p) == LARGER
File: /usr/local/include/CGAL/Arr_segment_traits_2.h Line: 706
据我从该文件中所见,函数 throwing 在给定分割点的情况下将一条曲线分成两条子曲线。前提是分割点在曲线上,而不是终点。显然,是否满足这个前提条件是调用者的问题,另一个 CGAL 例程。换句话说,您的输入没有任何问题。问题出在 CGAL 实现上,或者更准确地说,是它处理所用数字表示不精确的方式。
根据 sloriot 的评论进行编辑。
如果我更换
typedef CGAL::Simple_cartesian<double> K;
与
typedef CGAL::Exact_predicates_exact_constructions_kernel K;
并为 Area
和 totalArea
使用适当的类型(我只是使用了 auto
和 decltype(Area)
,分别),代码编译(你必须将它链接到 libgmp
和 libmpfr
)并且运行没有崩溃,报告
totalArea::0.000876944
关于c++ - 在 CGAL 中与多边形相交时的前提条件异常,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/42009418/
SQLite、Content provider 和 Shared Preference 之间的所有已知区别。 但我想知道什么时候需要根据情况使用 SQLite 或 Content Provider 或
警告:我正在使用一个我无法完全控制的后端,所以我正在努力解决 Backbone 中的一些注意事项,这些注意事项可能在其他地方更好地解决......不幸的是,我别无选择,只能在这里处理它们! 所以,我的
我一整天都在挣扎。我的预输入搜索表达式与远程 json 数据完美配合。但是当我尝试使用相同的 json 数据作为预取数据时,建议为空。点击第一个标志后,我收到预定义消息“无法找到任何内容...”,结果
我正在制作一个模拟 NHL 选秀彩票的程序,其中屏幕右侧应该有一个 JTextField,并且在左侧绘制弹跳的选秀球。我创建了一个名为 Ball 的类,它实现了 Runnable,并在我的主 Draf
这个问题已经有答案了: How can I calculate a time span in Java and format the output? (18 个回答) 已关闭 9 年前。 这是我的代码
我有一个 ASP.NET Web API 应用程序在我的本地 IIS 实例上运行。 Web 应用程序配置有 CORS。我调用的 Web API 方法类似于: [POST("/API/{foo}/{ba
我将用户输入的时间和日期作为: DatePicker dp = (DatePicker) findViewById(R.id.datePicker); TimePicker tp = (TimePic
放宽“邻居”的标准是否足够,或者是否有其他标准行动可以采取? 最佳答案 如果所有相邻解决方案都是 Tabu,则听起来您的 Tabu 列表的大小太长或您的释放策略太严格。一个好的 Tabu 列表长度是
我正在阅读来自 cppreference 的代码示例: #include #include #include #include template void print_queue(T& q)
我快疯了,我试图理解工具提示的行为,但没有成功。 1. 第一个问题是当我尝试通过插件(按钮 1)在点击事件中使用它时 -> 如果您转到 Fiddle,您会在“内容”内看到该函数' 每次点击都会调用该属
我在功能组件中有以下代码: const [ folder, setFolder ] = useState([]); const folderData = useContext(FolderContex
我在使用预签名网址和 AFNetworking 3.0 从 S3 获取图像时遇到问题。我可以使用 NSMutableURLRequest 和 NSURLSession 获取图像,但是当我使用 AFHT
我正在使用 Oracle ojdbc 12 和 Java 8 处理 Oracle UCP 管理器的问题。当 UCP 池启动失败时,我希望关闭它创建的连接。 当池初始化期间遇到 ORA-02391:超过
关闭。此题需要details or clarity 。目前不接受答案。 想要改进这个问题吗?通过 editing this post 添加详细信息并澄清问题. 已关闭 9 年前。 Improve
引用这个plunker: https://plnkr.co/edit/GWsbdDWVvBYNMqyxzlLY?p=preview 我在 styles.css 文件和 src/app.ts 文件中指定
为什么我的条形这么细?我尝试将宽度设置为 1,它们变得非常厚。我不知道还能尝试什么。默认厚度为 0.8,这是应该的样子吗? import matplotlib.pyplot as plt import
当我编写时,查询按预期执行: SELECT id, day2.count - day1.count AS diff FROM day1 NATURAL JOIN day2; 但我真正想要的是右连接。当
我有以下时间数据: 0 08/01/16 13:07:46,335437 1 18/02/16 08:40:40,565575 2 14/01/16 22:2
一些背景知识 -我的 NodeJS 服务器在端口 3001 上运行,我的 React 应用程序在端口 3000 上运行。我在 React 应用程序 package.json 中设置了一个代理来代理对端
我面临着一个愚蠢的问题。我试图在我的 Angular 应用程序中延迟加载我的图像,我已经尝试过这个2: 但是他们都设置了 src attr 而不是 data-src,我在这里遗漏了什么吗?保留 d
我是一名优秀的程序员,十分优秀!