qgspointxy.h

Go to the documentation of this file.

/***************************************************************************
                          qgspoint.h  -  description
                             -------------------
    begin                : Sat Jun 22 2002
    copyright            : (C) 2002 by Gary E.Sherman
    email                : sherman at mrcc.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 QGSPOINTXY_H
#define QGSPOINTXY_H
 
#include "qgis_core.h"
#include "qgsvector.h"
 
#include "qgis.h"
 
#include <iostream>
#include <QString>
#include <QPoint>
#include <QObject>
 
class QgsPoint;
 
class CORE_EXPORT QgsPointXY
{
    Q_GADGET
 
    Q_PROPERTY( double x READ x WRITE setX )
    Q_PROPERTY( double y READ y WRITE setY )
 
  public:
    QgsPointXY() = default;
 
    QgsPointXY( const QgsPointXY &p );
 
    QgsPointXY( double x, double y )
      : mX( x )
      , mY( y )
      , mIsEmpty( false )
    {}
 
    QgsPointXY( QPointF point )
      : mX( point.x() )
      , mY( point.y() )
      , mIsEmpty( false )
    {}
 
    QgsPointXY( QPoint point )
      : mX( point.x() )
      , mY( point.y() )
      , mIsEmpty( false )
    {}
 
    QgsPointXY( const QgsPoint &point );
 
    // IMPORTANT - while QgsPointXY is inherited by QgsReferencedPointXY, we do NOT want a virtual destructor here
    // because this class MUST be lightweight and we don't want the cost of the vtable here.
    // see https://github.com/qgis/QGIS/pull/4720#issuecomment-308652392
    ~QgsPointXY() = default;
 
    void setX( double x )
    {
      mX = x;
      mIsEmpty = false;
    }
 
    void setY( double y )
    {
      mY = y;
      mIsEmpty = false;
    }
 
    void set( double x, double y )
    {
      mX = x;
      mY = y;
      mIsEmpty = false;
    }
 
    double x() const
    {
      return mX;
    }
 
    double y() const
    {
      return mY;
    }
 
    QPointF toQPointF() const
    {
      return QPointF( mX, mY );
    }
 
    QString toString( int precision = -1 ) const;
 
    QString asWkt() const;
 
    double sqrDist( double x, double y ) const
    {
      return ( mX - x ) * ( mX - x ) + ( mY - y ) * ( mY - y );
    }
 
    double sqrDist( const QgsPointXY &other ) const
    {
      return sqrDist( other.x(), other.y() );
    }
 
    double distance( double x, double y ) const
    {
      return std::sqrt( sqrDist( x, y ) );
    }
 
    double distance( const QgsPointXY &other ) const
    {
      return std::sqrt( sqrDist( other ) );
    }
 
    double sqrDistToSegment( double x1, double y1, double x2, double y2, QgsPointXY &minDistPoint SIP_OUT, double epsilon = DEFAULT_SEGMENT_EPSILON ) const;
 
    double azimuth( const QgsPointXY &other ) const;
 
    QgsPointXY project( double distance, double bearing ) const;
 
    bool isEmpty() const { return mIsEmpty; }
 
    bool compare( const QgsPointXY &other, double epsilon = 4 * std::numeric_limits<double>::epsilon() ) const
    {
      return ( qgsDoubleNear( mX, other.x(), epsilon ) && qgsDoubleNear( mY, other.y(), epsilon ) );
    }
 
    bool operator==( const QgsPointXY &other )
    {
      if ( isEmpty() && other.isEmpty() )
        return true;
      if ( isEmpty() && !other.isEmpty() )
        return false;
      if ( ! isEmpty() && other.isEmpty() )
        return false;
 
      bool equal = true;
      equal &= qgsDoubleNear( other.x(), mX, 1E-8 );
      equal &= qgsDoubleNear( other.y(), mY, 1E-8 );
 
      return equal;
    }
 
    bool operator!=( const QgsPointXY &other ) const
    {
      if ( isEmpty() && other.isEmpty() )
        return false;
      if ( isEmpty() && !other.isEmpty() )
        return true;
      if ( ! isEmpty() && other.isEmpty() )
        return true;
 
      bool equal = true;
      equal &= qgsDoubleNear( other.x(), mX, 1E-8 );
      equal &= qgsDoubleNear( other.y(), mY, 1E-8 );
 
      return !equal;
    }
 
    void multiply( double scalar )
    {
      mX *= scalar;
      mY *= scalar;
    }
 
    QgsPointXY &operator=( const QgsPointXY &other )
    {
      if ( &other != this )
      {
        mX = other.x();
        mY = other.y();
        mIsEmpty = other.isEmpty();
      }
 
      return *this;
    }
 
    QgsVector operator-( const QgsPointXY &p ) const { return QgsVector( mX - p.mX, mY - p.mY ); }
 
    QgsPointXY &operator+=( QgsVector v ) { *this = *this + v; return *this; }
 
    QgsPointXY &operator-=( QgsVector v ) { *this = *this - v; return *this; }
 
    QgsPointXY operator+( QgsVector v ) const { return QgsPointXY( mX + v.x(), mY + v.y() ); }
 
    QgsPointXY operator-( QgsVector v ) const { return QgsPointXY( mX - v.x(), mY - v.y() ); }
 
    QgsPointXY operator*( double scalar ) const { return QgsPointXY( mX * scalar, mY * scalar ); }
 
    QgsPointXY operator/( double scalar ) const { return QgsPointXY( mX / scalar, mY / scalar ); }
 
    QgsPointXY &operator*=( double scalar ) { mX *= scalar; mY *= scalar; return *this; }
 
    QgsPointXY &operator/=( double scalar ) { mX /= scalar; mY /= scalar; return *this; }
 
    operator QVariant() const
    {
      return QVariant::fromValue( *this );
    }
 
#ifdef SIP_RUN
    SIP_PYOBJECT __repr__();
    % MethodCode
    QString str = QStringLiteral( "<QgsPointXY: %1>" ).arg( sipCpp->asWkt() );
    sipRes = PyUnicode_FromString( str.toUtf8().constData() );
    % End
 
    int __len__();
    % MethodCode
    sipRes = 2;
    % End
 
 
    SIP_PYOBJECT __getitem__( int );
    % MethodCode
    if ( a0 == 0 )
    {
      sipRes = Py_BuildValue( "d", sipCpp->x() );
    }
    else if ( a0 == 1 )
    {
      sipRes = Py_BuildValue( "d", sipCpp->y() );
    }
    else
    {
      QString msg = QString( "Bad index: %1" ).arg( a0 );
      PyErr_SetString( PyExc_IndexError, msg.toLatin1().constData() );
    }
    % End
 
    long __hash__() const;
    % MethodCode
    sipRes = qHash( *sipCpp );
    % End
#endif
 
  private:
 
    double mX = 0; //std::numeric_limits<double>::quiet_NaN();
 
    double mY = 0; //std::numeric_limits<double>::quiet_NaN();
 
    bool mIsEmpty = true;
 
    friend uint qHash( const QgsPointXY &pnt );
 
}; // class QgsPointXY
 
Q_DECLARE_METATYPE( QgsPointXY )
 
inline bool operator==( const QgsPointXY &p1, const QgsPointXY &p2 ) SIP_SKIP
{
  const bool nan1X = std::isnan( p1.x() );
  const bool nan2X = std::isnan( p2.x() );
  if ( nan1X != nan2X )
    return false;
  if ( !nan1X && !qgsDoubleNear( p1.x(), p2.x(), 1E-8 ) )
    return false;
 
  const bool nan1Y = std::isnan( p1.y() );
  const bool nan2Y = std::isnan( p2.y() );
  if ( nan1Y != nan2Y )
    return false;
 
  if ( !nan1Y && !qgsDoubleNear( p1.y(), p2.y(), 1E-8 ) )
    return false;
 
  return true;
}
 
inline std::ostream &operator << ( std::ostream &os, const QgsPointXY &p ) SIP_SKIP
{
  // Use Local8Bit for printouts
  os << p.toString().toLocal8Bit().data();
  return os;
}
 
inline uint qHash( const QgsPointXY &p ) SIP_SKIP
{
  uint hash;
  uint h1 = qHash( static_cast< quint64 >( p.mX ) );
  uint h2 = qHash( static_cast< quint64 >( p.mY ) );
  hash = h1 ^ ( h2 << 1 );
  return hash;
}
 
 
#endif //QGSPOINTXY_H