osgEarth2.10根据输入的经纬度坐标获取高程的C++代码

以下是使用osgEarth 2.10获取给定经纬度坐标的高程信息的C++代码示例：

#include <osgEarth/Map>
#include <osgEarth/MapNode>
#include <osgEarthUtil/ElevationQuery>
#include <osgEarthUtil/EarthManipulator>

osg::ref_ptr<osgEarth::MapNode> createMapNode()
{
    // 创建一个空的地图
    osg::ref_ptr<osgEarth::Map> map = new osgEarth::Map();

    // 创建一个地形图层
    osgEarth::Drivers::SimpleOceanLayerFactory oceanFactory;
    osg::ref_ptr<osgEarth::Layer> oceanLayer = oceanFactory.create(map.get(), nullptr);
    map->addLayer(oceanLayer.get());

    // 创建一个工具层
    osgEarth::Drivers::SimpleElevationLayerFactory elevationFactory;
    osg::ref_ptr<osgEarth::Layer> elevationLayer = elevationFactory.create(map.get(), nullptr);
    map->addLayer(elevationLayer.get());

    // 创建地图节点
    osgEarth::MapNodeOptions options;
    options.map() = map.get();
    return new osgEarth::MapNode(options);
}

double getElevation(double latitude, double longitude)
{
    // 创建地图节点
    osg::ref_ptr<osgEarth::MapNode> mapNode = createMapNode();

    // 创建高程查询对象
    osg::ref_ptr<osgEarth::Util::ElevationQuery> elevationQuery = new osgEarth::Util::ElevationQuery(mapNode.get());

    // 设置查询的经纬度坐标
    osgEarth::GeoPoint point(
        mapNode->getMapSRS(),
        osgEarth::Units::DEGREES,
        longitude,
        latitude);
    elevationQuery->setPoint(point);

    // 执行高程查询
    osgEarth::Util::ElevationQuery::QueryResult result;
    if (elevationQuery->getQuery(result))
    {
        return result.elevation().as(osgEarth::Units::METERS);
    }
    else
    {
        // 查询失败
        return 0.0;
    }
}

int main()
{
    double latitude = 40.7128;   // 纬度
    double longitude = -74.0060; // 经度

    double elevation = getElevation(latitude, longitude);
    std::cout << "Elevation at (" << latitude << ", " << longitude << "): " << elevation << " meters" << std::endl;

    return 0;
}

以上代码利用osgEarth库获取给定经纬度坐标的高程信息。需要注意的是，你需要正确配置osgEarth库及其依赖，并将相关头文件和库文件链接到你的项目中。此外，你也需要使用合适的地形图层数据来获得准确的高程信息。