+ outputGradientFilter = gradientFilter->GetOutput();
+ }
+ else {
+ typedef itk::GradientMagnitudeRecursiveGaussianImageFilter<InputImageType, OutputImageType> GradientMagnitudeImageFilterType;
+ typename GradientMagnitudeImageFilterType::Pointer gradientFilter=GradientMagnitudeImageFilterType::New();
+ gradientFilter->SetInput(input);
+ gradientFilter->Update();
+ //std::cout<<"sigma value="<<gradientFilter->GetSigma()<<std::endl;
+ outputGradientFilter = gradientFilter->GetOutput();
+ }