Eigen库要点「建议收藏」

Eigen库要点「建议收藏」旋转矩阵,旋转向量,四元数关系如下:另外,初始化(赋值)变换矩阵T的方式为:Isometry3dTcw(rotation);//rotation可以是旋转矩阵,可以是四元数,可以是旋转向量Tcw.pretranslate(t);//添加平移向量//或者:Isometry3dTcw=Isometry3d::Identity()//如果没有直接初始化,先设为单位阵Tcw.pre…

大家好,又见面了,我是你们的朋友全栈君。如果您正在找激活码,请点击查看最新教程,关注关注公众号 “全栈程序员社区” 获取激活教程,可能之前旧版本教程已经失效.最新Idea2022.1教程亲测有效,一键激活。

Jetbrains全系列IDE稳定放心使用

 旋转矩阵,旋转向量,四元数关系如下:

Eigen库要点「建议收藏」另外,初始化(赋值)变换矩阵T的方式为:

Isometry3d Tcw(rotation);//rotation可以是旋转矩阵,可以是四元数,可以是旋转向量
Tcw.pretranslate(t);//添加平移向量
//或者:
Isometry3d Tcw=Isometry3d::Identity()//如果没有直接初始化,先设为单位阵
Tcw.prerotate(rotation1);//然后添加旋转矩阵,或者向量,或者四元数
Tcw1.pretranslate(t1);//添加平移向量

 从变换矩阵获取旋转矩阵和平移矩阵:

Tcw.rotation();//返回旋转矩阵
Tcw.translation();//返回旋转矩阵

实例

设有小萝卜一号和小萝卜二号位于世界坐标系中,小萝卜一号的位姿为:q1=[0.35,0.2,0.3,0.1],
t2=[0.3,0.1,0.1]^T (q的第一项为实部。请你把q归一化后在进行计算)。这里的q和t的表达的是Tcw,也就是世界到相机的变换关系。小萝卜二号的位姿为q2=[-0.5,0.4,-0.1,0.2],t=[-0.1,0.5,0.3]^T.现在,小萝卜一号看到某个点在自身的坐标系下,坐标为p=[0.5,0,0.2]^T ,求该向量在小萝卜二号坐标系下的坐标,请编程实现此事。

方法一:直接用四元数

    Quaterniond q1={0.35,0.2,0.3,0.1},q2={-0.5,0.4,-0.1,0.2};
    Vector3d t1={0.3,0.1,0.1},t2={-0.1,0.5,0.3},p={0.5,0,0.2};
    Vector3d world,cam2;
    q1=q1.normalized();
    q2=q2.normalized();
    world=q1.inverse()*(p-t1);
    cam2=q2*world+t2;
    cout<<cam2<<endl;

方法二:使用旋转矩阵

    Matrix3d rotation1=q1.toRotationMatrix();
    Matrix3d rotation2=q2.toRotationMatrix();
    Matrix3d R1,R1_inv,R2;
    world=rotation1.inverse()*(p-t1);
    cam2=rotation2*world+t2;
    cout<<cam2<<endl;

方法三:使用变换矩阵(直接)

  Isometry3d Tcw1(rotation1),Twc1,Tcw2(rotation2),Twc2;//直接初始化,不需要设定单位阵
    Tcw1.pretranslate(t1);
    Tcw2.pretranslate(t2);
    //cout<<Tcw1.matrix()<<endl;
    Twc1=Tcw1.inverse();
    //cout<<Twc1.matrix()<<endl;
    world=Twc1*p;    cout<<rotation1<<endl;
    cam2=Tcw2*world;
    cout<<cam2<<endl;
    return 0;

方法三:使用变换矩阵(直接)

    Isometry3d Tcw1=Isometry3d::Identity(),Twc1,Tcw2=Isometry3d::Identity();//如果没有直接初始化,先设为单位阵
    Tcw1.prerotate(rotation1);
    Tcw2.prerotate(rotation2);
    Tcw1.pretranslate(t1);
    Tcw2.pretranslate(t2);
    //cout<<Tcw1.matrix()<<endl;
    Twc1=Tcw1.inverse();
    //cout<<Twc1.matrix()<<endl;
    world=Twc1*p;
    cam2=Tcw2*world;
    cout<<cam2<<endl;
    return 0;

方法四:间接使用变换矩阵(间接)

    Isometry3d Tcw1=Isometry3d::Identity(),Twc1,Tcw2=Isometry3d::Identity();//如果没有直接初始化,先设为单位阵
    Tcw1.prerotate(rotation1);
    Tcw2.prerotate(rotation2);
    Tcw1.pretranslate(t1);
    Tcw2.pretranslate(t2);
    world=Tcw1.rotation().inverse()*(p-Tcw1.translation());
    cam2=Tcw2.rotation()*world+Tcw2.translation();
    cout<<cam2<<endl;
    return 0;

C++数组和矩阵转换

使用Map函数,可以实现Eigen的矩阵和c++中的数组直接转换:

Map类用于通过C++中普通的连续指针或者数组 (raw C/C++ arrays)来构造Eigen里的Matrix类,这就好比Eigen里的Matrix类的数据和raw C++array 共享了一片地址,也就是引用。

数组转静态矩阵:

int i;
double *aMat = new double[20];
for(i =0;i<20;i++)
{
  aMat[i] = rand()%11;
}
Eigen:Map<Matrix<double,4,5> > staMat(aMat);
virtual void oplusImpl(const double* update) override
{
     Eigen::VectorXd::ConstMapType v(update,VectexCameraBAL::Dimension);
     _estimate+=v;
}

数组转动态矩阵:

int i;
double *aMat = new double[20];
for(i =0;i<20;i++)
{
  aMat[i] = rand()%11;
}
Eigen::Map<MatrixXd> dymMat(aMat,4,5);
virtual void oplusImpl(const double* update) override
{
     Eigen::Vector3d::ConstMapType v(update);
     _estimate+=v;
}

 

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请联系我们举报,一经查实,本站将立刻删除。

发布者:全栈程序员-站长,转载请注明出处:https://javaforall.net/180433.html原文链接:https://javaforall.net

(0)
全栈程序员-站长的头像全栈程序员-站长


相关推荐

  • COleVariant的转换

    COleVariant的转换1.转为int型var.ChangeType(VT_I4);  int i=var.intVal; 2.COleVariant是数据库常用到的数据类型。它可以是字串,整型值,日期等。知道怎样将它转换为CString很有用处。设有CStringA;COleVariantB;来看看怎样将COleVariant转换为CString:switch(B.vt){

    2022年7月18日
    20
  • HTML中document的作用,html中的document对象是什么?一篇文章让你了解document对象

    HTML中document的作用,html中的document对象是什么?一篇文章让你了解document对象本篇文章主要的介绍了关于HTMLdocument对象的解释,还有关于HTMLdocument对象的使用实例解析,接下来让我们一起来看这篇文章吧首先我们来介绍一下HTML中的document对象:文档对象(document)代表浏览器窗口中的文档,该对象是window对象的子对象,由于window对象是DOM对象模型中的默认对象,因此window对象中的方法和子对象不需要使用window来引用。…

    2022年7月19日
    19
  • JAVA布局模式:GridBagConstraints终极技巧参数详解「建议收藏」

    JAVA布局模式:GridBagConstraints终极技巧参数详解「建议收藏」布局模式:GridBagConstraints布局,先发一个实例:gridx=2;//X2gridy=0;//Y0gridwidth=1;//横占一个单元格gridheight=1;//列占一个单元格weightx=0.0;//当窗口放大时,长度不变weighty=0.0;//当窗口放大时,高度不变anchor=Gr

    2025年10月14日
    5
  • decimal 整数 mysql_MySQL decimal类型

    decimal 整数 mysql_MySQL decimal类型在本教程中,我们将向您介绍MySQLDECIMAL数据类型以及如何在数据库表中有效地使用它。MySQLDECIMAL数据类型简介MySQLDECIMAL数据类型用于在数据库中存储精确的数值。我们经常将DECIMAL数据类型用于保留准确精确度的列,例如会计系统中的货币数据。要定义数据类型为DECIMAL的列,请使用以下语法:column_nameDECIMAL(P,D);在上面的语法中:P…

    2022年7月17日
    19
  • Java初学者常用的开发工具都有哪些?[通俗易懂]

    Java初学者常用的开发工具都有哪些?[通俗易懂]俗话说的好“工欲善其事必先利其器”,要想学好java这门语言,选择一款好用顺手的开发工具是必不可少的。java初学者常用源码编辑工具介绍java源代码本质上其实就是普通的文本文件,所以理论上来说任何可以编辑文本文件的编辑器都可以作为我们的java代码编辑工具。比如:WINDOWS记事本、写字板、word等。但是这些简单工具没有语法的高亮提示、自动完成等功能,这些功能的缺失会大大降低代码的编写效率。学习开发时一般我们不会选用这些简…

    2022年7月7日
    17
  • pageinfo 分页实现_分页显示

    pageinfo 分页实现_分页显示项目中有时候遇到list多种来源,不能使用PageInfo在查询数据库时分页,需要查询全部之后,手动分页。 //手动分页的分割起始下标IntegerfromIndex=0;//手动分页的分割结尾下标IntegertoIndex=0;Integertotal=arrayList.size();…

    2025年7月28日
    2

发表回复

您的邮箱地址不会被公开。 必填项已用 * 标注

关注全栈程序员社区公众号