<div dir="ltr"><div>Hi Elvis,</div><div><br></div><div>I'd recommend generating a stencil with vtkPolyDataToImageStencil, then iterating over over all the voxels inside the stencil with vtkImagePointDataIterator:</div><div><br></div><div><a href="http://www.vtk.org/doc/nightly/html/classvtkPolyDataToImageStencil.html" target="_blank">http://www.vtk.org/doc/<wbr>nightly/html/<wbr>classvtkPolyDataToImageStencil<wbr>.html</a><br></div><div><a href="http://www.vtk.org/doc/nightly/html/classvtkImagePointDataIterator.html" target="_blank">http://www.vtk.org/doc/<wbr>nightly/html/<wbr>classvtkImagePointDataIterator<wbr>.html</a><br></div><div><br></div><div>The iterator can efficiently give you the pointId for every voxel enclosed by the polyhedral surface (works for concave as well as convex surfaces). Once you have the pointId, you can set that voxel to whatever you want. There's a brief and vaguely relevant description on gitlab: <a href="https://gitlab.kitware.com/vtk/vtk/merge_requests/1118" target="_blank">https://gitlab.<wbr>kitware.com/vtk/vtk/merge_<wbr>requests/1118</a></div><div><br></div><div> - David</div><br><div class="gmail_extra"><br><div class="gmail_quote">On Sun, Jul 23, 2017 at 6:45 PM, Elvis Stansvik <span dir="ltr"><<a href="mailto:elvis.stansvik@orexplore.com" target="_blank">elvis.stansvik@orexplore.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hi all,<br>
<br>
I have a set of points describing multiple convex polyhedra. The<br>
polyhedra are for the most part disjoint in space, but they may<br>
intersect slightly in some places.<br>
<br>
I need to create an empty voxel volume (all scalars 0), and then<br>
"fill" the voxels that are inside one (or several) of the polyhedrons<br>
with some Perlin noise. I also have access to normals for the faces<br>
constituting the polyhedra, if necessary.<br>
<br>
Finally, I also have a single separate polyhedron (non-convex this<br>
time), and as a finishing step, I'd like to fill the voxels that are<br>
inside this polyhedron with a constant scalar value.<br>
<br>
Can VTK help me here?<br>
<br>
As this is new territory for me as far as VTK goes, any<br>
pointers/advice on how to go about it would be much appreciated.<br>
<br>
Cheers,<br>
Elvis<br></blockquote></div></div></div>