void mitk::CLUtil::itkErodeBinary(TImageType * sourceImage, mitk::Image::Pointer& resultImage, int factor, MorphologicalDimensions d)
{
typedef TImageType ImageType;
typedef itk::BinaryBallStructuringElement<typename TImageType::PixelType, 3> BallType;
typedef typename itk::BinaryErodeImageFilter<TImageType, TImageType, BallType> BallErodeFilterType;
BallType strElem;
itkFitStructuringElement(strElem,d,factor);
typename BallErodeFilterType::Pointer erodeFilter = BallErodeFilterType::New();
erodeFilter->SetKernel(strElem);
erodeFilter->SetInput(sourceImage);
erodeFilter->SetErodeValue(1);
// erodeFilter->UpdateLargestPossibleRegion();
erodeFilter->Update();
mitk::CastToMitkImage(erodeFilter->GetOutput(), resultImage);
}
void mitk::MorphologicalOperations::itkErode(
itk::Image<TPixel, VDimension> *sourceImage,
mitk::Image::Pointer &resultImage,
int factor,
mitk::MorphologicalOperations::StructuralElementType structuralElementFlags)
{
typedef itk::Image<TPixel, VDimension> ImageType;
typedef itk::BinaryBallStructuringElement<TPixel, VDimension> BallType;
typedef itk::BinaryCrossStructuringElement<TPixel, VDimension> CrossType;
typedef typename itk::BinaryErodeImageFilter<ImageType, ImageType, BallType> BallErodeFilterType;
typedef typename itk::BinaryErodeImageFilter<ImageType, ImageType, CrossType> CrossErodeFilterType;
if (structuralElementFlags & (Ball_Axial | Ball_Coronal | Ball_Sagital))
{
BallType ball = CreateStructuringElement<BallType>(structuralElementFlags, factor);
typename BallErodeFilterType::Pointer erodeFilter = BallErodeFilterType::New();
erodeFilter->SetKernel(ball);
erodeFilter->SetInput(sourceImage);
erodeFilter->SetErodeValue(1);
erodeFilter->UpdateLargestPossibleRegion();
mitk::CastToMitkImage(erodeFilter->GetOutput(), resultImage);
}
else
{
CrossType cross = CreateStructuringElement<CrossType>(structuralElementFlags, factor);
typename CrossErodeFilterType::Pointer erodeFilter = CrossErodeFilterType::New();
erodeFilter->SetKernel(cross);
erodeFilter->SetInput(sourceImage);
erodeFilter->SetErodeValue(1);
erodeFilter->UpdateLargestPossibleRegion();
mitk::CastToMitkImage(erodeFilter->GetOutput(), resultImage);
}
}
