<html><body>
<p><tt>"Tilo Junge" <pult34@gmx.de> wrote on 27/07/2007 20:32:48:<br>
<br>
> Hello.<br>
> <br>
> I tried to cut a cylinder out of a cube, but it does not work right.<br>
> I use python and my code looks like this:<br>
> <br>
> import vtkpython<br>
> <br>
> cylinder = vtkpython.vtkCylinder()<br>
> cylinder.SetRadius(0.25)<br>
> cylinder.SetCenter(0.0, 0.0, 0.0)<br>
> vertPlane = vtkpython.vtkPlane()<br>
> vertPlane.SetOrigin(0, -2, 0)<br>
> vertPlane.SetNormal(0, -1, 0)<br>
> basePlane = vtkpython.vtkPlane()<br>
> basePlane.SetOrigin(0, 2, 0)<br>
> basePlane.SetNormal(0, 1, 0)<br>
> <br>
> cube = vtkpython.vtkCubeSource()<br>
> cube.SetXLength(2.0)<br>
> cube.SetYLength(2.0)<br>
> cube.SetZLength(2.0)<br>
> cube.SetCenter(0.0, 0.0, 0.0)<br>
> <br>
> theCylinder = vtkpython.vtkImplicitBoolean()<br>
> theCylinder.SetOperationTypeToIntersection()<br>
> theCylinder.AddFunction(cylinder)<br>
> theCylinder.AddFunction(vertPlane)<br>
> theCylinder.AddFunction(basePlane)<br>
> <br>
> clipper = vtkpython.vtkClipPolyData()<br>
> clipper.SetInput(cube.GetOutput())<br>
> clipper.SetClipFunction(theCylinder)<br>
> clipper.InsideOutOff()<br>
> <br>
> mapper = vtkpython.vtkPolyDataMapper()<br>
> mapper.SetInput(clipper.GetOutput())<br>
> <br>
> actor = vtkpython.vtkActor()<br>
> actor.SetMapper(mapper)<br>
> actor.GetProperty().SetColor(0, 1, 0)<br>
> <br>
> ren1= vtkpython.vtkRenderer()<br>
> ren1.AddActor(actor)<br>
> ren1.SetBackground(1, 1, 1)<br>
> <br>
> renWin = vtkpython.vtkRenderWindow()<br>
> renWin.AddRenderer( ren1 )<br>
> renWin.SetSize( 500, 500 )<br>
> iren = vtkpython.vtkRenderWindowInteractor()<br>
> iren.SetRenderWindow(renWin)<br>
> <br>
> iren.Initialize()<br>
> renWin.Render()<br>
> iren.Start()<br>
> <br>
> I also tried to use vtkCutter, but the result isnt what I wanted either.<br>
> Maybe someone can help me?</tt><br>
<br>
<tt>vtkCubeSource probably generates just 6 polygons (I've never actually used it), which is exactly what it should do. But using this as the input to clipping will not give what you expect. Clipping removes any cells completely inside the implicit surface and it clips any cells that have some points inside the surface and some points outside the surface. The points defining your cube are ALL outside the implicit surface so nothing is clipped! </tt><br>
<br>
<tt>What you need are cube surfaces that are made of many cells, some inside, some outside and some crossing the implicit surface, then you will see clipping. With say just 4 cells per cube surface you will get a very ugly square hole. As the number of cells (and points) increases you will get a better approximation to the circular hole you are expecting.</tt><br>
<br>
<tt>To get more cells/points on the cube surfaces you could use: </tt><br>
<tt>vtkPlaneSource (one per cube surface, thats 6 of them) and use SetResolution to increase the number of cells.</tt><br>
<tt>vtkLinearSubdivisionFilter to subdivide the cells in your vtkCubeSource, try 2 or 3 iterations for a start.</tt><br>
<br>
<tt> how do I know this.... yep I have been exactly where you are</tt><br>
<br>
<tt>good luck, Dave P</tt><br>
<tt><br>
> <br>
> Thanks<br>
> -- <br>
> Psssst! Schon vom neuen GMX MultiMessenger gehört?<br>
> Der kanns mit allen: <a href="http://www.gmx.net/de/go/multimessenger">http://www.gmx.net/de/go/multimessenger</a><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><br>
</tt></body></html>