osgEarth設置模型旋轉角度
阿新 • • 發佈：2017-11-01
axis ima 軌跡 tar child test orm 獲取 art
#include<windows.h> #include <osgViewer/Viewer> #include <osgEarthDrivers/gdal/GDALOptions> #include <osg/ShapeDrawable> #include <osgEarthUtil/EarthManipulator> #include <osg/MatrixTransform> #include <osgEarthFeatures/ConvertTypeFilter> #include <osgEarthDrivers/model_simple/SimpleModelOptions> #include<osgViewer/ViewerEventHandlers> #include <osgGA/StateSetManipulator> #ifdef _DEBUG #pragma comment(lib, "osgd.lib") //#pragma comment(lib, "osgDBd.lib") #pragma comment(lib, "osgViewerd.lib") #pragma comment(lib, "osgGAd.lib") #pragma comment(lib, "osgEarthd.lib") #pragma comment(lib, "osgEarthUtild.lib") #else#pragma comment(lib, "osg.lib") #pragma comment(lib, "osgDB.lib") #pragma comment(lib, "osgViewer.lib") #pragma comment(lib, "osgGA.lib") #endif // DEBUG // lonlat1:雷達波圓錐頂點 // lonlat2:軌跡點 void rotateCone(osg::MatrixTransform* mt, const osgEarth::SpatialReference* sr, osg::Vec3d lonlat1, osg::Vec3d lonlat2) {// 雷達波模型所在位置 osgEarth::GeoPoint geoPoint1( sr, lonlat1.x(), lonlat1.y(), lonlat1.z(), osgEarth::ALTMODE_ABSOLUTE); osg::Matrixd matrix1; // 獲取雷達波模型從原點變換到lonlat1位置的變換矩陣 geoPoint1.createLocalToWorld(matrix1); // 經緯度高程到xyz的變換 osg::Vec3d world1, world2; // geoPoint1.toWorld(world1);//等同於 sr->transformToWorld(lonlat1,world1); sr->transformToWorld(lonlat2, world2); // 計算軌跡點在雷達波模型坐標系下的位置 osg::Vec3 point2InRadarCoordinateSystem = world2*osg::Matrix::inverse(matrix1); // 在雷達波模型坐標系下，對Z軸進行旋轉，與連接原點指向軌跡點方向的矢量重合，計算出此旋轉矩陣 osg::Matrixd rotMat = osg::Matrixd::rotate(osg::Z_AXIS, point2InRadarCoordinateSystem-osg::Vec3(0,0,0)); // 將計算出的旋轉矩陣賦給雷達波模型所在的mt mt->setMatrix(rotMat); } int main(int argc, char** argv) { osgViewer::Viewer viewer; std::string world_tif = "data/world.tif"; osgEarth::Map* map = new osgEarth::Map(); // Start with a basemap imagery layer; we‘ll be using the GDAL driver // to load a local GeoTIFF file: osgEarth::Drivers::GDALOptions basemapOpt; basemapOpt.url() = world_tif; map->addImageLayer(new osgEarth::ImageLayer(osgEarth::ImageLayerOptions("basemap", basemapOpt))); osgEarth::MapNodeOptions mapNodeOptions; mapNodeOptions.enableLighting() = false; osgEarth::MapNode* mapNode = new osgEarth::MapNode(map, mapNodeOptions); osgEarth::Drivers::SimpleModelOptions opt; opt.location() = osg::Vec3(118, 40, 10000); //opt.url() = "cow.osg.1000,1000,1000.scale"; osg::Geode* geode = new osg::Geode; osg::ShapeDrawable* cone = new osg::ShapeDrawable(new osg::Cone(osg::Vec3(), 10000, 50000)); //osg::ShapeDrawable* cone = new osg::ShapeDrawable(new osg::Box(osg::Vec3(), 50000)); geode->addDrawable(cone); osg::MatrixTransform* mtCone = new osg::MatrixTransform; mtCone->addChild(geode); opt.node() = mtCone; map->addModelLayer(new osgEarth::ModelLayer("", opt)); rotateCone(mtCone, map->getProfile()->getSRS(), osg::Vec3(118, 40, 100), osg::Vec3(120, 40, 100)); osg::Group* root = new osg::Group(); root->addChild(mapNode); viewer.setSceneData(root); viewer.setCameraManipulator(new osgEarth::Util::EarthManipulator()); // Process cmdline args //MapNodeHelper().parse(mapNode, arguments, &viewer, root, new LabelControl("Features Demo")); //視點定位模型所在位置 osgEarth::Viewpoint vp("", 118, 40, 1000.0, -2.50, -90.0, 1.5e6); (dynamic_cast<osgEarth::Util::EarthManipulator*>(viewer.getCameraManipulator()))->setViewpoint(vp); // add some stock OSG handlers: viewer.addEventHandler(new osgViewer::StatsHandler()); viewer.addEventHandler(new osgViewer::WindowSizeHandler()); viewer.addEventHandler(new osgGA::StateSetManipulator(viewer.getCamera()->getOrCreateStateSet())); viewer.setUpViewInWindow(100, 500, 1024, 768); return viewer.run(); }
osgEarth設置模型旋轉角度