qgsgcptransformer.h

Go to the documentation of this file.

/***************************************************************************
    qgsgcptransformer.h
     --------------------------------------
    Date                 : 18-Feb-2009
    Copyright            : (c) 2009 by Manuel Massing
    Email                : m.massing at warped-space.de
 ***************************************************************************
 *                                                                         *
 *   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 QGSGCPTRANSFORMER_H
#define QGSGCPTRANSFORMER_H
 
#include <gdal_alg.h>
#include "qgspoint.h"
#include "qgis_analysis.h"
#include "qgis_sip.h"
 
class ANALYSIS_EXPORT QgsGcpTransformerInterface SIP_ABSTRACT
{
    Q_GADGET
 
  public:
 
    enum class TransformMethod : int
    {
      Linear, 
      Helmert, 
      PolynomialOrder1, 
      PolynomialOrder2, 
      PolynomialOrder3, 
      ThinPlateSpline, 
      Projective, 
      InvalidTransform = 65535 
    };
    Q_ENUM( TransformMethod )
 
    
    QgsGcpTransformerInterface() = default;
 
    virtual ~QgsGcpTransformerInterface() = default;
 
    QgsGcpTransformerInterface( const QgsGcpTransformerInterface &other ) = delete;
 
    QgsGcpTransformerInterface &operator=( const QgsGcpTransformerInterface &other ) = delete;
 
    virtual QgsGcpTransformerInterface *clone() const = 0 SIP_FACTORY;
 
    virtual bool updateParametersFromGcps( const QVector<QgsPointXY> &sourceCoordinates, const QVector<QgsPointXY> &destinationCoordinates, bool invertYAxis = false ) SIP_THROW( QgsNotSupportedException ) = 0;
 
    virtual int minimumGcpCount() const = 0;
 
    virtual TransformMethod method() const = 0;
 
    bool transform( double &x SIP_INOUT, double &y SIP_INOUT, bool inverseTransform = false ) const;
 
    static QString methodToString( TransformMethod method );
 
    static QgsGcpTransformerInterface *create( TransformMethod method ) SIP_FACTORY;
 
    static QgsGcpTransformerInterface *createFromParameters( TransformMethod method, const QVector<QgsPointXY> &sourceCoordinates, const QVector<QgsPointXY> &destinationCoordinates ) SIP_THROW( QgsNotSupportedException ) SIP_FACTORY;
 
#ifndef SIP_RUN
 
    virtual GDALTransformerFunc GDALTransformer() const = 0;
 
    virtual void *GDALTransformerArgs() const = 0;
#endif
 
  private:
 
#ifdef SIP_RUN
    QgsGcpTransformerInterface( const QgsGcpTransformerInterface &other )
#endif
};
 
class ANALYSIS_EXPORT QgsLinearGeorefTransform : public QgsGcpTransformerInterface SIP_SKIP
{
  public:
 
    QgsLinearGeorefTransform() = default;
 
    bool getOriginScale( QgsPointXY &origin, double &scaleX, double &scaleY ) const;
 
    QgsGcpTransformerInterface *clone() const override;
    bool updateParametersFromGcps( const QVector<QgsPointXY> &sourceCoordinates, const QVector<QgsPointXY> &destinationCoordinates, bool invertYAxis = false ) override;
    int minimumGcpCount() const override;
    GDALTransformerFunc GDALTransformer() const override;
    void *GDALTransformerArgs() const override;
    TransformMethod method() const override;
 
  private:
    struct LinearParameters
    {
      QgsPointXY origin;
      double scaleX = 1;
      double scaleY = 1;
      bool invertYAxis = false;
    } mParameters;
 
    static int linearTransform( void *pTransformerArg, int bDstToSrc, int nPointCount,
                                double *x, double *y, double *z, int *panSuccess );
 
};
 
class ANALYSIS_EXPORT QgsHelmertGeorefTransform : public QgsGcpTransformerInterface SIP_SKIP
{
  public:
 
    QgsHelmertGeorefTransform() = default;
 
    bool getOriginScaleRotation( QgsPointXY &origin, double &scale, double &rotation ) const;
 
    QgsGcpTransformerInterface *clone() const override;
    bool updateParametersFromGcps( const QVector<QgsPointXY> &sourceCoordinates, const QVector<QgsPointXY> &destinationCoordinates, bool invertYAxis = false ) override;
    int minimumGcpCount() const override;
    GDALTransformerFunc GDALTransformer() const override;
    void *GDALTransformerArgs() const override;
    TransformMethod method() const override;
 
  private:
 
    struct HelmertParameters
    {
      QgsPointXY origin;
      double scale = 0;
      double angle = 0;
      bool invertYAxis = false;
    };
    HelmertParameters mHelmertParameters;
 
    static int helmertTransform( void *pTransformerArg, int bDstToSrc, int nPointCount,
                                 double *x, double *y, double *z, int *panSuccess );
 
};
 
class ANALYSIS_EXPORT QgsGDALGeorefTransform : public QgsGcpTransformerInterface SIP_SKIP
{
  public:
 
    QgsGDALGeorefTransform( bool useTPS, unsigned int polynomialOrder );
    ~QgsGDALGeorefTransform() override;
 
    QgsGcpTransformerInterface *clone() const override;
    bool updateParametersFromGcps( const QVector<QgsPointXY> &sourceCoordinates, const QVector<QgsPointXY> &destinationCoordinates, bool invertYAxis = false ) override;
    int minimumGcpCount() const override;
    GDALTransformerFunc GDALTransformer() const override;
    void *GDALTransformerArgs() const override;
    TransformMethod method() const override;
 
  private:
    void destroyGdalArgs();
 
    QVector<QgsPointXY> mSourceCoords;
    QVector<QgsPointXY> mDestCoordinates;
    bool mInvertYAxis = false;
 
    const int mPolynomialOrder;
    const bool mIsTPSTransform;
 
    void *mGDALTransformerArgs = nullptr;
 
};
 
class ANALYSIS_EXPORT QgsProjectiveGeorefTransform : public QgsGcpTransformerInterface SIP_SKIP
{
  public:
 
    QgsProjectiveGeorefTransform();
 
    QgsGcpTransformerInterface *clone() const override;
    bool updateParametersFromGcps( const QVector<QgsPointXY> &sourceCoordinates, const QVector<QgsPointXY> &destinationCoordinates, bool invertYAxis = false ) override;
    int minimumGcpCount() const override;
    GDALTransformerFunc GDALTransformer() const override;
    void *GDALTransformerArgs() const override;
    TransformMethod method() const override;
 
  private:
    struct ProjectiveParameters
    {
      double H[9];        // Homography
      double Hinv[9];     // Inverted homography
      bool hasInverse;  // Could the inverted homography be calculated?
    } mParameters;
 
    static int projectiveTransform( void *pTransformerArg, int bDstToSrc, int nPointCount,
                                    double *x, double *y, double *z, int *panSuccess );
 
};
 
#endif //QGSGCPTRANSFORMER_H