BUG: Disallow zero-valued determinant directions

Disallow zero-valued determinant directions: the previous implementation
was not checking the actual method parameter values.

Modules/Core/Common/include/itkImageBase.hxx

@@ -141,6 +141,12 @@ ImageBase<VImageDimension>::SetDirection(const DirectionType & direction)
 {
   bool modified = false;
 
+  if (vnl_determinant(direction.GetVnlMatrix()) == 0.0)
+  {
+    itkExceptionMacro("Bad direction, determinant is 0. Refusing to change direction from " << this->m_Direction
+                                                                                            << " to " << direction);
+  }
+
   for (unsigned int r = 0; r < VImageDimension; ++r)
   {
     for (unsigned int c = 0; c < VImageDimension; ++c)
@@ -172,11 +178,6 @@ ImageBase<VImageDimension>::ComputeIndexToPhysicalPointMatrices()
     scale[i][i] = this->m_Spacing[i];
   }
 
-  if (vnl_determinant(this->m_Direction.GetVnlMatrix()) == 0.0)
-  {
-    itkExceptionMacro(<< "Bad direction, determinant is 0. Direction is " << this->m_Direction);
-  }
-
   this->m_IndexToPhysicalPoint = this->m_Direction * scale;
   this->m_PhysicalPointToIndex = m_IndexToPhysicalPoint.GetInverse();