About the rotation of the Earth #2397

luohuiiqng · 2023-12-09T13:19:13Z

luohuiiqng
Dec 9, 2023

When I set the rotation of the Earth, the render appears with white stripes. Is this problem caused by my use of GDALImageLayer to load elevation data? Because I found relevant information that GDALImageLayer is static. If it is caused by this problem, what kind of driver should I use to load tif elevation data? Thank you!

osgEarth::initialize();
m_viewer = this->getOsgViewer();
_terrainLayer = new osgEarth::GDALImageLayer();
_terrainLayer->setURL("D:/OSG/osgearth-master/data/world.tif");
_featureSource = new OGRFeatureSource;
_featureSource->setURL("D:/hqluo/osg/osgearth-3.2/data/boston_buildings_utm19.shp");
if (!_featureSource->open().isOK()) {
qDebug() << "OGRFeatureSource...";
return;
}
_featureModelLayer = new FeatureModelLayer;
_featureModelLayer->setFeatureSource(_featureSource);
_map = new osgEarth::Map;
_map->addLayer(_featureModelLayer);
_map->addLayer(_terrainLayer);
ref_ptr earthMt = new MatrixTransform;
_mapEarthNode = new osgEarth::MapNode(_map);
earthMt->addCullCallback(new EarthRotateCallBack(earthMt));
_modelGroupNodel = new osg::Group;
_modelGroupNodel->addChild(_mapEarthNode);
earthMt->addChild(_modelGroupNodel);
_mapNode = new osgEarth::MapNode;
_mapNode->addChild(earthMt);
Viewpoint vp("home", -90.335, 19, -1376, 0, -90, 10000000);
m_mainCamera->setViewpoint(vp);
m_viewer->setCameraManipulator(m_mainCamera);
m_viewer->setSceneData(_mapNode.get());

Answered by gwaldron

Dec 12, 2023

I still see two calls to new osgEarth::MapNode in there.

View full answer

luohuiiqng · 2023-12-09T13:53:37Z

luohuiiqng
Dec 9, 2023
Author

This is code of rotateCallback.

If it were rendered properly, it would look something like this

0 replies

gwaldron · 2023-12-09T17:08:05Z

gwaldron
Dec 9, 2023
Maintainer

Can you post a compilable .cpp file that we can use to reproduce the problem?
Also, please post the output of osgearth_version --caps.
Thanks.

10 replies

gwaldron Dec 11, 2023
Maintainer

You also appear to two MapNode's in your scene. Was that your intention?
(_mapNode -> earthMt -> modelGroupNodel -> _mapEarthNode)

luohuiiqng Dec 12, 2023
Author

No, only one!
if so,it's worng.

void ViewerWidget::initWindow()
{
    osgEarth::initialize();    
    m_viewer = this->getOsgViewer(); 

    _terrainLayer = new osgEarth::GDALImageLayer();
    _terrainLayer->setURL("D:/OSG/osgearth-master/data/world.tif");
    _featureSource = new OGRFeatureSource;
    _featureSource->setURL("D:/hqluo/osg/osgearth-3.2/data/boston_buildings_utm19.shp");
    if (!_featureSource->open().isOK()) {
        qDebug() << "OGRFeatureSource...";
        return;
    }
    _featureModelLayer = new FeatureModelLayer;
    _featureModelLayer->setFeatureSource(_featureSource);
    _map = new osgEarth::Map;      
    _map->addLayer(_featureModelLayer);
    _map->addLayer(_terrainLayer);
    
    MatrixTransform* earthMt = new MatrixTransform;      
    _mapEarthNode = new osgEarth::MapNode(_map);   
    
    earthMt->addCullCallback(new EarthRotateCallBack(earthMt));                            
    earthMt->addChild(_mapEarthNode);

    _mapNode = new osgEarth::MapNode;
    _mapNode->addChild(earthMt);

    Viewpoint vp("home", -90.335, 19, -1376, 0, -90, 10000000);
    m_mainCamera->setViewpoint(vp);
    m_mainCamera->setAutoComputeHomePosition(true);
    m_viewer->setCameraManipulator(m_mainCamera);
    m_viewer->setSceneData(_mapNode);
}

If so, it is correct, but you cannot use MapNode as data.

void ViewerWidget::initWindow()
{
    osgEarth::initialize();    
    m_viewer = this->getOsgViewer(); 
   
    _terrainLayer = new osgEarth::GDALImageLayer();
    _terrainLayer->setURL("D:/OSG/osgearth-master/data/world.tif");
    _featureSource = new OGRFeatureSource;
    _featureSource->setURL("D:/hqluo/osg/osgearth-3.2/data/boston_buildings_utm19.shp");
    if (!_featureSource->open().isOK()) {
        qDebug() << "OGRFeatureSource...";
        return;
    }
    _featureModelLayer = new FeatureModelLayer;
    _featureModelLayer->setFeatureSource(_featureSource);

    _map = new osgEarth::Map;   
    _map->addLayer(_featureModelLayer);
    _map->addLayer(_terrainLayer);
    
    MatrixTransform* earthMt = new MatrixTransform;              
    earthMt->addCullCallback(new EarthRotateCallBack(earthMt));                
    _modelGroupNodel = new osg::Group;
       
    _mapEarthNode = new osgEarth::MapNode(_map);
    earthMt->addChild(_mapEarthNode);

    /*
        If you annotate the following two lines, the earth stays still.
    */
    _mapNode = new osgEarth::MapNode;
    _mapNode->addChild(earthMt);

    Viewpoint vp("home", -90.335, 19, -1376, 0, -90, 10000000);
    m_mainCamera->setViewpoint(vp);
    m_mainCamera->setAutoComputeHomePosition(true);
    m_viewer->setCameraManipulator(m_mainCamera);

    m_viewer->setSceneData(earthMt);

}

luohuiiqng Dec 12, 2023
Author

My goal is to achieve the rotation of the Earth and manage the scene with MapNode. Do you have a better implementation? thank you.

gwaldron Dec 12, 2023
Maintainer

I still see two calls to new osgEarth::MapNode in there.

Answer selected by gwaldron

luohuiiqng Dec 13, 2023
Author

Yes, I'm sorry. I was wrong.
By looking at the code examples and your tips, I have already implemented the rotation function. However, what is confusing to me now is that the sun will revolve around the Earth and I will not be able to see the starry sky, the same code that shows the starry sky normally in programs that do not include qt. Is that normal?

class ECIReferenceFrame : public osg::MatrixTransform
{
public:
    ECIReferenceFrame()
    {
        Lighting::set(getOrCreateStateSet(), osg::StateAttribute::OFF);
    }

    void traverse(osg::NodeVisitor& nv)
    {
        osgUtil::CullVisitor* cv = Culling::asCullVisitor(nv);
        if (cv)
        {
            const osg::Camera* cam = cv->getRenderStage()->getCamera();
            cv->pushModelViewMatrix(new osg::RefMatrix(cam->getViewMatrix()), osg::Transform::ABSOLUTE_RF);
            osg::Group::traverse(nv);
            cv->popModelViewMatrix();
        }
        else osg::Group::traverse(nv);
    }
};


ref_ptr<SkyNode> ViewerWidget::addSky()
{
    ref_ptr<SkyNode> sky_node = SkyNode::create();
    SkyOptions skyOptions;
    skyOptions.coordinateSystem() = SkyOptions::COORDSYS_ECI;
    sky_node->attach(m_viewer);
    sky_node->setName("sky");   
    sky_node->getSunLight()->setAmbient(osg::Vec4(0.5, 0.5, 0.5, 1.0));

    return sky_node;
}

void ViewerWidget::initWindow()
{
    osgEarth::initialize();    
    m_viewer = this->getOsgViewer(); 
   
    _terrainLayer = new osgEarth::GDALImageLayer();
    _terrainLayer->setURL("D:/OSG/osgearth-master/data/world.tif");
    _featureSource = new OGRFeatureSource;
    _featureSource->setURL("D:/hqluo/osg/osgearth-3.2/data/boston_buildings_utm19.shp");
    if (!_featureSource->open().isOK()) {
        qDebug() << "OGRFeatureSource...";
        return;
    }
    _featureModelLayer = new FeatureModelLayer;
    _featureModelLayer->setFeatureSource(_featureSource);

    _map = new osgEarth::Map;   
    m_GRoot = new osg::Group;
    
    _map->addLayer(_featureModelLayer);
    _map->addLayer(_terrainLayer);
    
    ECIReferenceFrame* earthMt = new ECIReferenceFrame;
    earthMt->addCullCallback(new EarthRotateCallBack(earthMt));                
    _modelGroupNodel = new osg::Group;
       
    _mapEarthNode = new osgEarth::MapNode(_map);
    earthMt->addChild(MapNode::get(_mapEarthNode));
    _modelGroupNodel->addChild(earthMt);

    ref_ptr<SkyNode> skyNode = addSky();
    m_GRoot->addChild(_modelGroupNodel);
    m_GRoot->addChild(skyNode);
    Viewpoint vp("home", -90.335, 19, -1376, 0, -90, 10000000);
    m_mainCamera->setViewpoint(vp);
    m_viewer->setCameraManipulator(m_mainCamera);

    m_viewer->setSceneData(m_GRoot);

}

qt not include:

#include <osgEarth/MapNode>
#include <osgEarth/DateTime>
#include <osgEarth/NodeUtils>
#include <osgEarth/PointDrawable>
#include <osgEarth/CullingUtils>
#include <osgEarth/LineDrawable>
#include <osgEarth/Lighting>
#include <osgEarth/EarthManipulator>
#include <osgEarth/ExampleResources>
#include <osgEarth/Sky>
#include <osgEarth/Color>
#include <osgEarth/LabelNode>
#include <osgEarth/GDAL>
#include <osgEarth/TMS>
#include <osgEarth/MBTiles>
#include <osgEarth/WMS>
#include <osgEarth/XYZ>
#include <osgEarth/Composite>
#include <osgEarth/ContourMap>
#include <osgEarth/Style>
#include <osgEarth/FeatureModelLayer>
#include <osgEarth/FeatureImageLayer>
#include <osgEarth/AnnotationLayer>
#include <osgEarth/OGRFeatureSource>
#include <osgViewer/Viewer>
#include <iostream>


using namespace osgEarth;
using namespace osgEarth::Util;
using namespace osg;

class EarthRotateCallBack : public osg::NodeCallback
{
public:
    EarthRotateCallBack(ref_ptr<MatrixTransform>);
    virtual void operator()(osg::Node* node, osg::NodeVisitor* nv);
private:
    double _lastTimeStamp;
    double _duration;
    double _rotation;
    ref_ptr<MatrixTransform> m_mt;
    double g_timeScale = 100000;//时间放大倍数，程序时间乘以该倍数模拟显示时间，越大卫星速度越快
    const double g_daySeconds = 24 * 3600;//一天的秒数
};


EarthRotateCallBack::EarthRotateCallBack(ref_ptr<MatrixTransform> mt)
{
    _lastTimeStamp = 0;
    _duration = 0;
    _rotation = 0;
    m_mt = mt;
}

void EarthRotateCallBack::operator()(osg::Node* node, osg::NodeVisitor* nv)
{
    double currTimeStamp = nv->getFrameStamp()->getReferenceTime();
    double dura = currTimeStamp - _lastTimeStamp;
    _lastTimeStamp = currTimeStamp;
    dura *= g_timeScale;
    _duration += dura;
    if (_duration > g_daySeconds) {
        _duration = _duration - g_daySeconds;
    }
    _rotation = 2 * 3.1415926 * (_duration / g_daySeconds);

    osg::Matrix mat = osg::Matrix::rotate(Quat(_rotation, Vec3(0, 0, 1)));
    m_mt->setMatrix(mat);

    traverse(node, nv);
}



class ECIReferenceFrame : public osg::MatrixTransform
{
public:
    ECIReferenceFrame()
    {
        Lighting::set(getOrCreateStateSet(), osg::StateAttribute::OFF);
    }

    void traverse(osg::NodeVisitor& nv)
    {
        osgUtil::CullVisitor* cv = Culling::asCullVisitor(nv);
        if (cv)
        {
            const osg::Camera* cam = cv->getRenderStage()->getCamera();
            cv->pushModelViewMatrix(new osg::RefMatrix(cam->getViewMatrix()), osg::Transform::ABSOLUTE_RF);
            osg::Group::traverse(nv);
            cv->popModelViewMatrix();
        }
        else osg::Group::traverse(nv);
    }
};


int main(int argc, char** argv)
{
    osgEarth::initialize();
    osgViewer::Viewer m_viewer;

    osgEarth::GDALImageLayer* _terrainLayer = new osgEarth::GDALImageLayer();
    _terrainLayer->setURL("D:/OSG/osgearth-master/data/world.tif");
    OGRFeatureSource* _featureSource = new OGRFeatureSource;
    _featureSource->setURL("D:/hqluo/osg/osgearth-3.2/data/boston_buildings_utm19.shp");
    FeatureModelLayer* _featureModelLayer = new FeatureModelLayer;
    _featureModelLayer->setFeatureSource(_featureSource);

    osgEarth::Map* _map = new osgEarth::Map;
    osg::Group* m_GRoot = new osg::Group;

    _map->addLayer(_featureModelLayer);
    _map->addLayer(_terrainLayer);

    ECIReferenceFrame* earthMt = new ECIReferenceFrame;
    earthMt->addCullCallback(new EarthRotateCallBack(earthMt));
    osg::Group* _modelGroupNodel = new osg::Group;

    MapNode* _mapEarthNode = new osgEarth::MapNode(_map);
    earthMt->addChild(MapNode::get(_mapEarthNode));
    _modelGroupNodel->addChild(earthMt);


    osgEarth::DateTime dateTime(2023, 12, 12, 0);
    ref_ptr<SkyNode> sky_node = SkyNode::create();
    SkyOptions skyOptions;
    skyOptions.coordinateSystem() = SkyOptions::COORDSYS_ECI;
    sky_node->attach(&m_viewer);
    sky_node->setName("sky");
    sky_node->getSunLight()->setAmbient(osg::Vec4(0.5, 0.5, 0.5, 1.0));

    m_GRoot->addChild(_modelGroupNodel);
    m_GRoot->addChild(sky_node);
    Viewpoint vp("home", -90.335, 19, -1376, 0, -90, 10000000);
    EarthManipulator* m_mainCamera = new EarthManipulator;
    m_mainCamera->setViewpoint(vp);

    m_viewer.setCameraManipulator(m_mainCamera);

    m_viewer.setSceneData(m_GRoot);

    return m_viewer.run();
}

gwaldron · 2023-12-13T15:28:06Z

gwaldron
Dec 13, 2023
Maintainer

Marking this as Answered since the original question about a display anomaly is resolved (2 MapNodes in the scene).

As for your Qt question, I have no idea, we do not have much Qt expertise.
osgEarth by default runs in an ECEF reference frame. You can refer to the osgeath_eci.cpp example for code that rotates the earth into an ECI reference frame.

4 replies

luohuiiqng Dec 14, 2023
Author

Ok, thank you very much!

ZzWDlyy Dec 16, 2023

请问你解决星空的问题了嘛，我也是用的qt看不到星空只有太阳。

luohuiiqng Dec 17, 2023
Author

没有呢，后面打算自己写一个天空盒。但目前还没开始，在学习其它的模块中。

luohuiiqng Aug 14, 2024
Author

    osgEarth::SimpleSky::SimpleSkyOptions skyOptions;
    skyOptions.coordinateSystem() = osgEarth::SkyOptions::COORDSYS_ECEF;
    //Star texture
    skyOptions.starImageURI() = QApplication::applicationDirPath().append("/Resources/Image/galaxy_starfield.jpg").toStdString();
    skyOptions.starSize() = 0.0f;
    osgEarth::SkyNode* skyNode = osgEarth::Util::SkyNode::create(skyOptions);

Provide feedback

Saved searches

Use saved searches to filter your results more quickly

About the rotation of the Earth #2397

{{title}}

Replies: 3 comments 14 replies

{{title}}

{{title}}

{{title}}

{{title}}

{{title}}

{{editor}}'s edit

{{editor}}'s edit

{{title}}

{{title}}

{{title}}

{{title}}

{{title}}

{{title}}

{{title}}

Select a reply

About the rotation of the Earth #2397

luohuiiqng Dec 9, 2023

Replies: 3 comments · 14 replies

luohuiiqng Dec 9, 2023 Author

gwaldron Dec 9, 2023 Maintainer

gwaldron Dec 11, 2023 Maintainer

luohuiiqng Dec 12, 2023 Author

luohuiiqng Dec 12, 2023 Author

gwaldron Dec 12, 2023 Maintainer

luohuiiqng Dec 13, 2023 Author

gwaldron Dec 13, 2023 Maintainer

luohuiiqng Dec 14, 2023 Author

ZzWDlyy Dec 16, 2023

luohuiiqng Dec 17, 2023 Author

luohuiiqng Aug 14, 2024 Author

luohuiiqng
Dec 9, 2023

Replies: 3 comments 14 replies

luohuiiqng
Dec 9, 2023
Author

gwaldron
Dec 9, 2023
Maintainer

gwaldron Dec 11, 2023
Maintainer

luohuiiqng Dec 12, 2023
Author

luohuiiqng Dec 12, 2023
Author

gwaldron Dec 12, 2023
Maintainer

luohuiiqng Dec 13, 2023
Author

gwaldron
Dec 13, 2023
Maintainer

luohuiiqng Dec 14, 2023
Author

luohuiiqng Dec 17, 2023
Author

luohuiiqng Aug 14, 2024
Author