利用PlaneIntersector裁剪实现圆环

rubby

又向各位大牛请教一下：

我创建了两个圆柱面，半径不同，然后利用PlaneIntersector裁剪，是否可以将这两个圆柱面中间的圆环显示出来？
代码如下。现在的结果是： 此圆所在的面获取了，但是圆环没有，见附图，请各位赐教。 或者有别的方法实现吗吗？？
比如：将两个面都截取出来，然后通过求两个面的差集，osg貌似没有这种功能。
小弟实在想不出方法。跪谢。

1. int main( int argc, char** argv ){
   
2.         //1. 创建两个圆柱面
   
3.         float height = 5.f;
   
4.         osg::ref_ptr<osg::TessellationHints> hints;
   
5.         hints = new osg::TessellationHints;       
   
6.         hints->setDetailRatio(0.5f);//　控制精细度
   
7.         hints->setCreateBody(true);
   
8.         //hints->setCreateBackFace(true);
   
9.         hints->setCreateBottom(false);
   
10.         hints->setCreateTop(false);
    
11.         osg::ref_ptr<osg::Geode> mynode = new osg::Geode;
    
12.         osg::ref_ptr<osg::ShapeDrawable> drawable = new osg::ShapeDrawable(new osg::Cylinder(osg::Vec3(0.0f,0.0f,0.0f),1.0f,height),hints);
    
13.         osg::ref_ptr<osg::ShapeDrawable> drawable2 = new osg::ShapeDrawable(new osg::Cylinder(osg::Vec3(0.0f,0.0f,0.0f),4.0f,height),hints);;
    
14.        
    
15.         mynode->addDrawable(drawable.get());       
    
16.         mynode->addDrawable(drawable2.get());
    
17.        
    
18.         osg::ref_ptr<osg::Group> root = new osg::Group;
    
19.         root->addChild(mynode.get() );
    
20. 
    
21. 
    
22.         //2 PlaneIntersector 实现
    
23.         osg::Plane plane( osg::Vec3(0,0,1), osg::Vec3(0,0,0) ); //(norm point)        The normal of the plane. point        A point of the plane.
    
24.         osg::Polytope polytope;
    
25.         /*polytope.add( osg::Plane(osg::Vec3(1,0,0), osg::Vec3(-1,0,0)) );
    
26.         polytope.add( osg::Plane(osg::Vec3(-1,0,0), osg::Vec3(1,0,0)) );*/
    
27.         osg::ref_ptr<osgUtil::PlaneIntersector> intersector = new osgUtil::PlaneIntersector( plane, polytope );
    
28. 
    
29.         osgUtil::IntersectionVisitor iv( intersector.get() );
    
30.         mynode->accept( iv );
    
31. 
    
32.         if ( intersector->containsIntersections() )
    
33.         {
    
34.                 osgUtil::PlaneIntersector::Intersections
    
35.                         &intersections = intersector->getIntersections();
    
36.                 std::cout <<"intersection size "<<intersections.size()<<std::endl;
    
37. 
    
38.                 osg::Geode * geode=new osg::Geode;
    
39.                 osg::Geometry * geom=new osg::Geometry;
    
40.                 osg::Vec3Array* v3a=new osg::Vec3Array;
    
41.                 osg::ref_ptr<osg::Vec4Array> colors = new osg::Vec4Array;
    
42.                 colors->push_back( osg::Vec4(1.0, 0.0, 0.0, 1.0) );
    
43. 
    
44.                 osgUtil::PlaneIntersector::Intersections::iterator itr;
    
45.                 for ( itr = intersections.begin(); itr != intersections.end(); ++itr )
    
46.                 {
    
47.                         for ( unsigned int i=0; i<(*itr).polyline.size(); ++i )
    
48.                                 std::cout <<" X" << (*itr).polyline[i].x()
    
49.                                 << " Y" << (*itr).polyline[i].y()
    
50.                                 << " Z" <<(*itr).polyline[i].z() << std::endl;
    
51. 
    
52.                         //交点坐标
    
53.                         std::vector<osg::Vec3d> intersectPoints = itr->polyline;
    
54.                         std::cout<<        intersectPoints.size()<<std::endl;
    
55.                        
    
56.                         for (unsigned int i=0;i<intersectPoints.size();i++)
    
57.                         {
    
58.                                 v3a->push_back(intersectPoints.at(i));
    
59.                                
    
60.                         }       
    
61.                 }
    
62. 
    
63.                 geom->setVertexArray(v3a);
    
64.                 geom->setColorArray(colors.get() );
    
65.                 geom->setColorBinding(osg::Geometry::BIND_OVERALL);
    
66.                 geom->addPrimitiveSet(new osg::DrawArrays(osg::PrimitiveSet::POLYGON,0,v3a->size()));
    
67.                 //geom->addPrimitiveSet(new osg::DrawArrays(osg::PrimitiveSet::LINE_LOOP,0,v3a->size()));
    
68.                 geode->addDrawable(geom);
    
69.                 osg::BoundingSphere bs=geode->getBound();
    
70.                 std::cout<<"radius"<<bs.radius()<<std::endl;
    
71.                 root->addChild(geode);
    
72.         }
    
73. 
    
74.         osgViewer::Viewer viewer;
    
75.         viewer.setSceneData(root.get());
    
76.         viewer.addEventHandler(new osgViewer::WindowSizeHandler);
    
77.         return viewer.run();
    
78. 
    
79. }

复制代码