Fix definition and calculation of source angle

SimonRit · Feb 9, 2015 · b0dd888 · b0dd888
1 parent 6108d05
commit b0dd888
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 1 deletion.
diff --git a/code/rtkThreeDCircularProjectionGeometry.cxx b/code/rtkThreeDCircularProjectionGeometry.cxx
@@ -93,6 +93,7 @@ void rtk::ThreeDCircularProjectionGeometry::AddProjectionInRadians(
   z.Fill(0.);
   z[2] = 1.;
   HomogeneousVectorType sph = GetSourcePosition( m_GantryAngles.size()-1 );
+  sph[1] = 0.; // Project position to central plane
   VectorType sp( &(sph[0]) );
   sp.Normalize();
   double a = acos(sp*z);

diff --git a/code/rtkThreeDCircularProjectionGeometry.h b/code/rtkThreeDCircularProjectionGeometry.h
@@ -109,7 +109,8 @@ class RTK_EXPORT ThreeDCircularProjectionGeometry : public ProjectionGeometry<3>
   }
 
   /** Get a vector containing the source angles in radians. The source angle is
-   * defined as the angle between the z-axis and the isocenter-source line. */
+   * defined as the angle between the z-axis and the isocenter-source line
+   * projected on the central plane. */
   const std::vector<double> GetSourceAngles() const {
     return this->m_SourceAngles;
   }