void mitk::ImageToContourFilter::Itk2DContourExtraction(const itk::Image<TPixel, VImageDimension> *sliceImage)
{
typedef itk::Image<TPixel, VImageDimension> ImageType;
typedef itk::ContourExtractor2DImageFilter<ImageType> ContourExtractor;
typedef itk::ConstantPadImageFilter<ImageType, ImageType> PadFilterType;
typename PadFilterType::Pointer padFilter = PadFilterType::New();
typename ImageType::SizeType lowerExtendRegion;
lowerExtendRegion[0] = 1;
lowerExtendRegion[1] = 1;
typename ImageType::SizeType upperExtendRegion;
upperExtendRegion[0] = 1;
upperExtendRegion[1] = 1;
/*
* We need to pad here, since the ITK contour extractor fails if the
* segmentation touches more than one image edge.
* By padding the image for one row at each edge we overcome this issue
*/
padFilter->SetInput(sliceImage);
padFilter->SetConstant(0);
padFilter->SetPadLowerBound(lowerExtendRegion);
padFilter->SetPadUpperBound(upperExtendRegion);
typename ContourExtractor::Pointer contourExtractor = ContourExtractor::New();
contourExtractor->SetInput(padFilter->GetOutput());
contourExtractor->SetContourValue(0.5);
contourExtractor->Update();
unsigned int foundPaths = contourExtractor->GetNumberOfOutputs();
vtkSmartPointer<vtkPolyData> contourSurface = vtkSmartPointer<vtkPolyData>::New();
vtkSmartPointer<vtkPoints> points = vtkSmartPointer<vtkPoints>::New();
vtkSmartPointer<vtkCellArray> polygons = vtkSmartPointer<vtkCellArray>::New();
unsigned int pointId(0);
for (unsigned int i = 0; i < foundPaths; i++)
{
const ContourPath *currentPath = contourExtractor->GetOutput(i)->GetVertexList();
vtkSmartPointer<vtkPolygon> polygon = vtkSmartPointer<vtkPolygon>::New();
polygon->GetPointIds()->SetNumberOfIds(currentPath->Size());
Point3D currentPoint;
Point3D currentWorldPoint;
for (unsigned int j = 0; j < currentPath->Size(); j++)
{
currentPoint[0] = currentPath->ElementAt(j)[0];
currentPoint[1] = currentPath->ElementAt(j)[1];
currentPoint[2] = 0;
m_SliceGeometry->IndexToWorld(currentPoint, currentWorldPoint);
points->InsertPoint(pointId, currentWorldPoint[0], currentWorldPoint[1], currentWorldPoint[2]);
polygon->GetPointIds()->SetId(j, pointId);
pointId++;
} // for2
polygons->InsertNextCell(polygon);
} // for1
contourSurface->SetPoints(points);
contourSurface->SetPolys(polygons);
contourSurface->BuildLinks();
Surface::Pointer finalSurface = this->GetOutput();
finalSurface->SetVtkPolyData(contourSurface);
}
void mitk::ImageToContourModelFilter::Itk2DContourExtraction (const itk::Image<TPixel, VImageDimension>* sliceImage)
{
typedef itk::Image<TPixel, VImageDimension> ImageType;
typedef itk::ContourExtractor2DImageFilter<ImageType> ContourExtractor;
typedef itk::PolyLineParametricPath<2> PolyLineParametricPath2D;
typedef PolyLineParametricPath2D::VertexListType ContourPath;
typedef itk::ConstantPadImageFilter<ImageType, ImageType> PadFilterType;
typename PadFilterType::Pointer padFilter = PadFilterType::New();
typename ImageType::SizeType lowerExtendRegion;
lowerExtendRegion[0] = 1;
lowerExtendRegion[1] = 1;
typename ImageType::SizeType upperExtendRegion;
upperExtendRegion[0] = 1;
upperExtendRegion[1] = 1;
/*
* We need to pad here, since the ITK contour extractor fails if the
* segmentation touches more than one image edge.
* By padding the image for one row at each edge we overcome this issue
*/
padFilter->SetInput(sliceImage);
padFilter->SetConstant(0);
padFilter->SetPadLowerBound(lowerExtendRegion);
padFilter->SetPadUpperBound(upperExtendRegion);
typename ContourExtractor::Pointer contourExtractor = ContourExtractor::New();
contourExtractor->SetInput(padFilter->GetOutput());
contourExtractor->SetContourValue(m_ContourValue);
contourExtractor->Update();
unsigned int foundPaths = contourExtractor->GetNumberOfOutputs();
this->SetNumberOfIndexedOutputs(foundPaths);
for (unsigned int i = 0; i < foundPaths; i++)
{
const ContourPath* currentPath = contourExtractor->GetOutput(i)->GetVertexList();
mitk::Point3D currentPoint;
mitk::Point3D currentWorldPoint;
mitk::ContourModel::Pointer contour = this->GetOutput(i);
if (contour.IsNull())
{
contour = mitk::ContourModel::New();
}
if (contour.IsNull())
contour = mitk::ContourModel::New();
for (unsigned int j = 0; j < currentPath->Size(); j++)
{
currentPoint[0] = currentPath->ElementAt(j)[0];
currentPoint[1] = currentPath->ElementAt(j)[1];
currentPoint[2] = 0;
m_SliceGeometry->IndexToWorld(currentPoint, currentWorldPoint);
contour->AddVertex(currentWorldPoint);
}//for2
contour->Close();
}//for1
}