CloughTocherInterpolator.h

Go to the documentation of this file.

/***************************************************************************
                          CloughTocherInterpolator.h  -  description
                             -------------------
    copyright            : (C) 2004 by Marco Hugentobler
    email                : mhugent@geo.unizh.ch
 ***************************************************************************/
 
/***************************************************************************
 *                                                                         *
 *   This program is free software; you can redistribute it and/or modify  *
 *   it under the terms of the GNU General Public License as published by  *
 *   the Free Software Foundation; either version 2 of the License, or     *
 *   (at your option) any later version.                                   *
 *                                                                         *
 ***************************************************************************/
 
#ifndef CLOUGHTOCHERINTERPOLATOR_H
#define CLOUGHTOCHERINTERPOLATOR_H
 
#include "TriangleInterpolator.h"
#include "qgspoint.h"
#include "qgis_analysis.h"
 
class NormVecDecorator;
 
#define SIP_NO_FILE
 
class ANALYSIS_EXPORT CloughTocherInterpolator : public TriangleInterpolator
{
  protected:
    NormVecDecorator *mTIN = nullptr;
    double mEdgeTolerance = 0.00001;
    QgsPoint point1 = QgsPoint( 0, 0, 0 );
    QgsPoint point2 = QgsPoint( 0, 0, 0 );
    QgsPoint point3 = QgsPoint( 0, 0, 0 );
    QgsPoint cp1 = QgsPoint( 0, 0, 0 );
    QgsPoint cp2 = QgsPoint( 0, 0, 0 );
    QgsPoint cp3 = QgsPoint( 0, 0, 0 );
    QgsPoint cp4 = QgsPoint( 0, 0, 0 );
    QgsPoint cp5 = QgsPoint( 0, 0, 0 );
    QgsPoint cp6 = QgsPoint( 0, 0, 0 );
    QgsPoint cp7 = QgsPoint( 0, 0, 0 );
    QgsPoint cp8 = QgsPoint( 0, 0, 0 );
    QgsPoint cp9 = QgsPoint( 0, 0, 0 );
    QgsPoint cp10 = QgsPoint( 0, 0, 0 );
    QgsPoint cp11 = QgsPoint( 0, 0, 0 );
    QgsPoint cp12 = QgsPoint( 0, 0, 0 );
    QgsPoint cp13 = QgsPoint( 0, 0, 0 );
    QgsPoint cp14 = QgsPoint( 0, 0, 0 );
    QgsPoint cp15 = QgsPoint( 0, 0, 0 );
    QgsPoint cp16 = QgsPoint( 0, 0, 0 );
    double der1X = 0.0;
    double der1Y = 0.0;
    double der2X = 0.0;
    double der2Y = 0.0;
    double der3X = 0.0;
    double der3Y = 0.0;
    QgsPoint lpoint1 = QgsPoint( 0, 0, 0 );
    QgsPoint lpoint2 = QgsPoint( 0, 0, 0 );
    QgsPoint lpoint3 = QgsPoint( 0, 0, 0 );
    void init( double x, double y );
    double calcBernsteinPoly( int n, int i, int j, int k, double u, double v, double w );
 
  public:
    CloughTocherInterpolator() = default;
 
    CloughTocherInterpolator( NormVecDecorator *tin );
 
    bool calcNormVec( double x, double y, QgsPoint &result SIP_OUT ) override;
    bool calcPoint( double x, double y, QgsPoint &result SIP_OUT ) override;
    virtual void setTriangulation( NormVecDecorator *tin );
};
 
#endif