<div dir="ltr"><div><div><span id="result_box" class="" lang="en"><span class="">hello</span> <span class="">itk</span> <span class="">users</span>,<br></span><span id="result_box" class="" lang="en"><br> <span class="">I am</span> <span class="">a beginner</span> <span class="">ITK users</span><span>.</span> <span class="">I am trying to</span> <span class="">segment a</span> <span class="">volume using</span> <span class="">region growing</span> <span class="">algorithm.</span> <span class="">I am able to</span> <span class="">do it in</span> <span class="">a 2D image</span> <span class="">(</span><span>dicom</span><span>) but not for</span> <span class="">a volume.</span><br> <span class="">1)</span> <span class="">I saw that</span> <span class="">for 3D</span> <span class="">region growing</span> <span class="">the</span> <span class="">seed</span> <span class="">point have</span> <span class="">3 coordinates</span> <span class="">(x,</span> <span class="">y, z)</span><span>.</span> <span class="">In my case,</span> <span class="">z is</span> <span class="">the</span> <span class="">slice</span> <span class="">number (</span><span>i</span> <span class="">am using</span> <span class="">.img</span> <span class="">format).</span> <span class="">Can I use a</span>n (x,y) <span class="">seed</span> <span class="">point</span><span></span><span>?</span><br> <span class="">2)</span> <span class="">I have another idea, </span><span class=""> read</span>ing <span class="">an image series </span><span class=""> and</span> <span class="">applied the filter</span> <span class="">for each image and</span> <span class="">then</span> <span class="">writing them</span> <span class="">in volume.</span> <span class="">Are there any</span> <span class="">examples ( applied a filter in image series) that</span> <span class="">one</span> <span class="">could help me?</span><br><br></span></div><span id="result_box" class="" lang="en">regards<br></span></div><span id="result_box" class="" lang="en">Samah.<br></span><div><div><span id="result_box" class="" lang="en"><br> <span class="">My</span> <span class="">source</span> <span class="">code</span> <span class="">is as follows<br>================================<br>#include "itkConnectedThresholdImageFilter.h"<br>#include "itkImage.h"<br>#include "itkImageFileReader.h"<br>#include "itkImageFileWriter.h"<br><br>int main( int argc, char *argv[])<br>{<br>  typedef   signed short         InternalPixelType;<br>  typedef   unsigned char        OutputPixelType;<br>  const     unsigned int           Dimension = 3;<br>  typedef itk::Image< InternalPixelType, Dimension > InternalImageType;<br><br>  typedef itk::Image< OutputPixelType, Dimension > OutputImageType;<br><br>  typedef  itk::ImageFileReader< InternalImageType > ReaderType;<br>  typedef  itk::ImageFileWriter<  OutputImageType  > WriterType;<br><br>  ReaderType::Pointer reader = ReaderType::New();<br>  WriterType::Pointer writer = WriterType::New();<br><br>    reader->SetFileName( "/home/doctorant/Téléchargements/Images/in.img" );<br>    writer->SetFileName( "/home/doctorant/Téléchargements/Images/out.img");<br><br><br>  typedef itk::ConnectedThresholdImageFilter< InternalImageType,<br>OutputImageType > ConnectedFilterType;<br><br>  ConnectedFilterType::Pointer connectedThreshold =<br>ConnectedFilterType::New();<br><br>  connectedThreshold->SetInput( reader->GetOutput() );<br>  writer->SetInput( connectedThreshold->GetOutput() );<br><br><br>  const InternalPixelType lowerThreshold =  -1050  ;<br>  const InternalPixelType upperThreshold = -600;<br><br>  connectedThreshold->SetLower(  lowerThreshold  );<br>  connectedThreshold->SetUpper(  upperThreshold  );<br><br>  connectedThreshold->SetReplaceValue( 255 );<br><br><br><br><br>  InternalImageType::IndexType  index;<br>  index[0] = 250;<br>  index[1] = 239;<br><br>  connectedThreshold->SetSeed( index );<br><br><br>  try<br>    {<br>    writer->Update();<br>    }<br>  catch( itk::ExceptionObject & excep )<br>    {<br>    std::cerr << "Exception caught !" << std::endl;<br>    std::cerr << excep << std::endl;<br>    }<br><br>  return 0;<br>}<br></span></span></div></div></div>