FittedCurve.hpp

Go to the documentation of this file.

#pragma once
 
#include <cassert>
#include <cstddef>
#include <vector>
 
#include "vc/segmentation/lrps/Common.hpp"
#include "vc/segmentation/lrps/Spline.hpp"
 
namespace volcart::segmentation
{
class FittedCurve
{
private:
    std::size_t npoints_{0};
    int zIndex_{0};
    std::vector<double> ts_;
    std::vector<Voxel> points_;
    CubicSpline<double> spline_;
 
public:
    FittedCurve() = default;
 
    FittedCurve(const std::vector<Voxel>& vs, int zIndex);
    std::size_t size() const { return npoints_; }
 
    const std::vector<Voxel>& points() const { return points_; }
 
    const decltype(spline_)& spline() const { return spline_; }
 
    Pixel eval(double t) const { return spline_(t); }
 
    std::vector<Voxel> evenlySpacePoints() { return resample(1.0); }
 
    std::vector<Voxel> resample(double resamplePerc = 1.0);
 
    std::vector<Voxel> sample(std::size_t numPoints) const;
 
    Voxel operator()(int index) const;
 
    std::vector<double> curvature(int hstep = 1) const;
 
    double arclength() const;
};
}  // namespace volcart::segmentation