乐与大家分享。合成的PIV软件性能测试图像文件

PIV

一般PIV的互相关计算直接计算的结果都是整数位的（直接计算或者通过傅利叶变换的结果），然后可以通过亚像素拟合，就可以得到小数位了（目前比较成熟的算法就有重心法，抛物线法和高斯拟合算法）。你要是做PIV算法分析，这些算法还是很有用的。

robin

多谢PIV 老兄。 说的对极了。我没有任何补充。哈哈：）

清风明月82

我也明白了,但是不知道 robin 这幅图像亚像素拟合是采用什么算法的那.
我了解了一下,好像高斯拟合算法是最好的.
亚像素拟合是不是所有的算法都要计算四个方向?
robin 有没有亚像素拟合方面的资料,可以上传一下,看一看,谢谢啦!

robin

清风还是有一点没有搞清楚。这个txt文本给出的已知速度分布是用数学公式人为生成的，供人们用各种PIV软件算法计算出流场速度后对比计算精度，误差等用的。txt文件给出的速度值是以像素为单位的，但它本身并不是通过任何像素拟合计算出来的结果。但是当你用和这个txt文件相对应的两个粒子图像计算速度场时却需要采用合适的亚像素拟合算法，才能得到精度达到亚像素的速度结果。

清风明月82

我用FFT编了一个PIV，用这个标准PIV图像进行了测试（结果未进行亚象素拟和），发现基本上能检测到流场的物理形态，但是在图像的右下脚检测到许多错误矢量，不知到是为什么？（不论查询窗口取32*32或是64*64，偏移都为16）
还有就是，请问robin你给出的已知速度分布文件和我算出的速度对不上号啊，虽然我没有进行亚象素拟和，但是整数位都对不上号，可能是定义的x，y方向不同吧！，我定义的x，y方向零点在图像的最左下角，x向右增大，y向上增大，不知道这个文件的x，y方向如何定义，请回复，谢谢！！！

清风明月82

FFT——PIV：查询窗口取32*32，偏移为16
位移矢量图如示：

robin

清风
在这个帖子的第二段我就曾经作了如下说明：
“ 对刚刚上传的文件稍微说明一下。
这是一对用数学方法合成的粒子图像。所以速度的分布是已知的。这对图像特别设计了可以用来显示PIV软件性能的结构。由于其数学合成的特性，图像是没有噪声的。这样就可以避免硬件性能差异对软件性能比较的干扰。
这对图像的运算难点有二：
1:涡旋结构尺度从左至右逐渐变小。
2:涡旋的流速从上至下逐渐变大
于是就有了一种直观的判据，一种软件所能给出的质量理想的速度场分布的区域越靠近右下角，则这种软件算法性能越好。因为这种软件能够分辨更细微的速度和空间漩涡结构。用这样一对合成图像来评测 PIV软件性能可以避免被用速度梯度很大的流场测试的结果所“蒙蔽”。从而能够客观而公正的评价软件的性能。”
   上述说明，换一种说法就是，越靠近右下角，运算的难度就越大。只用普通简单的FFT变换来计算，得不到满意的结果也是意料之中的。应该采用更为先进的计算方法才可以。即使是目前比较成熟的商用PIV软件，其先前的早期版本计算这个测试流场的结果和最新版本的计算结果的差别也是很大的。我会发一个PowerPoint文件给您。供您参考。

robin

关于这个图像的坐标。原点在左上角。见附图。您把坐标在变换一下，应该能够对应上了。

清风明月82

谢谢robin， 你的回复速度可真够快的。我明白了你的意思了。希望早日看到您的PowerPoint，如果不好上传的话，可以发到我的邮箱：lxiaohui82@163.com
还有我看前面大家用商用PIV软件算的结果都还不错，不知道这些软件主要采用什么算法啊？
谢谢

foxes

楼主能否将您的文件再发一下呢，我这里下载说文件被损坏了。或者发到我的邮箱：cxzhao@aoxue.org

robin

foxes
没问题。这里重新上传一下。同时已将文件发给您，请您查收。