NormVecDecorator.h

Go to the documentation of this file.

/***************************************************************************
                          NormVecDecorator.h  -  description
                             -------------------
    copyright            : (C) 2004 by Marco Hugentobler
    email                : [email protected]
 ***************************************************************************/
 
/***************************************************************************
 *                                                                         *
 *   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 NORMVECDECORATOR_H
#define NORMVECDECORATOR_H
 
#include "TriDecorator.h"
#include "qgis_sip.h"
#include "TriangleInterpolator.h"
#include "MathUtils.h"
#include "qgslogger.h"
#include "qgis_analysis.h"
 
#define SIP_NO_FILE
 
class QgsFeedback;
 
class ANALYSIS_EXPORT NormVecDecorator: public TriDecorator
{
  public:
    enum PointState {Normal, BreakLine, EndPoint};
    NormVecDecorator();
    NormVecDecorator( QgsTriangulation *tin );
    ~NormVecDecorator() override;
    int addPoint( const QgsPoint &p ) override;
    bool calcNormal( double x, double y, QgsPoint &result SIP_OUT ) override;
    bool calcNormalForPoint( double x, double y, int pointIndex, Vector3D *result SIP_OUT );
    bool calcPoint( double x, double y, QgsPoint &result SIP_OUT ) override;
    void eliminateHorizontalTriangles() override;
    bool estimateFirstDerivative( int pointno );
    bool estimateFirstDerivatives( QgsFeedback *feedback = nullptr );
    Vector3D *getNormal( int n ) const;
    bool getTriangle( double x, double y, QgsPoint &p1 SIP_OUT, Vector3D *v1 SIP_OUT, QgsPoint &p2 SIP_OUT, Vector3D *v2 SIP_OUT, QgsPoint &p3 SIP_OUT, Vector3D *v3 SIP_OUT )  SIP_PYNAME( getTriangleVertices );
    bool getTriangle( double x, double y, QgsPoint &p1 SIP_OUT, int &ptn1 SIP_OUT, Vector3D *v1 SIP_OUT, PointState *state1 SIP_OUT, QgsPoint &p2 SIP_OUT, int &ptn2 SIP_OUT, Vector3D *v2 SIP_OUT, PointState *state2 SIP_OUT, QgsPoint &p3 SIP_OUT, int &ptn3 SIP_OUT, Vector3D *v3 SIP_OUT, PointState *state3 SIP_OUT );
    PointState getState( int pointno ) const;
    void setTriangleInterpolator( TriangleInterpolator *inter ) override;
    bool swapEdge( double x, double y ) override;
 
    bool saveTriangulation( QgsFeatureSink *sink, QgsFeedback *feedback = nullptr ) const override;
 
    QgsMesh triangulationToMesh( QgsFeedback *feedback = nullptr ) const override;
 
  protected:
    bool alreadyestimated;
    static const unsigned int DEFAULT_STORAGE_FOR_NORMALS = 100000;
    TriangleInterpolator *mInterpolator = nullptr;
    QVector<Vector3D *> *mNormVec;
    QVector<PointState> *mPointState;
    void setState( int pointno, PointState s );
 
  private:
    NormVecDecorator( const NormVecDecorator & ) = delete;
    NormVecDecorator &operator=( const NormVecDecorator & ) = delete;
};
 
#ifndef SIP_RUN
 
inline NormVecDecorator::NormVecDecorator()
  : mNormVec( new QVector<Vector3D*>( DEFAULT_STORAGE_FOR_NORMALS ) )
  , mPointState( new QVector<PointState>( DEFAULT_STORAGE_FOR_NORMALS ) )
{
  alreadyestimated = false;
}
 
inline NormVecDecorator::NormVecDecorator( QgsTriangulation *tin )
  : TriDecorator( tin )
  , mNormVec( new QVector<Vector3D*>( DEFAULT_STORAGE_FOR_NORMALS ) )
  , mPointState( new QVector<PointState>( DEFAULT_STORAGE_FOR_NORMALS ) )
{
  alreadyestimated = false;
}
 
inline void NormVecDecorator::setTriangleInterpolator( TriangleInterpolator *inter )
{
  mInterpolator = inter;
}
 
inline Vector3D *NormVecDecorator::getNormal( int n ) const
{
  if ( mNormVec )
  {
    return mNormVec->at( n );
  }
  else
  {
    QgsDebugMsg( QStringLiteral( "warning, null pointer" ) );
    return nullptr;
  }
}
 
#endif
 
#endif