
<html><body><div style="color:#000; background-color:#fff; font-family:HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif;font-size:10pt"><div>Hi,</div><div>I'm trying to find the angle between two vectors from three points, like in the vtkAngleWidget. I've checked vtkAngleRepresentation2D.cxx and made the functions below. Is there any simpler way to get the angle?</div><div><br></div><div style="background-color: transparent;">&nbsp;void get_angle_old(double p1[3], double &nbsp;c[3], double p2[3])</div><div style="background-color: transparent;">&nbsp;{</div><div style="background-color: transparent;">&nbsp; &nbsp; double theta1 = atan2(p1[1]-c[1],p1[0]-c[0]);</div><div style="background-color: transparent;">&nbsp; &nbsp; double theta2 = atan2(p2[1]-c[1],p2[0]-c[0]);</div><div style="background-color: transparent;">&nbsp; &nbsp; if ( p1[0]-c[0] == 0.0 || p2[0]-c[0] == 0.0 )</div><div style="background-color: transparent;">&nbsp;

 &nbsp; &nbsp; &nbsp; &nbsp;{</div><div style="background-color: transparent;">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;return;</div><div style="background-color: transparent;">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;}</div><div style="background-color: transparent;">&nbsp; &nbsp; if ( (theta1 &gt;= 0.0 &amp;&amp; theta1 &lt;= vtkMath::Pi() &amp;&amp;</div><div style="background-color: transparent;">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; theta2 &gt;= 0.0 &amp;&amp; theta2 &lt;= vtkMath::Pi()) ||</div><div style="background-color: transparent;">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;(theta1 &lt;= 0.0 &amp;&amp; theta1 &gt;= -vtkMath::Pi() &amp;&amp;</div><div style="background-color: transparent;">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; theta2 &lt;= 0.0 &amp;&amp; theta2 &gt;= -vtkMath::Pi()) )</div><div style="background-color: transparent;">&nbsp; &nbsp; &nbsp; &nbsp; {</div><div style="background-color: transparent;">&nbsp; &nbsp; &nbsp; &nbsp; ; //do nothin

 angles are fine</div><div style="background-color: transparent;">&nbsp; &nbsp; &nbsp; &nbsp; }</div><div style="background-color: transparent;">&nbsp; &nbsp; else if ( theta1 &gt;= 0.0 &amp;&amp; theta2 &lt;= 0.0 )</div><div style="background-color: transparent;">&nbsp; &nbsp; &nbsp; &nbsp; {</div><div style="background-color: transparent;">&nbsp; &nbsp; &nbsp; &nbsp; if ( (theta1 - theta2) &gt;= vtkMath::Pi() )</div><div style="background-color: transparent;">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; {</div><div style="background-color: transparent;">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; theta2 = theta2 + 2.0*vtkMath::Pi();</div><div style="background-color: transparent;">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; }</div><div style="background-color: transparent;">&nbsp; &nbsp; &nbsp; &nbsp; }</div><div style="background-color: transparent;">&nbsp; &nbsp; else //if ( theta1 &lt;= 0.0 &amp;&amp; theta2 &gt;= 0.0 )</div><div style="background-color: transparent;">&nbsp;

 &nbsp; &nbsp; &nbsp; {</div><div style="background-color: transparent;">&nbsp; &nbsp; &nbsp; &nbsp; if ( (theta2 - theta1) &gt;= vtkMath::Pi() )</div><div style="background-color: transparent;">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; {</div><div style="background-color: transparent;">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; theta1 = theta1 + 2.0*vtkMath::Pi();</div><div style="background-color: transparent;">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; }</div><div style="background-color: transparent;">&nbsp; &nbsp; &nbsp; &nbsp; }</div><div style="background-color: transparent;"><span style="background-color: transparent; font-size: 10pt;">&nbsp; &nbsp;std::cout &lt;&lt; "Angle = " &lt;&lt; vtkMath::DegreesFromRadians(abs(theta1-theta2)) &lt;&lt; std::endl;</span><br></div><div style="background-color: transparent;"><span style="background-color: transparent; font-size: 10pt;">&nbsp;}</span><br></div><div style="background-color: transparent;"><br></div><div

 style="background-color: transparent;">&nbsp;void get_angle_new(double p1[3], double c[3], double p2[3])</div><div style="background-color: transparent;">&nbsp;{</div><div style="background-color: transparent;">&nbsp; &nbsp; double vector2[3], vector1[3];</div><div style="background-color: transparent;"><span style="background-color: transparent; font-size: 10pt;">&nbsp; &nbsp; vector1[0] = p1[0] - c[0];</span><br></div><div style="background-color: transparent;">&nbsp; &nbsp; vector1[1] = p1[1] - c[1];</div><div style="background-color: transparent;">&nbsp; &nbsp; vector1[2] = p1[2] - c[2];</div><div style="background-color: transparent;">&nbsp; &nbsp; vector2[0] = p2[0] - c[0];</div><div style="background-color: transparent;">&nbsp; &nbsp; vector2[1] = p2[1] - c[1];</div><div style="background-color: transparent;">&nbsp; &nbsp; vector2[2] = p2[2] - c[2];</div><div style="background-color: transparent;">&nbsp; &nbsp; vtkMath::Normalize( vector1

 );</div><div style="background-color: transparent;">&nbsp; &nbsp; vtkMath::Normalize( vector2 );</div><div style="background-color: transparent;">&nbsp; &nbsp; double angle = acos( vtkMath::Dot( vector1, vector2 ) );</div><div style="background-color: transparent;"><span style="background-color: transparent; font-size: 10pt;">&nbsp; &nbsp; std::cout &lt;&lt; "angle: " &lt;&lt; vtkMath::DegreesFromRadians(angle) &lt;&lt; std::endl;</span><br></div><div style="background-color: transparent;"><span style="background-color: transparent; font-size: 10pt;">&nbsp;}</span><br></div></div></body></html>