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 "qgis_analysis.h"
#include "qgis_sip.h"
#include "qgspointxy.h"
 
#include <QObject>

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