C ++ OpenCV：迭代Mat中的像素，这是另一个Mat的ROI

Question

我有一个非常大的Mat，它实际上是另一个Mat的ROI（由otherMat获得（cv :: Rect（x，y，w，h）））。我想通过Mat的所有像素，做一些像素计算，并使用指针将结果写入另一个Mat。

通过所有像素，包括ROI之外的像素到目前为止工作正常，但我想知道在ROI之外跳过像素的最快方法是什么。我希望尽可能少地缓存未命中，并且我也不想进行低效的分支预测。什么是最好的方法呢？

编辑：我对获得感兴趣的特定区域的子矩阵不感兴趣。我感兴趣的是以最有效的方式通过指针迭代像素，而无需访问子矩阵之外的数据。区域。

Answer 1

只需使用子矩阵：

cv::Mat largeMat
cv::Rect roi(yourROI);
cv::Mat submatrix = largeMat(roi);

// now iterate over all the pixels of submatrix

每行末尾都会有缓存未命中

这是实际的代码示例，它显示了子像素之外的像素被跳过（你会在每一行的末尾获得额外的缓存未命中但是应该是全部）。

int main(int argc, char* argv[])
{
    cv::Mat input = cv::imread("C:/StackOverflow/Input/Lenna.png");

    cv::Rect roi(128, 128, 256, 256);
    cv::Mat submat = input(roi);

    cv::MatIterator_<cv::Vec3b> it; // = src_it.begin<cv::Vec3b>();
    for (it = submat.begin<cv::Vec3b>(); it != submat.end<cv::Vec3b>(); ++it)
    {
        (*it)[0] = 0;
        (*it)[1] = 0;
    }

    cv::imshow("input", input);
    cv::imwrite("C:/StackOverflow/Output/submatIter.png", input);
    cv::waitKey(0);
    return 0;
}

给出这个结果：

如果你想要它快一点，你可以使用行指针：http://docs.opencv.org/2.4/doc/tutorials/core/how_to_scan_images/how_to_scan_images.html

请注意，在链接中他们比较了调试模式运行时速度，这就是随机访问速度太慢的原因。在释放模式下，它应该比迭代器verson快（或者更快）。 但是这里是行-Ptr版本（它可以计算每个像素访问的行偏移量），它可以提供相同的结果并且应该是最快的方法（如果openCV的LUT函数不能用于你的任务）：

int main(int argc, char* argv[])
{
    cv::Mat input = cv::imread("C:/StackOverflow/Input/Lenna.png");

    cv::Rect roi(128, 128, 256, 256);
    cv::Mat submat = input(roi);

    cv::Vec3b * currentRow;

    for (int j = 0; j < submat.rows; ++j)
    {
        currentRow = submat.ptr<cv::Vec3b>(j);
        for (int i = 0; i < submat.cols; ++i)
        {
            currentRow[i][0] = 0;
            currentRow[i][1] = 0;
        }
    }

    cv::imshow("input", input);
    cv::imwrite("C:/StackOverflow/Output/submatIter.png", input);
    cv::waitKey(0);
    return 0;
}

Answer 2

由于OtherMat是原始垫子的一个子集，您希望对原始垫子进行操作，但只能在otherMat区域内进行操作

//As otherMat(cv::Rect(x,y,w,h)));
for(int j=x;j<x+w;j++) 
{
  for (int i=y;i<y+h;i++)
  {  
       original.at<uchar>(j,i) = 255; 
  }
}