连接球体表面的两点

Connect two points on the surface of a sphere

我需要连接球体上的两个点,这样线(边)就停留在球体的表面上而不穿过它。

现在我有:

  1. 这个领域:Evenly distributing n points on a sphere

  2. 绘制了边,但它们穿过球体。

  3. 想要的结果:

下面是 spherical linear interpolation or slerp proposed in this answer 的实现:

import numpy as np
import matplotlib.pylab as plt

def slerp(p1, p2, t):
    omega = np.arccos( p1.dot(p2) )
    sin_omega = np.sin(omega)    
    t = t[:, np.newaxis]
    return ( np.sin( (1-t)*omega )*p1 + np.sin( t*omega )*p2 )/sin_omega

p1 = np.array([1, 0, 0])
p2 = np.array([0, 1, 0])
t = np.linspace(0, 1, 30)

arc = slerp(p1, p2, t)

plt.plot( arc[:, 0], arc[:, 1] );
plt.axis('square');

在 2D 中给出: