osg model

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();
}