qgsvector3d.h

Go to the documentation of this file.

/***************************************************************************
  qgsvector3d.h
  --------------------------------------
  Date                 : November 2017
  Copyright            : (C) 2017 by Martin Dobias
  Email                : wonder dot sk at gmail dot com
 ***************************************************************************
 *                                                                         *
 *   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 QGSVECTOR3D_H
#define QGSVECTOR3D_H
 
#include "qgis.h"
#include "qgis_core.h"
 
#include <QString>
#include <QVector3D>
 
using namespace Qt::StringLiterals;

class CORE_EXPORT QgsVector3D
{
  public:
    QgsVector3D() = default;

    QgsVector3D( double x, double y, double z )
      : mX( x ), mY( y ), mZ( z ) {}

    QgsVector3D( const QVector3D &v )
      : mX( v.x() ), mY( v.y() ), mZ( v.z() ) {}

    bool isNull() const SIP_HOLDGIL { return mX == 0 && mY == 0 && mZ == 0; }

    double x() const SIP_HOLDGIL { return mX; }
    double y() const SIP_HOLDGIL { return mY; }
    double z() const SIP_HOLDGIL { return mZ; }

    void setX( double x ) SIP_HOLDGIL { mX = x; }

    void setY( double y ) SIP_HOLDGIL { mY = y; }

    void setZ( double z ) SIP_HOLDGIL { mZ = z; }

    void set( double x, double y, double z ) SIP_HOLDGIL
    {
      mX = x;
      mY = y;
      mZ = z;
    }
 
    // TODO c++20 - replace with = default
    bool operator==( const QgsVector3D &other ) const SIP_HOLDGIL
    {
      return mX == other.mX && mY == other.mY && mZ == other.mZ;
    }
    bool operator!=( const QgsVector3D &other ) const SIP_HOLDGIL
    {
      return !operator==( other );
    }

    QgsVector3D operator+( const QgsVector3D &other ) const SIP_HOLDGIL
    {
      return QgsVector3D( mX + other.mX, mY + other.mY, mZ + other.mZ );
    }

    QgsVector3D &operator+=( const QgsVector3D &other ) SIP_HOLDGIL
    {
      mX += other.mX;
      mY += other.mY;
      mZ += other.mZ;
      return *this;
    }

    QgsVector3D operator-( const QgsVector3D &other ) const SIP_HOLDGIL
    {
      return QgsVector3D( mX - other.mX, mY - other.mY, mZ - other.mZ );
    }

    const QgsVector3D operator-() const SIP_HOLDGIL
    {
      return QgsVector3D( -mX, -mY, -mZ );
    }

    QgsVector3D operator *( const double factor ) const SIP_HOLDGIL
    {
 
      return QgsVector3D( mX * factor, mY * factor, mZ * factor );
    }

    QgsVector3D operator /( const double factor ) const SIP_HOLDGIL
    {
      return QgsVector3D( mX / factor, mY / factor, mZ / factor );
    }

    static double dotProduct( const QgsVector3D &v1, const QgsVector3D &v2 ) SIP_HOLDGIL
    {
      return v1.x() * v2.x() + v1.y() * v2.y() + v1.z() * v2.z();
    }

    static QgsVector3D crossProduct( const QgsVector3D &v1, const QgsVector3D &v2 ) SIP_HOLDGIL
    {
      return QgsVector3D( v1.y() * v2.z() - v1.z() * v2.y(),
                          v1.z() * v2.x() - v1.x() * v2.z(),
                          v1.x() * v2.y() - v1.y() * v2.x() );
    }

    double length() const SIP_HOLDGIL
    {
      return sqrt( mX * mX + mY * mY + mZ * mZ );
    }

    double lengthSquared() const SIP_HOLDGIL
    {
      return mX * mX + mY * mY + mZ * mZ;
    }

    void normalize() SIP_HOLDGIL
    {
      const double len = length();
      if ( !qgsDoubleNear( len, 0.0 ) )
      {
        mX /= len;
        mY /= len;
        mZ /= len;
      }
    }

    double distance( const QgsVector3D &other ) const SIP_HOLDGIL
    {
      return std::sqrt( ( mX - other.x() ) * ( mX - other.x() ) +
                        ( mY - other.y() ) * ( mY - other.y() ) +
                        ( mZ - other.z() ) * ( mZ - other.z() ) );
    }

    static QgsVector3D perpendicularPoint( const QgsVector3D &v1, const QgsVector3D &v2, const QgsVector3D &vp ) SIP_HOLDGIL
    {
      const QgsVector3D d = ( v2 - v1 ) / v2.distance( v1 );
      const QgsVector3D v = vp - v2;
      const double t = dotProduct( v, d );
      QgsVector3D P = v2 + ( d * t );
      return P;
    }

    QString toString( int precision = 17 ) const SIP_HOLDGIL
    {
      QString str = "Vector3D (";
      str += qgsDoubleToString( mX, precision );
      str += ", ";
      str += qgsDoubleToString( mY, precision );
      str += ", ";
      str += qgsDoubleToString( mZ, precision );
      str += ')';
      return str;
    }

    QVector3D toVector3D() const SIP_HOLDGIL { return QVector3D( static_cast< float >( mX ), static_cast< float >( mY ), static_cast< float >( mZ ) ); }
 
#ifdef SIP_RUN
    SIP_PYOBJECT __repr__();
    % MethodCode
    QString str = u"<QgsVector3D: %1>"_s.arg( sipCpp->toString() );
    sipRes = PyUnicode_FromString( str.toUtf8().constData() );
    % End
#endif
  private:
    double mX = 0, mY = 0, mZ = 0;
};
 
#endif // QGSVECTOR3D_H