2015年3月21日
16:19
osgEarth是基于三维引擎osg开发的三维数字地球引擎库,在osg基础上实现了瓦片调度插件,可选的四叉树调度插件,更多的地理数据加载插件(包括GDAL,ogr,WMS,TMS,VPB,filesystem等),再结合一套地理投影转换插件,这样就能够实现高效处理加载调度地理数据在三维地球上的显示,实现三维虚拟地球。
想要实现一个简单的基于osgEarth的三维地球,有两种方式,这两种方式是互通的。一种基于XML标签的earth文件加载,另外一种是采用C++代码,本质是一样的,osgEarth内部支持直接解析XML标签文件,转换为代码实现,具体参考tests文件夹例子,代码则参考application下面例子。但是大多数情况下,你需要更灵活掌控性更强的代码来实现功能,因为只有一个加载了基础数据的三维地球是只能看,不能解决实际问题的,需要界面通常采用QT,更多的三维渲染和仿真业务则由osg来完成。因此学好osg是做这一切的基础。
osgEarth的特点:支持加载常见的栅格数据(影像和DEM),但是大数据必须建立金字塔,设置为地理投影,想要高效率最好处理为瓦片,这样也便于部署在服务端。矢量数据,最好尽可能的简化,因为大的矢量会十分影响渲染速度,当然也可以对矢量栅格化处理加快速度,对于模型的话,大数据量一定要做LOD或者pageLod。
osgEarth程序的常规流程:
创建osgViewer---->创建MapNode---->设置Earth操作器---->设置场景参数----->run
MapNode是继承自osg的Node,是osgEarth中地球节点,你所添加的影像,DEM,模型都包含在MapNode中,因为它们都加入到Map中,Map则类似二维中的Map可以添加各种图层。剩余的不管是模型节点Node,或者是标注Node,还是其他的都是可以直接添加到MapNode中或者另外的Group中。
Earth操作器则和其他osg操作器一样,只不过专门为三维地球浏览定制,具体参数可以设置。
场景参数则主要有自动地形裁剪,最小裁剪像素等其他优化场景的参数。
下面就简单阐述一个小例子说明:
代码功能主要实现了查询实时高程,并显示XYZ坐标的功能。
使用命令app.exe test.earth即可得到下面的效果。
//引入osg和osgEarth的头文件和命名空间 #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 <osgEarthUtil/EarthManipulator> #include <osgEarthUtil/Controls> #include <osgEarthUtil/LatLongFormatter> #include <iomanip> using namespace osgEarth; using namespace osgEarth::Util; using namespace osgEarth::Util::Controls; 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_mapLabel = 0L; static LabelControl* s_resLabel = 0L; // An event handler that will print out the elevation at the clicked point //查询高程的一个事件回调,在场景有事件更新触发时调用,详细参考osg或者osgGA::GUIEventHandler struct QueryElevationHandler : public osgGA::GUIEventHandler { //构造函数 QueryElevationHandler() : _mouseDown( false ), _terrain ( s_mapNode->getTerrain() ), _query ( s_mapNode->getMap() ) { _map = s_mapNode->getMap(); //初始化最大查询LOD级别 _query.setMaxTilesToCache(10); _path.push_back( s_mapNode->getTerrainEngine() ); } //更新回调,具体的内容可以参考父类,传进来的参数是屏幕坐标xy,和osgViewer void update( float x, float y, osgViewer::View* view ) { bool yes = false; // look under the mouse: //采用线去对地球做碰撞检测,根据鼠标点击点去检测,得到交点,就是当前点的xyz osg::Vec3d world; osgUtil::LineSegmentIntersector::Intersections hits; //判断求交结果是否为空 if ( view->computeIntersections(x, y, hits) ) { //得到世界坐标系下面的坐标,就是osg的xyz坐标 world = hits.begin()->getWorldIntersectPoint(); // convert to map coords: //将其转换为地球的地理坐标,转换方法都照抄即可 GeoPoint mapPoint; mapPoint.fromWorld( _terrain->getSRS(), world ); // do an elevation query: double query_resolution = 0; // 1/10th of a degree 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()->getGeodeticSRS(), mapPoint ); //经纬度坐标的格式化工具,也可以自己用字符串去拼接xyz数字 static LatLongFormatter s_f; //更新显示的xyz值,label是传入的控件 s_posLabel->setText( Stringify() << std::fixed << std::setprecision(2) << s_f.format(mapPointGeodetic.y()) << ", " << s_f.format(mapPointGeodetic.x()) ); //还可以输出分辨率,椭球体信息等 s_mslLabel->setText( Stringify() << out_hamsl ); s_haeLabel->setText( Stringify() << mapPointGeodetic.z() ); s_resLabel->setText( Stringify() << out_resolution ); yes = true; } // finally, get a normal ISECT HAE point. GeoPoint isectPoint; isectPoint.fromWorld( _terrain->getSRS()->getGeodeticSRS(), world ); s_mapLabel->setText( Stringify() << isectPoint.alt() ); } //如果查询不到高程的话 if (!yes) { s_posLabel->setText( "-" ); s_mslLabel->setText( "-" ); s_haeLabel->setText( "-" ); s_resLabel->setText( "-" ); } } //参数一个是事件的动作,另外一个是对应的操作 bool handle( const osgGA::GUIEventAdapter& ea, osgGA::GUIActionAdapter& aa ) { //判断如果是移动鼠标事件才进行更新当前的坐标显示 if (ea.getEventType() == osgGA::GUIEventAdapter::MOVE && aa.asView()->getFrameStamp()->getFrameNumber() % 10 == 0) { osgViewer::View* view = static_cast<osgViewer::View*>(aa.asView()); update( ea.getX(), ea.getY(), view ); } return false; } //Map对象 const Map* _map; //地形对象 const Terrain* _terrain; bool _mouseDown; //查询高程使用的对象 ElevationQuery _query; osg::NodePath _path; }; //main函数, int main(int argc, char** argv) { //这儿两个参数,第一个是命令参数的个数为,后面是字符串数组输入earth文件的路径osg::ArgumentParser arguments(&argc,argv); //osg的场景 osgViewer::Viewer viewer(arguments); //构造MapNode,arguments里面有earth文件的路径,命令行输入 s_mapNode = MapNode::load(arguments); //如果路径不正确或者earth文件错误,没有构造好MapNode if ( !s_mapNode ) { OE_WARN << "Unable to load earth file." << std::endl; return -1; } //建立一个组节点 osg::Group* root = new osg::Group(); //将组节点设置为场景节点 viewer.setSceneData( root ); // install the programmable manipulator. //设置earth操作器 viewer.setCameraManipulator( new osgEarth::Util::EarthManipulator() ); // The MapNode will render the Map object in the scene graph. //将MapNode添加到组节点中去 root->addChild( s_mapNode ); //下面是设置一个控件,grid的意思是用格网去布局里面的小控件 // Make the readout: Grid* grid = new Grid(); //设置几个Label文字控件显示在场景中的第行 grid->setControl(0,0,new LabelControl("Coords (Lat, Long):")); grid->setControl(0,1,new LabelControl("Vertical Datum:")); grid->setControl(0,2,new LabelControl("Height (MSL):")); grid->setControl(0,3,new LabelControl("Height (HAE):")); grid->setControl(0,4,new LabelControl("Isect (HAE):")); grid->setControl(0,5,new LabelControl("Resolution:")); //设置几个Label文字控件显示在场景中的第行 s_posLabel = grid->setControl(1,0,new LabelControl("")); s_vdaLabel = grid->setControl(1,1,new LabelControl("")); s_mslLabel = grid->setControl(1,2,new LabelControl("")); s_haeLabel = grid->setControl(1,3,new LabelControl("")); s_mapLabel = grid->setControl(1,4,new LabelControl("")); s_resLabel = grid->setControl(1,5,new LabelControl("")); //得到空间参考,椭球面信息,并显示对应上面的label const SpatialReference* mapSRS = s_mapNode->getMapSRS(); s_vdaLabel->setText( mapSRS->getVerticalDatum() ? mapSRS->getVerticalDatum()->getName() : Stringify() << "geodetic (" << mapSRS->getEllipsoid()->getName() << ")" ); //控件绘制容器 ControlCanvas* canvas = new ControlCanvas(); //将要显示的控件加入到root组节点中去 root->addChild(canvas); canvas->addControl( grid ); //添加刚刚自定义的查询高程的事件回调 // An event handler that will respond to mouse clicks: viewer.addEventHandler( new QueryElevationHandler() ); //添加状态显示,窗口改变等事件回调 // add some stock OSG handlers: viewer.addEventHandler(new osgViewer::StatsHandler()); viewer.addEventHandler(new osgViewer::WindowSizeHandler()); viewer.addEventHandler(new osgGA::StateSetManipulator(viewer.getCamera()->getOrCreateStateSet())); //run return viewer.run(); }