Hi All,<br><br>I would like to create a vtkImageData from a 3D surface mesh.  The mesh is closed.  It represents a 3D model in terms of the model volumetric outline.  You can assume this mesh is smooth and triangular, but it can be concave.  I want to have the output in vtkImageData having 0 voxel intensity if the whole voxel is outside the model, 1 if the whole voxel is lay inside the model, otherwise a value between 0 and 1 indicating the percentage of voxel volume that is inside the model.<br>
<br>I know how to use vtkImageStencil to generate a binary vtkImageData from vtkPolyData, but this approach does not take into account the partial volume.<br><br>May vtkProbeFilter help?  It seems that I need to have a volumetric mesh in order to do so.<br>
<br>Is there a combination of VTK filters that I can use to generate such vtkImageData?  ITK filters or a combination of both are also fine.<br><br>I believe this is a very fundamental problem in the medical imaging community.  May someone share his/her brilliant idea in tackling this basic problem?<br>
<br>Best regards,<br>Wilbur<br>