[Insight-users] getting an flat picture from dicom volume

[Hadi Rezaee]

Dear Mr. Luis,

Thanks alot for your time and reply.
i think what i am looking for is something called MPR
(multi-planar reconstruction), or resampling an image
in an arbitrary orientation as you mentioned.

I just googled the web and found one of your posts
about introducing the suitable classses for
implementing MPR:

a) ResampleImageFilter
b) VersorRigid3DTransform
c) LinearInterpolator (or BSpline is you have
   a powerful computer and what a higher quality
   resampling. 


Let me study these classes first and i will ask you my
questions :)

Thanks alot.

Regards,



--- Luis Ibanez <luis.ibanez at kitware.com> wrote:

> 
> 
> Hi Hadi,
> 
> Thanks for your clarifications and for posting the
> pictures.
> 
> I'm afraid however that your description is still
> confusing.
> 
> What you call a "flat picture":
> 
>      is it a slice ?
>      or a projection image ?
> 
> Here is what I can suggest,
> 
> Please download Slicer  (www.slicer.org) and play
> with its
> visualization functioalities (slicing and volume
> rendering)
> and let us know which one of them is closer to what
> you
> want to do.
> 
> With that information we will be able to suggest you
> way
> of doing this by using ITK and/or VTK.
> 
> 
> --
> 
> If the diagram in the second picture is what you
> want to
> get, then you could do this
> 
> in ITK by using the itk::ResampleImageFilter,
>
http://www.itk.org/Insight/Doxygen/html/classitk_1_1ResampleImageFilter.html
> 
> or in VTK by using the vtkImageReslice class
>
http://www.vtk.org/doc/nightly/html/classvtkImageReslice.html
> 
> Note that in either case you have to pass to the
> filter
> the information about the orientation of the line.
> 
> 
> 
>     Regards,
> 
> 
>         Luis
> 
> 
> -------------------
> Hadi Rezaee wrote:
> > Dear Luis,
> > 
> > Thanks alot for your reply.
> > 
> > I'm so sorry i am very newbie on this field.
> > I've checked the samples but i couldn't use the
> code
> > for my purpose, because i simply donno about
> concepts.
> > 
> > I'm also worry that i couldnt say my question
> clearly.
> > So, please check out the picture:
> > 
> >
>
http://www.freeimagehosting.net/uploads/c9c85bac6f.jpg
> > 
> > As you can see in the picture, the user specify a
> line
> > in first picture of a dicom serie.
> > So i should make an flat picture from where the
> user
> > wants to get. this flat picture should be
> generated
> > from the volume i have created for this dicom
> serie.
> > 
> > about Vtk i have a question, i used
> > vtkImagePlaneWidget object to store and display a
> > selected slice :
> > 
> >
>
http://www.freeimagehosting.net/image.php?b3495f819c.jpg
> > 
> > but i donno by which method i can save the
> > vtkImagePlaneWidget object output to jepg image.
> > 
> > Please help me.
> > 
> > Regards,
> > 
> > 
> > --- Luis Ibanez <luis.ibanez at kitware.com> wrote:
> > 
> > 
> >>Hi Hadi,
> >>
> >>You may want to use the ProjectionFilters:
> >>
> >>See the derived classes of:
> >>
> > 
> >
>
http://www.itk.org/Insight/Doxygen/html/classitk_1_1ProjectionImageFilter.html
> > 
> >>A common choice is to use the
> >>
> >>      itkMaximumProjectionImageFilter.h
> >>
> >>
> >>Please note that these filters will compute
> >>projections
> >>for lines that are parallel to the image axis.
> >>
> >>If you have a line that is not parallel then you
> >>should consider the methods used for Digitally
> >>Reconstructed
> >>Radiographs (DRR) illustrated in
> >>
> >>
> >>   InsightApplications/
> >>       IntensityBased2D3DRegistration/
> >>                 GenerateProjection.cxx
> >>
> >>
> >>That being said...
> >>If you are generating this flat image for the sole
> >>purpose of visualization, then you will have a lot
> >>more flexibility by doing Volume Rendering with
> VTK.
> >>
> >>
> >>
> >>    Regards,
> >>
> >>
> >>
> >>        Luis
> >>
> >>
> >>
> >>----------------------
> >>Hadi Rezaee wrote:
> >>
> >>>Dear all,
> >>>
> >>>i've made an volume from a series of dicom files
> >>
> >>via
> >>
> >>>itk and i can simply show the volume with vtk by
> >>
> >>using
> >>
> >>>itk::ImageToVTKImageFilter without any problem.
> >>>
> >>>What i want to do is to specify a line on the
> >>
> >>volume
> >>
> >>>and save the flat picture.
> >>>
> >>>let me clear my self. let's say i have 50 dicom
> >>>images. after i made the volume from these
> images,
> >>
> >>i
> >>
> >>>want user to give me two points on the top image
> >>
> >>that
> >>
> >>>i will show him ,for example: X= 10,10 Y=40,40 .
> >>
> >>by
> >>
> >>>specifying a line on a picture i want to get the
> >>
> >>flat
> >>
> >>>picture from volume and save it to a file. 
> >>>this flat picture should be produce from 3d
> volume
> >>
> >>by
> >>
> >>>scanning all below layers.
> >>>
> >>>please help me, it is so urgent.
> >>>
> >>>
> >>>      
> >>>_______________________________________________
> >>>Insight-users mailing list
> >>>Insight-users at itk.org
> 
=== message truncated ===


-----
Hadi Rezaee
GalaxyRoad.com