Function catmullRomSpline

jni/venus/region_operation.cpp:213–223

Cubic interpolation http://www.paulinternet.nl/?page=bicubic https://en.wikipedia.org/wiki/Cubic_Hermite_spline#Interpolation_on_the_unit_interval_without_exact_derivatives [0 2 0 0] [p0] f(n + t) = 0.5 * [1, t, t^2, t^3] * [-1 0 1 0] * [p1] [2 -5 4 -1] [p2] [-1 3 -3 1] [p3] */

										[-1  3 -3  1]   [p3]
*/
cv::Point2f catmullRomSpline(float t, const cv::Point2f& p0, const cv::Point2f& p1, const cv::Point2f& p2, const cv::Point2f& p3)
{
//	assert(0 <= t && t <= 1);
	float tt = t*t, ttt = tt*t;
	float b0 = -ttt + 2*tt - t;
	float b1 = 3*ttt - 5*tt + 2;
	float b2 = -3*ttt + 4*tt + t;
	float b3 = ttt - tt;
	
	return 0.5f * (b0*p0 + b1*p1 + b2*p2 + b3*p3);
}

cv::Mat susan(const cv::Mat& image, int radius, int tolerance)
{