[Insight-users] the sigmoid image is entirely white for GeodesicActiveContourLevelSet

[Goo]

Hi All:

This question is so confused to me.
The output of sigmoid image is shown white entirely but gradient image fine.

The following are my programs:
if the ExternalImageType used as float, the image of result will be white
entirely in smooth,gradient and sigmoid.

Please give me some guidance.
Regards.

//----------------------------------CODE-------------------------------------------

    typedef itk::Image< signed short, 3 >  InternalImageType;
    typedef itk::Image< signed short, 3 >  ExternalImageType;

    typedef itk::RescaleIntensityImageFilter<ExternalImageType,
InternalImageType> InputCastType;
    typedef itk::RescaleIntensityImageFilter<InternalImageType,
ExternalImageType> OutputCastType;

    OutputCastType::Pointer            m_SmoothObserver;
    OutputCastType::Pointer            m_GradientObserver;
    OutputCastType::Pointer            m_SigmoidObserver;

        typedef itk::ResampleImageFilter<
                            InternalImageType , InternalImageType >
FilterType;
        FilterType::Pointer SubSample = FilterType::New();

        typedef itk::AffineTransform< double, 3 >  TransformType;
        TransformType::Pointer transform = TransformType::New();

        typedef itk::NearestNeighborInterpolateImageFunction<
                           InternalImageType, double >  InterpolatorType;
        InterpolatorType::Pointer interpolator = InterpolatorType::New();

        SubSample->SetInterpolator( interpolator );
        SubSample->SetDefaultPixelValue( 50 );

        double spacing[ 3 ];
        spacing[0] = InputCaster->GetOutput()->GetSpacing()[0];
        spacing[1] = InputCaster->GetOutput()->GetSpacing()[1];
        spacing[2] = InputCaster->GetOutput()->GetSpacing()[2];
        SubSample->SetOutputSpacing( spacing );

        double origin [3];
        origin[0] = 80*spacing[0];
        origin[1] = 80*spacing[1];
        origin[2] = 15*spacing[2];
        SubSample->SetOutputOrigin(origin);

        InternalImageType::SizeType   size;

        size[0] = 100;
        size[1] = 100;
        size[2] = 5;

        SubSample->SetSize( size );
        SubSample->SetInput( InputCaster->GetOutput() );
        SubSample->SetTransform( transform );
        transform->SetIdentity();

        SubSample->Update();


        typedef itk::CurvatureAnisotropicDiffusionImageFilter
                    < InternalImageType,InternalImageType >
SmoothingFilterType;

        typedef itk::GradientMagnitudeRecursiveGaussianImageFilter
                    < InternalImageType,InternalImageType >
GradientFilterType;

        typedef itk::SigmoidImageFilter
                    < InternalImageType,InternalImageType >
SigmoidFilterType;

        SmoothingFilterType::Pointer    m_smoothing =
SmoothingFilterType::New();
        GradientFilterType::Pointer        m_gradientMagnitude =
GradientFilterType::New();
        SigmoidFilterType::Pointer        m_sigmoid =
SigmoidFilterType::New();

        m_smoothing->SetNumberOfIterations(  m_SmoothIter );
        m_smoothing->SetConductanceParameter( m_SmoothConductance );
        m_smoothing->SetTimeStep( m_TimeStep );

        m_gradientMagnitude->SetSigma(  m_sigma  );

        m_sigmoid->SetOutputMinimum(  m_SigmoidOutMin  );
        m_sigmoid->SetOutputMaximum(  m_SigmoidOutMax  );
        m_sigmoid->SetAlpha( m_alpha );
        m_sigmoid->SetBeta(  m_beta  );

        m_smoothing->SetInput(     SubSample->GetOutput() );
        m_gradientMagnitude->SetInput( m_smoothing->GetOutput() );
        m_sigmoid->SetInput( m_gradientMagnitude->GetOutput() );

        m_SmoothObserver->SetInput(m_smoothing->GetOutput());
        m_GradientObserver->SetInput(m_gradientMagnitude->GetOutput());
        m_SigmoidObserver->SetInput(m_sigmoid->GetOutput());

        m_SmoothObserver->Update();
        m_GradientObserver->Update();
        m_SigmoidObserver->Update();
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://public.kitware.com/pipermail/insight-users/attachments/20070327/2efaed9f/attachment.htm