如何实现地形对包络体的遮挡效果？

口口广大

osg有个例子如下图，实现了包络体SphereSegment.与地形模型相切。




想实现一个更复杂的效果，即是在生成包络体时，判断地形有否遮挡住该包络体，最后只生成该包络体未被地形遮蔽的部分。

想想原理很简单，查看SphereSegment.cpp，在画surfac的PowerSphereSegment::Surface_drawImplementation源码如下

1. for(int i=0; i+1<=_density; i++)
   
2.             {
   
3.                 // Because we're drawing quad strips, we need to work out
   
4.                 // two azimuth values, to form each edge of the (z-vertical)
   
5.                 // strips
   
6.                 float az1 = _azMin + (i*azIncr);
   
7.                 float az2 = _azMin + ((i+1)*azIncr);
   
8. 
   
9.                 glBegin(GL_QUAD_STRIP);
   
10.                 for (int j=0; j<=_edensity; j++)
    
11.                 {
    
12.                     float elev = _elevMin + (j*elevIncr);
    
13. 
    
14.                     // QuadStrip Edge formed at az1
    
15.                     // ----------------------------
    
16. 
    
17.                     // Work out the sphere normal
    
18.                     float x = cos(elev)*sin(az1);
    
19.                     float y = cos(elev)*cos(az1);
    
20.                     float z = sin(elev);
    
21. 
    
22.                     glNormal3f(-x, -y, -z);               
    
23.                     glVertex3f(_centre.x() + _radius*x,
    
24.                             _centre.y() + _radius*y,
    
25.                             _centre.z() + _radius*z);
    
26.                           
    
27.                     // QuadStrip Edge formed at az2
    
28.                     // ----------------------------
    
29. 
    
30.                     // Work out the sphere normal
    
31.                     x = cos(elev)*sin(az2);
    
32.                     y = cos(elev)*cos(az2);
    
33.                     // z = sin(elev);   z doesn't change
    
34. 
    
35.                     glNormal3f(-x, -y, -z);
    
36.                     glVertex3f(_centre.x() + _radius*x,
    
37.                             _centre.y() + _radius*y,
    
38.                             _centre.z() + _radius*z);
    
39.                 }
    
40.                 glEnd();
    
41.             }
    
42.           }
    

复制代码

无非将方位角与俯仰角划分成i*j份，然后每一份根据某个数学公式算出两个点，然后连起来即可。

那么我想有个简单的办法，即判断每一份两点的连线与地形模型是否有交，如果有交就退出当层循环，那么最后画出的片面应该就是未被地形遮挡的片面。


问题就是不知道如何写线段与地形模型相交的函数，呵呵。写了试试如下

1. 
   
2. int PowerSphereSegment::computeIntersectionTerrain(osg::Vec3d& start,osg::Vec3d& end)
   
3. {
   
4.   int result = 0;
   
5. 
   
6.    osgUtil::IntersectVisitor ivXY;
   
7.    osg::ref_ptr<osg::LineSegment> lineXY = new osg::LineSegment(start,end);
   
8.    ivXY.addLineSegment(lineXY.get());
   
9.       _terrainNode->accept(ivXY);
   
10. 
    
11.     if (ivXY.hits())
    
12.    {
    
13.                 osgUtil::IntersectVisitor::HitList& hitList = ivXY.getHitList(lineXY.get());
    
14.                 std::cout<<"  hitList SIZE = "<<hitList.size()<<std::endl;
    
15.    
    
16.                 result =  hitList.size();
    
17.     }
    
18. 
    
19.          return 0;
    
20. }

复制代码

执行时报警告如下
Warning: invalid line segment passed to IntersectVisitor::addLineSegment(..)
         0 0 0 0 0 0 segment ignored..

不知有何问题?而且我感觉这个osg::Vec3d& start,osg::Vec3d& end是该包络体的局部坐标系，而地形模型则是另一个局部坐标系，两个不同坐标系下的几何体能否直接相交？是不是要矩阵变化什么的的？谢谢！

liushui

写线段与地形模型相交的函数

用osgUtil:ineSegmentIntersector

口口广大

liushui 发表于 2011-11-8 15:15
写线段与地形模型相交的函数

用osgUtil:ineSegmentIntersector

谢谢，利用osgUtil:ineSegmentIntersector改写了求交函数

1. 
   
2. int PowerSphereSegment::computeIntersectionTerrain(osg::Vec3d& start,osg::Vec3d& end)
   
3. {
   
4. osgUtil::LineSegmentIntersector::Intersections _intersections;
   
5. osg::ref_ptr< osgUtil::LineSegmentIntersector > _lineSegmentIntersector = new osgUtil::LineSegmentIntersector(start,end);
   
6. osgUtil::IntersectionVisitor _iv(_lineSegmentIntersector.get());
   
7. _terrainNode->accept(_iv);
   
8. _intersections=_lineSegmentIntersector->getIntersections();
   
9. return_intersections.size();        
   
10. }

复制代码

最后返回的结果始终是0个交点，我想问题可能是传入的osg::Vec3d& start,osg::Vec3d& end点是该包罗体的局部坐标，是不是需要先求出这两个点的世界坐标才能用ineSegmentIntersector求交点？


如果利用osg::computeLocalToWorld函数求世界坐标，一般来说要用一个NodeVisitor。该访问器类用于返回某个节点的世界坐标， 它从起始节点开始向父节点遍历，并随时将历经的节点记录到nodePath中。
，第一次到达根节点之后，它将记录起始节点的世界坐标。连结起始节点到根节点路径上的 ，所有矩阵之后，即可获得节点的世界坐标。


这时仿佛有个矛盾。矛盾就是一般的思路是我先生成这个包罗体以后，然后通过旋转、缩放、平移矩阵将它加载到世界坐标系下，此时通过NodeVisitor能获得根节点路径。而我此刻在PowerSphereSegment::Surface_drawImplementation中调用 computeIntersectionTerrain函数时，我的包罗体尚未加载任何旋转矩阵，也就是还未有任何到父节点的路径，那么我此时该如何求包罗体两点的世界坐标？谢谢！

liushui

不矛盾，比如碰撞检测，你向前移动一个步长，是要首先判断前面是否有障碍物，如果有计算可以移动的距离等等，算好了预期位置之后，再把自己的位置放到那里。

我想是一回事