[Rtk-users] forward and back projection - MITK

[Tobias Stein]

Thanks for the information.

 

To keep it simple I want to project forward a single slice which I load as a DICOM file.

Now I want to compute the sinogram of that slice and write it in a new file.

I got a little confused about the ConstantImageSource which is used in the test class. Do I need it to compute a sinogram or is just a replacement of a itk::Image?

 

Here is my code so far:

 

                const unsigned int Dimension = 2;

       typedef float PixelType;

       typedef itk::Image< PixelType, Dimension >      ImageType;

                …

                (reading image by itk::ImageFileReader which works)

                …

       ImageType::Pointer inputImage = reader->GetOutput();

 

       const unsigned int NumberOfProjectionImages = 1;

 

       ImageType::Pointer outputImage = ImageType::New();

 

       typedef  rtk::JosephForwardProjectionImageFilter<ImageType, ImageType> ForwardType;

       ForwardType::Pointer forward = ForwardType::New();

       forward->SetInput(1, inputImage);

       forward->SetInput(0, outputImage);

       forward->Update();

 

                …

                (writing output image into a file)

                …

 

The problem is that I get some errors when I instantiate the forward-Filter.

The errors are like that one:

error C2784: "vnl_matrix<T> operator *(const vnl_diag_matrix<T> &,const vnl_matrix<T> &)": template-Argument für "const vnl_diag_matrix<T> &" konnte nicht von "vnl_matrix_fixed<T,4,4>" hergeleitet werden.

 

Maybe there is something wrong with the PixelType that i used?

 

If the forward projection works, which filter should i use to achieve a filtered back projection without aliasing?

 

Many thanks in advance.

Tobias

 

Von: simon.rit at gmail.com [mailto:simon.rit at gmail.com] Im Auftrag von Simon Rit
Gesendet: Mittwoch, 4. Mai 2016 09:02
An: Tobias Stein <tstein at stud.hs-heilbronn.de>
Cc: rtk-users at public.kitware.com
Betreff: Re: [Rtk-users] forward and back projection - MITK

 

Dear Tobias,

The forward projection has 2 inputs:

- input 0: the stack of projections in which you wish to forward project,

- input 1: the volume you wish to forward project.

For the backprojection, it's exactly the same:

- input 0: the volume in which you wish to backproject,

- input 1: the stack of projections you wish to backproject.

Be aware that JosephBackProjectionImageFilter is the transpose of the forward projection and will have some aliasing (see, e.g., De Man and Basu <http://iopscience.iop.org/article/10.1088/0031-9155/49/11/024/meta;jsessionid=87B598ABFDC2AA07D1DED2DEE607F0E7.c2.iopscience.cld.iop.org>  to see what I mean).

I have personally never used MITK but don't hesitate to share your experience on the mailing list.

Simon

 

On Wed, May 4, 2016 at 12:59 AM, Tobias Stein <tstein at stud.hs-heilbronn.de <mailto:tstein at stud.hs-heilbronn.de> > wrote:

Hi all,

 

i want to use the forward and backward projection to reduce metal artefacts in ct images. I’ve seen so far that  I may use the JosephForwardProjectionImageFilter to perform the forward projection. I’ve also seen the test for this class but I don’t get it, where should i put my 2D slice as input to execute the transformation. About the back transformation with the JosephBackProjectionImageFilter I also need more information how I can use it to transform a sinogram back to a ct slice. There is a test at the documentation, but the link is missing there.

 

I’ve got another independent question.

Are some of you familiar with MITK and know how to get the superbuild up and running with RTK? I am writing a MITK-Plugin and want to reduce the metal artifacts before a segmentation. So if some of you have experience with the combination of RTK and MITK it would be nice if you will share it ;)

 

Best regards,

Tobias

 


_______________________________________________
Rtk-users mailing list
Rtk-users at public.kitware.com <mailto:Rtk-users at public.kitware.com> 
http://public.kitware.com/mailman/listinfo/rtk-users

 

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://public.kitware.com/pipermail/rtk-users/attachments/20160513/9f1e3827/attachment-0010.html>