- if( visual_debug )
- {
- // Configure observer
- fpa::VTK::Image3DObserver< fpa::Image::Dijkstra< TDistanceMap, TScalar >, vtkRenderWindow >::Pointer obs =
- fpa::VTK::Image3DObserver< fpa::Image::Dijkstra< TDistanceMap, TScalar >, vtkRenderWindow >::New( );
- obs->SetImage( distanceMap->GetOutput( ), view.GetWindow( ) );
- paths->AddObserver( itk::AnyEvent( ), obs );
- paths->ThrowEventsOn( );
- }
- else
- paths->ThrowEventsOff( );
- start = std::clock( );
- paths->Update( );
- end = std::clock( );
- seconds = double( end - start ) / double( CLOCKS_PER_SEC );
- std::cout << "Paths extraction time = " << seconds << std::endl;
+ TImageWriter::Pointer segmentation_writer =
+ TImageWriter::New( );
+ segmentation_writer->SetInput( segmentor->GetOutput( ) );
+ segmentation_writer->SetFileName( output_segmentation_fn );
+ segmentation_writer->Update( );