射线求交获得地球表面高度

fqyd1987

各位大大好,我在用射线求交获得地球表面的经纬度高，获得的经纬度是正确的而高度不正确，测得大陆的高度却是负数，不知道为什么，使用的求交代码如下：
osg::ref_ptr<osg:ineSegment> lineSegment = new osg::LineSegment();
osg::Vec3 terrainHeight;
lineSegment->set(
  osg::Vec3(x, y, z) ,
  osg::Vec3(x1, y1, z1) );
osgUtil::IntersectVisitor intersectVisitor;
intersectVisitor.addLineSegment(lineSegment);
node->accept(intersectVisitor);
osgUtil::IntersectVisitor::HitList hitList;
hitList =intersectVisitor.getHitList(lineSegment);
osgUtil::Hit heightTestResults;
if ( hitList.empty() )
{
  std::cout << " couldn't place the model on terrain" << std::endl;
  AfxMessageBox(L"empty");
  }
heightTestResults = hitList.front();
osg::Vec3 terrainHeight=heightTestResults.getLocalIntersectPoint();
return terrainHeight;

array

1、请使用getWorldIntersectPoint()获取世界坐标而非局部坐标
2、IntersectVisitor有些太古老了，建议您使用IntersectionVisitor和LinesegmentIntersector

fqyd1987

哦，谢谢array老大回复，那个getWorldIntersectPoint()是写错了，已经改了，高度还是有问题，我去看了下osgPick， osg::ref_ptr< osgUtil:ineSegmentIntersector > picker = new osgUtil::LineSegmentIntersector(osgUtil::Intersector::MODEL, vecStart, vecEnd);
        osgUtil::IntersectionVisitor iv(picker.get());
        node->accept(iv);
if (picker->containsIntersections())
        {
                intersections = picker->getIntersections();
                for(osgUtil::LineSegmentIntersector::Intersections::iterator hitr = intersections.begin();
                        hitr != intersections.end();
                        ++hitr)
                {
                       
                        terrainHeight=hitr->getWorldIntersectPoint().z();

               
                }

        }
        改成这样，高度也不对。。。不知道写的哪个地方有问题

array

我不知道您所说的不对是怎么讲。既然您要求取经纬度高度的话，应该有相应的转换过程才对，因为这里getWorldIntersectPoint直接得到的是笛卡儿坐标系下的Z坐标。但是我并没有看到您之后的经纬度转换过程

fqyd1987

osg::EllipsoidModel::convertXYZToLatLongHeight我采用的这个把  hitr->getWorldIntersectPoint()  转换成经纬高的：
osg::ref_ptr<osg::EllipsoidModel> Em=new osg::EllipsoidModel;
double x,y,z;
Em->convertXYZToLatLongHeight(hitr->getWorldIntersectPoint().x(),hitr->getWorldIntersectPoint().y(),hitr->getWorldIntersectPoint().z(),lat,lon,hei)

我的方法就是把当前的经纬高通过convertLatLongHeightToXYZ转换成世界坐标，作为射线求交的起点，然后把（0.0，0.0，0.0）作为射线求交的终点，与地球节点进行求交得到交点，再转换通过convertXYZToLatLongHeight成经纬度高，高度负的不对着。。。

array

您是否考虑了地球中心位置的影响？如果地球中心不在零点，要把getWorldIntersectPoint减去一个偏移量才行