void mitk::CLUtil::itkDilateBinary(TImageType * sourceImage, mitk::Image::Pointer& resultImage, int factor, MorphologicalDimensions d)
{
typedef TImageType ImageType;
typedef itk::BinaryBallStructuringElement<typename TImageType::PixelType, 3> BallType;
typedef typename itk::BinaryDilateImageFilter<TImageType, TImageType, BallType> BallDilateFilterType;
BallType strElem;
itkFitStructuringElement(strElem,d,factor);
typename BallDilateFilterType::Pointer erodeFilter = BallDilateFilterType::New();
erodeFilter->SetKernel(strElem);
erodeFilter->SetInput(sourceImage);
erodeFilter->SetDilateValue(1);
erodeFilter->Update();
mitk::CastToMitkImage(erodeFilter->GetOutput(), resultImage);
}
void mitk::MorphologicalOperations::itkDilate(
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::BinaryDilateImageFilter<ImageType, ImageType, BallType> BallDilateFilterType;
typedef typename itk::BinaryDilateImageFilter<ImageType, ImageType, CrossType> CrossDilateFilterType;
if (structuralElementFlags & (Ball_Axial | Ball_Coronal | Ball_Sagital))
{
BallType ball = CreateStructuringElement<BallType>(structuralElementFlags, factor);
typename BallDilateFilterType::Pointer dilateFilter = BallDilateFilterType::New();
dilateFilter->SetKernel(ball);
dilateFilter->SetInput(sourceImage);
dilateFilter->SetDilateValue(1);
dilateFilter->UpdateLargestPossibleRegion();
mitk::CastToMitkImage(dilateFilter->GetOutput(), resultImage);
}
else
{
CrossType cross = CreateStructuringElement<CrossType>(structuralElementFlags, factor);
typename CrossDilateFilterType::Pointer dilateFilter = CrossDilateFilterType::New();
dilateFilter->SetKernel(cross);
dilateFilter->SetInput(sourceImage);
dilateFilter->SetDilateValue(1);
dilateFilter->UpdateLargestPossibleRegion();
mitk::CastToMitkImage(dilateFilter->GetOutput(), resultImage);
}
}
