[vtkusers] 3D image as N cubes

[Francisco Lopez de la Franca]

Hi.
I've created an example (see below) but everything is OK except that the
output image is not created. I'd be very pleased if someone can help me
(maybe with an explanatory example or whatever).

    char* filename = new char[255];
    strcpy(filename, "nifti_image.nii");

    vtkSmartPointer<vtkNIFTIImageReader> reader =
vtkSmartPointer<vtkNIFTIImageReader>::New();

    reader->SetFileName(filename);
    reader->Update();

    //Create the octree
    vtkSmartPointer<vtkOctreePointLocator> octree =
vtkSmartPointer<vtkOctreePointLocator>::New();
    octree->SetDataSet((vtkDataSet*) reader->GetOutput());
    octree->BuildLocator();

    std::cout << "Number of cells in tree: " <<
octree->GetDataSet()->GetNumberOfCells() << std::endl;
    std::cout << "Number of points in tree: " <<
octree->GetDataSet()->GetNumberOfPoints() << std::endl;

    char* outfile = new char[255];
    strcpy(outfile, "OcTree.vtp");

    vtkSmartPointer<vtkXMLPolyDataWriter> writer =
vtkSmartPointer<vtkXMLPolyDataWriter>::New();
    writer->SetInputData((vtkPolyData*) octree->GetDataSet());
    writer->SetFileName(outfile);
    writer->Write();


2015-05-29 12:40 GMT+02:00 Francisco Lopez de la Franca <
franciscolopezdelafranca at gmail.com>:

> Hi everybody,
> I'd like to represent a 3D image as a fixed number of cubes, I mean,
> something like a spatial occupancy matrix in 3D.
>
> Could you, please, recommend me a fast way to do it or does it exist  a
> data structure in VTK to achieve this?
>
> I've found the OcTree structure but I donot know if it is valid for this.
>
> Thank you very much.
> My best regards.
> /Francisco
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://public.kitware.com/pipermail/vtkusers/attachments/20150603/c596c698/attachment.html>