提升multi_array视图

Question

我有一个3维boost::multi_array表示2d瓦片地图的图层.我希望能够清除一个图层 - 也就是遍历一层上的所有图块并将其值设置为0,但我无法弄清楚如何做到这一点 - 我相信我必须使用视图,但文档并没有对我来说很有意义.例如,在查看此示例代码时:

// myarray = 2 x 3 x 4

//
// array_view dims: [base,bound) (dimension striding default = 1)
// dim 0: [0,2) 
// dim 1: [1,3) 
// dim 2: [0,4) (strided by 2), 
// 

typedef boost::multi_array_types::index_range range;
// OR typedef array_type::index_range range;
array_type::array_view<3>::type myview = myarray[ boost::indices[range(0,2)][range(1,3)][range(0,4,2)] ];

for (array_type::index i = 0; i != 2; ++i)
    for (array_type::index j = 0; j != 2; ++j)
        for (array_type::index k = 0; k != 2; ++k) 
            assert(myview[i][j][k] == myarray[i][j+1][k*2]);

boost::indices工作怎么样？这些范围甚至代表什么？对不起,如果这对某些人来说是显而易见的,我就是无法掌握它.

在此先感谢,ell.

Answer 1

该示例创建了一个3D视图,可以访问较大的3D数组的一部分,并且它有点复杂,因为它们在有界的同时显示跨步视图.

如果他们使用它会更容易理解[range(0,2)][range(1,3)][range(0,4)],这会切出一个坚固的3D块(沿第一维0..2,沿第二维1..3,沿第三维0..4)

你的目标是2D视图(切片),所以让我们看一下,我会简单地说:

你有一个3x3x3阵列,你想要将一个穿过它的中心的平面归零.有三个这样的平面,所以让我们为三个平面创建一个2D视图:

#include <iostream>
#include <iomanip>
#include <boost/multi_array.hpp>

int main()
{
    typedef boost::multi_array<int, 3> array_t;

    array_t a(boost::extents[3][3][3]); // 3D array of 3x3x3

    int values = 0;
    for(array_t::index i = 0; i != 3; ++i)
        for(array_t::index j = 0; j != 3; ++j)
            for(array_t::index k = 0; k != 3; ++k)
                a[i][j][k] = values++;

    for(array_t::index i = 0; i != 3; ++i) {
        for(array_t::index j = 0; j != 3; ++j) {
            for(array_t::index k = 0; k != 3; ++k)
                std::cout << std::setw(2) << a[i][j][k] << ' ';
            std::cout << '\n';
        }
        std::cout << '\n';
    }

到目前为止这么好..现在让我们剪掉一些2D飞机:

typedef boost::multi_array_types::index_range range_t;
array_t::index_gen indices;

// 2D plane that's constant in the 1st dimension
array_t::array_view<2>::type plane =
    a[ indices[1][range_t(0,3)][range_t(0,3)] ];

for (array_t::index i = 0; i != 3; ++i)
    for (array_t::index j = 0; j != 3; ++j)
        plane[i][j] = 0;

现在尝试打印这个立方体,你会看到三个平面中的一个被归零.

然后更换初始化为plane用

a[ indices[range_t(0,3)][1][range_t(0,3)] ];

与

a[ indices[range_t(0,3)][range_t(0,3)][1] ];

现在打印立方体..你会明白的!