基于VS与OpenCV的模板匹配学习（3）：OpenCV Mat快速遍历

基于VS与OpenCV的模板匹配学习（3）：

OpenCV Mat 快速遍历

---

---

前言

在前文的模板匹配中，基于边缘的模板匹配的速度并不友好，尝试对算法进行时间的优化。其中，对OpenCV的Mat遍历操作十分普遍，本文结合OpenCV的Mat的存储结构，试图找到一种的高效遍历Mat的方法，为后续进行时间优化作铺垫。
网上有类似文章可供参考：
opencv Mat遍历速度实验 https://blog.csdn.net/fanghtao/article/details/80995997
opencv高效遍历图像 https://blog.csdn.net/jacke121/article/details/94432519

---

一、OpenCV图像类型

OpenCV的图像类型的格式：CV_[bite]（U|S|F）C [channels]
bite : 比特数，位数，对应在Mat中每个像素的所占的空间大小。
U|S|F ：
U : unsigned int , 无符号整型
S : signed int , 有符号整型
F : float , 单精度浮点型
channels：图像的通道数

本场景中，使用的图像类型为CV_8UC1，即单通道的uchar格式，与CV_32FC1的单通道float格式在运算速度上差异较大。

cv::Mat resultImage;
resultImage.create(5000, 5000, CV_8UC1);

---

二、OpenCV指针ptr

OpenCV提供了方便访问Mat元素的指针ptr，正常的访问格式如下：

resultImage.create(5000, 5000, CV_8UC1);
int resultImageRow = 0;
int resultImageCol = 0;
//先访问行指针，再访问相应元素
uchar * resultImagePtr = resultImage.ptr<uchar>(resultImageRow);
resultImagePtr[resultImageCol] = 1;
//直接访问相应元素
resultImage.ptr<uchar>(resultImageRow)[resultImageCol]= 1;

通过上述代码，我们注意到：
resultImage.ptr(resultImageRow)返回的是一个行指针头；
resultImage.ptr(resultImageRow)[resultImageCol]返回的是相应的元素；
这点在后续时间的讨论中，两者的差异性会愈发明显。

注意事项：uchar指针对应CV_8UC1，int指针对应CV_32SC1，float指针对应CV_32FC1。

---

三、遍历时间对比

在场景设计中，我们使用VS的release模式，与debug模式时间差异性较大。

（1）直接遍历目标元素：

	for (int resultImageRow = 0; resultImageRow < 5000; resultImageRow++)
		for (int resultImageCol = 0; resultImageCol < 5000; resultImageCol++)
			resultImage.ptr<uchar>(resultImageRow)[resultImageCol]= 1;

代码精简，利用ptr直接遍历目标元素，效率最低。
Run Time:0.035s

（2）行指针逐行遍历目标元素：

	uchar* resultImagePtr;
	for (int resultImageRow = 0; resultImageRow < 5000; resultImageRow++)
	{
		resultImagePtr = resultImage.ptr<uchar>(resultImageRow);
		for (int resultImageCol = 0; resultImageCol < 5000; resultImageCol++)
			resultImagePtr[resultImageCol] = 1;
	}

先访问行指针ptr，再依次遍历目标元素，效率中等。
Run Time:0.004s

（3）指针按行列遍历目标元素：
除ptr指针之外，OpenCV还提供Image.data指向数据的首部。在大部分场景中，OpenCV的数据是连续存储的，可以通过isContinuous()判断矩阵元素在每行末尾是否连续存储而没有间隙。经过一些场景测试，目前未发现网上提到的行结尾不满4Byte会自动补零的情况。

	uchar* resultImageDataStart = resultImage.data;
	for (int resultImageRow = 0; resultImageRow < 5000; resultImageRow++)
		for (int resultImageCol = 0; resultImageCol < 5000; resultImageCol++)
			*(resultImageDataStart++) = 1;

利用元素存储的连续性，按行列滑动指针位置遍历目标元素，效率较高。
Run Time:0.003s

注意事项：resultImage.data默认uchar格式，倘若图像格式为CV_32FC1，需要强制格式转化。

float* resultImageDataStart = (float*) resultImage.data;

（4）指针单循环遍历目标元素：

	uchar* resultImageDataStart = resultImage.data;
	for (int resultImageSeq = 0; resultImageSeq < 25000000; resultImageSeq++)
		*(resultImageDataStart++) = 1;

指针连续遍历目标元素，代码精简，效率最高。
Run Time:0.003s

---

总结

综上所述：
遍历速度：ptr(row)[col] < ptr(row) < data