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 "qgis_core.h"
#include "qgis_sip.h"
#include <QVector>
#include <QList>
 
#include "qgsrectangle.h"
#include "qgscoordinatereferencesystem.h"
#include "qgscoordinatetransform.h"
#include "qgsrasterinterface.h"
 
#include <cmath>
 
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 );
 
    int matrixRow( int destRow );
 
    int matrixCol( int destCol );
 
    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 );
 
    bool checkRows( const QgsCoordinateTransform &ct );
 
    void calcHelper( int matrixRow, QgsPointXY *points );
 
    void nextHelper();
 
    QString cpToString() const;
 
    bool mApproximate;
 
    QgsCoordinateTransform mInverseCt;
 
    QgsRectangle mDestExtent;
 
    QgsRectangle mSrcExtent;
 
    QgsRectangle mExtent;
 
    int mDestRows;
 
    int mDestCols;
 
    double mDestXRes;
 
    double mDestYRes;
 
    int mSrcRows;
 
    int mSrcCols;
 
    double mSrcXRes;
 
    double mSrcYRes;
 
    double mDestRowsPerMatrixRow;
 
    double mDestColsPerMatrixCol;
 
    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;
 
    int mCPCols;
    int mCPRows;
 
    double mSqrTolerance;
 
    double mMaxSrcXRes;
    double mMaxSrcYRes;
 
};
 
#endif
 
#endif