Function BuildBaseSeg

Modules/Multilabel/Testing/mitkDICOMSegmentationIOWriterTest.cpp:97–146

Common scaffolding shared by every writer-test seg: 3-slice geometry, a single label with the metadata Validate requires, foreground voxels so dcmqi sees non-empty frames, and optionally the strict-mode identity-tag set. No source-image relation is attached here; callers layer that on top per the case they exercise. withTracking / withSegProperty toggle the per-label DICOM Type 3 / Type 1 metadat

  // withTracking / withSegProperty toggle the per-label DICOM Type 3 / Type 1
  // metadata so individual tests can exercise the writer's absence paths.
  mitk::MultiLabelSegmentation::Pointer BuildBaseSeg(bool stampIdentity = true,
                                                     unsigned char labelValue = 1,
                                                     bool withTracking = true,
                                                     bool withSegProperty = true)
  {
    auto geometryImage = mitk::Image::New();
    unsigned int dim[3] = {4u, 4u, 3u};
    geometryImage->Initialize(mitk::MakeScalarPixelType<mitk::Label::PixelType>(), 3, dim);

    auto seg = mitk::MultiLabelSegmentation::New();
    seg->Initialize(geometryImage);

    auto label = mitk::Label::New();
    label->SetName("L");
    label->SetValue(labelValue);
    label->SetAlgorithmType(mitk::Label::AlgorithmType::MANUAL);
    label->SetAlgorithmName("WriterTest");
    if (withSegProperty)
    {
      label->SetSegmentedPropertyCategory(mitk::DICOMCodeSequence("T-D0050", "SRT", "Tissue"));
      label->SetSegmentedPropertyType(mitk::DICOMCodeSequenceWithModifiers("T-D0050", "SRT", "Tissue"));
    }
    if (withTracking)
    {
      label->SetTrackingID("track-id");
      label->SetTrackingUID("track-uid");
    }
    seg->AddLabel(label, 0, true, true);

    // dcmqi skips empty slices and produces no DICOM SEG if every slice is
    // empty. Stamp a small foreground voxel block so each slice carries
    // the declared label value and the resulting SEG has frames.
    {
      auto groupImage = seg->GetGroupImage(0);
      mitk::ImageWriteAccessor writeAccessor(groupImage);
      auto *pixels = static_cast<mitk::Label::PixelType *>(writeAccessor.GetData());
      const auto dims = groupImage->GetDimensions();
      const auto sliceSize = static_cast<std::size_t>(dims[0]) * dims[1];
      for (unsigned int z = 0; z < dims[2]; ++z)
      {
        pixels[z * sliceSize + 0] = labelValue;
        pixels[z * sliceSize + 1] = labelValue;
      }
    }

    if (stampIdentity)
      StampIdentityTags(seg);

    return seg;
  }

  // Attach a full DICOM-flavoured rule connection to seg with per-slice
  // source SOPInstance + class UIDs. Matches the source-image-relation