找到一個三角形

的內角我有三個頂點（X，Y，Z）找到一個三角形

triangle = [[ 0. 0. 0.], [ 1. 0. 0.], [ 0.5 0.866 0. ]]

一個三角形我使用的答案從link但不知何故，我沒有得到正確的角度。

DEF角（三角形）：

def unit_vector(vector): 
    print(vector) 
    """ Returns the unit vector of the vector. """ 
    return vector/np.linalg.norm(vector) 

def angle_between(v1, v2): 
    """ Returns the angle in radians between vectors 'v1' and 'v2':: 

      >>> angle_between((1, 0, 0), (0, 1, 0)) 
      1.5707963267948966 
      >>> angle_between((1, 0, 0), (1, 0, 0)) 
      0.0 
      >>> angle_between((1, 0, 0), (-1, 0, 0)) 
      3.141592653589793 
    """ 
    v1_u = unit_vector(v1) 
    v2_u = unit_vector(v2) 
    return np.arccos(np.clip(np.dot(v1_u, v2_u), -1.0, 1.0)) 

v1, v2, v3 = triangle[0], triangle[1], triangle[2] 
vec1, vec2, vec3 = v2-v1, v3-v2, v3-v1 
a1, a2, a3 = angle_between(vec1,vec2),angle_between(vec2,vec3),angle_between(vec3,vec1) 


print(a1,a2,a3) 
2.09440780623 1.04722295888 1.04718484736

的角度的總和高於3.1416rad什麼是不可能的。

問題在哪裏？

來源

2017-06-16 Tobias Stein

它可能是整型問題。嘗試將.0添加到您的整數。 – Dan

對於第一個角度，您需要'angle_between（v1-v2，v3-v2）'，而不是'angle_between（v2-v1，v3-v2）'。矢量的方向很重要。 –

你需要得到您的載體的正確方向：

angle_between(vec1,-vec2) # note the sign here

如果向量的符號是錯的，你不要在這個角落讓你的三角形的內角（但「圓周率 - 的角度'）。畫草圖一定會澄清這一點。

來源

2017-06-16 19:53:05

找到一個三角形

回答

相關問題