Method InitializeByLabeledImageProcessing

Modules/Multilabel/src/mitkLabelSetImage.cpp:1177–1246

template <typename MultiLabelSegmentationType, typename ImageType>
void mitk::MultiLabelSegmentation::InitializeByLabeledImageProcessing(MultiLabelSegmentationType *labelSetImage, const ImageType *image)
{
  typedef itk::ImageRegionConstIteratorWithIndex<ImageType> SourceIteratorType;
  typedef itk::ImageRegionIterator<MultiLabelSegmentationType> TargetIteratorType;

  TargetIteratorType targetIter(labelSetImage, labelSetImage->GetRequestedRegion());
  targetIter.GoToBegin();

  SourceIteratorType sourceIter(image, image->GetRequestedRegion());
  sourceIter.GoToBegin();

  mitk::MultiLabelSegmentation::LabelValueType maxLabelValue = 0;
  std::vector<mitk::Label*> addedLabels;

  while (!sourceIter.IsAtEnd())
  {
    const auto originalSourceValue = sourceIter.Get();
    const auto sourceValue = static_cast<LabelValueType>(originalSourceValue);

    if (originalSourceValue > mitk::Label::MAX_LABEL_VALUE)
    {
      mitkThrow() << "Cannot initialize MultiLabelSegmentation by image. Image contains a pixel "
                  << "value that exceeds the label value range.Invalid pixel value : "
                  << originalSourceValue;
    }

    targetIter.Set(sourceValue);

    if (MultiLabelSegmentation::UNLABELED_VALUE!=sourceValue && !this->ExistLabel(sourceValue))
    {
      if (this->GetTotalNumberOfLabels() >= mitk::Label::MAX_LABEL_VALUE)
      {
        mitkThrow() << "Cannot initialize MultiLabelSegmentation by image. "
                    << "Image contains to many labels.";
      }

      double rgba[4];
      this->GetLookupTable()->GetTableValue(sourceValue, rgba);

      mitk::Color color;
      color.SetRed(rgba[0]);
      color.SetGreen(rgba[1]);
      color.SetBlue(rgba[2]);

      auto label = mitk::Label::New();
      label->SetColor(color);
      label->SetOpacity(rgba[3]);
      label->SetValue(sourceValue);

      addedLabels.push_back(label); // this is used for faster access in next loop.
      //we need to do the adding here already to ensure correct checking and correction of labels
      this->AddLabel(label,0,false);

      maxLabelValue = std::max(maxLabelValue, sourceValue);
    }

    ++sourceIter;
    ++targetIter;