qgsrasterprojector.h

Go to the documentation of this file.

/***************************************************************************
    qgsrasterprojector.h - Raster projector
     --------------------------------------
    Date                 : Jan 16, 2011
    Copyright            : (C) 2005 by Radim Blazek
    email                : radim dot blazek 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.                                   *
 *                                                                         *
 ***************************************************************************/
 
/* This code takes ideas from WarpBuilder in Geotools.
 * Thank to Ing. Andrea Aime, Ing. Simone Giannecchini and GeoSolutions S.A.S.
 * See : http://geo-solutions.blogspot.com/2011/01/developers-corner-improving.html
 */
 
#ifndef QGSRASTERPROJECTOR_H
#define QGSRASTERPROJECTOR_H
 
#include <cmath>
 
#include "qgis_core.h"
#include "qgis_sip.h"
#include "qgscoordinatereferencesystem.h"
#include "qgscoordinatetransform.h"
#include "qgsrasterinterface.h"
#include "qgsrectangle.h"
 
#include <QList>
#include <QVector>
 
class QgsPointXY;

class CORE_EXPORT QgsRasterProjector : public QgsRasterInterface
{
    Q_GADGET
 
  public:
    enum Precision
    {
      Approximate = 0, 
      Exact = 1,       
    };
    Q_ENUM( Precision )
 
    QgsRasterProjector();
 
    QgsRasterProjector *clone() const override SIP_FACTORY;
 
    int bandCount() const override;
 
    Qgis::DataType dataType( int bandNo ) const override;

    Q_DECL_DEPRECATED void setCrs( const QgsCoordinateReferenceSystem &srcCRS, const QgsCoordinateReferenceSystem &destCRS, int srcDatumTransform = -1, int destDatumTransform = -1 ) SIP_DEPRECATED;

    void setCrs( const QgsCoordinateReferenceSystem &srcCRS, const QgsCoordinateReferenceSystem &destCRS, QgsCoordinateTransformContext transformContext );

    QgsCoordinateReferenceSystem sourceCrs() const { return mSrcCRS; }

    QgsCoordinateReferenceSystem destinationCrs() const { return mDestCRS; }
 
    Precision precision() const { return mPrecision; }
    void setPrecision( Precision precision ) { mPrecision = precision; }
    // Translated precision mode, for use in ComboBox etc.
    static QString precisionLabel( Precision precision );
 
    QgsRasterBlock *block( int bandNo, const QgsRectangle &extent, int width, int height, QgsRasterBlockFeedback *feedback = nullptr ) override SIP_FACTORY;

    bool destExtentSize( const QgsRectangle &srcExtent, int srcXSize, int srcYSize, QgsRectangle &destExtent SIP_OUT, int &destXSize SIP_OUT, int &destYSize SIP_OUT );

    static bool extentSize( const QgsCoordinateTransform &ct, const QgsRectangle &srcExtent, int srcXSize, int srcYSize, QgsRectangle &destExtent SIP_OUT, int &destXSize SIP_OUT, int &destYSize SIP_OUT );
 
  private:
    QgsCoordinateReferenceSystem mSrcCRS;

    QgsCoordinateReferenceSystem mDestCRS;

    Q_DECL_DEPRECATED int mSrcDatumTransform = -1;

    Q_DECL_DEPRECATED int mDestDatumTransform = -1;

    Precision mPrecision = Approximate;
 
    QgsCoordinateTransformContext mTransformContext;
};
 
 
#ifndef SIP_RUN

class ProjectorData
{
  public:
    ProjectorData(
      const QgsRectangle &extent, int width, int height, QgsRasterInterface *input, const QgsCoordinateTransform &inverseCt, QgsRasterProjector::Precision precision, QgsRasterBlockFeedback *feedback = nullptr
    );
    ~ProjectorData();
 
    ProjectorData( const ProjectorData &other ) = delete;
    ProjectorData &operator=( const ProjectorData &other ) = delete;

    bool srcRowCol( int destRow, int destCol, int *srcRow, int *srcCol );
 
    QgsRectangle srcExtent() const { return mSrcExtent; }
    int srcRows() const { return mSrcRows; }
    int srcCols() const { return mSrcCols; }
 
  private:
    void destPointOnCPMatrix( int row, int col, double *theX, double *theY ) const;

    int matrixRow( int destRow ) const;

    int matrixCol( int destCol ) const;

    inline bool preciseSrcRowCol( int destRow, int destCol, int *srcRow, int *srcCol );

    inline bool approximateSrcRowCol( int destRow, int destCol, int *srcRow, int *srcCol );

    void insertRows( const QgsCoordinateTransform &ct );

    void insertCols( const QgsCoordinateTransform &ct );

    void calcCP( int row, int col, const QgsCoordinateTransform &ct );

    bool calcRow( int row, const QgsCoordinateTransform &ct );

    bool calcCol( int col, const QgsCoordinateTransform &ct );

    void calcSrcExtent();

    void calcSrcRowsCols();

    bool checkCols( const QgsCoordinateTransform &ct ) const;

    bool checkRows( const QgsCoordinateTransform &ct ) const;

    void calcHelper( int matrixRow, QgsPointXY *points );

    void nextHelper();

    QString cpToString() const;

    bool mApproximate = false;

    QgsCoordinateTransform mInverseCt;

    QgsRectangle mDestExtent;

    QgsRectangle mSrcExtent;

    QgsRectangle mExtent;

    int mDestRows;

    int mDestCols;

    double mDestXRes = 0.0;

    double mDestYRes = 0.0;

    int mSrcRows = 0;

    int mSrcCols = 0;

    double mSrcXRes = 0.0;

    double mSrcYRes = 0.0;

    double mDestRowsPerMatrixRow = 0.0;

    double mDestColsPerMatrixCol = 0.0;

    QList< QList<QgsPointXY> > mCPMatrix;

    /* Same size as mCPMatrix */
    QList< QList<bool> > mCPLegalMatrix;

    /* Warning: using QList is slow on access */
    QgsPointXY *pHelperTop = nullptr;

    /* Warning: using QList is slow on access */
    QgsPointXY *pHelperBottom = nullptr;

    int mHelperTopRow = 0;

    int mCPCols = 0;
    int mCPRows = 0;

    double mSqrTolerance = 0.0;

    double mMaxSrcXRes = 0;
    double mMaxSrcYRes = 0;
};

#endif
 
#endif