...

[FrontAlgorithms.git] / appli / CTBronchi / MoriSegmentation.cxx

diff --git a/appli/CTBronchi/MoriSegmentation.cxx b/appli/CTBronchi/MoriSegmentation.cxx
index 0b266492488f3f84727696061f8a977e061f49a6..39b9e0dcb3bcec9978554fa19e2c6bdec5a8f35d 100644 (file)
--- a/appli/CTBronchi/MoriSegmentation.cxx
+++ b/appli/CTBronchi/MoriSegmentation.cxx
@@ -17,17 +17,17 @@ typedef fpa::Image::Mori< TInputImage, TLabelImage > TFilter;
 int main( int argc, char* argv[] )
 {
   // Get arguments
-  if( argc < 17 )
+  if( argc < 16 )
   {
     std::cerr
       << "Usage: " << argv[ 0 ] << std::endl
       << "   input_image output_image output_signal" << std::endl
-      << "   init_threshold(-1024) end_threshold(0) delta(1)" << std::endl
-      << "   minimum_threshold(-850)" << std::endl
-      << "   inside_value(255) outside_value(0)" << std::endl
-      << "   signal_kernel_size(20) signal_threshold(500) signal_influence(0.5)"
+      << "   init_threshold end_threshold delta" << std::endl
+      << "   minimum_threshold" << std::endl
+      << "   inside_value outside_value" << std::endl
+      << "   signal_kernel_size signal_threshold signal_influence"
       << std::endl
-      << "   [index/point] seed_x seed_y seed_z"
+      << "   seed_x seed_y seed_z"
       << std::endl;
     return( 1 );
 
@@ -44,17 +44,9 @@ int main( int argc, char* argv[] )
   unsigned long signal_kernel_size = std::atoi( argv[ 10 ] );
   double signal_threshold = std::atof( argv[ 11 ] );
   double signal_influence = std::atof( argv[ 12 ] );
-  std::string seed_type = argv[ 13 ];
-  TInputImage::IndexType iseed;
   TInputImage::PointType pseed;
-  for( unsigned int i = 0; i < Dim; ++i )
-  {
-    if( seed_type == "index" )
-      iseed[ i ] = std::atoi( argv[ 14 + i ] );
-    else
-      pseed[ i ] = std::atof( argv[ 14 + i ] );
-
-  } // rof
+  for( int i = 0; i < Dim; ++i )
+    pseed[ i ] = std::atof( argv[ 13 + i ] );
 
   // Read image
   itk::ImageFileReader< TInputImage >::Pointer input_image_reader =
@@ -64,10 +56,7 @@ int main( int argc, char* argv[] )
   // Prepare filter
   TFilter::Pointer filter = TFilter::New( );
   filter->SetInput( input_image_reader->GetOutput( ) );
-  if( seed_type == "index" )
-    filter->SetSeed( iseed );
-  else
-    filter->SetSeed( pseed );
+  filter->SetSeed( pseed );
   filter->SetThresholds( init_threshold, end_threshold, delta );
   filter->SetMinimumThreshold( minimum_threshold );
   filter->SetInsideValue( inside_value );