void mitk::FastMarchingTool::ConfirmSegmentation()
{
// combine preview image with current working segmentation
if (dynamic_cast<mitk::Image*>(m_ResultImageNode->GetData()))
{
//logical or combination of preview and segmentation slice
OutputImageType::Pointer workingImageSliceInITK = OutputImageType::New();
mitk::Image::Pointer workingImageSlice;
mitk::Image::Pointer workingImage = dynamic_cast<mitk::Image*>(this->m_ToolManager->GetWorkingData(0)->GetData());
if(workingImage->GetTimeGeometry()->CountTimeSteps() > 1)
{
mitk::ImageTimeSelector::Pointer timeSelector = mitk::ImageTimeSelector::New();
timeSelector->SetInput( workingImage );
timeSelector->SetTimeNr( m_CurrentTimeStep );
timeSelector->UpdateLargestPossibleRegion();
// todo: make GetAffectedWorkingSlice dependant of current time step
workingImageSlice = GetAffectedWorkingSlice( m_PositionEvent );
CastToItkImage( workingImageSlice, workingImageSliceInITK );
}
else
{
workingImageSlice = GetAffectedWorkingSlice( m_PositionEvent );
CastToItkImage( workingImageSlice, workingImageSliceInITK );
}
typedef itk::OrImageFilter<OutputImageType, OutputImageType> OrImageFilterType;
OrImageFilterType::Pointer orFilter = OrImageFilterType::New();
orFilter->SetInput(0, m_ThresholdFilter->GetOutput());
orFilter->SetInput(1, workingImageSliceInITK);
orFilter->Update();
mitk::Image::Pointer segmentationResult = mitk::Image::New();
mitk::CastToMitkImage(orFilter->GetOutput(), segmentationResult);
segmentationResult->GetGeometry()->SetOrigin(workingImageSlice->GetGeometry()->GetOrigin());
segmentationResult->GetGeometry()->SetIndexToWorldTransform(workingImageSlice->GetGeometry()->GetIndexToWorldTransform());
//write to segmentation volume and hide preview image
// again, current time step is not considered
this->WriteBackSegmentationResult(m_PositionEvent, segmentationResult );
this->m_ResultImageNode->SetVisibility(false);
this->ClearSeeds();
}
mitk::RenderingManager::GetInstance()->RequestUpdateAll();
m_ToolManager->ActivateTool(-1);
}
void mitk::FastMarchingTool3D::ConfirmSegmentation()
{
// combine preview image with current working segmentation
if (dynamic_cast<mitk::Image*>(m_ResultImageNode->GetData()))
{
//logical or combination of preview and segmentation slice
OutputImageType::Pointer segmentationImageInITK = OutputImageType::New();
mitk::Image::Pointer workingImage = dynamic_cast<mitk::Image*>(GetTargetSegmentationNode()->GetData());
if(workingImage->GetTimeSlicedGeometry()->GetTimeSteps() > 1)
{
mitk::ImageTimeSelector::Pointer timeSelector = mitk::ImageTimeSelector::New();
timeSelector->SetInput( workingImage );
timeSelector->SetTimeNr( m_CurrentTimeStep );
timeSelector->UpdateLargestPossibleRegion();
CastToItkImage( timeSelector->GetOutput(), segmentationImageInITK );
}
else
{
CastToItkImage( workingImage, segmentationImageInITK );
}
typedef itk::OrImageFilter<OutputImageType, OutputImageType> OrImageFilterType;
OrImageFilterType::Pointer orFilter = OrImageFilterType::New();
orFilter->SetInput(0, m_ThresholdFilter->GetOutput());
orFilter->SetInput(1, segmentationImageInITK);
orFilter->Update();
//set image volume in current time step from itk image
workingImage->SetVolume( (void*)(m_ThresholdFilter->GetOutput()->GetPixelContainer()->GetBufferPointer()), m_CurrentTimeStep);
this->m_ResultImageNode->SetVisibility(false);
this->ClearSeeds();
workingImage->Modified();
}
mitk::RenderingManager::GetInstance()->RequestUpdateAll();
}
void mitk::FastMarchingTool::ConfirmSegmentation()
{
// combine preview image with current working segmentation
if (dynamic_cast<mitk::Image *>(m_ResultImageNode->GetData()))
{
// logical or combination of preview and segmentation slice
mitk::Image::Pointer workingImageSlice;
mitk::Image::Pointer workingImage = dynamic_cast<mitk::Image *>(this->m_ToolManager->GetWorkingData(0)->GetData());
workingImageSlice = GetAffectedImageSliceAs2DImage(m_WorkingPlane, workingImage, m_CurrentTimeStep);
mitk::Image::Pointer segmentationResult = mitk::Image::New();
bool isDeprecatedUnsignedCharSegmentation =
(workingImage->GetPixelType().GetComponentType() == itk::ImageIOBase::UCHAR);
if (isDeprecatedUnsignedCharSegmentation)
{
typedef itk::Image<unsigned char, 2> OutputUCharImageType;
OutputUCharImageType::Pointer workingImageSliceInITK = OutputUCharImageType::New();
CastToItkImage(workingImageSlice, workingImageSliceInITK);
typedef itk::OrImageFilter<OutputImageType, OutputUCharImageType, OutputUCharImageType> OrImageFilterType;
OrImageFilterType::Pointer orFilter = OrImageFilterType::New();
orFilter->SetInput1(m_ThresholdFilter->GetOutput());
orFilter->SetInput2(workingImageSliceInITK);
orFilter->Update();
mitk::CastToMitkImage(orFilter->GetOutput(), segmentationResult);
}
else
{
OutputImageType::Pointer workingImageSliceInITK = OutputImageType::New();
CastToItkImage(workingImageSlice, workingImageSliceInITK);
typedef itk::OrImageFilter<OutputImageType, OutputImageType> OrImageFilterType;
OrImageFilterType::Pointer orFilter = OrImageFilterType::New();
orFilter->SetInput(0, m_ThresholdFilter->GetOutput());
orFilter->SetInput(1, workingImageSliceInITK);
orFilter->Update();
mitk::CastToMitkImage(orFilter->GetOutput(), segmentationResult);
}
segmentationResult->GetGeometry()->SetOrigin(workingImageSlice->GetGeometry()->GetOrigin());
segmentationResult->GetGeometry()->SetIndexToWorldTransform(
workingImageSlice->GetGeometry()->GetIndexToWorldTransform());
// write to segmentation volume and hide preview image
// again, current time step is not considered
this->WriteBackSegmentationResult(m_WorkingPlane, segmentationResult, m_CurrentTimeStep);
this->m_ResultImageNode->SetVisibility(false);
this->ClearSeeds();
}
mitk::RenderingManager::GetInstance()->RequestUpdateAll();
m_ToolManager->ActivateTool(-1);
}
