Function TEST

test/unit/filters/CovarianceFeaturesTest.cpp:46–90

namespace pdal {

TEST(DimensionalityTest, Linearity)
{
    using namespace Dimension;

    PointTable table;
    table.layout()->registerDims({Id::X, Id::Y, Id::Z});

    BufferReader bufferReader;
    CovarianceFeaturesFilter filter;
    Options ops;
    ops.add("knn", 3);
    ops.add("feature_set", "Dimensionality,Omnivariance,Anisotropy,Eigenentropy,EigenvalueSum,SurfaceVariation,DemantkeVerticality");
    filter.setInput(bufferReader);
    filter.setOptions(ops);
    filter.prepare(table);

    PointViewPtr view(new PointView(table));
    view->setField(Id::X, 0,0);
    view->setField(Id::X, 1, 0);
    view->setField(Id::X, 2, 0);
    view->setField(Id::Y, 0, 0);
    view->setField(Id::Y, 1, 0);
    view->setField(Id::Y, 2, 0);
    view->setField(Id::Z, 0, 0);
    view->setField(Id::Z, 1, 1);
    view->setField(Id::Z, 2, 2);
    bufferReader.addView(view);

    PointViewSet viewSet = filter.execute(table);
    PointViewPtr outView = *viewSet.begin();

    for (point_count_t i =0; i < outView->size(); i++)
    {
        ASSERT_FLOAT_EQ(outView->getFieldAs<float>(Id::Linearity, i) ,1);
        ASSERT_FLOAT_EQ(outView->getFieldAs<float>(Id::Verticality, i) ,1);
        ASSERT_FLOAT_EQ(outView->getFieldAs<float>(Id::Planarity, i) ,0);
        ASSERT_FLOAT_EQ(outView->getFieldAs<float>(Id::Scattering, i) ,0);
        ASSERT_FLOAT_EQ(outView->getFieldAs<float>(Id::Omnivariance, i), 0);
        ASSERT_FLOAT_EQ(outView->getFieldAs<float>(Id::Anisotropy, i), 1);
        ASSERT_TRUE(std::isnan(outView->getFieldAs<float>(Id::Eigenentropy, i)));
        ASSERT_FLOAT_EQ(outView->getFieldAs<float>(Id::EigenvalueSum, i), 1);
        ASSERT_FLOAT_EQ(outView->getFieldAs<float>(Id::SurfaceVariation, i), 0);
        ASSERT_FLOAT_EQ(outView->getFieldAs<float>(Id::DemantkeVerticality, i), 1);
    }
}

TEST(DimensionalityTest, Planarity)
{