<div>hi Italo,</div>
<div>Thanks for your help.</div>
<div>i find the vtkContourFilter class only generate the isosurface volume.</div>
<div>the attachment is i do the test(use you give me code) result.</div>
<div>i want to generated 2D images of each 6° degree around the <strong><font style="BACKGROUND-COLOR: #ffff00">MIP </font></strong>volume.</div>
<div>i don't know if VTK has MIP filter like vtkContourFiler. so i can use than to generate MIP volume .</div>
<div> </div>
<div>Thanks.</div>
<div>best regards.</div>
<div>Yinghua.<br><br></div>
<div class="gmail_quote">2010/9/19 Italo Ribeiro <span dir="ltr"><<a href="mailto:italo.ribeiro@gmail.com">italo.ribeiro@gmail.com</a>></span><br>
<blockquote style="BORDER-LEFT: #ccc 1px solid; MARGIN: 0px 0px 0px 0.8ex; PADDING-LEFT: 1ex" class="gmail_quote">If i understand you want create 3D model of slice? Try this:<br><br>vtkDICOMImageReader *reader = vtkDICOMImageReader::New();<br>
reader->SetDirectoryName("C:/cranio");<br> reader->Update(); <br><br> vtkContourFilter *skinExtractor = vtkContourFilter::New();<br> skinExtractor->SetInputConnection(reader->GetOutputPort());<br>
skinExtractor->SetValue(0, 500); //skin human<br><br> vtkPolyDataNormals *skinNormals = vtkPolyDataNormals::New();<br> skinNormals->SetInputConnection(skinExtractor->GetOutputPort());<br> skinNormals->SetFeatureAngle(60.0);<br>
<br> vtkStripper *skinStripper = vtkStripper::New();<br> skinStripper->SetInputConnection(skinNormals->GetOutputPort());<br><br> vtkPolyDataMapper *skinMapper = vtkPolyDataMapper::New();<br> skinMapper->SetInputConnection(skinStripper->GetOutputPort());<br>
skinMapper->ScalarVisibilityOff();<br><br> vtkActor *skin = vtkActor::New();<br> skin->SetMapper(skinMapper);<br><br>vtkRenderer *renderer2 = vtkRenderer::New();<br> renderer2->SetViewport(0.5, 0.0, 1.0, 0.5);<br>
renderer2->SetBackground(0.0, 0.0, 0.0);<br> renderWindow->AddRenderer(renderer2);<br><br> renderer2->AddActor(skin);<br><br>bye.<br><br>
<div class="gmail_quote">2010/9/19 陆阳 <span dir="ltr"><<a href="mailto:thuluyang@gmail.com" target="_blank">thuluyang@gmail.com</a>></span>
<div>
<div></div>
<div class="h5"><br>
<blockquote style="BORDER-LEFT: rgb(204,204,204) 1px solid; MARGIN: 0pt 0pt 0pt 0.8ex; PADDING-LEFT: 1ex" class="gmail_quote">
<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" target="_blank">zhuyinghua01@gmail.com</a>></span>
<div>
<div></div>
<div><br>
<blockquote style="BORDER-LEFT: rgb(204,204,204) 1px solid; MARGIN: 0px 0px 0px 0.8ex; PADDING-LEFT: 1ex" class="gmail_quote"><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></div></div><font color="#888888"><br><br clear="all"><br>-- <br>陆阳<br>清华大学医学院生物医学工程系<br>Yang Lu<br>Department of Biomedical Engineering,<br>School of Medicine, Tsinghua University<br><br></font><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>
<br></blockquote></div></div></div><font color="#888888"><br><br clear="all"><br>-- <br>Ítalo Mendes<br><a href="http://ogre.irados.org/" target="_blank">ogre.irados.org</a><br></font><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>
<br></blockquote></div><br>