在c ++中将值分配给多维特征矩阵?

时间:2017-03-21 15:05:32

标签: c++ matrix eigen eigen3

我有一个特征矩阵,矢量作为标量类型,定义如下

Eigen::Matrix<Eigen::VectorXd, Eigen::Dynamic, Eigen::Dynamic> res = Eigen::Matrix<Eigen::VectorXd, Eigen::Dynamic, Eigen::Dynamic>::
        Constant(10,
            10,
            Eigen::VectorXd::Zero(3));

我可以访问并指定值,如下所示

res(0,0)(0) = 1;

我正在尝试将矩阵分配给每个第三个索引,如下所示

res(all_rows,all_cols)(0) = matrix

以下是我的代码

int main()
{
    Eigen::Matrix<Eigen::VectorXd, Eigen::Dynamic, Eigen::Dynamic> res = Eigen::Matrix<Eigen::VectorXd, Eigen::Dynamic, Eigen::Dynamic>::
        Constant(10,
            10,
            Eigen::VectorXd::Zero(3));

    Eigen::MatrixXf x = Eigen::MatrixXf::Ones(10, 1)*2;
    Eigen::MatrixXf y = Eigen::MatrixXf::Ones(10, 1);
    Eigen::MatrixXf product = Eigen::MatrixXf::Zero(10, 10);
    MatrixXd m(2, 2);
    res(0, 0)(0) = 1;
    m(0, 0) = 3;
    m(1, 0) = 2.5;

    m(0, 1) = -1;
    m(1, 1) = m(1, 0) + m(0, 1);
    product = x*y.transpose();
    res.block(0,0,10,10)(0) = x*y.transpose();
}

但是当我尝试编译代码时,这一行给出了错误

res.block(0,0,10,10)(0) = x*y.transpose();

我收到以下错误。

Error   C2338   THIS_COEFFICIENT_ACCESSOR_TAKING_ONE_ACCESS_IS_ONLY_FOR_EXPRESSIONS_ALLOWING_LINEAR_ACCESS

有没有办法克服这个问题,还是应该改变我的做法?

1 个答案:

答案 0 :(得分:3)

您正在尝试使用Eigen::VectorXd作为标量类型。因此,您应该确保按https://eigen.tuxfamily.org/dox-devel/TopicCustomizing_CustomScalar.html

中的说明扩展Eigen::NumTraits

看起来您正在尝试重新创建多通道阵列或图像。所以你通过上述方法获得了很多收益。更好的想法是使用等级3的Eigen::Tensor或尺寸(3,宽*高)的扁平矩阵。然后你得到所有的Eigen好处(矢量化标量操作,没有奇怪的访问模式,如(row,col)(channel)等)

Eigen::Tensor<double, 3> rgb_image(height, width, 3);
Eigen::MatrixXd rgb_image(3, width*height); // Each column is a pixel color