<html>
<head>
<style>
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
FONT-SIZE: 10pt;
FONT-FAMILY:Tahoma
}
</style>
</head>
<body class='hmmessage'>
Hi <br><br>I'm triying to convert an Image that I've created with itk and it pass a process to apply confidenceconected filter to vtk.<br><br>The problem is that when i see the image after the conversion it is draw the planes xz and yz in the limit of image<br><br>if we use axial view we can see the connected image ( X is supposed to be the image) and somwething like that: &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp; <br>_____&nbsp; &nbsp;&nbsp; &nbsp;&nbsp;  &nbsp;&nbsp;  &nbsp;&nbsp;  &nbsp;&nbsp;  &nbsp;&nbsp;  &nbsp;&nbsp;  &nbsp;&nbsp;  &nbsp;&nbsp;  &nbsp;&nbsp;  &nbsp;&nbsp;  &nbsp;&nbsp;  &nbsp;&nbsp;  &nbsp;&nbsp;  &nbsp;&nbsp;  &nbsp;&nbsp;  &nbsp;&nbsp;  &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;  &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp; &nbsp; &nbsp;&nbsp;&nbsp; &nbsp; |&nbsp;&nbsp;&nbsp;  &nbsp;&nbsp;  &nbsp;&nbsp;  &nbsp;&nbsp;  &nbsp;&nbsp;  &nbsp;&nbsp;  &nbsp;&nbsp;  &nbsp;&nbsp;  &nbsp;&nbsp;  &nbsp;&nbsp;  &nbsp;&nbsp;  &nbsp;&nbsp;  &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp; &nbsp;  <br>&nbsp; X&nbsp;&nbsp;&nbsp; |<br>&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;  &nbsp;&nbsp;  &nbsp;&nbsp;  &nbsp;&nbsp;  &nbsp;&nbsp;  &nbsp;&nbsp;  &nbsp;&nbsp;  &nbsp;&nbsp;  &nbsp;&nbsp;  &nbsp;&nbsp;  &nbsp;&nbsp;  &nbsp;&nbsp;  &nbsp;&nbsp;  &nbsp;&nbsp;  &nbsp;&nbsp;  &nbsp;&nbsp;  &nbsp;&nbsp;  &nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; <br><br>and i need to see the image without planes xz yz <br><br><br>here is the code if it can serve as aid<br><br>&nbsp;&nbsp;&nbsp;  ///The application is developed in visual studio an using a GUI <br><br>&nbsp;&nbsp;&nbsp; const&nbsp;&nbsp;&nbsp;&nbsp; unsigned int&nbsp;&nbsp;&nbsp; Dimension = 3;<br><br>&nbsp;&nbsp;&nbsp; typedef&nbsp;&nbsp; float&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; PixelType;<br>&nbsp;&nbsp;&nbsp; typedef itk::Image&lt; PixelType, Dimension &gt;&nbsp; ImageType;<br>&nbsp;&nbsp;&nbsp; typedef unsigned short OutputPixelType;<br>&nbsp;&nbsp;&nbsp; typedef itk::Image&lt; OutputPixelType, Dimension &gt; OutputImageType;<br>&nbsp;&nbsp;&nbsp; typedef itk::CastImageFilter&lt; ImageType, OutputImageType &gt;&nbsp;&nbsp;&nbsp; CastingFilterType;<br>&nbsp;&nbsp;&nbsp; typedef itkMeshTovtkPolyData::TriangleMeshType MeshType;<br>&nbsp;&nbsp;&nbsp; typedef itk::BinaryMask3DMeshSource&lt; OutputImageType, MeshType &gt; MeshSourceType;<br><br>&nbsp;&nbsp;&nbsp; typedef itk::CurvatureFlowImageFilter&lt; ImageType, ImageType &gt;&nbsp; CurvatureFlowImageFilterType;<br>&nbsp;&nbsp;&nbsp; typedef itk::VTKPolyDataWriter&lt;MeshType&gt;&nbsp;&nbsp; WriterType;<br>&nbsp;&nbsp;&nbsp; typedef itk::ConfidenceConnectedImageFilter&lt;OutputImageType, OutputImageType&gt; ConnectedFilterType; /////COMENTAR P4<br><br>&nbsp;&nbsp;&nbsp; MeshSourceType::Pointer meshSource = MeshSourceType::New();<br>&nbsp;&nbsp;&nbsp; CastingFilterType::Pointer caster = CastingFilterType::New();<br>&nbsp;&nbsp;&nbsp; CurvatureFlowImageFilterType::Pointer smoothing = CurvatureFlowImageFilterType::New();<br>&nbsp;&nbsp;&nbsp; ConnectedFilterType::Pointer confidenceConnected = ConnectedFilterType::New();/////COMENTAR P4<br><br><br><br>&nbsp;&nbsp; ////// I have omitted the reading code&nbsp; of&nbsp; ITK <br><br><br>&nbsp;&nbsp;&nbsp; caster-&gt;SetInput(inputImage);&nbsp;&nbsp;&nbsp; <br>&nbsp;&nbsp;&nbsp; caster-&gt;ReleaseDataFlagOn();<br>&nbsp;&nbsp;&nbsp;&nbsp; <br>&nbsp;&nbsp;&nbsp; confidenceConnected-&gt;SetInput( caster-&gt;GetOutput());&nbsp;&nbsp;  <br>&nbsp;&nbsp;&nbsp; confidenceConnected-&gt;SetMultiplier( 2.5 ); // confMultiplier<br>&nbsp;&nbsp;&nbsp; confidenceConnected-&gt;SetNumberOfIterations( 1 ); //confNumIterations<br>&nbsp;&nbsp;&nbsp; confidenceConnected-&gt;SetReplaceValue( 180 );&nbsp; //confReplaceValue<br>&nbsp;&nbsp;&nbsp; <br>&nbsp;&nbsp;&nbsp; OutputImageType::IndexType&nbsp; indexseed;<br>&nbsp;&nbsp;&nbsp; confidenceConnected-&gt;ClearSeeds();<br><br>&nbsp;&nbsp;&nbsp; int semillaZ=(pVO-&gt;sizeZ()-valsemz);<br>&nbsp;&nbsp;&nbsp; indexseed[0] = valseedx;// x<br>&nbsp;&nbsp;&nbsp; indexseed[1] = valseedy;// y<br>&nbsp;&nbsp;&nbsp; indexseed[2] = valseedz;// z<br>&nbsp;&nbsp;&nbsp; confidenceConnected-&gt;SetSeed( indexseed );<br>&nbsp;&nbsp;&nbsp; confidenceConnected-&gt;SetInitialNeighborhoodRadius( 2 ); //confRadius<br>&nbsp;&nbsp;&nbsp; confidenceConnected-&gt;ReleaseDataFlagOn();<br>&nbsp;&nbsp;&nbsp; confidenceConnected-&gt;Update();<br>&nbsp;&nbsp;&nbsp; <br>&nbsp;&nbsp;&nbsp; OutputImageType::Pointer maskImage=confidenceConnected-&gt;GetOutput(); <br>&nbsp;&nbsp;&nbsp; <br>&nbsp;&nbsp;&nbsp; ///here it is draw the planes (xz,yz) of the image limit<br>&nbsp;&nbsp;&nbsp; maskImage-&gt;Update();<br>&nbsp;&nbsp;&nbsp; meshSource-&gt;ReleaseDataFlagOn();&nbsp;&nbsp;&nbsp; <br>&nbsp;&nbsp;&nbsp; meshSource-&gt;SetInput(maskImage);<br>&nbsp;&nbsp;&nbsp; meshSource-&gt;Update();<br><br>&nbsp; ///here we convert to vtkpolydata <br><br>&nbsp; vtkPolyData *vtkPolyDataOut = NULL;<br><br><br>&nbsp;&nbsp;&nbsp; itkMeshTovtkPolyData * itkMesh2PolyData = new itkMeshTovtkPolyData();<br>&nbsp;&nbsp;&nbsp; itkMesh2PolyData-&gt;SetInput (meshSource-&gt;GetOutput());<br>&nbsp;&nbsp;&nbsp; vtkPolyDataOut = itkMesh2PolyData-&gt;GetOutput();<br>&nbsp;<br>&nbsp;&nbsp;&nbsp; vtkCleanPolyData *pClean = vtkCleanPolyData::New();<br>&nbsp;&nbsp;&nbsp; pClean-&gt;SetInput(vtkPolyDataOut);<br><br>&nbsp;&nbsp;&nbsp; vtkDecimatePro *pDecimador = vtkDecimatePro::New();<br>&nbsp;&nbsp;&nbsp; pDecimador-&gt;SetInput(pClean-&gt;GetOutput());<br>&nbsp;&nbsp;&nbsp; pDecimador-&gt;PreserveTopologyOn ();<br>&nbsp;&nbsp;&nbsp; pDecimador-&gt;SetTargetReduction(.30); <br>&nbsp;&nbsp;&nbsp; <br>&nbsp;&nbsp;&nbsp; vtkSmoothPolyDataFilter *pSmoother = vtkSmoothPolyDataFilter::New();<br>&nbsp;&nbsp;&nbsp; pSmoother-&gt;SetInput (pDecimador-&gt;GetOutput());<br>&nbsp;&nbsp;&nbsp; pSmoother-&gt;SetNumberOfIterations (3);<br>&nbsp;&nbsp;&nbsp; pSmoother-&gt;BoundarySmoothingOff ();<br>&nbsp;&nbsp;&nbsp; pSmoother-&gt;FeatureEdgeSmoothingOff ();<br><br>&nbsp;&nbsp;&nbsp; vtkStripper * stripper = vtkStripper::New();<br>&nbsp;&nbsp;&nbsp; stripper-&gt;SetInput(pSmoother-&gt;GetOutput());<br>&nbsp;&nbsp;&nbsp; stripper-&gt;Update();<br><br>&nbsp;&nbsp;&nbsp; vtkPolyData *pPolyData = stripper-&gt;GetOutput();<br><br>&nbsp;&nbsp; //Here convert the vtkPolydata to my own type of visualization<br>&nbsp;&nbsp;&nbsp; triObject_t *obj = build3dObject (pPolyData, bInvertirCaras);<br><br><br>greetings<br>Manuel<br><br><br><br /><hr />Discover the new Windows Vista <a href='http://search.msn.com/results.aspx?q=windows+vista&mkt=en-US&form=QBRE' target='_new'>Learn more!</a></body>
</html>