[ITK-users] Surface extraction problem

[Francisco López-Franca]

Hi everybody.
I've created a 3D cylinder by using the class CylinderSpatialObject. After
that, I pass the cylinder object to an object of class
SpatialObjectToImageFilter in order to prepare the spatial object to be
written to file.
I've set both the cylinder object and the image filter a value of 0 for
background and 255 for foreground.
Finally, I pass, as input, the image filter output to an image file writer
(ImageFileWriter) object to save it to file (Nifti format, for instance).
The result is the attached "cilindro.jpg" figure.

Up to here, all is ok. The problem comes next.
What I want is to extract the surface of the 3D cylinder. I search for the
web and found an example, from which I followed the implementation.
I copy below my code and attach the result of the surface extraction
("cilindro_superficie.jpg"), hoping that anyone see if I am committing any
error. I expected just the boundary of the 3D cylinder (with no inside
points), and it seems like the bounding box of the cylinder what I've got
including the inside points.

Thank you very much and my best regards.

My code is (I've erased the obvious code):


    typedef itk::Image< PixelType, Dimension > ImageType;
    typedef itk::Mesh<double, 3> MeshType;

    typedef itk::ImageFileReader< ImageType > ReaderType;
    ReaderType::Pointer reader = ReaderType::New();
    reader->SetFileName(
inputFileName);

        reader->Update();

    typedef itk::BinaryMask3DMeshSource< ImageType, MeshType >
MeshSourceType;
    MeshSourceType::Pointer meshSource = MeshSourceType::New();

    meshSource->SetObjectValue( 255);

    meshSource->SetInput( reader->GetOutput() );

        meshSource->Update();

    typedef itk::MeshSpatialObject<MeshType> SpatialObjectType;
    SpatialObjectType::Pointer mesh = SpatialObjectType::New();
    mesh->SetMesh(meshSource->GetOutput());

    typedef itk::SpatialObjectToImageFilter<SpatialObjectType, ImageType>
MeshFilterType;
    MeshFilterType::Pointer meshFilter = MeshFilterType::New();

    meshFilter->SetSize(size);
    meshFilter->SetOutsideValue(0);
    meshFilter->SetInsideValue(255);
    meshFilter->SetInput(mesh);

        meshFilter->Update();

    typedef itk::ImageFileWriter<ImageType> ImageWriterType;
    ImageWriterType::Pointer imageWriter = ImageWriterType::New();

    imageWriter->SetFileName(outputFileName);
    imageWriter->SetInput(meshFilter->GetOutput());

        imageWriter->Update();

<http://itk-insight-users.2283740.n2.nabble.com/file/n7585840/cilindro.jpg> 
<http://itk-insight-users.2283740.n2.nabble.com/file/n7585840/cilindro_superficie.jpg> 



--
View this message in context: http://itk-insight-users.2283740.n2.nabble.com/Surface-extraction-problem-tp7585840.html
Sent from the ITK Insight Users mailing list archive at Nabble.com.