Dear itk users,

 I outputed my intermediate embedding function every n iterations (I got it inside one of the levelset terms), but  I compare that with the one  that I got after evolution  (level_set->GetLabelMap()->GetPixel(idx) ), they are pretty much different.

Should itk embedding image in the last iteration be the same as the output embedding image?

Looking forward to your reply.

Best, Hui

Dear itk users,

 I outputed my intermediate embedding function every n iterations (I got it inside one of the levelset terms), but  I compare that with the one  that I got after evolution  (level_set->GetLabelMap()->GetPixel(idx) ), they are pretty much different.

Should itk embedding image in the last iteration be the same as the output embedding image?

Looking forward to your reply.

Best, Hui




_____________________________________ Powered by www.kitware.com Visit other Kitware open-source projects at http://www.kitware.com/opensource/opensource.html Kitware offers ITK Training Courses, for more information visit: http://www.kitware.com/products/protraining.html Please keep messages on-topic and check the ITK FAQ at: http://www.itk.org/Wiki/ITK_FAQ Follow this link to subscribe/unsubscribe: http://www.itk.org/mailman/listinfo/insight-users

ha Arnaud, thank you for reply.

I am using ITK4.....embedding image is the signed distance image whose zero value shows the current surface.
I runned a program that combines several itk levelset terms.
I would like to monitor the surface evolution, so I saved the signed distance image in every n  iteration.
but I expect the last signed distance image should be the same as the final segmentation, i.e. that one that I got outside ITK classes using the code from ITK example.

Best,

HUi

On Wed, Feb 1, 2012 at 2:41 PM, Arnaud Gelas <arnaudgelas-Re5JQEeQqe8AvxtiuMwx3w@public.gmane.org> wrote:

Hi Hui,

It's not clear to me what you actually mean by "embedding function" / "embedding image" ?
This function / image are embedded with respect to which class?
What kind of level-set representation are you using?

Thanks,
Arnaud

On 02/01/2012 08:32 AM, Hui Tang wrote:

Dear itk users,

 I outputed my intermediate embedding function every n iterations (I got it inside one of the levelset terms), but  I compare that with the one  that I got after evolution  (level_set->GetLabelMap()->GetPixel(idx) ), they are pretty much different.

Should itk embedding image in the last iteration be the same as the output embedding image?

Looking forward to your reply.

Best, Hui

_____________________________________ Powered by www.kitware.com Visit other Kitware open-source projects at http://www.kitware.com/opensource/opensource.html Kitware offers ITK Training Courses, for more information visit: http://www.kitware.com/products/protraining.html Please keep messages on-topic and check the ITK FAQ at: http://www.itk.org/Wiki/ITK_FAQ Follow this link to subscribe/unsubscribe: http://www.itk.org/mailman/listinfo/insight-users

Ha Arnaud
Thak you for your reply.

I used sparse Whitaker
I wrote a piece of code inside itkLevelsetXXTerm
I did it via       
      const LevelSetOutputRealType value =static_cast< LevelSetOutputRealType >( this->m_CurrentLevelSetPointer->Evaluate( it.GetIndex()  ) );
        //LevelSetPointer levelSet = this->m_LevelSetContainer->GetLevelSet( this->m_CurrentLevelSetId );
        const LevelSetOutputRealType d_val1 = this->m_Heaviside->Evaluate( -value );

I have an iterator "it" loop over an image "m_CurreneLevelSet" which is initialized.

I did save these images........

Best,

Hui

On Wed, Feb 1, 2012 at 2:57 PM, Arnaud Gelas <arnaudgelas-Re5JQEeQqe8AvxtiuMwx3w@public.gmane.org> wrote:

Hui,

With all level set discrete representations (dense, sparse Whitaker, sparse Shi or sparse Malcolm), the level set function at the last iteration (after the evolution) should be the same as the output.

Did you use an observer and a command to save these images? or did you modify some piece of code?

Do you save your image at the end of one iteration (after the evolution gets done) ?
How do you evaluate the level set function?

Arnaud

On 02/01/2012 08:50 AM, Hui Tang wrote:

ha Arnaud, thank you for reply.

I am using ITK4.....embedding image is the signed distance image whose zero value shows the current surface.
I runned a program that combines several itk levelset terms.
I would like to monitor the surface evolution, so I saved the signed distance image in every n  iteration.
but I expect the last signed distance image should be the same as the final segmentation, i.e. that one that I got outside ITK classes using the code from ITK example.

Best,

HUi

On Wed, Feb 1, 2012 at 2:41 PM, Arnaud Gelas <arnaudgelas-Re5JQEeQqe8AvxtiuMwx3w@public.gmane.org> wrote:

Hi Hui,

It's not clear to me what you actually mean by "embedding function" / "embedding image" ?
This function / image are embedded with respect to which class?
What kind of level-set representation are you using?

Thanks,
Arnaud

On 02/01/2012 08:32 AM, Hui Tang wrote:

Dear itk users,

 I outputed my intermediate embedding function every n iterations (I got it inside one of the levelset terms), but  I compare that with the one  that I got after evolution  (level_set->GetLabelMap()->GetPixel(idx) ), they are pretty much different.

Should itk embedding image in the last iteration be the same as the output embedding image?

Looking forward to your reply.

Best, Hui

_____________________________________ Powered by www.kitware.com Visit other Kitware open-source projects at http://www.kitware.com/opensource/opensource.html Kitware offers ITK Training Courses, for more information visit: http://www.kitware.com/products/protraining.html Please keep messages on-topic and check the ITK FAQ at: http://www.itk.org/Wiki/ITK_FAQ Follow this link to subscribe/unsubscribe: http://www.itk.org/mailman/listinfo/insight-users

Hello,

I am trying to use itk::VTKImageToImageFilter to to convert a VTK image to an ITK image.
The call vtkImageToImageFilter->Update(); never returns but there is no error message and program simply stops.

vtkImageData* vtkInput;


        typedef itk::Image<float, 3> ImageType;

        typedef itk::VTKImageToImageFilter<ImageType> VTKImageToImageType;

        typename VTKImageToImageType::Pointer vtkImageToImageFilter = VTKImageToImageType::New();
        QLOG_DEBUG() << "Prepared Filter";
        vtkImageToImageFilter->SetInput(vtkInput);
        QLOG_DEBUG() << "Set input to filter";
        vtkImageToImageFilter->Update();
        QLOG_DEBUG() << "Updated filter";
        typename ImageType::ConstPointer itkOutput = vtkImageToImageFilter->GetOutput();
        QLOG_DEBUG() << "Got output";
       



--
Karl Diedrich, Ph.D.
Principal Engineer
----------------------------------------
AZE Research and Development
1 Broadway, 14th floor
Cambridge, MA 02142
Tel: +1 (617) 475-1502
Cell: +1 (801) 350-1251
www.azeofamerica.com

Hi Felix,

I am using ITK 3.20.0 and itk::VTKImageToImageFilter was from the latest Slicer 4.0 source code.

I was just able to convert a VTK image to ITK with this template function without using the  itk::VTKImageToImageFilter .

IMAGEDATATYPE = short and DIMENSIONS = 3

typename itk::Image<IMAGEDATATYPE, DIMENSIONS>::Pointer convertVtkToItk(vtkImageData* vtkInput)
    {
        typedef itk::Image<IMAGEDATATYPE, DIMENSIONS> InputImageType;
        typedef itk::VTKImageImport<InputImageType> InputImageImportType;
        vtkSmartPointer<vtkImageExport> inputImageExporter = vtkImageExport::New();
        inputImageExporter->SetInput(vtkInput);
        typename InputImageImportType::Pointer inputImageImporter = InputImageImportType::New();
        ConnectPipelines(inputImageExporter.GetPointer(), inputImageImporter);
        typename InputImageType::Pointer outputImage = const_cast<InputImageType*>(inputImageImporter->GetOutput());
        outputImage->Update();
        if(outputImage.IsNull())
        {
            QLOG_WARN() << "ITK image is invalid.";
        }
        return outputImage;
    }

Cheers,