[vtkusers] Create an approximate cross-section of a RT ROI

Thu May 6 11:24:28 EDT 2010

Hello Xiaofeng,

vtkPolyDataTransformation is new to me.  I'll give it a try.

I was trying to position and orient the vtkActor to get it in the correct spot.

Thanks

                                -- Scott

From: Xiaofeng Z [mailto:xf10036 at hotmail.com] 
Sent: Thursday, May 06, 2010 9:32 AM
To: Scott Johnson; lars-friedrich at gmx.net; VTK
Subject: RE: [vtkusers] Create an approximate cross-section of a RT ROI

Scott,

Sorry for not responding earlier.  I was out on a trip.

Have you resolved your issue yet? Did you get a vtkPolyData from cutting the sphere?  It seems your cutting plane is not orientated the same way as your reslicer.  What I did is first cut the structure with the orientaiton and the origin corresponding to the reslicer, then use the vtkPodyDataTransformation to transform the plydata from the patient coordiation to the world coordination of vtkImageViewer.

Xiaofeng

________________________________

Subject: RE: [vtkusers] Create an approximate cross-section of a RT ROI
Date: Tue, 4 May 2010 14:35:09 -0500
From: Scott.Johnson at neuwave.com
To: xf10036 at hotmail.com; lars-friedrich at gmx.net; vtkusers at vtk.org

Hello Xiaofeng,

I'll attach a picture of what I'm seeing.  The image shows a coronal, sagittal, axial, and 3D MIP rendering.

The plane connected to the cutter has the same normal as the vtkImageReslice that I'm using to extract the images.  The vtkImageReslice is the input to the vtkImageViewer2.  I've added the vtkActor showing the cut ROI to the vtkRenderer managed by the vtkImageViewer2.  The little green lines on the coronal and sagittal images should be the resulting circles from cutting the sphere and they should appear over the image at the little bright spots.  They do show up when I get to the appropriate slice.

I've tried changing the orientation of the vtkActor to match the ResliceAxes and the inverse of the ResliceAxes and get the same result.  This works on the axial slices with no extra processing, but not for the coronal and sagittal images.

I'm sure I'm missing something fundamental, I just haven't been able to figure out what it is.

Thanks

                                -- Scott

From: vtkusers-bounces at vtk.org [mailto:vtkusers-bounces at vtk.org] On Behalf Of Xiaofeng Z
Sent: Tuesday, May 04, 2010 10:20 AM
To: Scott Johnson; lars-friedrich at gmx.net; VTK
Subject: Re: [vtkusers] Create an approximate cross-section of a RT ROI

Scott,

What I did is to transform the contour using the inverse of the reslice axies as the transform matrix, works pretty well in transversal!  Should work in the other orientation too.

Xiaofeng

> Subject: RE: [vtkusers] Create an approximate cross-section of a RT ROI
> Date: Tue, 4 May 2010 09:52:55 -0500
> From: Scott.Johnson at neuwave.com
> To: lars-friedrich at gmx.net; vtkusers at vtk.org; xf10036 at hotmail.com
> 
> Thanks Lars.
> 
> I got around that problem, but now I've moved to using vtkImageReslice for creating the other slice orientations and the ROI slice is floating out in space away from the images. It works fine for the axial but sagittal, coronal, and oblique slices of the ROI are in the wrong location.
> 
> I've been playing with origins, centers, and positions of the vtkActor connected to the vtkCutter without any luck.
> 
> -- Scott
> 
> -----Original Message-----
> From: Lars Friedrich Lars [mailto:lars-friedrich at gmx.net] 
> Sent: Tuesday, May 04, 2010 9:36 AM
> To: Scott Johnson; vtkusers at vtk.org; xf10036 at hotmail.com
> Subject: Re: [vtkusers] Create an approximate cross-section of a RT ROI
> 
> Scott,
> 
> I had a similar problem with 'free planes' at the beginning, but I usually define the cut plane 'a bit' (0.001 mm or so) in front of the image slice (towards view position).
> 
> regards,
> 
> lars
> 
> -------- Original-Nachricht --------
> > Datum: Tue, 4 May 2010 09:12:33 -0500
> > Von: "Scott Johnson" <Scott.Johnson at neuwave.com>
> > An: "Xiaofeng Z" <xf10036 at hotmail.com>, "VTK" <vtkusers at vtk.org>
> > Betreff: Re: [vtkusers] Create an approximate cross-section of a RT ROI
> 
> > Hello Xiaofeng,
> > 
> > 
> > 
> > I've been working on something similar, but am not quite there yet. I
> > can get you started.
> > 
> > 
> > 
> > The general idea is that you will need to create a polygonal shape from
> > the contours in the RT structure set, then create a vtkPlane which
> > represents the plane you wish to cut. You can then create a pipeline
> > which looks something like:
> > 
> > 
> > 
> > vtkCutter -> vtkPolyDataMapper -> vtkActor
> > 
> > 
> > 
> > Set the properties of the actor to have EdgeVisibilityOn and an
> > EdgeColor or the cut shape will not be visible. The actor will show a
> > cross section of the ROI intersected with the plane.
> > 
> > 
> > 
> > The issue I'm having is properly displaying the vtkActor with the
> > images. I've extended beyond axial sagittal and coronal and I seem to
> > be missing something in the coordinate conversions.
> > 
> > 
> > 
> > Good luck.
> > 
> > 
> > 
> > -- Scott
> > 
> > 
> > 
> > From: vtkusers-bounces at vtk.org [mailto:vtkusers-bounces at vtk.org] On
> > Behalf Of Xiaofeng Z
> > Sent: Monday, May 03, 2010 7:12 PM
> > To: VTK
> > Subject: [vtkusers] Create an approximate cross-section of a RT ROI
> > 
> > 
> > 
> > Hi vtkUsers,
> > 
> > I have a set of contours in axial orientation representing a ROI. I
> > need to sketch the ROI in coronal and sagittal plane, by connection the
> > points made from the intersection of the contours and the cutting plane.
> > Is there a function to do that in vtk or in ITK?
> > 
> > I know the result is not "precise" in a sense that the contours make
> > this way is not always consistent between themselves. But it seemed to
> > be the only viable way to sketch a ROI in a multi-planar rendering in
> > real time.
> > 
> > Many thanks!
> > Xiaofeng
> > 
> > 
> > 
> > 
> > ________________________________
> > 
> > The New Busy is not the old busy. Search, chat and e-mail from your
> > inbox. Get started.
> > <http://www.windowslive.com/campaign/thenewbusy?ocid=PID28326::T:WLMTAGL
> > :ON:WL:en-US:WM_HMP:042010_3> 
> > 
> 
> -- 
> GRATIS für alle GMX-Mitglieder: Die maxdome Movie-FLAT!
> Jetzt freischalten unter http://portal.gmx.net/de/go/maxdome01

________________________________

Hotmail is redefining busy with tools for the New Busy. Get more from your inbox. See how. <http://www.windowslive.com/campaign/thenewbusy?ocid=PID28326::T:WLMTAGL:ON:WL:en-US:WM_HMP:042010_2> 

________________________________

Hotmail is redefining busy with tools for the New Busy. Get more from your inbox. See how. <http://www.windowslive.com/campaign/thenewbusy?ocid=PID28326::T:WLMTAGL:ON:WL:en-US:WM_HMP:042010_2> 

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.vtk.org/pipermail/vtkusers/attachments/20100506/7e3dc99c/attachment.htm>