时间戳

首页 > 解决方案 > 如何通过取 n 行/列然后跳过 m 行/列来对 Eigen 中的矩阵进行切片

c++ - 如何通过取 n 行/列然后跳过 m 行/列来对 Eigen 中的矩阵进行切片

问题描述

我想问一个关于在 Eigen 中以有效方式切片的问题。我想以一种允许我取n行然后跳过m行并重复该过程直到结束索引的方式对矩阵进行切片。例如，如果 A 是1000 x 20矩阵，我想形成矩阵 B，即800 x 20，其中包含 A 的每个连续 5 行的前 4 行。

我希望我的问题很清楚。

非常感谢你，萨达姆

标签： c++eigen

解决方案

取自Eigen 文档

假设您的计算机上有库，您可以编译以下代码。测试如何将一列超过三列。

#include <eigen3/Eigen/Dense>
#include <iostream>
 
using namespace std;
using namespace Eigen;
int main()
{

MatrixXf M1 = MatrixXf::Random(3,8);
cout << "Column major input:" << endl << M1 << "\n";
cout << M1.outerStride() << endl;
cout << (M1.cols()+2)/3 << endl;
Map<MatrixXf,0,OuterStride<> > M2(M1.data(), M1.rows(), (M1.cols()+2)/3, OuterStride<>(M1.outerStride()*3));
cout << "1 column over 3:" << endl << M2 << "\n";
 
typedef Matrix<float,Dynamic,Dynamic,RowMajor> RowMajorMatrixXf;
RowMajorMatrixXf M3(M1);
cout << "Row major input:" << endl << M3 << "\n";
Map<RowMajorMatrixXf,0,Stride<Dynamic,3> > M4(M3.data(), M3.rows(), (M3.cols()+2)/3,
                                              Stride<Dynamic,3>(M3.outerStride(),3));
cout << "1 column over 3:" << endl << M4 << "\n";
}

推荐阅读