<html><head><style type="text/css"><!-- DIV {margin:0px;} --></style></head><body><div style="font-family:"times new roman", "new york", "times", serif;font-size:12pt"><DIV>
<P class="MsoNormal">
<P class="MsoNormal">Hello,<BR><BR></P><P class="MsoNormal"><o:p> </o:p><SPAN lang="EN-US" style="mso-ansi-language:EN-US">I have a
problem with cut (cropping) and visualization a 3D image, which I produced by
the segmentation-method “watershed and region growing”. In the beginning I read
the DICOM images with itkImageSeriesReader, then I run watershed or region
growing algorithm, in order to connect itk with vtk and a 3D image produce and
for the displaying of segmentation I use vtkContourFilter. After segmentation
for cropping the images I use vtkBoxClipDataSet and vtkDataSetSurfaceFilter,
but the result is a black screen.<SPAN style="mso-spacerun:yes"> </SPAN>When I
read the images with vtkDICOMImageReader or itkImageSeriesReader and then segmented
with volume rendering, marching cubes watershed and region growing, then it
works completely well. The problem arises only when cut for watershed and
region growing algorithm. What am I doing wrong?<BR><BR> </SPAN></P>
<P class="MsoNormal"><SPAN lang="EN-US" style="mso-ansi-language:EN-US">To
understand my problem better, please see my screenshots on the following link.</SPAN></P>
<P class="MsoNormal"><SPAN lang="EN-US" style="mso-ansi-language:EN-US"><o:p> </o:p></SPAN></P>
<P class="MsoNormal"><SPAN lang="EN-US" style="mso-ansi-language:EN-US"><A href="http://www.wopsys.com/bilder/segmentation.html">http://www.wopsys.com/bilder/segmentation.html</A></SPAN></P>
<P class="MsoNormal"><SPAN lang="EN-US" style="mso-ansi-language:EN-US"><o:p> </o:p></SPAN></P>
<P class="MsoNormal"><SPAN lang="EN-US" style="mso-ansi-language:EN-US">To
determine whether there missing of any geometrical informations. I had sent the
following outputs of tests.</SPAN></P>
<P class="MsoNormal"><SPAN lang="EN-US" style="mso-ansi-language:EN-US"><A href="http://www.wopsys.com/Outputs/Output_MarchingCubes.txt">http://www.wopsys.com/Outputs/Output_MarchingCubes.txt</A><BR>
<A href="http://www.wopsys.com/Outputs/Output_VolumeRendering.txt">http://www.wopsys.com/Outputs/Output_VolumeRendering.txt</A><BR>
<A href="http://www.wopsys.com/Outputs/Output_Watershed.txt">http://www.wopsys.com/Outputs/Output_Watershed.txt</A><BR>
<A href="http://www.wopsys.com/Outputs/Output_RegionGrowing.txt">http://www.wopsys.com/Outputs/Output_RegionGrowing.txt</A></SPAN></P>
<SPAN lang="EN-US" style="mso-ansi-language:EN-US"></SPAN></P><P class="MsoNormal"><SPAN lang="EN-US" style="mso-ansi-language:EN-US">Thanks.<BR></SPAN></P><P class="MsoNormal"><SPAN lang="EN-US" style="mso-ansi-language:EN-US"><BR></SPAN></P><BR><DIV style="font-family:times new roman, new york, times, serif;font-size:12pt;"><FONT size="2" face="Tahoma"><HR size="1"><B><SPAN style="font-weight:bold;">Von:</SPAN></B> Luis Ibanez <luis.ibanez@kitware.com><BR><B><SPAN style="font-weight:bold;">An:</SPAN></B> Yusuf OEZBEK <nasil122002@yahoo.de><BR><B><SPAN style="font-weight:bold;">CC:</SPAN></B> ITK Mailinglist <insight-users@itk.org>; VTK Mailinglist <vtkusers@vtk.org>; nasil122002@gmail.com<BR><B><SPAN style="font-weight:bold;">Gesendet:</SPAN></B> Sonntag, den 18. April 2010, 1:12:33 Uhr<BR><B><SPAN style="font-weight:bold;">Betreff:</SPAN></B> Re: [Insight-users] Cropping Problem<BR></FONT><BR><BR>Hi Yusuf,<BR><BR><BR>You should
start by separating the segmentation<BR>problem from the visualization problem.<BR><BR><BR>Please do the following:<BR><BR>0) Just after you call:<BR><BR> reader->SetFileNames( fileNames );<BR>
reader->Update();<BR><BR> add<BR><BR> reader->GetOutput()->Print( std::cout );<BR><BR> then<BR><BR>1) Run the ITK segmentation pipeline until you <BR> get the segmented image, and save it as a<BR> MetaImage file (with extension .mhd)<BR>
<BR><BR>2) Post the content of the resulting .mhd file <BR> back to the list, as well as the print out of <BR> the statement:<BR><BR> reader->GetOutput()->Print( std::cout );<BR>
<BR><BR>3) Use Paraview (<A rel="nofollow" target="_blank" href="http://www.paraview.org">www.paraview.org</A>) to load the<BR> resulting image and extracting an iso-surface.<BR> Compare that iso-surface with the one that <BR> you get by directly from the DICOM image.<BR>
(and let us know what you find).<BR><BR><BR>With the information from (2) we should be <BR>able to identify if any geometrical information<BR>is missing from the ITK processing pipeline.<BR><BR><BR> Thanks<BR><BR><BR>
Luis<BR><BR><BR></DIV></DIV><DIV style="font-family:times new roman, new york, times, serif;font-size:12pt"><DIV style="">
</DIV><BR></DIV>
</div><br></body></html>