<html><body>
<p><tt>"Steve Chall" <stevec@renci.org> wrote on 09/10/2007 06:35:29:<br>
<br>
> David,</tt><br>
<tt>> I just ran into a problem with your approach that I’m hoping you can<br>
> help me with. Each of the spheres needs to have a different color <br>
> assigned to it. If I add each vtkSphereSource to a <br>
> vtkAppendPolyData object, is there any way I maintain the <br>
> association between each sphere and its assigned color? The only way<br>
> I know of to do that is with vtkActor::GetProperty()::SetColor(…) <br>
> and with the vtkAppendPolyData approach, the whole vtkAppendPolyData<br>
> object is associated with a single actor, which means it gets just <br>
> one color assignment. Is that correct? If so, can you suggest a <br>
> workaround? I can still display all the spheres individually with <br>
> their colors, I think,, but I can’t figure out how to keep the color<br>
> associations in the 2D slice. Thanks.</tt><br>
<tt>> </tt><br>
<tt>> -Steve</tt><br>
<tt>> </tt><br>
<br>
<tt>Hi Steve, this is getting a bit outside my experience, ie I always use vtk to represent point scalar data and have never had to 'manually' assign colours. Typical work with vtk is to use point (or sometimes cell) scalar data to control colour or to assign actor colours, you are in between these cases. </tt><br>
<br>
<tt>You could create a data set of just points with scalars and use glyphs (vtkGlyph3D) scaled and coloured according to scalar values, then cut them... But what are your ultimate goals here? there may be a way to do what you are asking now: spheres of differing sizes and colours but that might be a dead end when you move to the next step (you mentioned more complex shapes for particles). </tt><br>
<br>
<tt>What kinds of particle shapes are you considering? How do you want colours to be assigned? Any other requirements? Are shapes, sizes and colours all independent? </tt><a href="http://www.vtk.org/doc/release/5.0/html/c2_vtk_t_9.html#c2_vtk_t_vtkProgrammableGlyphFilter"><u><font size="4" color="#0000FF">vtkProgrammableGlyphFilter</font></u></a><font size="4"> </font><tt> could be worth a look (I've never used it) but it sounds like you could create points with one or more attribute values associated, then use those attributes to generate glyphs for your particles with desired features: size, colour... You might need to find a way to associate (new?) scalars with each glyph so when it is cut it will have the desired colour....</tt><br>
<br>
<tt> regards, Dave P</tt><br>
<br>
<tt>> <br>
> From: Steve Chall [<a href="mailto:steve.chall@earthlink.net">mailto:steve.chall@earthlink.net</a>] <br>
> Sent: Monday, October 08, 2007 12:12 AM<br>
> To: David.Pont@ensisjv.com; Steve Chall<br>
> Cc: 'Theresa-Marie Rhyne'; vtkusers@vtk.org<br>
> Subject: Re: [vtkusers] Displaying the intersection of a plane with <br>
> a bunch of spheres</tt><br>
<tt>> </tt><br>
<tt>> Thanks, David, This looks eminently doable and makes plenty of <br>
> sense. I'm going to try it today.<br>
> <br>
> -Steve</tt><br>
<tt>> -----Original Message----- <br>
> From: David.Pont@ensisjv.com <br>
> Sent: Oct 7, 2007 4:32 PM <br>
> To: Steve Chall <br>
> Cc: 'Theresa-Marie Rhyne' , vtkusers@vtk.org <br>
> Subject: Re: [vtkusers] Displaying the intersection of a plane with <br>
> a bunch of spheres <br>
> <br>
</tt><br>
<tt>> "Steve Chall" <stevec@renci.org> wrote on 06/10/2007 07:35:07:<br>
> <br>
> > I want to generate a 3D rectangular volume of ~38,000 closely packed,<br>
> > sometimes overlapping "particles" - initially spheres simply defined as a<br>
> > list of centroids and a parallel list of radii, more complex objects later<br>
> > on - and interactively look at arbitrary slices through the volume, with the<br>
> > particle intersections represented as disks or circles in the 2D slice (down<br>
> > to a single pixel, of course, if the slice is tangent to the surface of a<br>
> > "particle"). Due to the number of "particles" I think performance scaling<br>
> > issues might also be a concern at some point.<br>
> > <br>
> > I’ve been working with VTK now for some time in pursuit of a solution and<br>
> > haven’t really come to a satisfactory conclusion. I’ve tried the following<br>
> > two approaches thus far:<br>
> > <br>
> > 1) Implicit booleans: I read in the sphere positions and radii and<br>
> > construct a vtkSphere for each one. As I do that I add each to a<br>
> > vtkImplicitBoolean object “sphereUnion” with the operation type set to<br>
> > union. I then create a vtkPlane object (right now I’m not worrying about<br>
> > interactively moving the plane around: I’ve just got it set through the<br>
> > middle of the cluster of spheres so that I know, because I can see it, that<br>
> > there’s a non-empty intersection between the plane and some of the<br>
> > spheres). Then I create another vtkImplicitBoolean called “slice” and add<br>
> > to it the sphereUnion and the plane with operation type set to<br>
> > intersection. I run slice into a vtkSampleFunction, that into a<br>
> > vtkContourFilter into a vtkPolyDataMapper into a vtkActor and into a<br>
> > vtkRenderer (the standard textbook sequence). I’ve tried looking at various<br>
> > subsets – just the sphereUnion, just the plane, etc. – and I see what I<br>
> > expect. But I’d think that the intersection of a bunch of spheres and a<br>
> > plane would look like calamari or sausage slices: a cluster of rings or<br>
> > disks on a 2D surface. However, what I’m seeing is a strange cluster of<br>
> > lumpy, irregular 3D objects scattered through the 3D sample space that<br>
> > perhaps roughly correspond to some or all of the original spheres. Thus far<br>
> > I haven't made any sense of it. So on to the second approach:<br>
> > <br>
> > 2) Polygons: I create the vtkSpheres, but this time convert each into a<br>
> > vtkPolyData object (vtkSphere -> vtkSampleFunction -> vtkContourFilter) and<br>
> > add the resulting polygonal sphere to a vtkAppendPolyData object. I create<br>
> > a vtkPlaneSource object (the polygonal equivalent to the vtkPlane implicit<br>
> > surface object in approach 1 above) and a vtkProbeFilter. I set the probe’s<br>
> > input connection to the vtkPlaneSource and its source connection to the<br>
> > vtkAppendPolydata object (the spheres) and then from there into the usual<br>
> > pipeline to the renderer. Once again I think the subunits are ok – I can<br>
> > display them and they look like I’d expect – but I don’t see any probe’d<br>
> > data at all. <br>
> > <br>
> > Does either of these approaches look promising, or is there an entirely<br>
> > different direction that’s more likely to produce useful results? I’d<br>
> > appreciate any insights you might have. I can provide source code if it<br>
> > would make it any clearer. Thanks for any help you can provide.<br>
> <br>
> You should try this: create your spheres as polydata objects with <br>
> vtkSphereSource (instead of implicit spheres with vtkSphere) and <br>
> merge into one data set with vtkAppendPolyData. Create an implicit <br>
> plane with vtkPlane (as before) and use vtkCutter to generate the 2D<br>
> slice through your polydata spheres. You should get 'circular' <br>
> polylines (and/or points), increase the resolution in <br>
> vtkSphereSource to get better approximation to circles (better calamari ;-).<br>
> <br>
> You can also cut more complex shapes this way, so long as you can <br>
> represent them as polydata. Performance with this approach should be<br>
> OK, unlike an approach based on vtkSampleFunction -> <br>
> vtkContourFilter which typically requires high sampling density to <br>
> get acceptable resolution. 38K particles is quite a few so vtkCutter<br>
> will be busy but using vtkSampleFunction on this problem would be an<br>
> order of magnitude (or 2) slower.<br>
> If you really need to specify the particles as implicit functions <br>
> you would use vtkSampleFunction -> vtkContour to turn them into <br>
> polydata and then apply vtkCutter on the polydata, that way if you <br>
> want to move the cutting plane around you will not have to re-<br>
> execute vtkSampleFunction. <br>
> <br>
> HTH, Dave P<br>
> <br>
> > <br>
> > Steve<br>
> > <br>
> > -Steve Chall<br>
> > Senior Research Software Developer<br>
> > Renaissance Computing Institute<br>
> > Phone: 919-515-0051<br>
> > Email: stevec@renci.org<br>
> > <br>
> > <br>
> > <br>
> > _______________________________________________<br>
> > This is the private VTK discussion list. <br>
> > Please keep messages on-topic. Check the FAQ at: <a href="http://www.vtk">http://www.vtk</a>.<br>
> > org/Wiki/VTK_FAQ<br>
> > Follow this link to subscribe/unsubscribe:<br>
> > <a href="http://www.vtk.org/mailman/listinfo/vtkusers">http://www.vtk.org/mailman/listinfo/vtkusers</a></tt></body></html>