You use an isovalue of 1 in your code, or is it meant to be 0.5 (half way between Foreground and background). <br><br>You could try also using vtkThresholdPoints. This will generate a polydata with vertex cells at each voxel. You can visualize the result with a vtkPolyDataMapper.<br>
<br>A more appropriate way to visualize is to construct a polydata out of the skeletonized data. Just visit every voxel and its neighbors add to a vtkCellArray, the centers of the voxels with any connected neighbors. You can then visualize the construcuted polydata. You get to set the line width etc... <br>
Or you get can pass the data through a tube filter and use a vtkPolyDataMapper...<br><br><div class="gmail_quote">On Mon, Jan 10, 2011 at 3:20 PM, Xiaopeng Yang <span dir="ltr"><<a href="mailto:yxp233@postech.ac.kr">yxp233@postech.ac.kr</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><div link="blue" vlink="purple" lang="EN-US"><div><p class="MsoNormal">Hi guys,</p><p class="MsoNormal">
</p><p class="MsoNormal">How are you doing?</p><p class="MsoNormal"> </p><p class="MsoNormal">I applied ITK Binary3D Thinning Image Filter to skeletonize portal vein of liver. Then I tried to visualize the result as Iso-surface generated by vtk. But in the screen, I can just see a few discrete and very small point (if you do not see carefully, you cannot recognize these points). Then I visualized my result as Iso-surface by MITK. The vessel skeleton is shown continuously and looks very good. Here are the two results (the former one is generated by my code; the later one is from MITK):</p>
<p class="MsoNormal"> </p><p class="MsoNormal"><img src="https://mail.google.com/a/kitware.com/?ui=2&ik=856cfef1d6&view=att&th=12d6f5686e9e9c6a&attid=0.0.1&disp=emb&zw" alt="cid:image006.jpg@01CB9EAA.770905A0" width="327" height="341"><img src="https://mail.google.com/a/kitware.com/?ui=2&ik=856cfef1d6&view=att&th=12d6f5686e9e9c6a&attid=0.0.2&disp=emb&zw" alt="cid:image007.jpg@01CB9EAA.770905A0" width="388" height="334"></p>
<p class="MsoNormal"> </p><p class="MsoNormal">Here is my code:</p><p class="MsoNormal"> </p><p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: DotumChe;"> reader->SetFileName( <span style="color: rgb(163, 21, 21);">"D:/Hybrid Method/Hybrid/BSP/PortalSkeleton/3D.dcm"</span> );</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: DotumChe;"> reader->Update();<span style="color: green;"></span></span></p><p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: DotumChe;"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: DotumChe;"> <span style="color: blue;">typedef</span> itk::VTKImageExport< InputImageType > ExportFilterType;</span></p><p class="MsoNormal" style="">
<span style="font-size: 10pt; font-family: DotumChe;"> </span></p><p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: DotumChe;"> ExportFilterType::Pointer itkExporter = ExportFilterType::New();</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: DotumChe;"> </span></p><p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: DotumChe;"> itkExporter->SetInput( reader->GetOutput() );</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: DotumChe;"> </span></p><p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: DotumChe;"> vtkImageImport* vtkImporter = vtkImageImport::New(); </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: DotumChe;"> ConnectPipelines(itkExporter, vtkImporter);</span></p><p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: DotumChe;"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: DotumChe;"> vtkRenderer *ren4 = vtkRenderer::New();</span></p><p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: DotumChe;"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: DotumChe;"> vtkRenderWindow *renWin = vtkRenderWindow::New();</span></p><p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: DotumChe;"> renWin->AddRenderer(ren4);</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: DotumChe;"> </span></p><p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: DotumChe;"> vtkRenderWindowInteractor *iren = vtkRenderWindowInteractor::New();</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: DotumChe;"> iren->SetRenderWindow(renWin);</span></p><p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: DotumChe;"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: DotumChe;"> <span style="color: green;">//Create an iso-surface using Marching Cubes</span></span></p><p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: DotumChe;"> vtkContourFilter *iso = vtkContourFilter::New();</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: DotumChe;"> iso->SetInput(vtkImporter->GetOutput());</span></p><p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: DotumChe;"> iso->ComputeScalarsOn();</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: DotumChe;"> iso->SetValue(1, 1);</span></p><p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: DotumChe; color: green;"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: DotumChe;"> <span style="color: green;">//vtkDataSetMapper *isoMapper = vtkDataSetMapper::New();</span></span></p><p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: DotumChe;"> vtkPolyDataMapper *isoMapper = vtkPolyDataMapper::New();</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: DotumChe;"> isoMapper->SetInput(iso->GetOutput());</span></p><p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: DotumChe;"> isoMapper->ScalarVisibilityOff();</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: DotumChe; color: green;"> </span></p><p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: DotumChe;"> vtkActor *isoActor = vtkActor::New();</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: DotumChe;"> isoActor->SetMapper(isoMapper);</span></p><p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: DotumChe;"> isoActor->GetProperty()->SetColor(1, 0, 0);</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: DotumChe;"> </span></p><p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: DotumChe;"> ren4->AddActor(isoActor);</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: DotumChe;"> </span></p><p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: DotumChe;"> renWin->SetSize(512,512);</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: DotumChe;"> ren4->ResetCamera();</span></p><p class="MsoNormal" style="text-indent: 10.5pt;"><span style="font-size: 10pt; font-family: DotumChe;">ren4->SetBackground(1,1,1);</span></p>
<p class="MsoNormal" style="text-indent: 10.5pt;"><span style="font-size: 10pt; font-family: DotumChe;"> </span></p><p class="MsoNormal" style="text-indent: 10.5pt;"><span style="font-size: 10pt; font-family: DotumChe;"> </span></p>
<p class="MsoNormal">Thank you,</p><p class="MsoNormal">Yang</p></div></div><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>