int
main(int argc, char * argv[])
{
if (argc < 5)
{
std::cerr << "Missing parameters." << std::endl;
std::cerr << "Usage: " << argv[0] << " inputImageFile outputImageFile threshold level" << std::endl;
return EXIT_FAILURE;
}
const auto input = itk::ReadImage<InputImageType>(argv[1]);
gradientMagnitudeImageFilter->SetInput(input);
gradientMagnitudeImageFilter->Update();
float threshold = std::stod(argv[3]);
float level = std::stod(argv[4]);
watershed->SetThreshold(threshold);
watershed->SetLevel(level);
watershed->SetInput(gradientMagnitudeImageFilter->GetOutput());
watershed->Update();
colormapImageFilter->SetInput(watershed->GetOutput());
colormapImageFilter->Update();
return EXIT_SUCCESS;
}