如何将泰勒级数展开应用于opencv中的图像？

Question

对于计算机视觉项目，我正在尝试实现D.Lowe paper(SIFT)。我坚持将这个（检查 ppt 图像）翻译成代码。你能告诉我如何将泰勒级数展开应用到代码中吗（公式如下）。你知道我可以用 numpy 中的任何函数来计算这个 D(x) 吗？你能在 python.

中给出 D(x) 的示例代码吗

Video Credits

Answer 1

要更好地理解，请阅读 SIFT Taylor Expansion working out subpixel locations，稍后在下面找到真正实现 SIFT 的 C++ 代码。

static void interp_step( IplImage*** dog_pyr, int octv, int intvl, int r, int c,
             double* xi, double* xr, double* xc )
{
  CvMat* dD, * H, * H_inv, X;
  double x[3] = { 0 };

  dD = deriv_3D( dog_pyr, octv, intvl, r, c );
  H = hessian_3D( dog_pyr, octv, intvl, r, c );
  H_inv = cvCreateMat( 3, 3, CV_64FC1 );
  cvInvert( H, H_inv, CV_SVD );
  cvInitMatHeader( &X, 3, 1, CV_64FC1, x, CV_AUTOSTEP );
  cvGEMM( H_inv, dD, -1, NULL, 0, &X, 0 );

  cvReleaseMat( &dD );
  cvReleaseMat( &H );
  cvReleaseMat( &H_inv );

  *xi = x[2];
  *xr = x[1];
  *xc = x[0];
}


static CvMat* hessian_3D( IplImage*** dog_pyr, int octv, int intvl, int r,
              int c )
{
  CvMat* H;
  double v, dxx, dyy, dss, dxy, dxs, dys;

  v = pixval32f( dog_pyr[octv][intvl], r, c );
  dxx = ( pixval32f( dog_pyr[octv][intvl], r, c+1 ) + 
      pixval32f( dog_pyr[octv][intvl], r, c-1 ) - 2 * v );
  dyy = ( pixval32f( dog_pyr[octv][intvl], r+1, c ) +
      pixval32f( dog_pyr[octv][intvl], r-1, c ) - 2 * v );
  dss = ( pixval32f( dog_pyr[octv][intvl+1], r, c ) +
      pixval32f( dog_pyr[octv][intvl-1], r, c ) - 2 * v );
  dxy = ( pixval32f( dog_pyr[octv][intvl], r+1, c+1 ) -
      pixval32f( dog_pyr[octv][intvl], r+1, c-1 ) -
      pixval32f( dog_pyr[octv][intvl], r-1, c+1 ) +
      pixval32f( dog_pyr[octv][intvl], r-1, c-1 ) ) / 4.0;
  dxs = ( pixval32f( dog_pyr[octv][intvl+1], r, c+1 ) -
      pixval32f( dog_pyr[octv][intvl+1], r, c-1 ) -
      pixval32f( dog_pyr[octv][intvl-1], r, c+1 ) +
      pixval32f( dog_pyr[octv][intvl-1], r, c-1 ) ) / 4.0;
  dys = ( pixval32f( dog_pyr[octv][intvl+1], r+1, c ) -
      pixval32f( dog_pyr[octv][intvl+1], r-1, c ) -
      pixval32f( dog_pyr[octv][intvl-1], r+1, c ) +
      pixval32f( dog_pyr[octv][intvl-1], r-1, c ) ) / 4.0;

  H = cvCreateMat( 3, 3, CV_64FC1 );
  cvmSet( H, 0, 0, dxx );
  cvmSet( H, 0, 1, dxy );
  cvmSet( H, 0, 2, dxs );
  cvmSet( H, 1, 0, dxy );
  cvmSet( H, 1, 1, dyy );
  cvmSet( H, 1, 2, dys );
  cvmSet( H, 2, 0, dxs );
  cvmSet( H, 2, 1, dys );
  cvmSet( H, 2, 2, dss );

  return H;
}

参考： https://github.com/robwhess/opensift/blob/master/src/sift.c

如何将泰勒级数展开应用于opencv中的图像？

How to apply taylor series expansion to an image in opencv?

python

opencv

numpy

taylor-series

sift