python - 如何从python中的二面角计算笛卡尔坐标

Question

我在笛卡尔空间中有一组点。我可以使用 python 和 numpy 计算由四个点（a、b、c、d）的给定子集合定义的二面角。以下是我的功能：

def getDihedral(a,b,c,d):
 v1 = getNormedVector(a, b)
 v2 = getNormedVector(b, c)
 v3 = getNormedVector(c, d)
 v1v2 = numpy.cross(v1,v2)
 v2v3 = numpy.cross(v2,v3)
 return getAngle(v1v2,v2v3)

def getNormedVector(a,b):
 return (b-a)/numpy.linalg.norm(b-a)

def getAngle(a,b):
 return numpy.rad2deg(numpy.arccos(numpy.dot(a/numpy.linalg.norm(a),b.T/numpy.linalg.norm(b))))[0,0]

我只想旋转一个二面角，如何使用带有 numpy 和 scipy 的 python 计算子集合的新坐标？

Answer 1

如果您可以计算二面角，我假设您可以获得要围绕其旋转点子集的轴。鉴于此，您可以通过在 vpython 中以您想要的角度围绕该轴旋转所有点来轻松做到这一点 - 请参阅此示例（转到“旋转矢量”）。否则，您需要编写适当的方程式（在此线程中详细说明）。