Hi all,
I am trying to interpolate some data with a catmull_rom spline available in
boost math (interpolation). However, the resulting spline has a strange
behavior; especially, at its ends (see the attached image).
Just out of curiosity, I also tried to normalise my data to [0-1] range.
However, the line has had even worse behavior (see the attached image).
Below, the c++ code segment of the interpolation is depicted:
std::vector> points_intr(16);
points_intr[0] = { 1.263270, 0.774614 };
points_intr[1] = { 1.876, 2.480 };
points_intr[2] = { 1.651110, 4.550 };
points_intr[3] = { 1.426340, 5.688 };
points_intr[4] = { 1.429, 7.054 };
points_intr[5] = { 2.073220, 8.377020 };
points_intr[6] = { 3.910, 9.140 };
points_intr[7] = { 6.430, 9.537 };
points_intr[8] = { 8.950, 9.859 };
points_intr[9] = { 11.470, 10.317 };
points_intr[10] = { 12.730, 10.6456 };
points_intr[11] = { 13.990, 11.0741 };
points_intr[12] = { 15.335, 11.6928 };
points_intr[13] = { 16.680, 12.5661 };
points_intr[14] = { 18.3538, 14.830 };
points_intr[15] = { 18.700, 16.056 };
bm::catmull_rom>
interpolator_cr(std::move(points_intr));
std::vector<double> x_linspace_cr(100);
std::vector> inter_points(100);
linspace(1.263270, 18.700, 100, x_linspace_cr); // evaluation points
along the line
for (std::size_t i{0}; i < x_linspace_cr.size(); ++i)
{
inter_points[i] = interpolator_cr(x_linspace_cr[i]);
}
Should I have to add additional (ghost) points at both ends? If yes, is
there any optimum way to calculate where exactly to poisition these points?
I understand if I get the derivatives at both end I can calculate a
straight line that the additional points should lie within. However, how
should I calculate the distance from the end points?
I hope I am not doing anything fundamentally wrong. Any help would be much
appreciated.
Many thanks!
