+
+ // Account for direction in transform. The offset is already accounted for
+ // in the VTK image coordinates, no need to put it in the transform.
+ for(unsigned int j=0; j<input->GetImageDimension(); j++)
+ for(unsigned int i=0; i<input->GetImageDimension(); i++)
+ mTransform->GetMatrix()->SetElement(i,j, input->GetDirection()[i][j]);
+
+ // Create the corresponding transformed image