Method processPoint

io/EptReader.cpp:808–890

This code runs in a single thread, so doesn't need locking.

// This code runs in a single thread, so doesn't need locking.
bool EptReader::processPoint(PointRef& dst, const ept::TileContents& tile)
{
    using namespace Dimension;

    BasePointTable& t = tile.table();

    // Save current point ID and increment so that we can return without
    // worrying about m_pointId being correct on exit.
    PointId pointId = m_pointId++;

    PointRef p(t, pointId);
    if (m_queryOriginId != -1 &&
        p.getFieldAs<int64_t>(Id::OriginId)!= m_queryOriginId)
    {
        return false;
    }

    auto passesBoundsFilter = [this](double x, double y, double z)
    {
        if (!m_p->bounds.box.valid())
            return true;
        m_p->bounds.xform.transform(x, y, z);
        return m_p->bounds.box.contains(x, y, z);
    };

    auto passesPolyFilter = [this](double xo, double yo, double zo)
    {
        if (m_p->polys.empty())
            return true;

        for (PolyXform& ps : m_p->polys)
        {
            double x = xo;
            double y = yo;
            double z = zo;

            ps.xform.transform(x, y, z);
            if (ps.poly.contains(x, y))
                return true;
        }
        return false;
    };

    double x = p.getFieldAs<double>(Id::X);
    double y = p.getFieldAs<double>(Id::Y);
    double z = p.getFieldAs<double>(Id::Z);

    // If there is a spatial filter, make sure it passes.
    if (m_p->hasSpatialFilter())
        if (!passesBoundsFilter(x, y, z) || !passesPolyFilter(x, y, z))
            return false;

    for (auto& el : m_p->info->dims())
    {
        DimType& dt = el.second;
        if (dt.m_id != Dimension::Id::X &&
                dt.m_id != Dimension::Id::Y &&
                dt.m_id != Dimension::Id::Z)