X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;f=appli%2FCTBronchi%2FMoriSegmentation.cxx;fp=appli%2FCTBronchi%2FMoriSegmentation.cxx;h=39b9e0dcb3bcec9978554fa19e2c6bdec5a8f35d;hb=e977d1bb57ca61759a1fb1a1d0d8937718ef8722;hp=0b266492488f3f84727696061f8a977e061f49a6;hpb=2047276c8f1a02432fbcc7014722d460d6c1e60f;p=FrontAlgorithms.git diff --git a/appli/CTBronchi/MoriSegmentation.cxx b/appli/CTBronchi/MoriSegmentation.cxx index 0b26649..39b9e0d 100644 --- 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 );