<div dir="ltr"><div dir="ltr"><div>If you already have the spheres as polydata, and if they are already in your volume coord system, then the simplest approach is probably to use these two filters:</div><div>- vtkCutter, to cut the sphere (this generates a polyline contour of the intersection)</div><div>- vtkContourTriangulator to fill the polyline contour with triangles</div><div><br></div><div>When rendering the result on the image, the polydata and your image slice will be coincident in the z-buffer.  It might be possible to fix this by calling SetResolveCoincidentTopologyToPolygonOffset() on the polydata mapper.</div><div><br></div><div>Rendering the filled spheres with voxels would be more complicated.  The vtkPolyDataToImageStencil class works volumetrically, i.e. it would be used to fill the entire sphere and produce a binary volume.  Then the desired slice of that volume would be rendered as an overlay.</div><div><br></div><div>   David</div><div><br></div><br><div class="gmail_quote"><div dir="ltr">On Fri, Jan 11, 2019 at 8:02 AM The Merper <<a href="mailto:msmerps@gmail.com">msmerps@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Thanks for the speedy reply, David.</div><div>The application is a GUI that allows users to mark points in a 3D MRI volume. The user browses the volume and marks the points in 2D slices of the volume, which are rendered via a vtkImageViewer2. Once a point is marked, I represent the location with a sphere. I use vtkSphereSource via vtkPolyDataMapper to render the spheres.</div><div><br></div><div>I'd be happy using polygons or voxels.</div><div>   thanks again,</div><div>      -Merps</div><div><br></div><div><div><br></div><div><br></div></div></div><br><div class="gmail_quote"><div dir="ltr">On Thu, Jan 10, 2019 at 5:08 PM David Gobbi <<a href="mailto:david.gobbi@gmail.com" target="_blank">david.gobbi@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div dir="ltr">On Thu, Jan 10, 2019 at 12:13 PM The Merper <<a href="mailto:msmerps@gmail.com" target="_blank">msmerps@gmail.com</a>> wrote:<br></div><div dir="ltr"><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div><br></div><div>  I have a set of spheres in a volume that I visualize in 2D slices. When the plane of a 2D slice intersects a sphere, I want the interior of the sphere to be filled.</div></div></div></div></div></div></blockquote><div><br></div><div>How are the spheres defined?  Do you already have them as vtkPolyData objects, or do you just have a list of "center and radius" values to describe the spheres?</div><div> </div><div>When you say "visualize in 2D slices" do you mean that you are already visualizing slices of something (e.g. of an image) and you want to display the sphere cuts as an overlay on the image?  For your application, does it matter whether the sphere cuts are displayed as voxels (i.e. a vtkImageData) vs polygons?</div><div><br></div><div>Also, for context, what is the general purpose of your application?</div><div><br></div><div>Cheers,</div><div>  David</div><div><br></div><div><br></div></div></div>
</div>
</blockquote></div>
</blockquote></div></div></div>