void mitk::ImageStatisticsHolder::ComputeImageStatistics(int t, unsigned int component)
{
// timestep valid?
if (!m_Image->IsValidTimeStep(t)) return;
// image modified?
if (this->m_Image->GetMTime() > m_LastRecomputeTimeStamp.GetMTime())
this->ResetImageStatistics();
Expand(t+1);
// do we have valid information already?
if( m_ScalarMin[t] != itk::NumericTraits<ScalarType>::max() ||
m_Scalar2ndMin[t] != itk::NumericTraits<ScalarType>::max() ) return; // Values already calculated before...
// used to avoid statistics calculation on qball images. property will be replaced as soons as bug 17928 is merged and the diffusion image refactoring is complete.
mitk::BoolProperty* isqball = dynamic_cast< mitk::BoolProperty* >( m_Image->GetProperty( "IsQballImage" ).GetPointer() );
const mitk::PixelType pType = m_Image->GetPixelType(0);
if(pType.GetNumberOfComponents() == 1 && (pType.GetPixelType() != itk::ImageIOBase::UNKNOWNPIXELTYPE) && (pType.GetPixelType() != itk::ImageIOBase::VECTOR) )
{
// recompute
mitk::ImageTimeSelector::Pointer timeSelector = this->GetTimeSelector();
if(timeSelector.IsNotNull())
{
timeSelector->SetTimeNr(t);
timeSelector->UpdateLargestPossibleRegion();
const mitk::Image* image = timeSelector->GetOutput();
AccessByItk_2( image, _ComputeExtremaInItkImage, this, t );
}
}
else if (pType.GetPixelType() == itk::ImageIOBase::VECTOR && (!isqball || !isqball->GetValue())) // we have a vector image
{
// recompute
mitk::ImageTimeSelector::Pointer timeSelector = this->GetTimeSelector();
if(timeSelector.IsNotNull())
{
timeSelector->SetTimeNr(t);
timeSelector->UpdateLargestPossibleRegion();
const mitk::Image* image = timeSelector->GetOutput();
AccessVectorPixelTypeByItk_n( image, _ComputeExtremaInItkVectorImage, (this, t, component) );
}
}
else
{
m_ScalarMin[t] = 0;
m_ScalarMax[t] = 255;
m_Scalar2ndMin[t] = 0;
m_Scalar2ndMax[t] = 255;
}
}
mitk::LabelSetImage::Pointer mitk::LabelSetImageConverter::ConvertImageToLabelSetImage(const mitk::Image::Pointer input)
{
mitk::LabelSetImage::Pointer output;
if (input->GetChannelDescriptor().GetPixelType().GetPixelType() == itk::ImageIOBase::VECTOR)
{
AccessVectorPixelTypeByItk_n(input, MitkImageToMitkLabelSetImage, (output));
}
else
{
output = mitk::LabelSetImage::New();
output->InitializeByLabeledImage( input );
}
return output;
}
