添加光照
osg::ref_ptr<osg::Node> MyOSGLoadEarth::CreateNode() { osg::ref_ptr<osg::Group> _root = new osg::Group; //定义并读取高程文件 //真实高程文件名称为:ASTGTM2_N34E110_dem.tif //属于特殊的tiff格式,GEOTiff //读取的时候使用osg的gdal插件进行读取,所以在路径后面加上了.gdal //.gdal后缀名只要在这里加就可以了,真实的高程文件后缀名不需要修改 osg::ref_ptr<osg::HeightField> heightMap = osgDB::readHeightFieldFile("G:\\ASTER.GDEM.V2-DEM_sn\\ASTGTM2_N34E110_dem.tif.gdal"); //创建一个叶结点对象 osg::ref_ptr<osg::Geode> geode = new osg::Geode; if (heightMap != nullptr) { //由于原始数据过大,创建三维对象会失败,所以重新构造一个对象 //相当于数据抽稀了一次。当然,可以直接把原图使用特殊工具裁了 //创建一个新的HeightField对象,用来拷贝heightMap osg::ref_ptr<osg::HeightField> heightMap1 = new osg::HeightField; //从原对象中拷贝一些熟悉过来 heightMap1->setOrigin(heightMap->getOrigin()); heightMap1->setRotation(heightMap->getRotation()); heightMap1->setSkirtHeight(heightMap->getSkirtHeight()); //XY方向的间隔设置为原来的两倍, heightMap1->setXInterval(heightMap->getXInterval() * 2); heightMap1->setYInterval(heightMap->getYInterval() * 2); //设置新的高程数据量的行列数目为原来的一半 heightMap1->allocate(heightMap->getNumColumns()/2,heightMap->getNumRows()/2); //把真实的数据值放进来 for (size_t r = 0; r < heightMap1->getNumRows(); ++r) { for (size_t c = 0; c < heightMap1->getNumColumns();++c) { //加载的数据中XY方向的间隔是0.0002左右(经纬度偏移),3600个格子,数量级太小,高程值动辄在千级别,如果没有进行坐标转换(GPS转换成米),显示出来之后结果会严重失常。所以此处简单的给高度值除以50000(这个是按照这个tif文件来试出来的,不同高程文件可能不同) heightMap1->setHeight(c, r, heightMap->getHeight(c * 2, r * 2)/50000); } } //添加到叶子节点中 geode->addDrawable(new osg::ShapeDrawable(heightMap1)); } _root->addChild(geode.get()); return _root.get(); }
参照:https://blog.csdn.net/woaicd/article/details/77163012