运行时间错误使用CV :: Mat :: AT分配CV :: MAT元素值

Question

在以下代码中，我想为循环中垫变量的元素分配一个值。我在下面遇到运行时错误。

pair<Mat, Mat> meshgrid(vector<int> x, vector<int> y) {

    int sx = (int)x.size();
    int sy = (int)y.size();
    
    Mat xmat = Mat::ones(sy, sx, CV_16U);
    Mat ymat = Mat::ones(sy, sy, CV_16U);

    for (int i = 0; i < sx; i++) {
        for (int j = 0; j < sy; j++) {
            xmat.at<int>(i, j) = j;  // <------- here is place of error.
            cout << j << "\t";
        }
        cout << endl;
    }

    for (int i = 0; i < sx; i++) {
        for (int j = 0; j < sy; j++) {
            ymat.at<int>(i, j) = i;  // <------- here is place of error.
            cout << j << "\t";
        }
        cout << endl;
    }

    return make_pair(xmat, ymat);
}

调试时的这张照片；

这是我收到的运行时间错误：

OpenCV(...) Error: Assertion failed
(((((sizeof(size_t)<<28)|0x8442211) >> ((traits::Depth<_Tp>::value) &
((1 << 3) - 1))*4) & 15) == elemSize1()) in cv::Mat::at, file
...\include\opencv2\core\mat.inl.hpp, line 1108

谢谢您的回答。

Answer 1

我假设您打算生成类似于 /code> 和 matlab matlab meshgrid meshgrid 。

您的代码中有几个错误：

1. cv :: mat使用类型CV_16U初始化（即16位无符号值），但是当您使用 at 访问元素时，您会使用int（为32位签名）。
   您应该将其更改为at＆lt; unsigned Short＆gt;（或将cv :: mat的类型更改为32位签名 - cv_32s）。
2. 您初始化了cv :: mat带错误的大小：xmat的大小为 （sy，sx） ，但是ymat的大小为 （SY，SY） 。
3. 您用来访问垫子元素的索引（行，col）不正确。为了使其更易于正确使用，我将装饰的名称更改为行，cols，
   循环索引IROW，ICOL。
4. 矩阵中的值应来自x和y向量（而不是索引）中的值。

请参阅下面的更新代码（以及有关更改的下面的注释）：

#include <opencv2/core/core.hpp>
#include <vector>
#include <iostream>

std::pair<cv::Mat, cv::Mat> meshgrid(std::vector<unsigned short> const & x, std::vector<unsigned short> const & y)
{
    int cols = static_cast<int>(x.size());
    int rows = static_cast<int>(y.size());
    cv::Mat xmat(rows, cols, CV_16U);
    cv::Mat ymat(rows, cols, CV_16U);

    for (int iRow = 0; iRow < rows; ++iRow) {
        auto * pRowData = xmat.ptr<unsigned short>(iRow);
        for (int iCol = 0; iCol < cols; ++iCol) {
            pRowData[iCol] = x[iCol];
            std::cout << pRowData[iCol] << "\t";
        }
        std::cout << std::endl;
    }

    std::cout << std::endl;

    for (int iRow = 0; iRow < rows; ++iRow) {
        auto * pRowData = ymat.ptr<unsigned short>(iRow);
        for (int iCol = 0; iCol < cols; ++iCol) {
            pRowData[iCol] = y[iRow];
            std::cout << pRowData[iCol] << "\t";
        }
        std::cout << std::endl;
    }

    return std::make_pair(std::move(xmat), std::move(ymat));
}

int main()
{
    std::vector<unsigned short> xxx{ 1,2 };
    std::vector<unsigned short> yyy{ 10,11,12 };
    auto p = meshgrid(xxx, yyy);
    return 0;
}

输出：

1       2
1       2
1       2

10      10
11      11
12      12

一些注释：

1. 我可能会误解您要在cv :: mat's中设置的值。但是至少现在您的代码不会崩溃。您可以根据需要更改分配的值。
2. 使用at访问cv :: mat一个接一个的元素效率非常低，因为at包含每个访问的一些验证。
   使用cv :: mat方法ptr，它为您提供了一排数据的指针。然后，您可以使用此指针来更有效地穿越该行 - 请参见上文
3. 以任何方法，遍历cv :: mat一排接头（而不是列列）更有效。这会导致您访问连续的内存，并减少缓存失误的数量。
4. 在您的实际代码中，最好将计算与I/O分开。因此，如果您的meshgrid函数只会创建矩阵，那就更好了。如果需要，将它们打印在外面。
5. 无需将cv :: mat s初始化，因为之后立即为所有元素设置值。
6. 在我的代码x和y中，通过const refernce传递给该函数。它更有效（避免复制），还迫使编译器验证向量未经修改。
7. 最好避免使用命名空间std - 请参阅此处。为什么“使用名称空间std”；被认为是不好的做法吗？。
   从类似的原因中，我也推荐使用命名空间CV 避免。