Method filterPoints

filters/M3C2Filter.cpp:305–337

}

std::vector<double> M3C2Filter::filterPoints(Eigen::Vector3d cylCenter, Eigen::Vector3d cylNormal,
    const PointView& view, const PointIdList& neighbors)
{
    std::vector<double> dists;

    if (!neighbors.size())
        return dists;
    dists.reserve(neighbors.size());

    Eigen::Vector3d point(view.getFieldAs<double>(Dimension::Id::X, 0),
        view.getFieldAs<double>(Dimension::Id::Y, 0),
        view.getFieldAs<double>(Dimension::Id::Z, 0));

    size_t start = 0;
    // If the first point is the test point, ignore it.
    if (Comparison::closeEnough(cylCenter(0), point(0)) &&
    Comparison::closeEnough(cylCenter(1), point(1)) &&
    Comparison::closeEnough(cylCenter(2), point(2)))
        start++;

    for (size_t i = start; i < neighbors.size(); i++)
    {
        PointId id = neighbors[i];
        Eigen::Vector3d point(view.getFieldAs<double>(Dimension::Id::X, id),
            view.getFieldAs<double>(Dimension::Id::Y, id),
            view.getFieldAs<double>(Dimension::Id::Z, id));

        double dist = pointPasses(point, cylCenter, cylNormal);
        if (!std::isnan(dist))
            dists.push_back(dist);
    }
    return dists;
}

double M3C2Filter::pointPasses(Eigen::Vector3d point, Eigen::Vector3d cylCenter,
    Eigen::Vector3d cylNormal)