diff -ru UnmodifiedInsightToolkit-3.2.0/Code/Algorithms/itkMattesMutualInformationImageToImageMetric.txx InsightToolkit-3.2.0/Code/Algorithms/itkMattesMutualInformationImageToImageMetric.txx
--- UnmodifiedInsightToolkit-3.2.0/Code/Algorithms/itkMattesMutualInformationImageToImageMetric.txx	2006-10-30 17:05:11.000000000 -0500
+++ InsightToolkit-3.2.0/Code/Algorithms/itkMattesMutualInformationImageToImageMetric.txx	2007-05-23 11:13:14.000000000 -0400
@@ -1194,15 +1194,18 @@
   ImageDerivativesType& gradient ) const
 {
   
+  typedef typename MovingImageType::DirectionType DirectionType;
+  const DirectionType & direction = this->GetMovingImage()->GetDirection();
+
   if( m_InterpolatorIsBSpline )
     {
     // Computed moving image gradient using derivative BSpline kernel.
-    gradient = m_BSplineInterpolator->EvaluateDerivative( mappedPoint );
+    gradient = direction * m_BSplineInterpolator->EvaluateDerivative( mappedPoint );
     }
   else
     {
     // For all generic interpolator use central differencing.
-   gradient = m_DerivativeCalculator->Evaluate( mappedPoint );
+   gradient = direction * m_DerivativeCalculator->Evaluate( mappedPoint );
     }
 
 }
diff -ru UnmodifiedInsightToolkit-3.2.0/Code/Algorithms/itkMeanSquaresImageToImageMetric.txx InsightToolkit-3.2.0/Code/Algorithms/itkMeanSquaresImageToImageMetric.txx
--- UnmodifiedInsightToolkit-3.2.0/Code/Algorithms/itkMeanSquaresImageToImageMetric.txx	2004-12-21 17:47:27.000000000 -0500
+++ InsightToolkit-3.2.0/Code/Algorithms/itkMeanSquaresImageToImageMetric.txx	2007-05-23 11:12:30.000000000 -0400
@@ -213,8 +213,11 @@
         mappedIndex[j] = static_cast<long>( vnl_math_rnd( tempIndex[j] ) );
         }
 
+      typedef typename MovingImageType::DirectionType DirectionType;
+      const DirectionType & direction = this->GetMovingImage()->GetDirection();
+
       const GradientPixelType gradient = 
-        this->GetGradientImage()->GetPixel( mappedIndex );
+        direction * this->GetGradientImage()->GetPixel( mappedIndex );
 
       for(unsigned int par=0; par<ParametersDimension; par++)
         {
@@ -347,8 +350,11 @@
         mappedIndex[j] = static_cast<long>( vnl_math_rnd( tempIndex[j] ) );
         }
 
+      typedef typename MovingImageType::DirectionType DirectionType;
+      const DirectionType & direction = this->GetMovingImage()->GetDirection();
+
       const GradientPixelType gradient = 
-        this->GetGradientImage()->GetPixel( mappedIndex );
+        direction * this->GetGradientImage()->GetPixel( mappedIndex );
 
       for(unsigned int par=0; par<ParametersDimension; par++)
         {
diff -ru UnmodifiedInsightToolkit-3.2.0/Code/Algorithms/itkMutualInformationImageToImageMetric.txx InsightToolkit-3.2.0/Code/Algorithms/itkMutualInformationImageToImageMetric.txx
--- UnmodifiedInsightToolkit-3.2.0/Code/Algorithms/itkMutualInformationImageToImageMetric.txx	2006-03-18 23:36:55.000000000 -0500
+++ InsightToolkit-3.2.0/Code/Algorithms/itkMutualInformationImageToImageMetric.txx	2007-05-23 11:13:14.000000000 -0400
@@ -472,9 +472,12 @@
   
   CovariantVector<double,MovingImageDimension> imageDerivatives;
 
+  typedef typename MovingImageType::DirectionType DirectionType;
+  const DirectionType & direction = this->GetMovingImage()->GetDirection();
+
   if ( m_DerivativeCalculator->IsInsideBuffer( mappedPoint ) )
     {
-    imageDerivatives = m_DerivativeCalculator->Evaluate( mappedPoint );
+    imageDerivatives = direction * m_DerivativeCalculator->Evaluate( mappedPoint );
     }
   else
     {