急啊！！有做亚像素拟合的吗？--高斯拟合

cnnoah

用MATLAB计算速度，后面有必要进行亚像素拟合，主要是利用高斯曲线拟合。
不知道我下面列的拟合公式正确与否？？？
用这些公式拟合后，日本标准视频协会上的PIV标准图进行验证时，有很大的差别。
======================================================
横向拟合：
  i1 = i + [R(i-1,j)-R(i+1,j)]/[2*(R(i-1,j)+R(i+1,j)-2*R(i,j)]
  斜向拟合1：
i2 = i1 + [R(i-1,j-1)-R(i+1,j+1)]/[2*sqrt(2)*(R(i-1,j-1)+R(i+1,j+1)-2*sqrt(2)*

R(i,j)]
  斜向拟合2：
i3 = i2 + [R(i-1,j+1)-R(i+1,j-1)]/[2*sqrt(2)*(R(i-1,j+1)+R(i+1,j-1)-2*sqrt(2)*

R(i,j)]
  纵向拟合：
  j1 = j + [R(i,j-1)-R(i,j+1)]/[2*(R(i,j-1)+R(i,j+1)-2*R(i,j)]
  斜向拟合1：
j2 = j1 + [R(i-1,j-1)-R(i+1,j+1)]/[2*sqrt(2)*(R(i-1,j-1)+R(i+1,j+1)-2*sqrt(2)*

R(i,j)]
  斜向拟合2：
j3 = j2 + [R(i-1,j+1)-R(i+1,j-1)]/[2*sqrt(2)*(R(i-1,j+1)+R(i+1,j-1)-2*sqrt(2)*

R(i,j)]

piv

不知道你的公式从哪里查的，有错误（里面还多了个左括号，这样的程序也能编译通过 ！？）：
[R(i-1,j)-R(i+1,j)]/[2*(R(i-1,j)+R(i+1,j)-2*R(i,j)]
应该是：
[R(i-1,j)-R(i+1,j)]/[R(i-1,j)+R(i+1,j)-2*R(i,j)]
前面的加减号记不清了，你自己试一下吧。

cnnoah

谢啦，发贴子的时候粗心写错了。
公式不是从程序中粘贴过来的，我编的程序运行通过了。