Method inspect

io/CopcReader.cpp:573–617

QuickInfo CopcReader::inspect()
{
    PointTable t;
    QuickInfo qi;

    initialize(t);

    const las::Header& h = m_p->header;
    qi.m_bounds = h.bounds;
    qi.m_srs = getSpatialReference();
    qi.m_pointCount = h.pointCount();

    PointLayout layout;
    addDimensions(&layout);
    for (Dimension::Id dim : layout.dims())
        qi.m_dimNames.push_back(layout.dimName(dim));

    // If there is a spatial filter from an explicit --bounds, an origin query,
    // or polygons, then we'll limit our number of points to be an upper bound,
    // and clip our bounds to the selected region.
    if (hasSpatialFilter())
    {
        loadHierarchy();

        qi.m_pointCount = m_p->hierarchy.pointCount();

        //ABELL - This is wrong since we're not transforming the tile bounds to the
        //  SRS of each clip region, but that seems like a lot of mess for
        //  little value. Wait until someone complains. (Note that's it's a bit
        //  different from queryOverlaps or we'd just call that.)
        // Clip the resulting bounds to the intersection of:
        //  - the query bounds (from an explicit bounds or an origin query)
        //  - the extents of the polygon selection
        BOX3D b;
        b.grow(m_p->clip.box);
        for (const auto& poly : m_args->polys)
            b.grow(poly.bounds());

        if (b.valid())
            qi.m_bounds.clip(b);
    }
    qi.m_valid = true;

    return qi;
}


void CopcReader::addDimensions(PointLayoutPtr layout)