<div>i have the exatly desire as finsher has...</div>
<div>i am waiting for the answeer to...</div>
<div>best regards<br><br></div>
<div class="gmail_quote">2010/9/19 finsher <span dir="ltr"><<a href="mailto:zhuyinghua01@gmail.com">zhuyinghua01@gmail.com</a>></span><br>
<blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid"><br>Hi all,<br>I need to get the Maximum intensity projection(MIP) slices.<br>so I have a series of MR images. I use these series to do MIP(or volume<br>
render etc...). when MIP rendering completed I need to get the rendered<br>slices. if i use vtkResliceImage dor imageplanewidget s the input data is<br>the voxel raw data, i got the slices is the original image slices, not cut<br>
from MIP volume slices.<br>so could you help me how to get the correct slices?<br><br>best regards.<br><br>and below is my code :<br><br>//step 1: read raw image series data<br>vtkVolume16Reader* v16 = vtkVolume16Reader::New();<br>
v16->SetFilePrefix("d:\\aaa\\new\\a");<br>...<br>v16->Update();<br><br>//step.2: do MIP<br>vtkPiecewiseFunction *opacityTransferFunc = vtkPiecewiseFunction::New();<br>opacityTransferFunc->AddPoint(0.0, 0.0);<br>
...<br>vtkColorTransferFunction *colorTransferFunc =<br>vtkColorTransferFunction::New();<br>colorTransferFunc->AddRGBPoint(0.0, 0.0, 0.0, 0.0);<br>...<br>volumeProperty = vtkVolumeProperty::New();<br>volumeProperty->SetColor(colorTransferFunc);<br>
volumeProperty->SetScalarOpacity(opacityTransferFunc);<br><br>//step 3: set mip function and create volume data<br>vtkVolumeRayCastMIPFunction *mipFunc = vtkVolumeRayCastMIPFunction::New();<br>volumeMapper = vtkVolumeRayCastMapper::New();<br>
volumeMapper->SetVolumeRayCastFunction(mipFunc);<br>volumeMapper->SetInput(v16->GetOutput());<br>volume = vtkVolume::New();<br>volume->SetMapper(volumeMapper);<br>volume->SetProperty(volumeProperty);<br>volume->Update();<br>
<br>//step.4: renderer volume show rendered window<br>render->AddVolume(volume);<br>renderWindow->Render();<br>...<br><br>//setp.5: get MIP slices<br>vtkImagePlaneWidget* planeWidget = vtkImagePlaneWidget::New();<br>
...<br>planeWidget->SetInput(v16->GetOutput()); <----i think this input data could<br>not be correct. but i don't know how ...<br><br>//step.6 display the slices<br>imageViewer2->SetInput(planeWidget->GetResliceOutput());<br>
renderWindowInteractor->Initialize();<br>imageViewer2->Render();<br>renderWindowInteractor->Start();<br>....<br>--<br>View this message in context: <a href="http://vtk.1045678.n5.nabble.com/How-to-Get-rendered-Volume-Slices-tp2845243p2845243.html" target="_blank">http://vtk.1045678.n5.nabble.com/How-to-Get-rendered-Volume-Slices-tp2845243p2845243.html</a><br>
Sent from the VTK - Users mailing list archive at Nabble.com.<br>_______________________________________________<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>
</blockquote></div><br><br clear="all"><br>-- <br>陆阳<br>清华大学医学院生物医学工程系<br>Yang Lu<br>Department of Biomedical Engineering,<br>School of Medicine, Tsinghua University<br><br>