[Insight-users] Rigid3DPerspectiveTransform

[Luis Ibanez]

Hi Peter,

You may find useful the instructions on the following email:
http://www.itk.org/pipermail/insight-users/2006-January/016377.html

regarding the use of the DigitallyReconstructedRadiograph1.cxx example


For a synthetic cube image, you can use the following command line:


    DigitallyReconstructedRadiograph1 -t -90 -90 -90 -o CubeDRR.img


For a real image you can use the following:


    DigitallyReconstructedRadiograph1 -t -90 -109 -90
                          -o brain.img  brainweb165a10f17.mhd


(in a single line)

The brainweb image can be downloaded from the Data page of www.itk.org:

              http://www.itk.org/HTML/Data.htm
              ftp://public.kitware.com/pub/itk/Data/BrainWeb



    Regards,



       Luis


-----------------------------
Boettcher, Dr. Peter wrote:
> Dear Bill
> 
>  
> 
> Yes I already used the 
> Insight/Examples/DigitallyReconstructedRadiograph1.cxx example. 
> Unfortunately I have problems defining the proper transformation. 
> Instead of transforming the object (e.g. bone volume) the virtual C-arm 
> is rotated as well as translated around the image volume. When I define 
> only a rotation or a translation every thing works fine. However when I 
> apply both rotation as well as translation the result is not what I 
> expected.
> 
>  
> 
> Using a parallel projection for the generation of DRR would be of less 
> value for me as I would like to match those DRRs to a real fluoroscopic 
> image. In fact I must admit that I miss a function like 
> itkRayCaster-SetViewAngle(). As the itk::ResampleImageFilter defining 
> the dimensions as well as the image resolution of the DRR is generated 
> after the itk::RayCaster has done its work suggests, that the projection 
> angle at which the RayCaster calculates the DRR is fixed.
> 
>  
> 
> However if anyone could explain me the proper usage of 
> DigitallyReconstructedRadiograph1.cxx, maybe that will solve my problem. 
> Anyway I thought using itk::Rigid3DPerspectiveTransform would be more 
> efficient as it has been especially designed to generate DRRs ;-)
> 
>  
> 
>  
> 
>  
> 
> Regards, Peter.
> 
>  
>  
> 
> ------------------------------------------------------------------------
> *Von:* Bill Lorensen [mailto:bill.lorensen at gmail.com]
> *Gesendet:* Mo 21.04.2008 18:42
> *An:* Boettcher, Dr. Peter
> *Cc:* insight-users @ itk. org
> *Betreff:* Re: [Insight-users] Rigid3DPerspectiveTransform
> 
> Peter,
> 
> I the past at GE, we generated DRR's using the ExpNegativeImageFilter
> and AccummulateImageFilter. These resulted in parallel projections
> along each of the major planes, not Perspective projections. I do
> recall that we also successfully generated arbitrary DRR's with
> Insight/Examples/DigitallyReconstructedRadiograph1.cxx. I have no
> experience with similar code in InsightApplications. Which code were
> you using in InsightApplications?
> 
> Bill
> 
> On Mon, Apr 21, 2008 at 11:13 AM, Boettcher, Dr. Peter
> <Boettcher at kleintierklinik.uni-leipzig.de> wrote:
>  >
>  > I try to use Rigid3DPerspectiveTransform for the generation of digitally
>  > reconstructed radiographs (DRR). The example within the 
> ITKApplication uses
>  > the RayCaster instead of Rigid3DPerspectiveTransform. Searching the 
> mailing
>  > list as well as the web I did not find any example on how to use that
>  > filter. The testing script does only transform points and not an image
>  > volume.
>  > I suppose I will have to declare Rigid3DPerspectiveTransform as my
>  > transformer. In additon to that I will need an image resample filter 
> and an
>  > interpolator. However pluggin those components together does not work.
>  >
>  > ImageResampleFilter->SetTransform(Rigid3DPerspectiveTransform) does ot
>  > compile.
>  >
>  > Regards, Peter.
>  >
>  >
>  > _______________________________________________
>  > Insight-users mailing list
>  > Insight-users at itk.org
>  > http://www.itk.org/mailman/listinfo/insight-users
>  >
>  >
> 
> 
> ------------------------------------------------------------------------
> 
> _______________________________________________
> Insight-users mailing list
> Insight-users at itk.org
> http://www.itk.org/mailman/listinfo/insight-users