好吧这个部分比较复杂所以本来没有打算多说的（事实上期中发表的时候这个地方也不是重点，所以时间有限决定跳过），8过既然conanshang问道，这里就简单说一说把w
其实对比相似度的算法有很多，比如PIV通常用的标准“直接相互相关法”的算法的话是这样：Rfg = A / (B * C)，其中：
A = sum[]sum[](f(x + i, y + j) – fbar) * (g(x + ζ + i, y + η + j) – gbar)
B = sqrt(sum[]sum[](f(x + i, y + j) – fbar)2)
C = sqrt(sum[]sum[](g(x + ζ + i, y + η + j) – gbar)2)
fbar = sum[]sum[](f(x + i, y + j)) / (n * m)
gbar = sum[]sum[](g(x + ζ + i, y + η + j)) / (n * m)
（sum[]sum[]是sum[j = 0 to n]sum[i = 0 to m]，相当于for (j = 0; j < n; j++){ for (i = 0; i < m; i++) }}）
Rfg是上一篇文章中的红色正方形区域和白色正方形区域的相关度，x, y分别是小区域的起点坐标，对应上一篇文章的红色正方形的左上角坐标；n, m分别是长宽，i, j分别是正方形内坐标；ζ, η分别是白色正方形相对红色正方形的相对坐标偏差，一点点移动ζ和η，得出每一个位置的R值，然后R值最大的位置的ζ和η就是上一篇文章红色正方形范围的流体的移动方向距离 Continue reading “PIV解析时用的算法”
