Method compute

filters/NormalFilter.cpp:140–182

}

void NormalFilter::compute(PointView& view, KD3Index& kdi)
{
    log()->get(LogLevel::Debug) << "Computing normal vectors\n";
    for (auto&& p : view)
    {
        math::NormalResult result;

        // Perform eigen decomposition of covariance matrix computed from
        // neighborhood composed of k-nearest neighbors, or within radius.
        if (m_radiusArg->set())
            result = math::findNormal(view, kdi.radius(p.pointId(), m_args->m_radius));
        else
            result = math::findNormal(view, kdi.neighbors(p.pointId(), m_args->m_knn));

        if (result.normal.isZero())
        {
            log()->get(LogLevel::Info)
                << "Skipping point " << p.pointId()
                << ": " << result.msg << "\n";
            continue;
        }

        if (m_viewpointArg->set())
        {
            // If a viewpoint has been specified, orient the normals to face the
            // viewpoint by taking the dot product of the vector connecting the
            // point with the viewpoint and the normal. Flip the normal, where
            // the dot product is negative.
            double dx = m_args->m_viewpoint.x() - p.getFieldAs<double>(Id::X);
            double dy = m_args->m_viewpoint.y() - p.getFieldAs<double>(Id::Y);
            double dz = m_args->m_viewpoint.z() - p.getFieldAs<double>(Id::Z);
            result.normal = math::orientToViewpoint({ dx, dy, dz }, result.normal);
        }
        else if (m_args->m_up)
            result.normal = math::orientUp(result.normal);

        // Set the computed normal and curvature dimensions.
        p.setField(Id::NormalX, result.normal[0]);
        p.setField(Id::NormalY, result.normal[1]);
        p.setField(Id::NormalZ, result.normal[2]);
        p.setField(Id::Curvature, result.curvature);
    }
}

void NormalFilter::update(
    PointView& view, KD3Index& kdi, std::vector<bool> inMST,