comp.lang.idl-pvwave archive
Messages from Usenet group comp.lang.idl-pvwave, compiled by Paulo Penteado

Home » Public Forums » archive » CT LUNG VISUALISATION
Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend 
Switch to threaded view of this topic Create a new topic Submit Reply
CT LUNG VISUALISATION [message #11541] Tue, 21 April 1998 00:00 Go to next message
Fraser Hatfield is currently offline  Fraser Hatfield
Messages: 2
Registered: April 1998
Junior Member
I was wondering if someone advise me on the following.
I am trying to create a 3-D visualisation of the lungs from CT data
within IDL. To do this I use the Shade_volume and polyshade commands.
The problem I have is that the voxels are anisotropic, and so
the visualisation looks squashed on storing the slices in an array
of size (512,512,26)
The dimensions are as follows.
x=512,y=512,z=26, xres = 0.5mm, yres=0.5mm, zres=10mm.
(NB the ct slices are contiguous)
Could someone suggest the best way to manipulate the data so
that the visualisation has the correct proportions in all directions

Cheers
Fraser

Dr Fraser Niles Hatfield
Centre for Industrial and Medical Informatics (CIMI)
University of Nottingham
E-mail:fnh@cs.nott.ac.uk ; http://www.cimi.nottingham.ac.uk
Re: CT LUNG VISUALISATION [message #11606 is a reply to message #11541] Mon, 27 April 1998 00:00 Go to previous messageGo to next message
Martin Schultz is currently offline  Martin Schultz
Messages: 515
Registered: August 1997
Senior Member
Fraser Hatfield wrote:
>
> I was wondering if someone advise me on the following.
> I am trying to create a 3-D visualisation of the lungs from CT data
> within IDL. To do this I use the Shade_volume and polyshade commands.
> The problem I have is that the voxels are anisotropic, and so
> the visualisation looks squashed on storing the slices in an array
> of size (512,512,26)
> The dimensions are as follows.
> x=512,y=512,z=26, xres = 0.5mm, yres=0.5mm, zres=10mm.
> (NB the ct slices are contiguous)
> Could someone suggest the best way to manipulate the data so
> that the visualisation has the correct proportions in all directions
>

Although I have never seen my loungs in IDL, you may want to try congrid
in order to interpolate your data in the z dimension:

newdata = congrid(data,512,512,520,/interpolate)


Hope this helps,
Martin.

------------------------------------------------------------ -------
Dr. Martin Schultz
Department for Earth&Planetary Sciences, Harvard University
186 Pierce Hall, 29 Oxford St., Cambridge, MA-02138, USA

phone: (617)-496-8318
fax : (617)-495-4551

e-mail: mgs@io.harvard.edu
IDL-homepage: http://www-as.harvard.edu/people/staff/mgs/idl/
------------------------------------------------------------ -------
Re: CT LUNG VISUALISATION [message #11727 is a reply to message #11541] Thu, 30 April 1998 00:00 Go to previous message
David Foster is currently offline  David Foster
Messages: 341
Registered: January 1996
Senior Member
Martin Schultz wrote:
>
> Fraser Hatfield wrote:
>>
>> I was wondering if someone advise me on the following.
>> I am trying to create a 3-D visualisation of the lungs from CT data
>> within IDL. To do this I use the Shade_volume and polyshade commands.
>> The problem I have is that the voxels are anisotropic, and so
>> the visualisation looks squashed on storing the slices in an array
>> of size (512,512,26)
>> The dimensions are as follows.
>> x=512,y=512,z=26, xres = 0.5mm, yres=0.5mm, zres=10mm.
>> (NB the ct slices are contiguous)
>> Could someone suggest the best way to manipulate the data so
>> that the visualisation has the correct proportions in all directions
>>
>
> Although I have never seen my loungs in IDL, you may want to try congrid
> in order to interpolate your data in the z dimension:
>
> newdata = congrid(data,512,512,520,/interpolate)
>

Since CT data is usually 16-bit integer data, this would create a
270MB array! I think a better approach would be to use IDL's system
variables !X.S, !Y.S, !Z.S, and !P.T together with the /T3D keyword
in the POLYSHADE() call, together with the Z buffer.

Here is a code snippet from one of our programs. Email me if you
have questions.

!P.T = d.isosurf_pt
!X.S = d.isosurf_x_s
!Y.S = d.isosurf_y_s
!Z.S = d.isosurf_z_s

set_plot,'Z'
erase
b = polyshade(verts, polys, /t3d) ; This brain is upside down

Dave

--

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ ~~~~~~
David S. Foster Univ. of California, San Diego
Programmer/Analyst Brain Image Analysis Laboratory
foster@bial1.ucsd.edu Department of Psychiatry
(619) 622-5892 8950 Via La Jolla Drive, Suite 2240
La Jolla, CA 92037
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ ~~~~~~
  Switch to threaded view of this topic Create a new topic Submit Reply
Previous Topic: Re: IDL5 License Manager on DEC UNIX 4.0d
Next Topic: Re: Help : I cannot poly_fit to 9 degree!!!

-=] Back to Top [=-
[ Syndicate this forum (XML) ] [ RSS ] [ PDF ]

Current Time: Fri Oct 10 05:01:22 PDT 2025

Total time taken to generate the page: 7.04266 seconds