在OpenCV 2.3中为Mat数组赋值的问题 - 似乎很简单

Question

使用OpenCV 2.3的新API,我无法将值分配给循环内的Mat数组(或说图像).这是我正在使用的代码片段;

    int paddedHeight = 256 + 2*padSize; 
    int paddedWidth = 256 + 2*padSize;  

    int n = 266; // padded height or width

    cv::Mat fx = cv::Mat(paddedHeight,paddedWidth,CV_64FC1);
    cv::Mat fy = cv::Mat(paddedHeight,paddedWidth,CV_64FC1);        
    float value = -n/2.0f;

    for(int i=0;i<n;i++)
    {
        for(int j=0;j<n;j++)
            fx.at<cv::Vec2d>(i,j) = value++;                    

        value = -n/2.0f;
    }

    meshElement = -n/2.0f;

    for(int i=0;i<n;i++)
    {
        for(int j=0;j<n;j++)
            fy.at<cv::Vec2d>(i,j) = value;
        value++;
    }

现在在第一个循环中,一旦j = 133,我得到一个异常,这似乎与图像的深度有关,我无法弄清楚我在这里做错了什么.

请指教!谢谢!

Answer 1

您正在以双组分双向量(使用.at<cv::Vec2d>())访问数据,但您创建的矩阵仅包含1个组件双精度(使用CV_64FC1).要么创建矩阵以包含每个元素的两个组件(带CV_64FC2),或者更适合您的代码,请使用简单的双精度访问值.at<double>().这恰好在j = 133处爆炸,因为它是图像大小的一半,并且当它仅包含1时被视为包含2分量向量,它只有一半宽.

或者也许您可以将这两个矩阵合并为一个,每个元素包含两个组件,但这取决于您将来使用这些矩阵的方式.在这种情况下,您还可以将两个循环合并在一起,并真正设置一个双组件向量:

cv::Mat f = cv::Mat(paddedHeight,paddedWidth,CV_64FC2);
float yValue = -n/2.0f;

for(int i=0;i<n;i++)
{
    float xValue = -n/2.0f;

    for(int j=0;j<n;j++)
    {
        f.at<cv::Vec2d>(i,j)[0] = xValue++;
        f.at<cv::Vec2d>(i,j)[1] = yValue;
    }

    ++yValue;
}

如果您始终需要两个值(fx中的值和来自fy的值),则可能会产生更好的内存访问方案.