Method done

filters/HexBinFilter.cpp:224–358

}

void HexBin::done(PointTableRef table)
{
    if (!createGrid())
        return;

    m_metadata.add("threshold", m_grid->denseLimit(),
        "Minimum number of points inside a hexagon to be considered full");
    m_metadata.add("sample_size", m_sampleSize, "Number of samples used for "
        "estimating hexagon edge size. Only used if 'edge_length' or "
        "'h3_resolution' is not set.");

    Utils::OStringStreamClassicLocale polygon;
    polygon.setf(std::ios_base::fixed, std::ios_base::floatfield);
    polygon.precision(m_precision);
    m_grid->toWKT(polygon);

    if (m_outputTesselation)
    {
        m_metadata.add("hex_boundary", polygon.str(),
            "Boundary MULTIPOLYGON of domain");
    }

    // density and boundary writing with OGR does not support polygon smoothing
    if (m_DensityOutput.size())
    {
        OGR writer(m_DensityOutput, m_srs.getWKT(), m_driver, "hexbins");
        writer.writeDensity(*m_grid);
    }
    if (m_boundaryOutput.size())
    {
        OGR writer(m_boundaryOutput, m_srs.getWKT(), m_driver, "hexbins");
        writer.writeBoundary(*m_grid);
    }

    pdal::Polygon p(polygon.str(), m_srs);

    // If the SRS was geographic, use relevant
    // UTM for area and density computation
    Polygon density_p(p);
    if (m_srs.isGeographic())
    {
        // Compute a UTM polygon
        BOX3D box = p.bounds();
        int zone = SpatialReference::calculateZone(box.minx, box.miny);
        if (!density_p.transform(SpatialReference::wgs84FromZone(zone)))
            density_p = Polygon();
    }

    double area = density_p.area();
    double density = m_count / area;
    if (std::isinf(density))
    {
        density = -1.0;
        area = -1.0;
    }

    m_metadata.add("density", density,
        "Number of points per square unit (total area)");