1
0
forked from mirrors/0ad
Files
0ad/binaries/data/mods/public/globalscripts/interpolation.js
T
elexis a482f7bc1e Extend the cubicInterpolation function to consume a tension argument allowing to modulate the smoothness of the interpolation.
Thereby unify the chordal Catmull-Rom spline interpolation of the
ClumpPlacer (C++ bd53b14f58, JS 0e0ed94926),
the copy of that in the PathPlacer (bc805bd357, refs #892) and
the centripetal Catmull-Rom spline of the bicubicInterpolation function
from 93aefe0787, refs #4218
and don't claim the latter to be a uniform Catmull-Rom spline.

Reviewed in part by fatherbushido, discussed in 93aefe0787.

This was SVN commit r20383.
2017-10-30 20:43:01 +00:00

43 lines
1.4 KiB
JavaScript

/**
* Catmull-Rom spline.
* Interpolates the value of a point that is located between four known equidistant points with given values by
* constructing a polynomial of degree three that goes through all points.
*
* @param {Number} tension - determines the geometry of the curve, between 0 (tight curve) and 1 (smooth curve).
* Depending on the tension, the spline is called uniform (0), centripetal (0.5) or chordal (1).
* @param {Number} x - Location of the point to interpolate, relative to p1
*/
function cubicInterpolation(tension, x, p0, p1, p2, p3)
{
let P = -tension * p0 + (2 - tension) * p1 + (tension - 2) * p2 + tension * p3;
let Q = 2 * tension * p0 + (tension - 3) * p1 + (3 - 2 * tension) * p2 - tension * p3;
let R = -tension * p0 + tension * p2;
let S = p1;
return ((P * x + Q) * x + R) * x + S;
}
/**
* Two dimensional interpolation within a square grid using a polynomial of degree three.
*
* @param {Number} x, y - Location of the point to interpolate, relative to p11
*/
function bicubicInterpolation
(
x, y,
p00, p01, p02, p03,
p10, p11, p12, p13,
p20, p21, p22, p23,
p30, p31, p32, p33
)
{
let tension = 0.5;
return cubicInterpolation(
tension,
x,
cubicInterpolation(tension, y, p00, p01, p02, p03),
cubicInterpolation(tension, y, p10, p11, p12, p13),
cubicInterpolation(tension, y, p20, p21, p22, p23),
cubicInterpolation(tension, y, p30, p31, p32, p33));
}