- /* TODO
- TDistanceMapWriter::Pointer distancemap_writer =
- TDistanceMapWriter::New( );
- distancemap_writer->SetInput( distanceMap->GetOutput( ) );
- distancemap_writer->SetFileName( output_distancemap_fn );
- distancemap_writer->Update( );
-
- TImageWriter::Pointer segmentation_writer =
- TImageWriter::New( );
- segmentation_writer->SetInput( segmentor->GetOutput( ) );
- segmentation_writer->SetFileName( output_segmentation_fn );
- segmentation_writer->Update( );
- */
-
- // Show result
- /*
- TVTKImage::Pointer output_image_vtk = TVTKImage::New( );
- output_image_vtk->SetInput( segmentor->GetOutput( ) );
- output_image_vtk->Update( );
-
- vtkSmartPointer< vtkImageMarchingCubes > mc =
- vtkSmartPointer< vtkImageMarchingCubes >::New( );
- mc->SetInputData( output_image_vtk->GetOutput( ) );
- mc->SetValue(
- 0,
- double( segmentor->GetInsideValue( ) ) * double( 0.95 )
- );
- mc->Update( );
-
- // Let some interaction and close program
- view.AddPolyData( mc->GetOutput( ), 0.1, 0.6, 0.8, 0.5 );
- view.Start( );
+ itk::ImageFileWriter< TImage >::Pointer segmentation_writer =
+ itk::ImageFileWriter< TImage >::New( );
+ segmentation_writer->SetInput( segmentor->GetOutput( ) );
+ segmentation_writer->SetFileName( output_segmentation_fn );