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

Scott Johnson Scott.Johnson at neuwave.com
Tue May 4 16:21:37 EDT 2010


Because I am going to extend the functionality to oblique angles of the plane.

My thought is that if I solve the problems with the coronal and sagittal slices the oblique slices will work as well.

		-- Scott

-----Original Message-----
From: Jothy [mailto:jothybasu at gmail.com] 
Sent: Tuesday, May 04, 2010 3:17 PM
To: Scott Johnson
Cc: Xiaofeng Z; lars-friedrich at gmx.net; VTK
Subject: Re: [vtkusers] Create an approximate cross-section of a RT ROI

If you are using vtkImageViewer2, why are you using vtkImageReslice
instead of setSlice and SetOrientation?

Jothy

On Tue, May 4, 2010 at 8:35 PM, Scott Johnson <Scott.Johnson at neuwave.com> wrote:
> 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.
>
> _______________________________________________
> 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 VTK FAQ at:
> http://www.vtk.org/Wiki/VTK_FAQ
>
> Follow this link to subscribe/unsubscribe:
> http://www.vtk.org/mailman/listinfo/vtkusers
>
>



More information about the vtkusers mailing list