PosVelAtt.hpp

Go to the documentation of this file.

// This file is part of INSTINCT, the INS Toolkit for Integrated
// Navigation Concepts and Training by the Institute of Navigation of
// the University of Stuttgart, Germany.
//
// This Source Code Form is subject to the terms of the Mozilla Public
// License, v. 2.0. If a copy of the MPL was not distributed with this
// file, You can obtain one at https://mozilla.org/MPL/2.0/.
 
 
#pragma once
 
#include "Navigation/Transformations/CoordinateFrames.hpp"
 
#include "Navigation/Transformations/Units.hpp"
#include "NodeData/State/PosVel.hpp"
 
namespace NAV
{
class PosVelAtt : public PosVel
{
  public:
    [[nodiscard]] static std::string type()
    {
        return "PosVelAtt";
    }
 
    [[nodiscard]] std::string getType() const override { return type(); }
 
    [[nodiscard]] static std::vector<std::string> parentTypes()
    {
        auto parent = PosVel::parentTypes();
        parent.push_back(PosVel::type());
        return parent;
    }
 
    [[nodiscard]] static std::vector<std::string> GetStaticDataDescriptors()
    {
        auto desc = PosVel::GetStaticDataDescriptors();
        desc.reserve(GetStaticDescriptorCount());
        desc.emplace_back("Roll [deg]");
        desc.emplace_back("Pitch [deg]");
        desc.emplace_back("Yaw [deg]");
        desc.emplace_back("Quaternion::w");
        desc.emplace_back("Quaternion::x");
        desc.emplace_back("Quaternion::y");
        desc.emplace_back("Quaternion::z");
        return desc;
    }
 
    [[nodiscard]] static constexpr size_t GetStaticDescriptorCount() { return PosVel::GetStaticDescriptorCount() + 7; }
 
    [[nodiscard]] std::vector<std::string> staticDataDescriptors() const override { return GetStaticDataDescriptors(); }
 
    [[nodiscard]] size_t staticDescriptorCount() const override { return GetStaticDescriptorCount(); }
 
    [[nodiscard]] std::optional<double> getValueAt(size_t idx) const override
    {
        INS_ASSERT(idx < GetStaticDescriptorCount());
        if (idx < PosVel::GetStaticDescriptorCount()) { return PosVel::getValueAt(idx); }
        switch (idx)
        {
        case PosVel::GetStaticDescriptorCount() + 0: // Roll [deg]
            if (_e_Quat_b.norm() != 0) { return rad2deg(rollPitchYaw().x()); }
            break;
        case PosVel::GetStaticDescriptorCount() + 1: // Pitch [deg]
            if (_e_Quat_b.norm() != 0) { return rad2deg(rollPitchYaw().y()); }
            break;
        case PosVel::GetStaticDescriptorCount() + 2: // Yaw [deg]
            if (_e_Quat_b.norm() != 0) { return rad2deg(rollPitchYaw().z()); }
            break;
        case PosVel::GetStaticDescriptorCount() + 3: // Quaternion::w
            if (_e_Quat_b.norm() != 0) { return n_Quat_b().w(); }
            break;
        case PosVel::GetStaticDescriptorCount() + 4: // Quaternion::x
            if (_e_Quat_b.norm() != 0) { return n_Quat_b().x(); }
            break;
        case PosVel::GetStaticDescriptorCount() + 5: // Quaternion::y
            if (_e_Quat_b.norm() != 0) { return n_Quat_b().y(); }
            break;
        case PosVel::GetStaticDescriptorCount() + 6: // Quaternion::z
            if (_e_Quat_b.norm() != 0) { return n_Quat_b().z(); }
            break;
        default:
            return std::nullopt;
        }
        return std::nullopt;
    }
 
    [[nodiscard]] bool setValueAt(size_t idx, double value) override
    {
        INS_ASSERT(idx < GetStaticDescriptorCount());
        if (idx < PosVel::GetStaticDescriptorCount()) { return PosVel::setValueAt(idx, value); }
        // switch (idx)
        // {
        // case PosVel::GetStaticDescriptorCount() + 0: // Roll [deg]
        // case PosVel::GetStaticDescriptorCount() + 1: // Pitch [deg]
        // case PosVel::GetStaticDescriptorCount() + 2: // Yaw [deg]
        // case PosVel::GetStaticDescriptorCount() + 3: // Quaternion::w
        // case PosVel::GetStaticDescriptorCount() + 4: // Quaternion::x
        // case PosVel::GetStaticDescriptorCount() + 5: // Quaternion::y
        // case PosVel::GetStaticDescriptorCount() + 6: // Quaternion::z
        // default:
        // }
 
        return false;
    }
 
    /* -------------------------------------------------------------------------------------------------------- */
    /*                                           Rotation Quaternions                                           */
    /* -------------------------------------------------------------------------------------------------------- */
 
    [[nodiscard]] const Eigen::Quaterniond& n_Quat_b() const
    {
        return _n_Quat_b;
    }
 
    [[nodiscard]] Eigen::Quaterniond b_Quat_n() const
    {
        return n_Quat_b().conjugate();
    }
 
    [[nodiscard]] const Eigen::Quaterniond& e_Quat_b() const
    {
        return _e_Quat_b;
    }
 
    [[nodiscard]] Eigen::Quaterniond b_Quat_e() const
    {
        return e_Quat_b().conjugate();
    }
 
    [[nodiscard]] Eigen::Vector3d rollPitchYaw() const
    {
        return trafo::quat2eulerZYX(n_Quat_b());
    }
 
    // ###########################################################################################################
    //                                                  Setter
    // ###########################################################################################################
 
    template<typename Derived>
    void setAttitude_e_Quat_b(const Eigen::QuaternionBase<Derived>& e_Quat_b)
    {
        _e_Quat_b = e_Quat_b;
        _n_Quat_b = n_Quat_e() * e_Quat_b;
    }
 
    template<typename Derived>
    void setAttitude_n_Quat_b(const Eigen::QuaternionBase<Derived>& n_Quat_b)
    {
        _e_Quat_b = e_Quat_n() * n_Quat_b;
        _n_Quat_b = n_Quat_b;
    }
 
    template<typename DerivedP, typename DerivedV, typename DerivedA>
    void setState_e(const Eigen::MatrixBase<DerivedP>& e_position, const Eigen::MatrixBase<DerivedV>& e_velocity, const Eigen::QuaternionBase<DerivedA>& e_Quat_b)
    {
        setPosition_e(e_position);
        setVelocity_e(e_velocity);
        setAttitude_e_Quat_b(e_Quat_b);
    }
 
    template<typename DerivedP, typename DerivedV, typename DerivedA>
    void setState_n(const Eigen::MatrixBase<DerivedP>& lla_position, const Eigen::MatrixBase<DerivedV>& n_velocity, const Eigen::QuaternionBase<DerivedA>& n_Quat_b)
    {
        setPosition_lla(lla_position);
        setVelocity_n(n_velocity);
        setAttitude_n_Quat_b(n_Quat_b);
    }
 
    template<typename DerivedP, typename DerivedP2, typename DerivedV, typename DerivedV2, typename DerivedA>
    void setStateAndStdDev_e(const Eigen::MatrixBase<DerivedP>& e_position, const Eigen::MatrixBase<DerivedP2>& e_positionCovarianceMatrix,
                             const Eigen::MatrixBase<DerivedV>& e_velocity, const Eigen::MatrixBase<DerivedV2>& e_velocityCovarianceMatrix,
                             const Eigen::QuaternionBase<DerivedA>& e_Quat_b)
    {
        setPositionAndStdDev_e(e_position, e_positionCovarianceMatrix);
        setVelocityAndStdDev_e(e_velocity, e_velocityCovarianceMatrix);
        setAttitude_e_Quat_b(e_Quat_b);
    }
 
    template<typename DerivedP, typename DerivedP2, typename DerivedV, typename DerivedV2, typename DerivedA>
    void setStateAndStdDev_n(const Eigen::MatrixBase<DerivedP>& lla_position, const Eigen::MatrixBase<DerivedP2>& n_positionCovarianceMatrix,
                             const Eigen::MatrixBase<DerivedV>& n_velocity, const Eigen::MatrixBase<DerivedV2>& n_velocityCovarianceMatrix,
                             const Eigen::QuaternionBase<DerivedA>& n_Quat_b)
    {
        setPositionAndStdDev_lla(lla_position, n_positionCovarianceMatrix);
        setVelocityAndStdDev_n(n_velocity, n_velocityCovarianceMatrix);
        setAttitude_n_Quat_b(n_Quat_b);
    }
 
    /* -------------------------------------------------------------------------------------------------------- */
    /*                                             Member variables                                             */
    /* -------------------------------------------------------------------------------------------------------- */
 
  private:
    Eigen::Quaterniond _e_Quat_b{ 0, 0, 0, 0 };
    Eigen::Quaterniond _n_Quat_b{ 0, 0, 0, 0 };
};
 
} // namespace NAV