osgEarth例子学习-feature_elevation.earth

最新推荐文章于 2024-07-10 09:39:44 发布

努力的包老板

最新推荐文章于 2024-07-10 09:39:44 发布

阅读量3.9k

点赞数 6

分类专栏： osgEarth2019 文章标签： osgEarth例子压平 elevation

本文链接：https://blog.csdn.net/qq_36881934/article/details/87933122

版权

本文记录了使用osgEarth进行地形处理的探索过程，包括地形数据的加载、高程显示、以及尝试创建下沉平台效果。通过解决坐标投影问题和数据转换，最终实现了地形的3D效果，但遇到的挑战还有待进一步解决。

摘要由CSDN通过智能技术生成

例子学习第一天

今天是2019年2.26，接下来的一段时间开始研究关于在osgEarth中地形整（压平）的事儿，最后的应用就是开路，在osgEarth地球的表面，比如山上能自动生成道路之类的。
首先看下例子的效果介绍
在这里插入图片描述
还是很酷炫的，羡慕啊，开始干吧。
好久没摸oe了，只能一点一点来，在vs2010中随便选了个例子cpp，这里选择的是

但是里面有好多用不上的代码（应该吧），先注销掉。上面一堆包含都没管，菜鸟一只，也分辨不出来，啥Handler都留着，没管。



#include <osgGA/StateSetManipulator>
#include <osgGA/GUIEventHandler>
#include <osgViewer/Viewer>
#include <osgViewer/ViewerEventHandlers>
#include <osgUtil/LineSegmentIntersector>
#include <osgEarth/MapNode>
#include <osgEarth/TerrainEngineNode>
#include <osgEarth/ElevationQuery>
#include <osgEarth/StringUtils>
#include <osgEarth/Terrain>
#include <osgEarth/VerticalDatum>
#include <osgEarthUtil/EarthManipulator>
#include <osgEarthUtil/Controls>
#include <osgEarthUtil/LatLongFormatter>
#include <osgEarthUtil/ExampleResources>
#include <osgEarthAnnotation/PlaceNode>
#include <iomanip>

using namespace osgEarth;
using namespace osgEarth::Util;
using namespace osgEarth::Util::Controls;
using namespace osgEarth::Annotation;

static MapNode*       s_mapNode     = 0L;
static LabelControl*  s_posLabel    = 0L;
static LabelControl*  s_vdaLabel    = 0L;
static LabelControl*  s_mslLabel    = 0L;
static LabelControl*  s_haeLabel    = 0L;
static LabelControl*  s_egm96Label  = 0L;
static LabelControl*  s_mapLabel    = 0L;
static LabelControl*  s_resLabel    = 0L;
static PlaceNode*     s_marker      = 0L;


// An event handler that will print out the elevation at the clicked point
struct QueryElevationHandler : public osgGA::GUIEventHandler 
{
    QueryElevationHandler()
        : _mouseDown( false ),
          _terrain  ( s_mapNode->getTerrain() ),
          _query    ( s_mapNode->getMap() )
    {
        _map = s_mapNode->getMap();
        _query.setMaxTilesToCache(10);
        _query.setFallBackOnNoData( false );
        _path.push_back( s_mapNode->getTerrainEngine() );
    }

    void update( float x, float y, osgViewer::View* view )
    {
        bool yes = false;

        // look under the mouse:
        osg::Vec3d world;
        osgUtil::LineSegmentIntersector::Intersections hits;
        if ( view->computeIntersections(x, y, hits) )
        {
            world = hits.begin()->getWorldIntersectPoint();

            // convert to map coords:
            GeoPoint mapPoint;
            mapPoint.fromWorld( _terrain->getSRS(), world );

            // do an elevation query:
            double query_resolution = 0; // max.
            double out_hamsl        = 0.0;
            double out_resolution   = 0.0;

            bool ok = _query.getElevation( 
                mapPoint,
                out_hamsl,
                query_resolution, 
                &out_resolution );

            if ( ok )
            {
                // convert to geodetic to get the HAE:
                mapPoint.z() = out_hamsl;
                GeoPoint mapPointGeodetic( s_mapNode->getMapSRS()->