<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">Hello Andrej,<div><br></div><div>Well, maybe someone else can give their opinion, too, but personally I would just use a standard vtkRenderWindow, and add all three actors to the vtkRenderer that you add to that render window. </div><div><br></div><div>Each of the images can be shown with a vtkImageActor, and those can be placed anywhere in 3D space with SetPosition(), and then you can just use a regular vtkActor for the contour filter output and place it "in front" of the image actors.</div><div><br></div><div>Here's a piece of Python code that does this, if you want to have a look. (Sorry the beginning is more complicated than it should be, I was already playing with VTK/Examples/ImageProcessing/Python/ImageSlicing.py and this was an easy way to get two grayscale images to superimpose...) It takes two images, gives one an orange and one a yellow color map with transparency, then overlays them. (I positioned them separately in space, but I think you can put these on the same Z plane.) Then, I place the contour output "over the top" of the images. (If you change the interactor style to TrackballCamera, you can roll them around and see the separation, which doesn't have to be this extreme.)</div><div><br></div><div>Again, maybe others have a better way to do this, but it works for keeping the spatial order as you want it, and they all move together in the window.</div><div><br></div><div>Good luck,</div><div>-Eric</div><div><br></div><div>===================</div><div><br></div><div><div>#!/usr/bin/env python</div><div><br></div><div>import vtk</div><div><br></div><div>VTK_DATA_ROOT = '/Users/emonson/Programming/VTK_cvs/VTKData'</div><div><br></div><div># Start by loading some data.</div><div>reader = vtk.vtkImageReader2()</div><div>reader.SetFilePrefix(VTK_DATA_ROOT + "/Data/headsq/quarter")</div><div>reader.SetDataExtent(0, 63, 0, 63, 1, 93)</div><div>reader.SetDataSpacing(1.0, 1.0, 1.0)</div><div>reader.SetDataOrigin(0.0, 0.0, 0.0)</div><div>reader.SetDataScalarTypeToUnsignedShort()</div><div>reader.UpdateWholeExtent()</div><div><br></div><div># Calculate the center of the volume</div><div>reader.GetOutput().UpdateInformation()</div><div>(xMin, xMax, yMin, yMax, zMin, zMax) = reader.GetOutput().GetWholeExtent()</div><div>(xSpacing, ySpacing, zSpacing) = reader.GetOutput().GetSpacing()</div><div>(x0, y0, z0) = reader.GetOutput().GetOrigin()</div><div><br></div><div>center = [x0 + xSpacing * 0.5 * (xMin + xMax),</div><div> y0 + ySpacing * 0.5 * (yMin + yMax),</div><div> z0 + zSpacing * 0.5 * (zMin + zMax)]</div><div><br></div><div>renderer = vtk.vtkRenderer()</div><div>renderer.SetBackground(1,1,1)</div><div><br></div><div>resliceList = []</div><div>tableList = []</div><div>colorList = []</div><div>actorList = []</div><div><br></div><div>imgOffset = 70</div><div><br></div><div>for ii in range(2):</div><div><span class="Apple-tab-span" style="white-space:pre"> </span></div><div><span class="Apple-tab-span" style="white-space:pre"> </span>zpos = z0 + zSpacing*(zMin+ii*imgOffset)</div><div><span class="Apple-tab-span" style="white-space:pre"> </span>axial = vtk.vtkMatrix4x4()</div><div><span class="Apple-tab-span" style="white-space:pre"> </span>axial.DeepCopy((1, 0, 0, center[0],</div><div><span class="Apple-tab-span" style="white-space:pre"> </span>0, 1, 0, center[1],</div><div><span class="Apple-tab-span" style="white-space:pre"> </span>0, 0, 1, zpos,</div><div><span class="Apple-tab-span" style="white-space:pre"> </span>0, 0, 0, 1))</div><div><span class="Apple-tab-span" style="white-space:pre"> </span></div><div><span class="Apple-tab-span" style="white-space:pre"> </span># Extract a slice in the desired orientation</div><div><span class="Apple-tab-span" style="white-space:pre"> </span>reslice = vtk.vtkImageReslice()</div><div><span class="Apple-tab-span" style="white-space:pre"> </span>reslice.SetInputConnection(reader.GetOutputPort())</div><div><span class="Apple-tab-span" style="white-space:pre"> </span>reslice.SetOutputDimensionality(2)</div><div><span class="Apple-tab-span" style="white-space:pre"> </span>reslice.SetResliceAxes(axial)</div><div><span class="Apple-tab-span" style="white-space:pre"> </span>reslice.SetInterpolationModeToNearestNeighbor()</div><div><span class="Apple-tab-span" style="white-space:pre"> </span>reslice.Update()</div><div><span class="Apple-tab-span" style="white-space:pre"> </span></div><div><span class="Apple-tab-span" style="white-space:pre"> </span>resliceList.append(reslice)</div><div><span class="Apple-tab-span" style="white-space:pre"> </span></div><div><span class="Apple-tab-span" style="white-space:pre"> </span># Create a greyscale lookup table</div><div><span class="Apple-tab-span" style="white-space:pre"> </span>table = vtk.vtkLookupTable()</div><div><span class="Apple-tab-span" style="white-space:pre"> </span>table.SetRange(0, 2000) # image intensity range</div><div><span class="Apple-tab-span" style="white-space:pre"> </span>table.SetHueRange(ii*0.1,0.1+ii*0.1)</div><div><span class="Apple-tab-span" style="white-space:pre"> </span>table.SetValueRange(0.0, 1.0) # from black to white</div><div><span class="Apple-tab-span" style="white-space:pre"> </span>table.SetSaturationRange(0.0, 0.8) # no color saturation</div><div><span class="Apple-tab-span" style="white-space:pre"> </span>table.SetAlphaRange(0,0.8)</div><div><span class="Apple-tab-span" style="white-space:pre"> </span>table.SetRampToLinear()</div><div><span class="Apple-tab-span" style="white-space:pre"> </span>table.Build()</div><div><span class="Apple-tab-span" style="white-space:pre"> </span></div><div><span class="Apple-tab-span" style="white-space:pre"> </span>tableList.append(table)</div><div><span class="Apple-tab-span" style="white-space:pre"> </span></div><div><span class="Apple-tab-span" style="white-space:pre"> </span># Map the image through the lookup table</div><div><span class="Apple-tab-span" style="white-space:pre"> </span>color = vtk.vtkImageMapToColors()</div><div><span class="Apple-tab-span" style="white-space:pre"> </span>color.SetLookupTable(tableList[ii])</div><div><span class="Apple-tab-span" style="white-space:pre"> </span>color.SetInputConnection(resliceList[ii].GetOutputPort())</div><div><span class="Apple-tab-span" style="white-space:pre"> </span></div><div><span class="Apple-tab-span" style="white-space:pre"> </span>colorList.append(color)</div><div><span class="Apple-tab-span" style="white-space:pre"> </span></div><div><span class="Apple-tab-span" style="white-space:pre"> </span># Display the image</div><div><span class="Apple-tab-span" style="white-space:pre"> </span>actor = vtk.vtkImageActor()</div><div><span class="Apple-tab-span" style="white-space:pre"> </span>actor.SetInput(colorList[ii].GetOutput())</div><div><span class="Apple-tab-span" style="white-space:pre"> </span>actor.SetPosition(0, 0, ii*1.0)</div><div><span class="Apple-tab-span" style="white-space:pre"> </span>actorList.append(actor)</div><div><span class="Apple-tab-span" style="white-space:pre"> </span></div><div><span class="Apple-tab-span" style="white-space:pre"> </span>renderer.AddActor(actorList[ii])</div><div><span class="Apple-tab-span" style="white-space:pre"> </span></div><div>contour = vtk.vtkContourFilter()</div><div>contour.SetInputConnection(resliceList[0].GetOutputPort(0))</div><div>contour.GenerateValues(5,0,2000)</div><div>contour.SetComputeScalars(True)</div><div><br></div><div>polyMapper = vtk.vtkPolyDataMapper()</div><div>polyMapper.SetInputConnection(contour.GetOutputPort(0))</div><div># Uncomment next two to color contour lines</div><div># polyMapper.SetScalarModeToUsePointData()</div><div># polyMapper.SetScalarRange(0,2000)</div><div>polyMapper.SetScalarVisibility(False)</div><div>polyActor = vtk.vtkActor()</div><div>polyActor.SetMapper(polyMapper)</div><div>polyActor.GetProperty().SetColor(0,0,0)</div><div>polyActor.GetProperty().SetLineWidth(1)</div><div>polyActor.SetPosition(0,0,2.0)</div><div><br></div><div>renderer.AddActor(polyActor)</div><div><br></div><div># Set up the interaction</div><div>interactorStyle = vtk.vtkInteractorStyleImage()</div><div>interactor = vtk.vtkRenderWindowInteractor()</div><div>interactor.SetInteractorStyle(interactorStyle)</div><div><br></div><div>window = vtk.vtkRenderWindow()</div><div>window.AddRenderer(renderer)</div><div>window.SetInteractor(interactor)</div><div>window.Render()</div><div><br></div><div># Start interaction</div><div>interactor.Start()</div><div><br></div><div><br></div></div><div><br><div><div>On Apr 8, 2010, at 9:49 AM, Andrej Gluhov wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><span id="result_box" class="long_text"><span style="background-color: rgb(255, 255, 255);" title="Hi, Eric.">Hi, Eric.<br><br></span><span style="background-color: rgb(255, 255, 255);" title="Спасибо за
отклик."Images from vtkActor2D gets under 2D images" это
ошибка.">Thanks for the response. "Images from vtkActor2D gets under 2D
images" is a mistake. </span><span style="background-color: rgb(255, 255, 255);" title="Я конечно имел ввиду vtkActor.">Of course, I meant
vtkActor. </span><span style="background-color: rgb(255, 255, 255);" title="Проблема в том, что контур не видно при отображении, но если
растянуть окно, то vtkActor2D остается в левом нижнем углу, а Contours
видно на заднем фоне(под vtkActor2D).">The problem is that the contour is
not visible when displaying, but if you stretch the window and then
vtkActor2D remains in the lower left corner, and the Contours can be
seen in the background (under vtkActor2D). </span><span style="background-color: rgb(255, 255, 255);" title="Я пробовал изменить
позицию по координате Z, но это не дало результата.">I tried to change
the position of the coordinate Z, but it does not yield results. </span><span style="background-color: rgb(255, 255, 255);" title="vtkActor2D не
имеет координаты Z.">vtkActor2D has only x,y coordinates.<br></span><span style="background-color: rgb(255, 255, 255);" title="А окно я создаю для
того, чтобы потом прикрепить vtkImageViewer2 in WindowsForm.">And the
window I create in order to then attach vtkImageViewer2 in WindowsForm.</span></span><br><br><div class="gmail_quote">8 апреля 2010 г. 17:16 пользователь Eric E. Monson <span dir="ltr"><<a href="mailto:emonson@cs.duke.edu">emonson@cs.duke.edu</a>></span> написал:<br>
<blockquote class="gmail_quote" style="margin-top: 0pt; margin-right: 0pt; margin-bottom: 0pt; margin-left: 0.8ex; border-left-width: 1px; border-left-style: solid; border-left-color: rgb(204, 204, 204); padding-left: 1ex; position: static; z-index: auto; ">Hey Andrej,<br>
<br>
I'm trying to understand your problem, but I'm not sure exactly what you mean by "Images from vtkActor2D gets under 2D images". If the problem is that your images show up in front of your contour lines so that the contours are hidden, you can use something like:<br>
<br>
actor.SetPosition(0, 0, 0.01);<br>
<br>
to move the contours out in front of the image, and not have to change the polydata into imagedata.<br>
<br>
Also, you seem to be trying to create a vtkRenderWindow at the beginning of your code but you are using the vtkImageViewer2 later. If you want to view all of your objects in one window, you can use just one or the other (ImageViewer2 creates a render window for you).<br>
<br>
If this doesn't address your problem, be sure to try again to clarify the trouble you're having and I'm sure someone can help.<br>
<br>
Talk to you later,<br>
-Eric<br>
<br>
------------------------------------------------------<br>
Eric E Monson<br>
Duke Visualization Technology Group<br>
<div><div></div><div class="h5"><br>
<br>
On Apr 8, 2010, at 3:17 AM, Andrej Gluhov wrote:<br>
<br>
> Hi.<br>
> I have this problem of combining vtkActor and vtkActor2D for vtkImageViewer2D. Images from vtkActor2D gets under 2D images. I understand that the problem is converting vtkPolyData to vtkImageData, I could use vtkImageData, but do not know how to create it on the contour points.<br>
><br>
> //Create the RenderWindow, Renderer and RenderWindowInteractor<br>
> vtkRenderer ren1 = new vtkRenderer();<br>
> vtkRenderWindow renWin = reWin.GetRenderWindow();<br>
> renWin.AddRenderer(ren1);<br>
> vtkRenderWindowInteractor iren = new vtkRenderWindowInteractor();<br>
> iren.SetRenderWindow(renWin);<br>
><br>
> vtkDICOMImageReader DicomReader = new vtkDICOMImageReader();<br>
> DicomReader.SetFileName(m_strFilePath);<br>
> DicomReader.Update();<br>
><br>
> vtkImageMapper DicomMapper = new vtkImageMapper();<br>
> DicomMapper.SetInputConnection(DicomReader.GetOutputPort());<br>
> DicomMapper.SetColorWindow(255.0);<br>
> DicomMapper.SetColorLevel(127.5);<br>
><br>
> vtkActor2D DicomActor = new vtkActor2D();<br>
> DicomActor.SetMapper(DicomMapper);<br>
><br>
> vtkImageData vol = new vtkImageData();<br>
> vol.SetDimensions(512, 512, 1);<br>
> vol.SetSpacing(1, 1, 1);<br>
> vol.SetOrigin(0, 0, 0);<br>
> vol.SetScalarTypeToFloat();<br>
> vol.AllocateScalars();<br>
><br>
> vtkFloatArray scalars = new vtkFloatArray();<br>
> for (int i = 0; i < 512; i++)<br>
> {<br>
> for (int j = 0; j < 512; j++)<br>
> {<br>
> scalars.InsertTuple1(i * 512 + j, m_DoseMap[i * 512 + j]);<br>
> }<br>
> }<br>
><br>
> vol.GetPointData().SetScalars(scalars);<br>
> vol.Update();<br>
><br>
> double[] DoseRange = vol.GetScalarRange();<br>
><br>
> vtkLookupTable LUT = new vtkLookupTable();<br>
> LUT.SetTableRange(DoseRange[0], DoseRange[1]);//min dose,max dose<br>
> LUT.SetNumberOfColors(1000);<br>
> LUT.SetSaturationRange(1, 1);<br>
> LUT.SetHueRange(0.67, 0);//blue to red<br>
> LUT.SetAlphaRange(0, 0.3);//opacity<br>
> LUT.Build();<br>
><br>
> vtkImageMapToColors mapToRGBA = new vtkImageMapToColors();<br>
> mapToRGBA.SetInput(vol);<br>
> mapToRGBA.SetOutputFormatToRGBA();<br>
> mapToRGBA.SetLookupTable(LUT);<br>
><br>
> vtkImageMapper ColorMapper = new vtkImageMapper();<br>
> ColorMapper.SetInputConnection(mapToRGBA.GetOutputPort());<br>
> ColorMapper.SetColorWindow(255.0);<br>
> ColorMapper.SetColorLevel(127.5);<br>
><br>
> vtkActor2D ColorActor = new vtkActor2D();<br>
> ColorActor.SetMapper(ColorMapper);<br>
><br>
> vtkContourFilter cf = new vtkContourFilter();<br>
> cf.SetInput(vol);<br>
> //cf.SetValue(0, Convert.ToDouble(dataGridViewIsodose.Rows[i].Cells[1].Value));<br>
> cf.SetValue(0, 100);<br>
> cf.Update();<br>
><br>
> vtkPolyDataMapper mapper = new vtkPolyDataMapper();<br>
> mapper.SetInputConnection(cf.GetOutputPort());<br>
> mapper.ScalarVisibilityOff();<br>
><br>
> vtkActor actor = new vtkActor();<br>
> actor.SetMapper(mapper);<br>
> actor.GetProperty().SetColor(1, 0, 0);<br>
><br>
> vtkImageViewer2 viewer = new vtkImageViewer2();<br>
> viewer.SetupInteractor(iren);<br>
> viewer.GetRenderer().AddActor(DicomActor);<br>
> viewer.GetRenderer().AddActor(ColorActor);<br>
> viewer.GetRenderer().AddActor(actor);<br>
> viewer.GetRenderer().ResetCamera();<br>
> viewer.Render();<br>
><br>
> Thanks in advance.<br>
><br>
> --<br>
> С Уважением,<br>
> Андрей.<br>
> Best regards, Andrew<br>
</div></div>> _______________________________________________<br>
> Powered by <a href="http://www.kitware.com/" target="_blank">www.kitware.com</a><br>
><br>
> Visit other Kitware open-source projects at <a href="http://www.kitware.com/opensource/opensource.html" target="_blank">http://www.kitware.com/opensource/opensource.html</a><br>
><br>
> Please keep messages on-topic and check the VTK FAQ at: <a href="http://www.vtk.org/Wiki/VTK_FAQ" target="_blank">http://www.vtk.org/Wiki/VTK_FAQ</a><br>
><br>
> Follow this link to subscribe/unsubscribe:<br>
> <a href="http://www.vtk.org/mailman/listinfo/vtkusers" target="_blank">http://www.vtk.org/mailman/listinfo/vtkusers</a><br>
<br>
</blockquote></div><br><br clear="all"><br>-- <br>С Уважением,<br>Андрей.<br>Best regards, Andrew<br>
</blockquote></div><br></div></body></html>