void DiffusionMultiShellQballReconstructionImageFilter<T,TG,TO,L,NODF>
::SetGradientImage(const GradientDirectionContainerType *gradientDirection
, const GradientImagesType *gradientImage
, float bvalue)
{
m_BValue = bvalue;
m_NumberOfBaselineImages = 0;
this->m_GradientDirectionContainer = GradientDirectionContainerType::New();
for(GradientDirectionContainerType::ConstIterator it = gradientDirection->Begin();
it != gradientDirection->End(); it++)
{
this->m_GradientDirectionContainer->push_back(it.Value());
}
if(m_BValueMap.size() == 0) {
itkWarningMacro(<< "DiffusionMultiShellQballReconstructionImageFilter.cpp : no GradientIndexMapAvalible");
GradientDirectionContainerType::ConstIterator gdcit;
for( gdcit = m_GradientDirectionContainer->Begin(); gdcit != m_GradientDirectionContainer->End(); ++gdcit)
{
double bValueKey = int(((m_BValue * gdcit.Value().two_norm() * gdcit.Value().two_norm())+7.5)/10)*10;
m_BValueMap[bValueKey].push_back(gdcit.Index());
}
}
void RadialMultishellToSingleshellImageFilter<TInputScalarType, TOutputScalarType>
::BeforeThreadedGenerateData()
{
// test whether BvalueMap contains all necessary information
if(m_BValueMap.size() == 0)
{
itkWarningMacro(<< "No BValueMap given: create one using GradientDirectionContainer");
GradientDirectionContainerType::ConstIterator gdcit;
for( gdcit = m_OriginalGradientDirections->Begin(); gdcit != m_OriginalGradientDirections->End(); ++gdcit)
{
double bValueKey = int(((m_OriginalBValue * gdcit.Value().two_norm() * gdcit.Value().two_norm())+7.5)/10)*10;
m_BValueMap[bValueKey].push_back(gdcit.Index());
}
}