同样遍历 Mat,为什么你的代码慢 10 倍?
文章目录前言一、什么是不连续Mat1.产生不连续内存的常见场景2.连续与不连续内存本质区别二、常见错误遍历方式踩坑分析1.错误一:at()逐像素访问速度慢2.错误二强行使用一维 data 指针高危崩溃3.错误三使用clone() 转连续内存高开销三、不连续 Mat 高性能优化方案无拷贝、零额外开销1.行指针遍历通用稳定2.parallel_for_ 多核并行遍历四、不连续 Mat遍历方式速度对比总结前言做 OpenCV 视觉开发几乎所有人都听过一句话不要用 at() 遍历大图要用连续指针提速。但现实工程里绝大多数慢代码、莫名卡顿、偶发崩溃并不是出在普通连续图像上而是出在不连续 Mat。当我们使用 ROI 抠图、submat 截取、图像局部裁剪、矩阵分片处理后Mat 的内存结构会发生改变行与行之间存在内存空白填充内存不再线性连续。很多开发者习惯性套用连续矩阵的一维指针写法轻则遍历耗时暴涨几十倍重则直接内存越界、程序闪退。更棘手的是大型 Mat 场景4K/8K高清影像、红外遥感大图、工业高精度图像一旦内存不连续错误的遍历方式会让 CPU 缓存命中率暴跌即使开启 Release 优化代码依旧卡顿严重。很多人为了省事直接使用 clone() 强制转为连续内存看似解决问题却引入了超大内存拷贝开销在高频处理、嵌入式、实时视觉项目中完全不可用。那么问题来了在不拷贝、不重构内存、保留不连续结构的前提下如何最大限度提高大型 Mat 的访问速度本文避开网上泛滥的基础连续 Mat 教程从内存原理、避坑禁忌、基础优化、多核并行优化逐层讲解给出生产环境最优遍历方案帮你彻底解决不连续矩阵遍历慢的难题。一、什么是不连续Mat1.产生不连续内存的常见场景OpenCV 中 Mat 默认创建为连续内存但是以下操作会生成非连续 Mat也是工程开发最高频场景ROI 感兴趣区域截取mat(Rect(x,y,w,h))submat 子矩阵提取mat.rowRange() / colRange()图像不规则裁剪、局部切片图像金字塔、边缘填充后的矩阵2.连续与不连续内存本质区别连续内存整张图像像素在物理内存中一字排开、无间隔data 指针可以直接线性遍历全部像素。不连续内存每行像素末尾存在内存对齐填充字节padding下一行起始地址不紧跟当前行末尾行与行之间存在内存断层。我们可以通过两个关键参数理解step单行真实内存字节长度包含填充字节cols * channels单行有效像素字节长度判定公式step cols * channels 即为连续反之不连续。std::string pathR(d:\detect.jpg);cv::Mat srccv::imread(path,0);//16K图像大小为16384*16000intchannelsrc.channels();inttotalPixelssrc.total();// 总像素数intstepsrc.step;intcontinusrc.cols*src.channels();if(stepcontinu)std::coutsrc is continuousstd::endl;elsestd::coutsrc is not continuousstd::endl;cv::Mat roisrc(cv::Rect(1344,1152,8000,14000));steproi.step;continuroi.cols*roi.channels();if(stepcontinu)std::coutroi is continuousstd::endl;elsestd::coutroi is not continuousstd::endl;上述代码在原图中截取一个ROI截取后Mat的step不等于continu此时会输出roi是不连续的。二、常见错误遍历方式踩坑分析1.错误一:at()逐像素访问速度慢cv::Matdst(roi.size(),CV_8UC1);autostartstd::chrono::high_resolution_clock::now();for(inti0;iroi.rows;i){for(intj0;jroi.cols;j)dst.atuchar(i,j)255-roi.atuchar(i,j);}autoendstd::chrono::high_resolution_clock::now();autodurationstd::chrono::duration_caststd::chrono::microseconds(end-start);底层每次都要进行行列地址计算边界检查函数调用开销极大在大图不连续场景下CPU 缓存完全失效速度很慢。2.错误二强行使用一维 data 指针高危崩溃// 高危代码不连续Mat直接崩溃uchar*proi.data;for(inti0;itotal;i){}上述代码对于不连续内存将跳过行填充字节读取非法内存大概率会出现闪退、段错误等严重错误因此需谨慎使用。3.错误三使用clone() 转连续内存高开销//利用clone拷贝一份至连续内存if(!roi.isContinuous())roiroi.clone();clone 会完整拷贝一份内存4K、8K 、16K大图拷贝耗时极高实时项目、嵌入式设备完全无法承受。三、不连续 Mat 高性能优化方案无拷贝、零额外开销1.行指针遍历通用稳定利用 ptr() 获取每行起始指针行内线性连续访问避开行尾填充字节兼容所有不连续矩阵编译器可自动向量化优化。cv::Matdst(roi.size(),CV_8UC1);autostartstd::chrono::high_resolution_clock::now();for(inti0;iroi.rows;i){uchar*pSrcroi.ptruchar(i);uchar*pDstdst.ptruchar(i);for(intj0;jroi.cols;j)pDst[j]255-pSrc[j];}autoendstd::chrono::high_resolution_clock::now();autodurationstd::chrono::duration_caststd::chrono::microseconds(end-start);2.parallel_for_ 多核并行遍历针对超大尺寸不连续 Mat4K/8K利用 OpenCV 原生并行接口自动拆分行任务、绑定 CPU 多核无需手动创建线程性能翻倍提升。classParallelTraverse:publicParallelLoopBody{private:Mat*src;Mat*dst;public:ParallelTraverse(Mat_src,Mat_dst):src(_src),dst(_dst){}voidoperator()(constRangerange)constoverride{intcolssrc-cols;intchannelssrc-channels();for(intirange.start;irange.end;i){constuchar*pSrcsrc-ptruchar(i);uchar*pDstdst-ptruchar(i);for(intj0;jcols*channels;j){pDst[j]255-pSrc[j];}}}};cv::Matdst(roi.size(),CV_8UC1);autostartstd::chrono::high_resolution_clock::now();parallel_for_(Range(0,roi.rows),ParallelTraverse(roi,dst));autoendstd::chrono::high_resolution_clock::now();autodurationstd::chrono::duration_caststd::chrono::microseconds(end-start);四、不连续 Mat遍历方式速度对比在同一台电脑上读取一张大小为16384x16000的图像并截取大小为8000x14000的ROI对ROI进行遍历反转测试at()ptr()和parallel_for_ 三种方法的速度对比结果如下遍历方式耗时(ms)安全性at() 逐像素遍历350~380安全ptr() 逐行遍历130~150安全parallel_for_遍历20~28安全运行结果如下图上面耗时是在优化已禁用(/Od)情况下测试的结果若启用最大优化速度(/O2)情况下耗时会更短。总结很多开发者把优化重心放在连续矩阵却忽略了 ROI、切片带来的不连续内存这也是工业视觉项目中隐性卡顿、偶发闪退的核心原因。OpenCV 的性能优化从来不是单一语法技巧而是对内存布局、CPU缓存、编译指令、多核调度的综合把控。本文摒弃冗余理论全部代码均可直接复制编译、投入生产。如果你正在做实时检测、工业分拣、遥感图像处理、嵌入式视觉开发这套不连续 Mat 优化方案可以帮你低成本拉高程序帧率、减少内存占用、规避莫名崩溃。后续我会持续更新 OpenCV 工业级优化干货如果本文对你有帮助欢迎点赞收藏关注有任何问题可以在评论区留言交流。
本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/2608201.html
如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!