src/analysis/interpolation/NormVecDecorator.h

Go to the documentation of this file.

/***************************************************************************
                          NormVecDecorator.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 NORMVECDECORATOR_H
#define NORMVECDECORATOR_H

#include "TriDecorator.h"
#include <TriangleInterpolator.h>
#include <MathUtils.h>
#include "qgslogger.h"
class QProgressDialog;

class ANALYSIS_EXPORT NormVecDecorator: public TriDecorator
{
  public:
    enum pointState {NORMAL, BREAKLINE, ENDPOINT};
    NormVecDecorator();
    NormVecDecorator( Triangulation* tin );
    virtual ~NormVecDecorator();
    int addPoint( Point3D* p );
    bool calcNormal( double x, double y, Vector3D* result );
    bool calcNormalForPoint( double x, double y, int point, Vector3D* result );
    bool calcPoint( double x, double y, Point3D* result );
    virtual void eliminateHorizontalTriangles();
    bool estimateFirstDerivative( int pointno );
    bool estimateFirstDerivatives( QProgressDialog* d = 0 );
    Vector3D* getNormal( int n ) const;
    bool getTriangle( double x, double y, Point3D* p1, Vector3D* v1, Point3D* p2, Vector3D* v2, Point3D* p3, Vector3D* v3 );
    bool getTriangle( double x, double y, Point3D* p1, int* ptn1, Vector3D* v1, pointState* state1, Point3D* p2, int* ptn2, Vector3D* v2, pointState* state2, Point3D* p3, int* ptn3, Vector3D* v3, pointState* state3 );
    pointState getState( int pointno ) const;
    void setTriangleInterpolator( TriangleInterpolator* inter );
    virtual bool swapEdge( double x, double y );
    virtual bool saveAsShapefile( const QString& fileName ) const;

  protected:
    bool alreadyestimated;
    const static unsigned int mDefaultStorageForNormals = 100000;
    TriangleInterpolator* mInterpolator;
    QVector<Vector3D*>* mNormVec;
    QVector<pointState>* mPointState;
    void setState( int pointno, pointState s );
};

inline NormVecDecorator::NormVecDecorator(): TriDecorator(), mInterpolator( 0 ), mNormVec( new QVector<Vector3D*>( mDefaultStorageForNormals ) ), mPointState( new QVector<pointState>( mDefaultStorageForNormals ) )
{
  alreadyestimated = false;
}

inline NormVecDecorator::NormVecDecorator( Triangulation* tin ): TriDecorator( tin ), mInterpolator( 0 ), mNormVec( new QVector<Vector3D*>( mDefaultStorageForNormals ) ), mPointState( new QVector<pointState>( mDefaultStorageForNormals ) )
{
  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( "warning, null pointer" );
    return 0;
  }
}

#endif