OpenCV cv::Mat类操作

分类： OpenCV 2011-12-29 11:47 2485人阅读 评论(17) 收藏 举报

首先还是要感谢箫鸣朋友在我《OpenCV学习笔记（四十）——再谈OpenCV数据结构Mat详解》的留言，告诉我M.at<float>(3, 3)在Debug模式下运行缓慢，推荐我使用M.ptr<float>(i)此类方法。这不禁勾起了我测试一下的冲动。下面就为大家奉上我的测试结果。

我这里测试了三种操作Mat数据的办法，套用流行词，普通青年，文艺青年，为啥第三种我不叫2b青年，大家慢慢往后看咯。

普通青年的操作的办法通常是M.at<float>(i, j)

文艺青年一般会走路线M.ptr<float>( i )[ j ]

暴力青年通常直接强制使用我第40讲提到的M.data这个指针

实验代码如下：

[cpp] view plain copy print ?

1. t = (double)getTickCount();
2. Mat img1(1000, 1000, CV_32F);
4. for (int i=0; i<1000; i++)
5. {
6. for (int j=0; j<1000; j++)
7. {
8. img1.at<float>(i,j) = 3.2f;
9. }
10. }
11. t = (double)getTickCount() - t;
12. printf("in %gms\n", t*1000/getTickFrequency());
13. /
14. t = (double)getTickCount();
15. Mat img5(1000, 1000, CV_32F);
16. float *pData1;
17. for (int i=0; i<1000; i++)
18. {
19. pData1=img5.ptr<float>(i);
20. for (int j=0; j<1000; j++)
21. {
22. pData1[j] = 3.2f;
23. }
24. }
25. t = (double)getTickCount() - t;
26. printf("in %gms\n", t*1000/getTickFrequency());
28. t = (double)getTickCount();
29. Mat img6(1000, 1000, CV_32F);
30. float *pData2;
31. Size size=img6.size();
32. if(img2.isContinuous())
33. {
34. size.width = size.width*size.height;
35. size.height = 1;
36. }
37. size.width*=img2.channels();
38. for(int i=0; i<size.height; i++)
39. {
40. pData2 = img6.ptr<float>(i);
41. for(int j=0; j<size.width; j++)
42. {
43. pData2[j] = saturate_cast<float>(3.2f);
44. }
45. }
46. t = (double)getTickCount() - t;
47. printf("in %gms\n", t*1000/getTickFrequency());

t = (double)getTickCount(); Mat img5(1000, 1000, CV_32F); float *pData1; for (int i=0; i<1000; i++) { pData1=img5.ptr<float>(i); for (int j=0; j<1000; j++) { pData1[j] = 3.2f; } } t = (double)getTickCount() - t; printf("in %gms\n", t*1000/getTickFrequency()); t = (double)getTickCount(); Mat img6(1000, 1000, CV_32F); float *pData2; Size size=img6.size(); if(img2.isContinuous()) { size.width = size.width*size.height; size.height = 1; } size.width*=img2.channels(); for(int i=0; i<size.height; i++) { pData2 = img6.ptr<float>(i); for(int j=0; j<size.width; j++) { pData2[j] = saturate_cast<float>(3.2f); } } t = (double)getTickCount() - t; printf("in %gms\n", t*1000/getTickFrequency());
测试结果：

	Debug	Release
加强版文艺青年	5.74ms	2.43ms
终极版文艺青年	40.12ms	2.34ms

我的测试结果感觉这两种方案只是锦上添花的效果，也使大家的操作有了更多的选择，但感觉在速度上并没有数量级的提升，再次感谢箫铭对我blog的支持。后来箫铭说saturate_cast才把速度降下来，我很同意，就不贴上去测试结果了。但我查看资料了解了一下 saturate_cast的作用。 可以看成是类型的强制转换，比如对于saturate_cast<uchar>来说，就是把数据转换成8bit的0~255区间，负值变成0，大于255的变成255。如果是浮点型的数据，变成round最近的整数，还是很有用处的函数，推荐大家在需要的时候尝试。