2013-11-23 51 views
2

调查this question,我将不得不处理边界由线段和圆弧组成的形状。看来CGAL应该可以在这里帮助我:根据this section of the User's ManualGeneral_polygon_set_2Gps_segment_traits_2作为其特征类应该能够表达大部分我需要的操作,特别是交集和差异。CGAL一般多边形:刚性运动和面积

我到目前为止在文档中没有找到的东西是将刚性运动应用于这些形状的方法,以及计算所得形状的面积的方法。

我想我可以解决这两个问题。对于刚性运动,我可以在转换原始定义对象之后重新创建形状。为了计算面积,我可以使用鞋带方法的差异来调整以适应圆弧。手册中的示例打印了有关支持圆的详细信息,并通过标题进行挖掘,发现我的多边形的每个curve确实有supporting_circle()方法,所以我想它实际上是Arr_circle_segment_traits_2<K>::X_monotone_curve_2。所以我应该能够获得足够的圈子信息来计算面积。在使用故意的编译器错误消息来了解手册简单描述为unspecified_type的某些对象的类型之后,我才在头文件中找到了它。

尽管如此,这两个操作都需要很多工作,而且我很惊讶似乎没有内置的方法来完成这些操作。另一方面,CGAL通过模板参数进行自定义的方式,我可能简单地错过了一种方法来完成这些对圆形段有效的事情,尽管它可能不适用于其他常规多边形。你知道我可能使用的任何快捷方式吗?

+0

嗨。 [我是CGAL开发者](http://www.cgal。组织/ people.html)。我不确定自己的答案,因为我不太清楚CGAL的布尔操作和安排包。我已将您的问题转发给开发人员的内部邮件列表。我希望我们中的一个能够帮助你。 – lrineau

+0

@lrineau:非常感谢!转换输入对象的解决方法可能会在性能方面产生问题。我写了一些代码来逐个变换一个多边形集合,但是在重构变形的多边形时我违反了一些先决条件。可能是因为我正在使用'Epick'内核,并且简单地将一个根点转换为double。目前我正在认真考虑从头开始自己实现这个布尔值。但也许有些开发人员可以帮助我,至少概述CAGL如何扩展以满足我的需求。 – MvG

回答

2

恐怕你不会喜欢我的答案。

我是CGAL开发人员,实际上是Reg的开发人员之一。布尔操作和安排包。

首先,您要求的操作不受支持。

关于面积计算,您的方法似乎是可行的。然而,对于我们来说,要在这个概念中要求这样的操作将是一个很大的努力,因为那样我们就需要为我们支持的所有特质类实施操作。我想,始终可以从一个开始,然后逐个添加它们。我会把它扔进我们的待办事项列表,但我不会把我的赌注放在快速交货...

关于转型,答案更涉及。正如您已经注意到的那样,对(确切的)几何形状(例如,排列,一般多边形集合,或甚至一个小的线性简单凸多边形)应用非精确变换可能是有害的。您必须拿出一个精确的变换,例如,一个变换矩阵,它包含一个确切类型的数字 - 用于表示几何元素的(确切)坐标的相同(或至少可以相互转换)的类型。问题自然就是旋转,因为通常从一个角度开始,并使用三角函数(例如sin()和cos())来计算旋转矩阵。假设你想旋转一个给定的角度,比如alpha。您需要计算α的近似值,使得sin(alpha)和cos(alpha)是有理​​数,因此可以用上述精确类型的数字表示。自由函数CGAL :: rational_rotation_approximation()可以提供帮助。正如在手动输入这个函数中提到的那样,近似值是基于Canny和Ressler在1992年第八届SoCG上提出的有理旋转方法中描述的Farey序列。

祝你好运!

+0

非常感谢您的回答!为了澄清事情,这是否意味着我必须将所有计算转换为合理的内核和有理数?毕竟,如果我将它们与不精确的双打相乘,计算精确的有理几乎旋转矩阵将无济于事,是吗?什么与双打有什么区别?什么不是?对于平等的检查就像我失败的前提条件一样,在某些情况下,它只会在双打中起作用,并可能在其他情况下打破。是这样吗? – MvG

+1

不完全是,您必须使用精确谓词精确构造内核。有支持代数的内核不一定是合理的。我建议你只有在确定合理的内核不足时才使用它们(例如,你需要计算一些数字的根)。但是,不存在支持超越号码的内核,例如PI。如你所说,双打的代码可能会破坏,并且在某些情况下它最终会破坏,而当涉及复杂的数据结构时,比如安排,它就会崩溃。 –

+0

比你想像的更多 –