LineSegmentIntersector 存在bug？

rubby

各位大哥，国庆节时候，发现LineSegmentIntersector 有问题：
代码如下： 程序运行后，本来应该有8条线，结果只有7条线，如附件中的图：

1. int main( int argc, char** argv ){
   
2. 
   
3.         float height = 5.f;
   
4. 
   
5.         osg::ref_ptr<osg::TessellationHints> hints;
   
6.         hints = new osg::TessellationHints;       
   
7.         hints->setDetailRatio(0.5f);//TODO　控制精细度
   
8.         hints->setCreateBody(true);
   
9.         //hints->setCreateBackFace(true); //必须去除这句话，否则贴图闪烁。
   
10.         hints->setCreateBottom(false);
    
11.         hints->setCreateTop(false);
    
12.         osg::ref_ptr<osg::Geode> mynode = new osg::Geode;
    
13.         osg::ref_ptr<osg::ShapeDrawable> drawable = new osg::ShapeDrawable(new osg::Cylinder(osg::Vec3(0.0f,0.0f,0.0f),1.0f,height),hints);
    
14.         osg::ref_ptr<osg::ShapeDrawable> drawable2 = new osg::ShapeDrawable(new osg::Cylinder(osg::Vec3(0.0f,0.0f,0.0f),4.0f,height),hints);
    
15. 
    
16. 
    
17.         mynode->addDrawable(drawable.get());        // geom.get() or gemo. both are ok.
    
18.         mynode->addDrawable(drawable2.get());
    
19. 
    
20. 
    
21. 
    
22.         osg::ref_ptr<osg::Group> root = new osg::Group;
    
23.         root->addChild(mynode.get() );
    
24. 
    
25. 
    
26.         float angle = 0.0;
    
27.         float angleDelta = 2.0f*osg::PI/8;       
    
28. 
    
29.         for(int i=0; i<8;i++){  
    
30. 
    
31.                 float c = cosf(angle);
    
32.                 float s = sinf(angle);
    
33.                 osg::ref_ptr<osgUtil::LineSegmentIntersector> intersector = new osgUtil::LineSegmentIntersector( osg::Vec3(0.0f,0.0f,0.0f),  osg::Vec3(10.0f*c,10.0f*s,0.0f) );
    
34. 
    
35.                 osgUtil::IntersectionVisitor iv( intersector.get() );
    
36.                 mynode->accept( iv );
    
37. 
    
38.                 if ( intersector->containsIntersections() )
    
39.                 {
    
40.                         osgUtil::LineSegmentIntersector::Intersections
    
41.                                 &intersections = intersector->getIntersections();
    
42.                         std::cout <<"intersection size "<<intersections.size()<<std::endl;
    
43. 
    
44.                         osg::ref_ptr<osg::Geode>  geode=new osg::Geode;
    
45.                         osg::ref_ptr<osg::Geometry> geom=new osg::Geometry;
    
46.                         osg::ref_ptr<osg::Vec3Array> v3a=new osg::Vec3Array;
    
47.                         osg::ref_ptr<osg::Vec4Array> colors = new osg::Vec4Array;
    
48.                         colors->push_back( osg::Vec4(1.0, 0.0, 0.0, 1.0) );
    
49. 
    
50. 
    
51. 
    
52.                         osgUtil::LineSegmentIntersector::Intersections::iterator itr;
    
53.                         for ( itr = intersections.begin(); itr != intersections.end(); ++itr ) //起始点：有四个交点，前两个交点一样，后两个交点一样（因为为了封闭）。其余点正常。
    
54.                         {
    
55.                                 std::cout <<" X" << itr->getWorldIntersectPoint().x()
    
56.                                         << " Y" <<  itr->getWorldIntersectPoint().y()
    
57.                                         << " Z" << itr->getWorldIntersectPoint().z();
    
58.                                 if(itr == intersections.begin())
    
59.                                         v3a->push_back( itr->getWorldIntersectPoint());
    
60.                                
    
61.                         }
    
62.                                 itr--;
    
63.                                 v3a->push_back( itr->getWorldIntersectPoint());
    
64.                         std::cout<<std::endl;
    
65.                         geom->setVertexArray(v3a);
    
66.                         geom->setColorArray(colors.get() );
    
67.                         geom->setColorBinding(osg::Geometry::BIND_OVERALL);
    
68.                         for(int k=0; k<v3a->size(); k++){
    
69.                                 std::cout<<"vectex: "<<v3a->at(k).x()<<" "<<v3a->at(k).y()<<" "<<v3a->at(k).z();
    
70.                         }
    
71.                         std::cout<<std::endl;
    
72.                         geom->addPrimitiveSet(new osg::DrawArrays(osg::PrimitiveSet::LINES,0,v3a->size()));
    
73.                         geode->addDrawable(geom);
    
74.                         root->addChild(geode);
    
75.                 }else{
    
76.                         std::cout<<"no intersection for i="<<i<<std::endl;
    
77.                         std::cout<<"x "<<intersector->getStart().x()<<" y "<<intersector->getStart().y()<<" z "<<intersector->getStart().z()<<std::endl;
    
78.                         std::cout<<"x "<<intersector->getEnd().x()<<" y "<<intersector->getEnd().y()<<" z "<<intersector->getEnd().z()<<std::endl;
    
79.                 }
    
80.                 angle+=angleDelta;
    
81.         }
    
82.         osgViewer::Viewer viewer;
    
83.         viewer.setSceneData(root.get());
    
84.         viewer.addEventHandler(new osgViewer::WindowSizeHandler);
    
85.         return viewer.run();
    
86. 
    
87. }

复制代码

rubby

此外，osg:rimitiveSet:INES绘制出来的线，弯弯曲曲，有方法设置，使之平滑吗？

多谢各位神仙大哥。

rubby

自己来结贴：
同样遇到此问题的兄弟可以如下解决：

intersector->setPrecisionHint(osgUtil::Intersector:recisionHint::USE_DOUBLE_CALCULATIONS);

注意：osg3.0.1 无此函数。