[Insight-users] Problems with voxelization using VTK

[Luis Ibanez]

Hi Shi,

This is the mailing list of the Insight Toolkit ITK.

Your question refers to the Visualization Toolkit VTK.

Please join the VTK users list at:

http://public.kitware.com/mailman/listinfo/vtkusers

and post your question there.


   Thanks,


      Luis



-------------------
shi yuan wrote:
> Hi all,
> I’m a beginner to VTK, I have requirements as follows;
> 1)Rasterizing a polygonal mesh object to a model represents as voxel, 
> and the voxelized model is not “hollow”;
> 2)Getting infomation from the voxelized model such as the position of 
> voxel-cube (x,y,z);
> 
>> From this mail list, vtkImplicitModeller and vtkVoxelModeller seems 
>> useful 
> 
> to me, and I made a test program, but it doesn’t match my need, because 
> the program produce a triangular mesh model which seems “hollow”, I mean 
> the model should be “accumulated” by many voxel which has the same size. 
> I am not very familiar with the VTK, maybe I use wrong method. Besides 
> the problem mentioned above, I have some problems as follows:
> 
> 1)Whether VTK has ability to Voxelize/Rasterize a geometry or not? (The 
> voxelized model must be “accumulated” by the same unit voxel). In the 
> other words, should I apply certain voxelization algorithm to program? 
> Whether VTK only plays the role of visualization or not?
> 2)If VTK has ability to Voxelize a geometry, which class(es) should I 
> use? vtkImplicitModeller and vtkVoxelModeller doesn’t work well for my 
> need (code segments listed as attachment). And which method(s) should I 
> use for Getting infomation from the voxelized model.
> 
> Dose anyone give me some hints? I would be grateful for any help!
> 
> Best wishes,
> Shiyuan
> 
> -------------------------------------------------------------------------------------------- 
> 
> 
> Code Segments:
> vtkSmartPointer<vtkCubeSource> cube = 
> vtkSmartPointer<vtkCubeSource>::New();
> cube->SetCenter(0,0,0);
> cube->SetXLength(20);
> cube->SetYLength(20);
> cube->SetZLength(20);
> 
> vtkSmartPointer<vtkVoxelModeller> voxel = 
> vtkSmartPointer<vtkVoxelModeller>::New();
> voxel->SetInputConnection(cube->GetOutputPort());
> voxel->SetSampleDimensions(8,8,8);
> voxel->SetModelBounds(-10,10,-10,10,-10,10);
> 
> vtkSmartPointer<vtkContourFilter> contour = 
> vtkSmartPointer<vtkContourFilter>::New();
> contour->SetInputConnection(voxel->GetOutputPort());
> contour->SetValue(0,1);
> 
> vtkSmartPointer<vtkPolyDataMapper> voxelMapper = 
> vtkSmartPointer<vtkPolyDataMapper>::New();
> voxelMapper ->SetInputConnection(contour->GetOutputPort());
> 
> vtkSmartPointer<vtkActor> voxelActor = vtkSmartPointer<vtkActor>::New();
> voxelActor ->SetMapper(voxelMapper);
> voxelActor ->GetProperty()->SetColor(0,0,1);
> voxelActor ->GetProperty()->SetOpacity(0.5);
> voxelActor ->GetProperty()->SetRepresentationToWireframe();
> 
> _________________________________________________________________
> 享用世界上最大的电子邮件系统― MSN Hotmail。  http://www.hotmail.com 
> _______________________________________________
> Insight-users mailing list
> Insight-users at itk.org
> http://www.itk.org/mailman/listinfo/insight-users
>