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 );
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 =
// 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 );