PIVlab软件使用(第二篇)

fly往事随风都随风

于 2024-09-24 10:53:39 发布

阅读量714

点赞数 20

文章标签：数码相机计算机视觉人工智能

本文链接：https://blog.csdn.net/qq_14865693/article/details/142479375

版权

我滴老天奶，快9个月没写这个系列了，因为被裁员了，悲了个崔的，想想还是把这个系列写完吧

前言

看了下源码，原作者大大今年又更新啦，这次界面整的更好看了些

最简单直接的就是logo好看了很多。。。

找到PIVLAB_GUI.m，运行，界面变成了这样

上一次导入数据图像讲过了，这个版本基本上都一样的。这次我们先看看怎么出矢量场吧

计算矢量场

通常情况下，第一步是需要标定的，即m/s级，但很多图像序列没有提供对应的标定图像，可以看pix/s级的情况。在工程项目的Examples文件夹下面，有一些示例图像，读取两张卡门涡的吧

左上角有个分析Analysis

PIV 参数设置

PIV settings，主要是设置相关参数，根据当前的图像和先验知识配置，然后再进行计算。通常情况下，可以认为是一个预览，先看看第一帧结果大概上长什么样子，然后再去计算几百几千长的数据

ANALYZE!就是直接运行，运行当前帧或者全部帧，没什么说的，有一点，当几百张几千张同时载入计算，有的时候可能卡壳了，分析所有帧的那个按钮下面有cancel，可以结束当前进程，省的电脑死机了。。。

PIV算法

主要包括PIV相关计算的算法

FFT变形窗算法

这个算法是该作者主推的算法，代码在piv_FFTmulti.m中体现。详细原理在作者论文中或者其他论文里都有介绍...多级网格迭代的图像变形算法，用粗网格的计算结果为图像变形提供第一次的估计位移，对细分网格采用图像变形提高位移计算精度。如果看英文原文比较吃力，也可以看看国内的一些相关算法的介绍，比如多级网格迭代的图像变形算法在PIV中的应用 (hanspub.org)

[xtable, ytable, utable, vtable, typevector, correlation_map,correlation_matrices] = piv_FFTmulti (image1,image2,interrogationarea, step, subpixfinder, mask_inpt, roi_inpt,passes,int2,int3,int4,imdeform,repeat,mask_auto,do_linear_correlation,do_correlation_matrices,repeat_last_pass,delta_diff_min)

输入参数

输入image1,image2就是待计算的相邻的两帧图像（单帧相机）或AB图像（跨帧相机）

interrogationarea，界面上设置的窗尺寸，也是相关计算时图像块的尺寸，注意PIVLAB中只用正方形块图像计算，没有长方形窗。

step,Pass1设置的步长，可以是0%，25%，50%，75%，也就是其他软件中需要设置的重叠率

因为底层算法逻辑不一样，所以会存在一些区别，主要体现在矢量网格的密度上

mask_inpt就是遮罩点，roi_inpt想要分析的区域的端点，passes表示界面设置的层数，只勾选了Pass1,passes=1，若勾选了Pass2,passes=2,勾选了Pass4，那么前面几个强制勾选，passes=4

PIVlab中2-4层的重叠率都是窗的50%,是受底层其重建算法有关，其他软件貌似不受这个限制，而且做矢量对比也不好做，因为PIVLAB算法会自动将网格划分到居中计算，跟其他软件也不一样。

其他几个参数就和鲁棒相关标准有关了，不怎么用。

输出参数

xtable, ytable, utable, vtable, typevector, correlation_map,correlation_matrices

xtable, ytable, utable, vtable表示了u,v速度矢量（像素级）在图像上划分的网格坐标位置x,y（像素级）上的情况，typevector表示当前的uv矢量标志，是遮罩特征或者其他，correlation_map表示了图像块相关计算收敛时所得的系数，注意只有相关计算才会有这个值