关于线段求交

huahuawei88

画一条线，和场景中的模型求交，如果有交点的话输出交点坐标，为什么明明不相交，还是会输出坐标呢？

liuzhiyu123

肯定是您的使用方法有问题

huahuawei88

liuzhiyu123 发表于 2013-3-13 07:56
肯定是您的使用方法有问题

我刚接触osg，还请大师指点，该怎么用啊

liuzhiyu123

使用LineSegmentIntersection 就行了，注意坐标系的设置

huahuawei88

liuzhiyu123 发表于 2013-3-13 12:32
使用LineSegmentIntersection 就行了，注意坐标系的设置

还有个问题就是intersectionvisitor 没有addLineSegment？  就是我看的书上是intersectvisitor，，现在不是改成intersectionvisitor 了怎么就没有addLineSegment这个方法了

liuzhiyu123

不需要了 直接作为构造函数的参数传进去就行了

huahuawei88

liuzhiyu123 发表于 2013-3-13 12:43
不需要了 直接作为构造函数的参数传进去就行了

能举个例子么，我这水平实在是不行，呵呵，给您添麻烦了

liuzhiyu123

看osgPick 这个例子就行了 基本都是固定用法

Berly

我也遇到同样的问题了，都说是我计算错误，我都没计算啊，碰撞的结果总是不对，鼠标点击碰撞，可能会碰撞到模型内部，都不知道这么点是怎么来的，点击用的下面的代码：

黄色部分在模型内部
x=ea.getX();
y=ea.getY();
if (ea.getButton()==1)
{
osgUtil:ineSegmentIntersector::Intersections intersections;

if (viewer->computeIntersections(x,y,intersections))
{
osgUtil::LineSegmentIntersector::Intersection intersection=*intersections.begin();
_point=osg::Vec3(intersection.getWorldIntersectPoint());
}
这段是用线求交算的，有的时候明明碰撞到很多点，结果就给我输出了起点和终点，不知道为什么呀？？？？
还有一个问题，同一个面片，起点和终点都在这个面片上居然都会产生交点，如图 红色黄色区分产生的交点
osg::ref_ptr< osgUtil::LineSegmentIntersector > _lineSegmentIntersector = new osgUtil::LineSegmentIntersector(p_start,p_end);  
        osgUtil::IntersectionVisitor _iv(_lineSegmentIntersector.get());
        model->accept(_iv);
        osg::Vec3Array *bt=new osg::Vec3Array;
        osgUtil::LineSegmentIntersector::Intersections _intersections = _lineSegmentIntersector->getIntersections();  
        int _intersectionNumber=_intersections.size();
        osgUtil::LineSegmentIntersector::Intersections::iterator hitr = _intersections.begin();
        for (; hitr != _intersections.end();hitr++)
        {
                bt->push_back(hitr->getWorldIntersectPoint());
        }

求大师传道授业解惑……

liuzhiyu123

intersections 记录了所有的交点，您区begin 当然就是一个

Berly

liuzhiyu123 发表于 2013-3-14 08:11
intersections 记录了所有的交点，您区begin 当然就是一个

非常感谢您的回复，begin()是取得鼠标点击到的第一个交点对么？之前做相交测试的时候，array说碰撞有误差，用double类型的也不行，因为碰撞内部还是用的float，让我把鼠标点击的线段两端都延伸一段，我试了一下，碰撞的结果，还是没有规律，如图：
这里面，线是通过鼠标分别点击模型表面形成的（线段的两端点都位于模型的同一个平面内），蓝色的球代表线段和模型的交点，可见，当相交用的线段全部在平面上时，出现的情况就是可能相交，可能不交，相交可能在端点，也可能在中间任何一个点。

osg::Vec3d n1=p_start-p_end;n1.normalize();
        osg::Vec3d n2=p_end-p_start;n2.normalize();
        osg::Vec3d pStart=p_start+n1;//延长后的起点
        osg::Vec3d pEnd=p_end+n2;//延长后的终点
        osg::ref_ptr< osgUtil:ineSegmentIntersector > _lineSegmentIntersector = new osgUtil::LineSegmentIntersector(pStart,pEnd);

不知道是什么原因，还没有找到源码是怎么写的，求助……

liuzhiyu123

没有找到源码是怎么写的？ 完完整整的都在LineSegmentIntersector.cpp 中

array

1. osg::Vec3d n1=p_start-p_end;n1.normalize();
   
2. osg::Vec3d n2=p_end-p_start;n2.normalize();
   
3. osg::Vec3d pStart=p_start+n1;//延长后的起点
   
4. osg::Vec3d pEnd=p_end+n2;//延长后的终点

复制代码

。。。您这就叫做延长了么，还有您可以提交一个足够短的，可以运行且说明问题的程序让别人帮助调试。否则单凭您每次提供一点点散碎的信息，没有办法确定问题所在

Berly

array 发表于 2013-3-14 10:20
。。。您这就叫做延长了么，还有您可以提交一个足够短的，可以运行且说明问题的程序让别人帮助调试。否则单 ...

那个不就是按照对应的方向向外延长单位向量么？难道不对？希望指点，拜谢……


下面是测试代码

1. #include <osgViewer/Viewer>
   
2. #include <osgViewer/ViewerEventHandlers>
   
3. 
   
4. #include <osg/Geode>
   
5. #include <osg/Geometry>
   
6. #include <osg/Node>
   
7. #include <osg/Group>
   
8. #include <osg/LineSegment>
   
9. #include <osg/ShapeDrawable>
   
10. #include <osg/BlendColor>
    
11. #include <osg/Material>
    
12. #include <osg/BoundingBox>
    
13. #include <osg/ComputeBoundsVisitor>
    
14. #include <osg/Shape>
    
15. #include <osg/PositionAttitudeTransform>
    
16. 
    
17. #include <osgDB/ReadFile>
    
18. #include <osgDB/WriteFile>
    
19. 
    
20. #include <osgGA/GUIEventHandler>
    
21. #include <osgGA/StateSetManipulator>
    
22. 
    
23. #include <osgUtil/Optimizer>
    
24. #include <osgUtil/IntersectionVisitor>
    
25. #include <osgUtil/IntersectVisitor>
    
26. 
    
27. #include <math.h>
    
28. 
    
29. using namespace std;
    
30. typedef std::vector<osg::Vec3dArray *> ArrayVector;
    
31. osg::ref_ptr<osg::Group>modelGroup=new osg::Group;
    
32. osg::Vec3Array *pointList=new osg::Vec3Array;
    
33. 
    
34. void getIntersections( osg::Vec3dArray *buttomPointList,osg::Group* sphereGroup, osg::Node* model, osg::Vec3d p_start,osg::Vec3d p_end)
    
35. {
    
36.         osg::Vec3d n1=p_start-p_end;n1.normalize();
    
37.         osg::Vec3d n2=p_end-p_start;n2.normalize();
    
38.         osg::Vec3d pStart=p_start+n1;
    
39.         osg::Vec3d pEnd=p_end+n2;
    
40.         osg::ref_ptr< osgUtil::LineSegmentIntersector > _lineSegmentIntersector = new osgUtil::LineSegmentIntersector(pStart,pEnd);  
    
41.         osgUtil::IntersectionVisitor _iv(_lineSegmentIntersector.get());
    
42.         model->accept(_iv);
    
43.         osgUtil::LineSegmentIntersector::Intersections _intersections = _lineSegmentIntersector->getIntersections();  
    
44.         int _intersectionNumber=_intersections.size();
    
45.         osgUtil::LineSegmentIntersector::Intersections::iterator hitr = _intersections.begin();
    
46. 
    
47.         osg::ref_ptr<osg::Geode>viewLine=new osg::Geode;
    
48.         osg::ref_ptr<osg::Geometry>geom=new osg::Geometry;
    
49.         osg::ref_ptr<osg::Vec3Array>v=new osg::Vec3Array;
    
50.         osg::Vec4 color=osg::Vec4(1.0f,0.0f,0.0f,1.0f);
    
51.         v->push_back(p_start);
    
52.         v->push_back(p_end);
    
53.         geom->setVertexArray(v.get());
    
54.         osg::ref_ptr<osg::Vec4Array>vc=new osg::Vec4Array;
    
55.         vc->push_back(color);
    
56.         geom->setColorArray(vc.get());
    
57.         geom->addPrimitiveSet(new osg::DrawArrays(osg::PrimitiveSet::LINES,0,2));
    
58.         viewLine->addDrawable(geom);
    
59.     sphereGroup->addChild(viewLine);
    
60.         for (; hitr != _intersections.end();hitr++)
    
61.         {   
    
62.                 osg::Vec3d point=hitr->getWorldIntersectPoint();
    
63.                 cout<<"    x:   "<<point.x()<<"    y:   "<<point.y()<<"    z:   "<<point.z()<<endl;
    
64.                 buttomPointList->push_back(point);
    
65.                                  osg::Geode* geode=new osg::Geode;
    
66.                                  osg::ShapeDrawable *viewPoint=new osg::ShapeDrawable(new osg::Sphere(point,0.1f));
    
67.                                  viewPoint->setColor(osg::Vec4(0.0f,0.0f,1.0f,1.0f));
    
68.                                  geode->addDrawable(viewPoint);
    
69.                                  sphereGroup->addChild(geode);
    
70.         }
    
71. }
    
72. class PolygonHandler:public osgGA::GUIEventHandler
    
73. {
    
74. public:
    
75.         PolygonHandler()
    
76.         {
    
77.         }
    
78.         ~PolygonHandler()
    
79.         {}
    
80.         bool PolygonHandler::handle(const osgGA::GUIEventAdapter &ea,osgGA::GUIActionAdapter &aa)
    
81.         {
    
82.                 osgViewer::View *viewer=dynamic_cast<osgViewer::View*>(&aa);
    
83.                 switch(ea.getEventType())
    
84.                 {
    
85.                 case (osgGA::GUIEventAdapter::PUSH):
    
86.                         {
    
87.                                 x=ea.getX();
    
88.                                 y=ea.getY();
    
89.                                 if (ea.getButton()==1)
    
90.                                 {
    
91.                                         osgUtil::LineSegmentIntersector::Intersections intersections;
    
92.                                         if (viewer->computeIntersections(x,y,intersections))
    
93.                                         {
    
94.                                                 osgUtil::LineSegmentIntersector::Intersection intersection=*intersections.begin();
    
95.                                                 _point=osg::Vec3d(intersection.getWorldIntersectPoint());
    
96.                                             pointList->push_back(_point);
    
97.                                                 if (pointList->size()>=2&&pointList->size()%2==0)
    
98.                                                 {
    
99.                                                         osg::Vec3d Point1=pointList->at(pointList->size()-2);
    
100.                                                         osg::Vec3d Point2=pointList->at(pointList->size()-1);
     
101.                                                         osg::ref_ptr< osg::Vec3dArray> buttomList  =  new osg::Vec3dArray;
     
102.                                                         osg::ref_ptr<osg::Group>group=new osg::Group;
     
103.                                                         getIntersections(buttomList,group,modelGroup,Point1,Point2);
     
104.                                                         viewer->getSceneData()->asGroup()->addChild(group);
     
105.                                                 }
     
106.                                         }
     
107.                                 }
     
108.                                 return false;
     
109.                         }break;
     
110.                 default:
     
111.                         return false;
     
112.                 }
     
113.         }
     
114. private:
     
115.         osg::Vec3d _point;
     
116.         double x;
     
117.         double y;
     
118. };
     
119. 
     
120. 
     
121. int main()
     
122. {
     
123.         osg::ref_ptr<osgViewer::Viewer>viewer=new osgViewer::Viewer;
     
124.         osg::ref_ptr<osg::Group>root=new osg::Group;
     
125.         osg::ref_ptr<osg::Node>model=osgDB::readNodeFile("moxing.osg");
     
126. 
     
127.         modelGroup->addChild(model.get());
     
128.         root->addChild(modelGroup.get());
     
129. 
     
130.         viewer->addEventHandler(new osgViewer::WindowSizeHandler);
     
131.         viewer->addEventHandler(new osgGA::StateSetManipulator(viewer->getCamera()->getOrCreateStateSet()));
     
132.         viewer->addEventHandler(new PolygonHandler);
     
133. 
     
134.         osgUtil::Optimizer optimizer;
     
135.         optimizer.optimize(root.get());
     
136. 
     
137.         viewer->setSceneData(root.get());
     
138.         viewer->realize();
     
139.         viewer->run();
     
140.         return 0;
     
141. }

复制代码

moxing.osg.gz (110.35 KB, 下载次数: 1)

2013-3-14 10:50 上传

点击文件名下载附件


moxing.osg在附件里面，传不上rar,加了个后缀，去掉“.gz”就行了,没有纹理，这样小点，点击同一个面的时候，有时相交，有时不交，有时有起始点，有时候没有……，希望善良好心的人帮我看看呗，谢谢哈

mxl12315

妹纸，我也最近做这个~对于碰撞检测也是不太清楚，总有些“灵异现象”（如你贴的图中）……还望此贴能早日解惑

huahuawei88

array 发表于 2013-3-14 10:20
。。。您这就叫做延长了么，还有您可以提交一个足够短的，可以运行且说明问题的程序让别人帮助调试。否则单 ...

大神出现了，我顺便再问个问题：从一个点出发做多条射线，做相交检测，怎么返回有多少条线是有相交多少条是没相交的。

liuzhiyu123

huahuawei88 发表于 2013-3-14 20:13
大神出现了，我顺便再问个问题：从一个点出发做多条射线，做相交检测，怎么返回有多少条线是有相交多少条 ...

使用 IntersectorGroup

liuzhiyu123

这是什么逻辑？
                                        if (viewer->computeIntersections(x,y,intersections))
                                        {
                                                osgUtil::LineSegmentIntersector::Intersection intersection=*intersections.begin();
                                                _point=osg::Vec3d(intersection.getWorldIntersectPoint());
                                            pointList->push_back(_point);
                                                if (pointList->size()>=2&&pointList->size()%2==0)
                                                {
                                                        osg::Vec3d Point1=pointList->at(pointList->size()-2);
                                                        osg::Vec3d Point2=pointList->at(pointList->size()-1);
                                                        osg::ref_ptr< osg::Vec3dArray> buttomList  =  new osg::Vec3dArray;
                                                        osg::ref_ptr<osg::Group>group=new osg::Group;
                                                        getIntersections(buttomList,group,modelGroup,Point1,Point2);
                                                        viewer->getSceneData()->asGroup()->addChild(group);
                                                }

Berly

liuzhiyu123 发表于 2013-3-15 08:14
这是什么逻辑？
                                        if (viewer->computeIntersections(x,y,interse ...

就是只有当单击的点数是2或者2的倍数的时候，就取最后点击的两个点作为线段相交的起点和终点

huahuawei88

liuzhiyu123 发表于 2013-3-15 08:10
使用 IntersectorGroup

从一点出发画很多射线，有什么简单的方法么

liuzhiyu123

huahuawei88 发表于 2013-3-18 13:47
从一点出发画很多射线，有什么简单的方法么

暂时没有，除非自己领想别的算法了

huahuawei88

liuzhiyu123 发表于 2013-3-18 14:01
暂时没有，除非自己领想别的算法了

那该怎么弄啊，一条一条的话吗？能详细讲讲吗？或者给点代码也行，新手，麻烦大神了！

liuzhiyu123

IntersectorGroup 将每个LinesegmentIntersector 添加到其中void addIntersector(Intersector* intersector)，然后获取的时候 对添加的每个LinesegmentIntersector 进行操作就行了

huahuawei88

liuzhiyu123 发表于 2013-3-18 15:01
IntersectorGroup 将每个LinesegmentIntersector 添加到其中void addIntersector(Intersector* intersector ...

就是现在我想把这些线也当做Geometry画出来，就是可以看见的那种，怎么办？
DrawArrays(osg:rimitiveSet:INES,0,2) ，每次只能画一条，要是线很多的话，上百条，该怎么办？（别笑哈，实在是太笨了）

liuzhiyu123

huahuawei88 发表于 2013-3-18 18:15
就是现在我想把这些线也当做Geometry画出来，就是可以看见的那种，怎么办？
DrawArrays(osg:rimitiveS ...

哦？您那么确定就只能画一条？看看下面的您或许会有答案

garyliyong

Berly 发表于 2013-3-14 09:50
非常感谢您的回复，begin()是取得鼠标点击到的第一个交点对么？之前做相交测试的时候，array说碰撞有误差 ...

其实可以把intesect函数里面的第三句话注释的 具体名字忘了好像是判断直线和几何体的包围盒是否相交（里面会将直线裁剪的）

huahuawei88

liuzhiyu123 发表于 2013-3-19 07:46
哦？您那么确定就只能画一条？看看下面的您或许会有答案

1. #include <Windows.h>
   
2. #include <osgViewer/Viewer>
   
3. #include <osgViewer/ViewerEventHandlers>
   
4. #include <osgGA/StateSetManipulator>
   
5. #include <osgDB/ReadFile>
   
6. #include <osg/LineSegment>
   
7. #include <osgUtil/LineSegmentIntersector>
   
8. #include <osgUtil/IntersectionVisitor>
   
9. #include <osgUtil/IntersectVisitor>
   
10. #include <osg/MatrixTransform>
    
11. #include <osg/Group>
    
12. #include <osg/Node>
    
13. #include "iostream"
    
14. using namespace std;
    
15. int main()
    
16. {
    
17.         osg::ref_ptr<osg::Node> cow = osgDB::readNodeFile("cow.osg");//第一个节点
    
18.         osg::ref_ptr<osg::Node> axes = osgDB::readNodeFile("axes.osgt");//第二个节点
    
19.         osg::ref_ptr<osg::MatrixTransform> mt_cow = new osg::MatrixTransform;
    
20.         mt_cow->setMatrix(osg::Matrix::scale(osg::Vec3(0.5,0.5,0.5))*osg::Matrix::translate(osg::Vec3(-10,0,0)));
    
21.         mt_cow->addChild(cow.get());
    
22.         osg::ref_ptr<osg::MatrixTransform> mt_axes = new osg::MatrixTransform;
    
23.         mt_axes->addChild(axes.get());
    
24.         osg::ref_ptr<osg::Group> root = new osg::Group;
    
25.         root->addChild(mt_cow.get());
    
26.         root->addChild(mt_axes.get());
    
27.         osgViewer::Viewer * viewer = new osgViewer::Viewer;
    
28.         viewer->addEventHandler(new osgGA::StateSetManipulator(viewer->getCamera()->getOrCreateStateSet()));
    
29.         osg::Vec4Array* colors = new osg::Vec4Array;
    
30.         colors->push_back(osg::Vec4(1.0f,0.0f,0.0f,1.0f));
    
31. 
    
32. // 碰撞检测
    
33.         osg::ref_ptr<osg::Group> group = new osg::Group;
    
34.         osg::Vec3f p_end;
    
35.         osg::Vec3 p_start(0,0,0);
    
36.         int j=0;
    
37.         for (int i=0;i<100;++i)
    
38.         {               
    
39.                 osg::Vec3Array* coords = new osg::Vec3Array;
    
40.                 coords->push_back(p_start);  
    
41.                 float dt = 5/100.0f;
    
42.                 coords->push_back(osg::Vec3f(-10,dt*i,0));
    
43.                 osg::Vec3Array* normals = new osg::Vec3Array;
    
44.                 normals->push_back(osg::Vec3(0.0f,0.0f,1.0f));                  
    
45.                 osg::ref_ptr<osg::Geometry> geom = new osg::Geometry;
    
46.                 geom->setVertexArray(coords);
    
47.                 geom->setNormalArray(normals);
    
48.                 geom->setColorArray(colors);
    
49.                 geom->setColorBinding(osg::Geometry::BIND_PER_PRIMITIVE);
    
50.                 geom->addPrimitiveSet(new osg::DrawArrays(osg::PrimitiveSet::LINES,0,2));
    
51.                 osg::ref_ptr<osg::Geode> geode = new osg::Geode;
    
52.                 geode->addDrawable(geom.get());
    
53.                 geode->setName("line");
    
54.                 group->addChild(geode.get());
    
55.                 root->addChild(group.get());
    
56. 
    
57.                 osgUtil::LineSegmentIntersector::Intersections _intersections;       
    
58.                 osg::ref_ptr< osgUtil::LineSegmentIntersector > _lineSegmentIntersector = new osgUtil::LineSegmentIntersector(p_start,osg::Vec3f(-10,dt*i,0));
    
59.                 osgUtil::IntersectionVisitor _iv(_lineSegmentIntersector.get());
    
60.                 //viewer->getSceneData()->accept(_iv);
    
61.                 cow->accept(_iv);
    
62.                 _intersections=_lineSegmentIntersector->getIntersections();
    
63.                 int _intersectionNumber=_intersections.size();  
    
64.                 if (_intersectionNumber!=0)
    
65.                 {   
    
66.                         osgUtil::LineSegmentIntersector::Intersections::iterator hitr = _intersections.begin();
    
67.                         osg::Vec3 vec3_p1 =hitr->getWorldIntersectPoint();
    
68.                         mt_axes->setMatrix(osg::Matrix::translate(vec3_p1));
    
69.                         hitr++;
    
70.                         colors->push_back(osg::Vec4(0.0f,1.0f,0.0f,1.0f));
    
71.                         cout<<vec3_p1.x()<<"  "<<vec3_p1.y()<<"   "<<vec3_p1.z()<<endl;
    
72.                         j++;                       
    
73.                 }                       
    
74.         }
    
75.         cout<<j;
    
76.         system("pause");       
    
77.         viewer->setSceneData(root.get());
    
78.         return viewer->run();
    
79. 
    
80. }

复制代码

javascript:;javascript:;

这是我做的一个相交检测的例子，输出交点的坐标，，这100条线里面有几条是个模型有交的，但是却输出了100个交点坐标，这些点的坐标看上去也不对，求指导啊，还是对相交检测这段代码不是很明白，要是能单独开个帖子把这一块代码好好讲一下就好了，方便新手啊，呵呵，，真是给您添麻烦了，感谢你一直给我解答

huahuawei88

garyliyong 发表于 2013-3-19 10:08
其实可以把intesect函数里面的第三句话注释的 具体名字忘了好像是判断直线和几何体的包围盒是否相交（里面 ...

您说的这个跟物体的包围盒有交点我也有点这个感觉，，，能不能说的详细一点啊，在下新手，急等，盼……

liuzhiyu123

mt_cow->setMatrix(osg::Matrix::scale(osg::Vec3(0.5,0.5,0.5))*osg::Matrix::translate(osg::Vec3(-10,0,0))); 牛已经有变换了，cow->accept(_iv); 已经不是在同一个父坐标下了，显示的东西估计是不对的

garyliyong

huahuawei88 发表于 2013-3-19 15:33
您说的这个跟物体的包围盒有交点我也有点这个感觉，，，能不能说的详细一点啊，在下新手，急等，盼……

osg::Vec3d s(_start), e(_end);
if ( !intersectAndClip( s, e, drawable->getBound() ) ) return;
这两句是将直线与包围盒求交，并将直线与包围盒相交的最小点和最大点作为直线的起始点，你将第二句注释掉，就是原始的直线段求交了，这个你试试看吧。我也遇到过求交失败的情况，就将这句注掉了。