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

Question

我在笛卡尔空间中有一组点。我可以使用带有 numpy 的 python 计算由给定的四个点 (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 中以您想要的角度围绕该轴旋转所有点来轻松地做到这一点 - 参见 this example (转到“旋转矢量”)。否则，您需要编程 the appropriate equation (在 this thread 中详细说明)。