pslgprep.h

Go to the documentation of this file.

 
#ifndef OPENSWMMVIS_MESH_PSLGPREP_H
#define OPENSWMMVIS_MESH_PSLGPREP_H
 
#include "mesh/meshgenerator.h"
 
#include <QPointF>
#include <QPolygonF>
#include <QRectF>
#include <QVector>
 
#include <functional>
 
namespace mesh {
namespace pslg {
 
[[nodiscard]] QVector<QPointF> simplifyPolyline(const QVector<QPointF> &pts,
                                                double epsilon);
 
[[nodiscard]] QVector<QPointF> simplifyRing(const QVector<QPointF> &ring,
                                            double epsilon);
 
[[nodiscard]] QVector<QPointF> densifyRing(const QVector<QPointF> &ring,
                                           double maxLen);
 
[[nodiscard]] double polylineLength(const QVector<QPointF> &pts);
 
[[nodiscard]] QVector<QPointF> resampleMinLength(const QVector<QPointF> &pts,
                                                 double minLen,
                                                 double maxDeviation,
                                                 int *flaggedOut = nullptr);
 
[[nodiscard]] QVector<QPointF> resampleRingMinLength(const QVector<QPointF> &ring,
                                                     double minLen,
                                                     double maxDeviation,
                                                     int *flaggedOut = nullptr);
 
[[nodiscard]] double ringSignedArea(const QVector<QPointF> &ring);
 
[[nodiscard]] double distSqToSegment(const QPointF &p, const QPointF &a,
                                     const QPointF &b);
 
void snapAndDedupe(QVector<mesh::SteinerPoint> &pts, double snapEps);
 
[[nodiscard]] QVector<bool> greedyMinSeparation(const QVector<QPointF> &pts,
                                                double minSep);
 
struct PreparedRing
{
    QVector<QPointF> ring;   
    QPointF          seed;   
    bool             valid = false;  
};
 
[[nodiscard]] PreparedRing prepareHoleRing(const QVector<QPointF> &raw,
                                           double simplifyEps,
                                           double maxEdgeLen);
 
bool prepareHoleRings(const QVector<QVector<QPointF>> &raw,
                      double simplifyEps, double maxEdgeLen,
                      QVector<PreparedRing> *out,
                      const std::function<bool()> &isCancelled = {},
                      const std::function<void(int, int)> &onChunk = {},
                      int *skippedOut = nullptr);
 
class PointInRingsIndex
{
public:
    void build(const QVector<QPolygonF>        &domains,
               const QVector<QVector<QPointF>> &holes  = {},
               int                              nBands = 1024);
 
    [[nodiscard]] bool   contains(const QPointF &p) const;
    [[nodiscard]] QRectF boundingBox() const;
    [[nodiscard]] bool   isEmpty() const { return m_edges.isEmpty(); }
 
private:
    void addRing(const QVector<QPointF> &ring);
 
    QVector<QPair<QPointF, QPointF>> m_edges;
    QVector<QVector<int>>            m_bandEdges;
    double m_x0 = 0.0, m_x1 = 0.0, m_y0 = 0.0, m_y1 = 0.0;
    double m_bandH  = 0.0;
    int    m_nBands = 0;
};
 
} // namespace pslg
} // namespace mesh
 
#endif // OPENSWMMVIS_MESH_PSLGPREP_H