Mat数据操作之普通青年、文艺青年、暴力青年

转自http://blog.csdn.net/yang_xian521/article/details/7161335

首先还是要感谢箫鸣朋友在我《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这个指针

实验代码如下：

  

   

    

     

    
    

     

      t = (double)getTickCount();
     
    
   

    

     

    
    

     

      Mat img1(1000, 1000, CV_32F);
     
    
   

    

     

    
    

     

      for (int i=0; i<1000; i++)
     
    
   

    

     

    
    

     

      { 
     
    
   

    

     

    
    

     
	for (int j=0; j<1000; j++)
     
    
   

    

     

    
    

     

      	{ 
     
    
   

    

     

    
    

     

      		img1.at<float>(i,j) = 3.2f;
     
    
   

    

     

    
    

     

      	}
     
    
   

    

     

    
    

     

      }
     
    
   

    

     

    
    

     

      t = (double)getTickCount() - t;
     
    
   

    

     

    
    

     

      printf("in %gms\n", t*1000/getTickFrequency());
     
    
   

    

     

    
    

     

      //***************************************************************
     
    
   

    

     

    
    

     

      t = (double)getTickCount();
     
    
   

    

     

    
    

     

      Mat img2(1000, 1000, CV_32F);
     
    
   

    

     

    
    

     

      for (int i=0; i<1000; i++) 
     
    
   

    

     

    
    

     

      { 
     
    
   

    

     

    
    

     
	for (int j=0; j<1000; j++)
     
    
   

    

     

    
    

     

      	{ 
     
    
   

    

     

    
    

     

      		img2.ptr<float>(i)[j] = 3.2f; 
     
    
   

    

     

    
    

     

      	}
     
    
   

    

     

    
    

     

      } 
     
    
   

    

     

    
    

     

      t = (double)getTickCount() - t;
     
    
   

    

     

    
    

     

      printf("in %gms\n", t*1000/getTickFrequency());
     
    
   

    

     

    
    

     

      //***************************************************************
     
    
   

    

     

    
    

     

      t = (double)getTickCount();
     
    
   

    

     

    
    

     

      Mat img3(1000, 1000, CV_32F);
     
    
   

    

     

    
    

     

      float* pData = (float*)img3.data;
     
    
   

    

     

    
    

     

      for (int i=0; i<1000; i++)
     
    
   

    

     

    
    

     

      { 
     
    
   

    

     

    
    

     
	for (int j=0; j<1000; j++)
     
    
   

    

     

    
    

     

      	{
     
    
   

    

     

    
    

     

      		*(pData) = 3.2f; pData++; 
     
    
   

    

     

    
    

     

      	}
     
    
   

    

     

    
    

     

      } 
     
    
   

    

     

    
    

     

      t = (double)getTickCount() - t;
     
    
   

    

     

    
    

     

      printf("in %gms\n", t*1000/getTickFrequency()); 
     
    
   

    

     

    
    

     

      //***************************************************************
     
    
   

    

     

    
    

     

      t = (double)getTickCount();
     
    
   

    

     

    
    

     

      Mat img4(1000, 1000, CV_32F);
     
    
   

    

     

    
    

     

      for (int i=0; i<1000; i++)
     
    
   

    

     

    
    

     

      { 
     
    
   

    

     

    
    

     
	for (int j=0; j<1000; j++) 
     
    
   

    

     

    
    

     

      	{ 
     
    
   

    

     

    
    

     

      		((float*)img3.data)[i*1000+j] = 3.2f; 
     
    
   

    

     

    
    

     

      	}
     
    
   

    

     

    
    

     

      }
     
    
   

    

     

    
    

     

      t = (double)getTickCount() - t;
     
    
   

    

     

    
    

     

      printf("in %gms\n", t*1000/getTickFrequency());
     
    
  
 

最后两招可以都看成是暴力青年的方法，因为反正都是指针的操作，局限了各暴力青年手段就不显得暴力了。

在Debug、Release模式下的测试结果分别为：

测试结果

根据测试结果，我觉得箫铭说的是很可信的，普通青年的操作在Debug模式下果然缓慢，他推荐的文艺青年的路线确实有提高。值得注意的是本来后两种办法确实是一种比较2b青年的做法，因为at操作符或者ptr操作符，其实都是有内存检查的，防止操作越界的，而直接使用data这个指针确实很危险。不过从速度上确实让人眼前一亮，所以我不敢称这样的青年为2b，尊称为暴力青年吧。

不过在Release版本下，几种办法的速度差别就不明显啦，都是很普通的青年。所以如果大家最后发行程序的时候，可以不在意这几种操作办法的，推荐前两种哦，都是很好的写法，操作指针的事还是留给大神们用吧。就到这里吧~~

ps:在OpenCV. 2. Computer Vision Application Programming Cookbook一书中还提到了一种比at方法简便的方法，如下

使用at的方法的话，我们在访问矩阵的每个元素时都要制定返回值的类型，这样显得过于麻烦，而如果通过形如

cv::Mat_<uchar> img = imread("lena.jpg");
 

img(50,100) = 0;
 

文章来源: panda1234lee.blog.csdn.net，作者：panda1234lee，版权归原作者所有，如需转载，请联系作者。

原文链接：panda1234lee.blog.csdn.net/article/details/8926186