FittedCurve.hpp

Go to the documentation of this file.

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