Function Math_linear_interpolate

Code/Source/PythonAPI/Math_PyModule.cxx:311–357

");

static PyObject *
Math_linear_interpolate(PyObject *self, PyObject *args)
{
  auto api = PyUtilApiFunction("Oi", PyRunTimeErr, __func__);
  PyObject *pointsArg;
  int numInterpPoints = 0;

  if (!PyArg_ParseTuple(args, api.format, &pointsArg, &numInterpPoints)) {
      return api.argsError();
  }

  // Get an array of points from the pointsArgs list.
  int dim = 2;
  auto pts = GetPointsFromList(api, pointsArg, dim, "points");
  if (pts == nullptr) {
      return nullptr;
  }
  int nlistpts = PyList_Size(pointsArg);

  // Calculate dt so that our time series will go from 0 to T.
  double t0 = pts[0][0];
  double dt = (pts[nlistpts-1][0]-t0)/(numInterpPoints-1);
  double **outPts = NULL;

  // Perform the linear interpolation.
  auto mathObj = cvMath();
  if (mathObj.linearInterpolate(pts, nlistpts, t0, dt, numInterpPoints, &outPts) == SV_ERROR) {
       mathObj.deleteArray(pts,nlistpts,dim);
       api.error("Error linear interplating points.");
       return nullptr;
  }

  // create result string
  char r[2048];
  PyObject *pylist = PyList_New(numInterpPoints);
  for (int i = 0; i < numInterpPoints; i++) {
    r[0] = '\0';
    sprintf(r,"%.6le %.6le",outPts[i][0],outPts[i][1]);
    PyObject *rr = PyBytes_FromString(r);
    PyList_SET_ITEM(pylist, i, rr);
  }

  mathObj.deleteArray(pts,nlistpts,dim);
  mathObj.deleteArray(outPts,numInterpPoints,dim);

  return pylist;
}

//---------------------
// Math_curve_length