[ITK-dev] itkImageToImageFilter issue with negative space from Analyze file

Seun Odutola seun at rogue-research.com
Thu Nov 26 16:04:04 EST 2015


Hi Everyone,

      I have an Anaylze 7.5 file with negative spacing in it and previously my unit test for this file passed in ITK 3.20.x which is now failing in ITK 4.9 (git master). I happen to have noticed there is a tolerance check in the file itkImageToImageFilter.hxx which in my opinion needs to be of absolute value. My question is, is the calculation of the tolerance deliberate?  Note that for this file, inputPtr1->GetSpacing()[0] gives -1.

const SpacePrecisionType coordinateTol
        = this->m_CoordinateTolerance * inputPtr1->GetSpacing()[0]; // use first dimension spacing
        
      if ( !inputPtr1->GetOrigin().GetVnlVector().is_equal(inputPtrN->GetOrigin().GetVnlVector(), coordinateTol) ||
           !inputPtr1->GetSpacing().GetVnlVector().is_equal(inputPtrN->GetSpacing().GetVnlVector(), coordinateTol) ||
           !inputPtr1->GetDirection().GetVnlMatrix().as_ref().is_equal(inputPtrN->GetDirection().GetVnlMatrix(), this->m_DirectionTolerance) )
        {
        std::ostringstream originString, spacingString, directionString;
        if ( !inputPtr1->GetOrigin().GetVnlVector().is_equal(inputPtrN->GetOrigin().GetVnlVector(), coordinateTol) )
          {
          originString.setf( std::ios::scientific );
          originString.precision( 7 );
          originString << "InputImage Origin: " << inputPtr1->GetOrigin()
                       << ", InputImage" << it.GetName() << " Origin: " << inputPtrN->GetOrigin() << std::endl;
          originString << "\tTolerance: " << coordinateTol << std::endl;
          }
        

………….


Regards,
  Seun.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://public.kitware.com/pipermail/insight-developers/attachments/20151126/074d9478/attachment.html>


More information about the Insight-developers mailing list