使用osg::Geometry类如何绘制圆柱体

Sailent

如题，因为shape类画出来的没办法进行一些复杂操作。我想出来的代码是这样的：

1. osg::ref_ptr<osg::Node> CreateCylinder()
   
2. {
   
3.         // Create an object to store geometry in.
   
4.         osg::ref_ptr<osg::Geometry> geom = new osg::Geometry;
   
5.         // Create an array of four vertices.
   
6.         osg::ref_ptr<osg::Vec3Array> v = new osg::Vec3Array;
   
7.         geom->setVertexArray( v.get() );
   
8. 
   
9.         //------------------关键代码，设置圆柱体顶点------------------
   
10.         float radius=2.0;//半径
    
11.         float fHeight=4.0;//高
    
12.         //上底
    
13.         for (float a=-1.0*radius;a<=1.0*radius+0.0001;a+=radius/25.0)
    
14.         {
    
15.                 v->push_back(osg::Vec3(a,0.0f,sqrt(fabs(pow(radius,2)-pow(a,2)))));
    
16.         }
    
17.         for (float a=1.0*radius;a>=-1.0*radius-0.0001;a-=radius/25.0)
    
18.         {
    
19.                 v->push_back(osg::Vec3(a,0.0f,-1*sqrt(fabs(pow(radius,2)-pow(a,2)))));
    
20.         }
    
21.         //下底
    
22.         for (float a=-1.0*radius;a<=1.0*radius+0.0001;a+=radius/25.0)
    
23.         {
    
24.                 v->push_back(osg::Vec3(a,0.0f,sqrt(fabs(pow(radius,2)-pow(a,2))))+osg::Vec3(0.0,fHeight,0.0));
    
25.         }
    
26.         for (float a=1.0*radius;a>=-1.0*radius-0.0001;a-=radius/25.0)
    
27.         {
    
28.                 v->push_back(osg::Vec3(a,0.0f,-1*sqrt(fabs(pow(radius,2)-pow(a,2))))+osg::Vec3(0.0,fHeight,0.0));
    
29.         }
    
30.         float fPoint=204;//总共使用的顶点数
    
31.         //----------------------------------------------------------
    
32. 
    
33. 
    
34.         // Create an array of four colors.
    
35.         osg::ref_ptr<osg::Vec4Array> c = new osg::Vec4Array;
    
36.         c->push_back( osg::Vec4( 1.f, 1.f, 1.f, 1.f ) );
    
37.         geom->setColorArray( c.get() );
    
38.         geom->setColorBinding( osg::Geometry::BIND_OVERALL );
    
39.         // Draw a four-vertex quad from the stored data.
    
40.         geom->addPrimitiveSet(new osg::DrawArrays( osg::PrimitiveSet::POLYGON, 0,fPoint) );
    
41.         // Add the Geometry (Drawable) to a Geode and
    
42.         // return the Geode.
    
43.         osg::ref_ptr<osg::Geode> geode = new osg::Geode;
    
44.         geode->addDrawable( geom.get() );
    
45. 
    
46.         return geode.get();
    
47. }

复制代码

很可惜，运行之后圆柱体始终出不来，不知各位有何良策？

[ 本帖最后由 Sailent 于 2008-9-18 17:45 编辑 ]

array

用osg:: PrimitiveSet:: POLYGON有点问题，应该用TRIANGLES或者QUADS，并且设置了正确的顶点索引才对。此外我没有细看，估计还要改一些东西

我以前做过一个Lathe库，用来生成各种旋转体的，算法什么的都不是很复杂，不过现在找不到了~~

tan_dunming

OSG官方网站上下载过一个LMBs_OSG_Examples.zip，里面有一个12_CustomDrawable.cpp，可以算是OSG与OPENGL的混合编程吧，如果OPENGL功底还可以的话可以考虑采用这种办法实现。

FlySky

估计图元方式，，，，索引有问题吧～～～你可以看看shape.cpp以及shapedrawable.cpp中的源代码，，里面有相关几何体的绘制～～～～～～～～～～～～

Sailent

原帖由 FlySky 于 2008-9-18 22:51 发表
估计图元方式，，，，索引有问题吧～～～你可以看看shape.cpp以及shapedrawable.cpp中的源代码，，里面有相关几何体的绘制～～～～～～～～～～～～

我看过，shape里面是用glBegin、glEnd实现的，参考价值不大吧……

FlySky

原帖由 Sailent 于 2008-9-19 09:26 发表

我看过，shape里面是用glBegin、glEnd实现的，参考价值不大吧……

OSG本来就是可以和OpenGL结合来绘制～～～～～那样绘制，，，你同样可以设置纹理等各种渲染～～～～～

Sailent

终于画出来了，不过我明明设置材质的漫反射和镜面发射为什么不起作用呢？要怎么设置光照呢？

1. osg::ref_ptr<osg::Node> CreateCylinder()
   
2. {
   
3.         // Create an object to store geometry in.
   
4.         osg::ref_ptr<osg::Geometry> geom = new osg::Geometry;
   
5.         // Create an array of four vertices.
   
6.         osg::ref_ptr<osg::Vec3Array> v = new osg::Vec3Array;
   
7.         geom->setVertexArray( v.get() );
   
8. 
   
9.         //------------------关键代码，设置圆柱体顶点------------------
   
10.         float radius=2.0;//半径
    
11.         float radius2=1.0;
    
12.         float fHeight=4.0;//高
    
13.         //上底
    
14.         for (float a=-1.0*radius;a<=1.0*radius+0.0001;a+=radius/25.0)
    
15.         {
    
16.                 v->push_back(osg::Vec3(a,0.0f,sqrt(fabs(pow(radius,2)-pow(a,2)))));
    
17.                 v->push_back(osg::Vec3(a*(radius2/radius),0.0f,sqrt(fabs(pow(radius2,2)-pow(a*(radius2/radius),2))))+osg::Vec3(0.0,fHeight,0.0));
    
18.         }
    
19.         for (float a=1.0*radius;a>=-1.0*radius-0.0001;a-=radius/25.0)
    
20.         {
    
21.                 v->push_back(osg::Vec3(a,0.0f,-1*sqrt(fabs(pow(radius,2)-pow(a,2)))));
    
22.                 v->push_back(osg::Vec3(a*(radius2/radius),0.0f,-1*sqrt(fabs(pow(radius2,2)-pow(a*(radius2/radius),2))))+osg::Vec3(0.0,fHeight,0.0));
    
23.         }
    
24. 
    
25.         float fPoint=204;//总共使用的顶点数
    
26.         //----------------------------------------------------------
    
27. 
    
28. 
    
29.         // Create an array of four colors.
    
30.         osg::ref_ptr<osg::Vec4Array> c = new osg::Vec4Array;
    
31.         c->push_back( osg::Vec4( 1.f, 1.f, 1.f, 1.f ) );
    
32.         geom->setColorArray( c.get() );
    
33.         geom->setColorBinding( osg::Geometry::BIND_OVERALL );
    
34.         // Draw a four-vertex quad from the stored data.
    
35.         geom->addPrimitiveSet(new osg::DrawArrays( osg::PrimitiveSet::TRIANGLE_STRIP, 0,fPoint) );
    
36.         // Add the Geometry (Drawable) to a Geode and
    
37.         // return the Geode.
    
38.         osg::ref_ptr<osg::Geode> geode = new osg::Geode;
    
39.         geode->addDrawable( geom.get() );
    
40. 
    
41.         osg::Matrix RoLast1;
    
42.         RoLast1.makeRotate(180.0*3.14159265/180.0,osg::Vec3f(1.0,0.0,0.0));
    
43. 
    
44.         osg::Matrix RoLast2;
    
45.         RoLast2.makeTranslate(osg::Vec3(0.0,2.75,0.0));
    
46. 
    
47.         osg::ref_ptr<osg::MatrixTransform> TransNode = new osg::MatrixTransform;
    
48.         TransNode->addChild(geode.get());
    
49.         TransNode->setMatrix(RoLast1*RoLast2);
    
50. 
    
51.         osg::StateSet* state4 = TransNode->getOrCreateStateSet();
    
52.         osg::ref_ptr<osg::Material> mat3 = new osg::Material;
    
53.         //漫反射光（决定物体颜色）
    
54.         mat3->setDiffuse( osg::Material::FRONT_AND_BACK,
    
55.                 osg::Vec4( 1.f, 0.f, 0.f, 1.f ) );
    
56.         //镜面反射光
    
57.         mat3->setSpecular( osg::Material::FRONT_AND_BACK,
    
58.                 osg::Vec4( 1.f, 1.f, 1.f, 1.f ) );
    
59.         //环境反射光
    
60.         mat3->setShininess( osg::Material::FRONT_AND_BACK, 128.f );
    
61.         state4->setAttribute( mat3.get() );
    
62.         state4->setRenderingHint( osg::StateSet::OPAQUE_BIN );
    
63. 
    
64.         return TransNode.get();
    
65. }

复制代码

[ 本帖最后由 Sailent 于 2008-9-19 15:51 编辑 ]

array

没有设置法线吧？setNormalArray

Sailent

可是应该往哪个方向设呢？要设多少条呢？要设成圆形的一圈吗？

Sailent

已经加了法线了，怎么它还是白色的呢？

1. osg::ref_ptr<osg::Node> CreateCylinder()
   
2. {
   
3.         // Create an object to store geometry in.
   
4.         osg::ref_ptr<osg::Geometry> geom = new osg::Geometry;
   
5.         // Create an array of four vertices.
   
6.         osg::ref_ptr<osg::Vec3Array> v = new osg::Vec3Array;
   
7.         geom->setVertexArray( v.get() );
   
8. 
   
9.         //------------------关键代码，设置圆柱体顶点------------------
   
10.         float radius=2.0;//半径
    
11.         float radius2=1.0;
    
12.         float fHeight=4.0;//高
    
13.         //上底
    
14.         for (float a=-1.0*radius;a<=1.0*radius+0.0001;a+=radius/25.0)
    
15.         {
    
16.                 v->push_back(osg::Vec3(a,0.0f,sqrt(fabs(pow(radius,2)-pow(a,2)))));
    
17.                 v->push_back(osg::Vec3(a*(radius2/radius),0.0f,sqrt(fabs(pow(radius2,2)-pow(a*(radius2/radius),2))))+osg::Vec3(0.0,fHeight,0.0));
    
18.         }
    
19.         for (float a=1.0*radius;a>=-1.0*radius-0.0001;a-=radius/25.0)
    
20.         {
    
21.                 v->push_back(osg::Vec3(a,0.0f,-1*sqrt(fabs(pow(radius,2)-pow(a,2)))));
    
22.                 v->push_back(osg::Vec3(a*(radius2/radius),0.0f,-1*sqrt(fabs(pow(radius2,2)-pow(a*(radius2/radius),2))))+osg::Vec3(0.0,fHeight,0.0));
    
23.         }
    
24. 
    
25.         float fPoint=204;//总共使用的顶点数       
    
26.        
    
27.         //设置法线
    
28.         osg::ref_ptr<osg::Vec3Array> n = new osg::Vec3Array;
    
29.         geom->setNormalArray( n.get() );
    
30.         geom->setNormalBinding( osg::Geometry::BIND_OVERALL );
    
31.         n->push_back( osg::Vec3( 0.0f,1.0f, 0.0f ) );
    
32.         //----------------------------------------------------------
    
33. 
    
34.         // Create an array of four colors.
    
35.         osg::ref_ptr<osg::Vec4Array> c = new osg::Vec4Array;
    
36.         c->push_back( osg::Vec4( 1.f, 1.f, 1.f, 1.f ) );
    
37.         geom->setColorArray( c.get() );
    
38.         geom->setColorBinding( osg::Geometry::BIND_OVERALL );
    
39.         // Draw a four-vertex quad from the stored data.
    
40.         geom->addPrimitiveSet(new osg::DrawArrays( osg::PrimitiveSet::TRIANGLE_STRIP, 0,fPoint) );
    
41.         // Add the Geometry (Drawable) to a Geode and
    
42.         // return the Geode.
    
43.         osg::ref_ptr<osg::Geode> geode = new osg::Geode;
    
44.         geode->addDrawable( geom.get() );
    
45. 
    
46.         osg::Matrix RoLast1;
    
47.         RoLast1.makeRotate(180.0*3.14159265/180.0,osg::Vec3f(1.0,0.0,0.0));
    
48. 
    
49.         osg::Matrix RoLast2;
    
50.         RoLast2.makeTranslate(osg::Vec3(0.0,2.75,0.0));
    
51. 
    
52.         osg::ref_ptr<osg::MatrixTransform> TransNode = new osg::MatrixTransform;
    
53.         TransNode->addChild(geode.get());
    
54.         TransNode->setMatrix(RoLast1*RoLast2);
    
55. 
    
56.         osg::StateSet* state4 = TransNode->getOrCreateStateSet();
    
57.         osg::ref_ptr<osg::Material> mat3 = new osg::Material;
    
58.         //漫反射光（决定物体颜色）
    
59.         mat3->setDiffuse( osg::Material::FRONT_AND_BACK,
    
60.                 osg::Vec4( 1.f, 0.f, 0.f, 1.f ) );
    
61.         //镜面反射光
    
62.         mat3->setSpecular( osg::Material::FRONT_AND_BACK,
    
63.                 osg::Vec4( 1.f, 1.f, 1.f, 1.f ) );
    
64.         //环境反射光
    
65.         mat3->setShininess( osg::Material::FRONT_AND_BACK, 128.f );
    
66.         state4->setAttribute( mat3.get() );
    
67.         state4->setRenderingHint( osg::StateSet::OPAQUE_BIN );
    
68. 
    
69.         return TransNode.get();
    
70. }

复制代码

array

要得到正确效果的话，法线的数量应该等同于顶点的数量，即每个顶点都要设置法线。

觉得麻烦的话也可以选用SmoothVisitor，不过它会自动生成一种模拟表面光滑效果的法线坐标，比较有意思。

Sailent

1. osg::ref_ptr<osg::Node> CreateCylinder()
   
2. {
   
3.         // Create an object to store geometry in.
   
4.         osg::ref_ptr<osg::Geometry> geom = new osg::Geometry;
   
5.         // Create an array of four vertices.
   
6.         osg::ref_ptr<osg::Vec3Array> v = new osg::Vec3Array;
   
7.         geom->setVertexArray( v.get() );
   
8. 
   
9.         //------------------关键代码，设置圆柱体顶点------------------
   
10.         float radius=2.0;//半径
    
11.         float radius2=1.0;
    
12.         float fHeight=4.0;//高
    
13.         //上底
    
14.         for (float a=-1.0*radius;a<=1.0*radius+0.0001;a+=radius/25.0)
    
15.         {
    
16.                 v->push_back(osg::Vec3(a,0.0f,sqrt(fabs(pow(radius,2)-pow(a,2)))));
    
17.                 v->push_back(osg::Vec3(a*(radius2/radius),0.0f,sqrt(fabs(pow(radius2,2)-pow(a*(radius2/radius),2))))+osg::Vec3(0.0,fHeight,0.0));
    
18.         }
    
19.         for (float a=1.0*radius;a>=-1.0*radius-0.0001;a-=radius/25.0)
    
20.         {
    
21.                 v->push_back(osg::Vec3(a,0.0f,-1*sqrt(fabs(pow(radius,2)-pow(a,2)))));
    
22.                 v->push_back(osg::Vec3(a*(radius2/radius),0.0f,-1*sqrt(fabs(pow(radius2,2)-pow(a*(radius2/radius),2))))+osg::Vec3(0.0,fHeight,0.0));
    
23.         }
    
24. 
    
25.         float fPoint=204;//总共使用的顶点数       
    
26.        
    
27.         //设置法线
    
28.         osg::ref_ptr<osg::Vec3Array> n = new osg::Vec3Array;
    
29.         geom->setNormalArray( n.get() );
    
30.         geom->setNormalBinding( osg::Geometry::BIND_PER_VERTEX );
    
31. //        n->push_back( osg::Vec3( 0.0f,1.0f, 0.0f ) );
    
32.         float x2=radius*fHeight*fHeight/(fHeight*fHeight+pow((radius-radius2),2));
    
33.         float x3=x2*(radius-radius2)/fHeight;
    
34.         for (float a=-1.0*radius;a<=1.0*radius+0.0001;a+=radius/25.0)
    
35.         {
    
36.                 n->push_back(osg::Vec3(a*(x2/radius),0.0f,sqrt(fabs(pow(x2,2)-pow(a*(x2/radius),2))))+osg::Vec3(0.0,x3,0.0));
    
37.         }
    
38.         for (float a=1.0*radius;a>=-1.0*radius-0.0001;a-=radius/25.0)
    
39.         {
    
40.                 n->push_back(osg::Vec3(a*(x2/radius),0.0f,-1*sqrt(fabs(pow(x2,2)-pow(a*(x2/radius),2))))+osg::Vec3(0.0,x3,0.0));
    
41.         }
    
42. 
    
43.        
    
44.         //----------------------------------------------------------
    
45. 
    
46.         // Create an array of four colors.
    
47.         osg::ref_ptr<osg::Vec4Array> c = new osg::Vec4Array;
    
48.         c->push_back( osg::Vec4( 1.f, 1.f, 1.f, 1.f ) );
    
49.         geom->setColorArray( c.get() );
    
50.         geom->setColorBinding( osg::Geometry::BIND_OVERALL );
    
51.         // Draw a four-vertex quad from the stored data.
    
52.         geom->addPrimitiveSet(new osg::DrawArrays( osg::PrimitiveSet::TRIANGLE_STRIP, 0,fPoint) );
    
53. 
    
54.         //使用smoothingVisitor
    
55.         osgUtil::SmoothingVisitor::smooth(*geom.get());
    
56. 
    
57.         // Add the Geometry (Drawable) to a Geode and
    
58.         // return the Geode.
    
59.         osg::ref_ptr<osg::Geode> geode = new osg::Geode;
    
60.         geode->addDrawable(geom.get());
    
61. 
    
62.         osg::Matrix RoLast1;
    
63.         RoLast1.makeRotate(180.0*3.14159265/180.0,osg::Vec3f(1.0,0.0,0.0));
    
64. 
    
65.         osg::Matrix RoLast2;
    
66.         RoLast2.makeTranslate(osg::Vec3(0.0,2.75,0.0));
    
67. 
    
68.         osg::ref_ptr<osg::MatrixTransform> TransNode = new osg::MatrixTransform;
    
69.         TransNode->addChild(geode.get());
    
70.         TransNode->setMatrix(RoLast1*RoLast2);
    
71. 
    
72.         osg::StateSet* state4 = TransNode->getOrCreateStateSet();
    
73.         osg::ref_ptr<osg::Material> mat3 = new osg::Material;
    
74.         //漫反射光（决定物体颜色）
    
75.         mat3->setDiffuse( osg::Material::FRONT_AND_BACK,
    
76.                 osg::Vec4( 1.f, 0.f, 0.f, 1.f ) );
    
77.         //镜面反射光
    
78.         mat3->setSpecular( osg::Material::FRONT_AND_BACK,
    
79.                 osg::Vec4( 1.f, 1.f, 1.f, 1.f ) );
    
80.         //环境反射光
    
81.         mat3->setShininess( osg::Material::FRONT_AND_BACK, 128.f );
    
82.         state4->setAttribute( mat3.get() );
    
83.         state4->setRenderingHint( osg::StateSet::OPAQUE_BIN );
    
84. 
    
85.         return TransNode.get();
    
86. }

复制代码

是这样吗？还是没有效果啊~~

array

您的法线作了归一化处理了吗？还有您的数学运算我不知道如何实现的，我印象中应该没有这么多sqrt和pow才对~~

下个月我打算启动一个新的osgModeling工程的计划，把所有3dsmax等软件中常用的基本体都用Geometry的派生类来实现，并计划编写Extrude（挤压），Lathe（旋转体），Loft（放样），BeizierSurface和NurbsSurface这几种常用建模手段类，远期还希望包括Subdivision（细化），BoolOperator（布尔）等高级建模功能。希望到时候能为大家提供一些方便。

Sailent

原帖由 array 于 2008-9-22 10:14 发表
您的法线作了归一化处理了吗？还有您的数学运算我不知道如何实现的，我印象中应该没有这么多sqrt和pow才对~~

下个月我打算启动一个新的osgModeling工程的计划，把所有3dsmax等软件中常用的基本体都用Geometry的派 ...

归一化处理？代码就这么多啊，最多加了下面一点：

1. osgViewer::Viewer viewer;
   
2. viewer.SetSceneData(CreateCyliner().get());
   
3. viewer.run();

复制代码

而且好像用了smoothingVisitor之后，自己写的法线设置就应该不起作用了，但运行的时候外面是黑的，里面是亮的，怎么反过来呢？

Sailent

我的数学运算其实就是实现：x*x+y*y=r*r这样圆的函数啊，设置每个顶点为下图所示。X就是循环控制的精度了。

[ 本帖最后由 Sailent 于 2008-9-22 10:55 编辑 ]

Sailent

将法线算法做了一些改进，取消了smoothingvisitor，但为什么会有2根黑黑的线呢？

1. osg::ref_ptr<osg::Node> CreateCylinder()
   
2. {
   
3.         // Create an object to store geometry in.
   
4.         osg::ref_ptr<osg::Geometry> geom = new osg::Geometry;
   
5.         // Create an array of four vertices.
   
6.         osg::ref_ptr<osg::Vec3Array> v = new osg::Vec3Array;
   
7.         geom->setVertexArray( v.get() );
   
8. 
   
9.         //------------------关键代码，设置圆柱体顶点------------------
   
10.         float radius=2.0;//半径
    
11.         float radius2=1.0;
    
12.         float fHeight=4.0;//高
    
13.         //上底
    
14.         for (float a=-1.0*radius;a<=1.0*radius+0.0001;a+=radius/25.0)
    
15.         {
    
16.                 v->push_back(osg::Vec3(a,0.0f,sqrt(fabs(pow(radius,2)-pow(a,2)))));
    
17.                 v->push_back(osg::Vec3(a*(radius2/radius),0.0f,sqrt(fabs(pow(radius2,2)-pow(a*(radius2/radius),2))))+osg::Vec3(0.0,fHeight,0.0));
    
18.         }
    
19.         for (float a=1.0*radius;a>=-1.0*radius-0.0001;a-=radius/25.0)
    
20.         {
    
21.                 v->push_back(osg::Vec3(a,0.0f,-1*sqrt(fabs(pow(radius,2)-pow(a,2)))));
    
22.                 v->push_back(osg::Vec3(a*(radius2/radius),0.0f,-1*sqrt(fabs(pow(radius2,2)-pow(a*(radius2/radius),2))))+osg::Vec3(0.0,fHeight,0.0));
    
23.         }
    
24. 
    
25.         float fPoint=204;//总共使用的顶点数       
    
26.        
    
27.         //设置法线
    
28.         osg::ref_ptr<osg::Vec3Array> n = new osg::Vec3Array;
    
29.         geom->setNormalArray( n.get() );
    
30.         geom->setNormalBinding( osg::Geometry::BIND_PER_VERTEX);
    
31. //        n->push_back( osg::Vec3( 0.0f,0.0f, 1.0f ) );
    
32.         float x2=radius*fHeight*fHeight/(fHeight*fHeight+pow((radius-radius2),2));
    
33.         float x3=x2*(radius-radius2)/fHeight;
    
34.         for (float a=-1.0*radius;a<=1.0*radius+0.0001;a+=radius/25.0)
    
35.         {
    
36.                 n->push_back(osg::Vec3(0.0,0.0,0.0)+(osg::Vec3(a*(x2/radius),0.0f,sqrt(fabs(pow(x2,2)-pow(a*(x2/radius),2))))+osg::Vec3(0.0,x3,0.0)));
    
37.                 n->push_back(osg::Vec3(0.0,0.0,0.0)+(osg::Vec3(a*(x2/radius),0.0f,sqrt(fabs(pow(x2,2)-pow(a*(x2/radius),2))))+osg::Vec3(0.0,x3,0.0)));
    
38.         }
    
39.         for (float a=1.0*radius;a>=-1.0*radius-0.0001;a-=radius/25.0)
    
40.         {
    
41.                 n->push_back(osg::Vec3(0.0,0.0,0.0)+(osg::Vec3(a*(x2/radius),0.0f,-1*sqrt(fabs(pow(x2,2)-pow(a*(x2/radius),2))))+osg::Vec3(0.0,x3,0.0)));
    
42.                 n->push_back(osg::Vec3(0.0,0.0,0.0)+(osg::Vec3(a*(x2/radius),0.0f,-1*sqrt(fabs(pow(x2,2)-pow(a*(x2/radius),2))))+osg::Vec3(0.0,x3,0.0)));
    
43.         }
    
44. 
    
45.        
    
46.         //----------------------------------------------------------
    
47. 
    
48.         // Create an array of four colors.
    
49.         osg::ref_ptr<osg::Vec4Array> c = new osg::Vec4Array;
    
50.         c->push_back( osg::Vec4( 1.f, 1.f, 1.f, 1.f ) );
    
51.         geom->setColorArray( c.get() );
    
52.         geom->setColorBinding( osg::Geometry::BIND_OVERALL );
    
53.         // Draw a four-vertex quad from the stored data.
    
54.         geom->addPrimitiveSet(new osg::DrawArrays( osg::PrimitiveSet::TRIANGLE_STRIP, 0,fPoint) );
    
55. 
    
56.         //使用smoothingVisitor
    
57. //        osgUtil::SmoothingVisitor::smooth(*geom.get());
    
58. 
    
59.         // Add the Geometry (Drawable) to a Geode and
    
60.         // return the Geode.
    
61.         osg::ref_ptr<osg::Geode> geode = new osg::Geode;
    
62.         geode->addDrawable(geom.get());
    
63. 
    
64.         osg::Matrix RoLast1;
    
65.         RoLast1.makeRotate(180.0*3.14159265/180.0,osg::Vec3f(1.0,0.0,0.0));
    
66. 
    
67.         osg::Matrix RoLast2;
    
68.         RoLast2.makeTranslate(osg::Vec3(0.0,2.75,0.0));
    
69. 
    
70.         osg::ref_ptr<osg::MatrixTransform> TransNode = new osg::MatrixTransform;
    
71.         TransNode->addChild(geode.get());
    
72.         TransNode->setMatrix(RoLast1*RoLast2);
    
73. 
    
74.         osg::StateSet* state4 = TransNode->getOrCreateStateSet();
    
75.         osg::ref_ptr<osg::Material> mat3 = new osg::Material;
    
76.         //漫反射光（决定物体颜色）
    
77.         mat3->setDiffuse( osg::Material::FRONT_AND_BACK,
    
78.                 osg::Vec4( 1.f, 0.f, 0.f, 1.f ) );
    
79.         //镜面反射光
    
80.         mat3->setSpecular( osg::Material::FRONT_AND_BACK,
    
81.                 osg::Vec4( 1.f, 1.f, 1.f, 1.f ) );
    
82.         //环境反射光
    
83.         mat3->setShininess( osg::Material::FRONT_AND_BACK, 128.f );
    
84.         state4->setAttribute( mat3.get() );
    
85.         state4->setRenderingHint( osg::StateSet::OPAQUE_BIN );
    
86. 
    
87.         return TransNode.get();
    
88. }

复制代码

array

法线是必须进行归一化的，也就是法线向量的模必须为1。否则OpenGL的光照计算会不正确。建议不要直接用setMode(GL_NORMALIZE)，会大大增加场景开销

使用SmoothVisitor的话，就不要再自己计算法线了，还有您的多边形方向可能反了，所以SmoothVisitor的法线方向也是反的。

用圆的方程x*x+y*y=r*r的话，精度损失恐怕不小，计算式看起来也有点恐怖，尤其sqrt和pow都是编程的时候应该尽量避免的；通常的方法是使用三角函数。

您所说的那两条线可能是圆锥体没有闭合的缘故，细看就知道了，事实上您的圆锥体是有缺口的。这恐怕也是浮点数以及精度的损失计算造成的，不要奢望C语言的浮点数计算准确，也永远都不要把两个浮点数作“==”的比较

Sailent

归一化的结果是精确为1还是近似为1（0.9999999或者1.00000001）呢？通过数学计算让模为1的话肯定都是近似于1的。
三角函数画圆的方法应该是（r*cos(a),0,r*sin(a)）吧（0<=a<=3.14）。
还有，我说的两条线是由光照引起的，像截图里的那样。

Sailent

归一化处理后那两条黑线没有了，设置法线代码改成这样就可以了：

1.         //设置法线
   
2.         osg::ref_ptr<osg::Vec3Array> n = new osg::Vec3Array;
   
3.         geom->setNormalArray( n.get() );
   
4.         geom->setNormalBinding( osg::Geometry::BIND_PER_VERTEX);
   
5. //        n->push_back( osg::Vec3( 0.0f,0.0f, 1.0f ) );
   
6.         float x2=radius*fHeight*fHeight/(fHeight*fHeight+pow((radius-radius2),2));
   
7.         float x3=x2*(radius-radius2)/fHeight;
   
8.         osg::Vec3 Vec3temp;
   
9.         for (float a=-1.0*radius;a<=1.0*radius+0.0001;a+=radius/25.0)
   
10.         {
    
11.                 Vec3temp=osg::Vec3(0.0,0.0,0.0)+(osg::Vec3(a*(x2/radius),0.0f,sqrt(fabs(pow(x2,2)-pow(a*(x2/radius),2))))+osg::Vec3(0.0,x3,0.0));
    
12.                 Vec3temp=Vec3temp*(1.0/sqrt(pow(Vec3temp.x(),2)+pow(Vec3temp.y(),2)+pow(Vec3temp.z(),2)));//归一化处理，让法线模为1.
    
13.                 n->push_back(Vec3temp);
    
14.                 n->push_back(Vec3temp);
    
15.         }
    
16.         for (float a=1.0*radius;a>=-1.0*radius-0.0001;a-=radius/25.0)
    
17.         {
    
18.                 Vec3temp=osg::Vec3(0.0,0.0,0.0)+(osg::Vec3(a*(x2/radius),0.0f,-1*sqrt(fabs(pow(x2,2)-pow(a*(x2/radius),2))))+osg::Vec3(0.0,x3,0.0));
    
19.                 Vec3temp=Vec3temp*(1.0/sqrt(pow(Vec3temp.x(),2)+pow(Vec3temp.y(),2)+pow(Vec3temp.z(),2)));
    
20.                 n->push_back(Vec3temp);
    
21.                 n->push_back(Vec3temp);
    
22.         }

复制代码