如何求一条线段和网格的所有交点

vincent

一条线段，给定起始点，如何把这条线段映射到地形表面，求出与地形网格的所有交点。
或者怎么贴着地面画线？
谢谢！

附：
附件中图片，红黄交界处是线段和地形的交界处

array

可以用LinesegmentIntersector来求，但是线段与地形网格的交点是不能得到贴地面的线的。
您要得到线段与地面的交线路径的话，可以用线段所在的垂直平面与地面求交，也就是PlaneIntersector

vincent

好的，我试试，谢谢array！

vincent

再次请教版主及各位高手。
用平面与地形求交，求出的交点有问题，应该是线段所在的垂直平面求错了。
不知道怎么解决？谢谢各位。

1. //起点startPt，终点endPt
   
2. //垂直地形的平面--,我选取的确定平面的第三点是不是有问题？？？
   
3. osg::Plane plane(startPt, endPt,osg::Vec3(startPt.x(),startPt.y(),
   
4.     (startPt.z()+endPt.z())/2.0));
   
5. //法线
   
6. osg::Vec3f norm(endPt-startPt);
   
7. norm.normalize();       
   
8. osg::Polytope polytope;
   
9. polytope.add( osg::Plane(norm,startPt));
   
10. polytope.add( osg::Plane(-norm, endPt));
    
11.                                                
    
12. osg::ref_ptr<osgUtil::PlaneIntersector>intersector = new osgUtil::PlaneIntersector( plane, polytope );
    
13. osgUtil::IntersectionVisitor iv( intersector.get() );
    
14. root->accept( iv );
    
15. osg::ref_ptr <osg::Vec3Array>  linepoint = new osg::Vec3Array;
    
16. 
    
17. if ( intersector->containsIntersections() )
    
18. {
    
19. osgUtil::PlaneIntersector::Intersections& intersections =
    
20. intersector->getIntersections();
    
21. osgUtil::PlaneIntersector::Intersections::iterator itr;
    
22. for ( itr = intersections.begin(); itr != intersections.end(); ++itr )
    
23. {
    
24.    for( unsigned int i=0; i<(*itr).polyline.size(); ++i)
    
25.     {
    
26.         fx1 = (*itr).polyline[i].x();       
    
27.         fy1 = (*itr).polyline[i].y()       
    
28.         fz1 = (*itr).polyline[i].z();       
    
29.         linepoint->push_back(osg::Vec3(fx1,fy1,fz1));
    
30.     }
    
31. }
    
32. osg::Geode* pyramidGeode = new osg::Geode();
    
33. osg::Geometry* pyramidGeometry = new osg::Geometry();
    
34. pyramidGeometry->setVertexArray(linepoint.get());
    
35. 
    
36. osg::Vec4Array* colors = new osg::Vec4Array;
    
37. colors->push_back(osg::Vec4(1.0f, 0.0f, 0.0f, 1.0f) ); //index 0 red
    
38. 
    
39. pyramidGeometry->setColorArray(colors);
    
40. pyramidGeometry->setColorBinding(osg::Geometry::BIND_PER_VERTEX);   
    
41. pyramidGeometry->addPrimitiveSet(new osg::DrawArrays(osg::PrimitiveSet::LINE_LOOP, 0, linepoint->size()));
    
42. 
    
43. pyramidGeometry->getOrCreateStateSet()->setMode(GL_LIGHTING, false);
    
44. pyramidGeode->addDrawable(pyramidGeometry);
    
45. root->addChild(pyramidGeode);

复制代码

array

您可以参看一下osgSim::ElevationSlice::computeIntersections()函数中的做法