gageassignment.h

Go to the documentation of this file.

 
#ifndef OPENSWMMVIS_CORE_GAGEASSIGNMENT_H
#define OPENSWMMVIS_CORE_GAGEASSIGNMENT_H
 
#include <QPair>
#include <QPointF>
#include <QString>
#include <QVector>
 
namespace GageAssignment
{
 
constexpr double kWeightEpsilon = 1e-4;
 
constexpr double kDefaultClusterTol = 0.01;
 
constexpr int kDefaultSampleTarget = 200;
 
// ===========================================================================
// Thiessen area-majority
// ===========================================================================
 
[[nodiscard]] QVector<QPointF> clipHalfPlane(const QVector<QPointF> &ring,
                                             const QPointF &keep,
                                             const QPointF &drop);
 
[[nodiscard]] QVector<double> thiessenAreaShares(const QVector<QPointF> &ring,
                                                 const QVector<QPointF> &gages);
 
[[nodiscard]] int areaMajorityGage(const QVector<double> &shares,
                                   double *fractionOut = nullptr);
 
// ===========================================================================
// Interpolation support
// ===========================================================================
 
[[nodiscard]] QVector<QPointF> samplePolygon(const QVector<QPointF> &ring,
                                             int target = kDefaultSampleTarget);
 
[[nodiscard]] QVector<QPair<int, double>> idwWeights(const QPointF &p,
                                                     const QVector<QPointF> &sites);
 
// ===========================================================================
// Weight-vector clustering
// ===========================================================================
 
struct ClusterKey
{
    QVector<QPair<int, long long>> terms;      
    QString                        serialized; 
 
    [[nodiscard]] bool isEmpty() const { return terms.isEmpty(); }
};
 
[[nodiscard]] ClusterKey quantizeWeights(const QVector<double> &w,
                                         double tol  = kDefaultClusterTol,
                                         double wEps = kWeightEpsilon);
 
[[nodiscard]] QVector<double> dequantizeWeights(const ClusterKey &key,
                                                double tol,
                                                int nGages);
 
} // namespace GageAssignment
 
#endif // OPENSWMMVIS_CORE_GAGEASSIGNMENT_H