X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;f=registration%2FclitkConvertBLUTCoeffsToVFFilter.txx;h=ae1f9d539bd70d9425efbab1a760e324ba93771a;hb=599d581f20dc1f11699c722d2cf40b958c517c74;hp=86659936a098aad33d09925c3425d6c6cea26867;hpb=0e4f77b4892597234d90ee974044aa9f5a81ce89;p=clitk.git diff --git a/registration/clitkConvertBLUTCoeffsToVFFilter.txx b/registration/clitkConvertBLUTCoeffsToVFFilter.txx index 8665993..ae1f9d5 100644 --- a/registration/clitkConvertBLUTCoeffsToVFFilter.txx +++ b/registration/clitkConvertBLUTCoeffsToVFFilter.txx @@ -162,8 +162,23 @@ namespace clitk component_filter[i]->Update(); coefficient_images[i] = component_filter[i]->GetOutput(); } + #if ITK_VERSION_MAJOR >= 4 - m_ITKTransform->SetCoefficientImages(coefficient_images); + // RP: 16/01/2013 + // ATTENTION: Apparently, there's a bug in the SetCoefficientImages function of ITK 4.x + // I needed to use SetParametersByValue instead. + // + typename ITKTransformType::ParametersType params(input->GetPixelContainer()->Size() * BLUTCoefficientImageType::ImageDimension); + for (unsigned int i=0; i < BLUTCoefficientImageType::ImageDimension; i++) { + for (unsigned int j=0; j < coefficient_images[i]->GetPixelContainer()->Size(); j++) + params[input->GetPixelContainer()->Size() * i + j] = coefficient_images[i]->GetPixelContainer()->GetBufferPointer()[j]; + } + + m_ITKTransform->SetGridOrigin(input->GetOrigin()); + m_ITKTransform->SetGridDirection(input->GetDirection()); + m_ITKTransform->SetGridRegion(input->GetLargestPossibleRegion()); + m_ITKTransform->SetGridSpacing(input->GetSpacing()); + m_ITKTransform->SetParametersByValue(params); #else m_ITKTransform->SetCoefficientImage(coefficient_images); #endif @@ -178,8 +193,8 @@ namespace clitk m_Filter->Update(); - SetNthOutput(0, m_Filter->GetOutput()); + this->SetNthOutput(0, m_Filter->GetOutput()); } -} \ No newline at end of file +}