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

[Dan Mueller]

Hi Hadi,

You might also be interested in the following Insight Journal article:
    http://www.insight-journal.org/browse/publication/168

This contribution provides a "look at" transform initializer, and
ObliqueSectionImageFilter which extracts a slice from a 3-D image at a
given point and direction.

Hope this helps.

Regards, Dan

2009/2/9 Hadi Rezaee <hadi_rezaie at yahoo.com>:
>
> 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
>
>
>
> _____________________________________
> Powered by www.kitware.com
>
> Visit other Kitware open-source projects at
> http://www.kitware.com/opensource/opensource.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
>