建议的包/函数来计算3单纯形(四面体)的顶点

时间:2013-09-26 23:33:41

标签: r

我想绘制三维单形,它包含3D中的一些随机点。例如:

pts <- rnorm(30)
pts <- matrix(pts, ncol = 3)

有了这些点,我想计算仅包围这些点的3单形(不规则四面体)的顶点。有人可以推荐一个可以做到这一点的包/功能吗?搜索与单纯形相关的材料的所有方式都由与将单纯形用于其他目的相关的答案所主导,其中有许多。我只是想计算一个并绘制它。看起来很简单,但我似乎并不知道我需要的相关关键字。

2 个答案:

答案 0 :(得分:1)

我认为你在convhulln包中寻找geometry,但我不是专家,所以也许这不是你想要的。

答案 1 :(得分:1)

如果没有人可以为此找到合适的包装,你必须自己解决这个问题,如果你不要求它是绝对最紧密的,那就不那么困难了。请参阅this question over at mathexchange

在这个问题中提出的最简单的方法在我看来是翻译原点,使得所有点都位于正的orthant(即,所有点维度都是正的),然后将点投影到每个单位表示的单一内容中向量。要在原始坐标系中获得此单纯形,您可以采用此单形图中的点的反投影和反向平移。

另一种方法建议找到包络球体(例如你可以使用Ritter的算法),然后找到球体的包络单形,这可能是一个更容易的任务,取决于你最熟悉的。< / p>