osgearth射线求交获取经纬度高程问题

gisbdf

我采用如下的方式来获取osgearth的经纬度和高程，osgearth模型只有一个world.tif的影像和某一块区域的DEM

1. 
   
2. 
   
3.   osg::Group* root=dynamic_cast<osg::Group*>(viewer->getSceneData());
   
4.    
   
5. osgUtil::LineSegmentIntersector::Intersections inter;
   
6. std::string mess;
   
7. std::stringstream ss;
   
8. float x=ea.getX();
   
9. float y=ea.getY();
   
10. 
    
11. if (viewer->computeIntersections(x,y,inter))
    
12. {
    
13.   
    
14.   
    
15.   for(osgUtil::LineSegmentIntersector::Intersections::iterator hitr = inter.begin();hitr != inter.end();++hitr)
    
16.   {
    
17.      //获取世界坐标系XYZ
    
18.       osg::Vec3f temp(hitr->getWorldIntersectPoint());
    
19.       double lat,lon,height;         
    
20.       osg::CoordinateSystemNode* csn=findTopMostNodeOfType<osg::CoordinateSystemNode>(root);
    
21.      //获取wgs84地理坐标
    
22.      csn->getEllipsoidModel()->convertXYZToLatLongHeight(temp.x(),temp.y(),temp.z(),lat,lon,height);
    
23.    
    
24.      ss<<"坐标： lat="<<lat*180/osg::PI<<"，lon="<<lon*180/osg::PI<<"，height="<<height<<std::endl;
    
25.   

复制代码

      



这样获取的坐标经纬度是正确的，但是高程却不正确，在远离视点的时候，高程大概都是负的几百，但是如果拉近视点，贴近表面的时候，高程大概变为正负零点几， 在加上某一块的DEM时，这一块区域的高程值通过上面的代码可以正确的获取到，这是什么原因？
      按我的理解，没有加DEM的区域，高程应该是0.
   感觉是射线求交的问题，请问这是什么原因呢？

gisbdf

通过查看一些资料，射线相交时，是与球体上的某个平面（三角行）相交，所以导致高程是负值，那怎样才能获取到正确的高程呢

FlySky

没有高程理论应该得到是无效值，上面的代码就没有问题，你这里可以加一个判断就可以设置为0

gisbdf

回复 2# gisbdf


    我查看了一些资料，好像是这样的，就是射线相交时，是与地球上的某个多边形（三角形）相交的，因为地球也是用三角形逼近的。 所以上面获取的结果的Z值，并不是地球表面上的，而是多边形的。 不知道我的理解正不正确。

gisbdf

回复 3# FlySky


   
    我查看了一些资料，好像是这样的，就是射线相交时，是与地球上的某个多边形（三角形）相交的，因为地球也是用三角形逼近的。 所以上面获取的结果的Z值，并不是地球表面上的，而是多边形的。 不知道我的理解正不正确。