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 <QVector>
#include <QList>

#include "qgsrectangle.h"
#include "qgscoordinatereferencesystem.h"
#include "qgscoordinatetransform.h"
#include "qgsrasterinterface.h"

#include <cmath>

class QgsPoint;

class CORE_EXPORT QgsRasterProjector : public QgsRasterInterface
{
  public:
    enum Precision
    {
      Approximate = 0, 
      Exact = 1,   
    };

    Q_DECL_DEPRECATED
    QgsRasterProjector( const QgsCoordinateReferenceSystem& theSrcCRS,
                        const QgsCoordinateReferenceSystem& theDestCRS,
                        int theSrcDatumTransform,
                        int theDestDatumTransform,
                        const QgsRectangle& theDestExtent,
                        int theDestRows, int theDestCols,
                        double theMaxSrcXRes, double theMaxSrcYRes,
                        const QgsRectangle& theExtent
                      );
    Q_DECL_DEPRECATED
    QgsRasterProjector( const QgsCoordinateReferenceSystem& theSrcCRS,
                        const QgsCoordinateReferenceSystem& theDestCRS,
                        const QgsRectangle& theDestExtent,
                        int theDestRows, int theDestCols,
                        double theMaxSrcXRes, double theMaxSrcYRes,
                        const QgsRectangle& theExtent
                      );
    Q_DECL_DEPRECATED
    QgsRasterProjector( const QgsCoordinateReferenceSystem& theSrcCRS,
                        const QgsCoordinateReferenceSystem& theDestCRS,
                        double theMaxSrcXRes, double theMaxSrcYRes,
                        const QgsRectangle& theExtent
                      );
    QgsRasterProjector();

    // To avoid synthesized which fails on copy of QgsCoordinateTransform
    // (QObject child) in Python bindings
    Q_DECL_DEPRECATED QgsRasterProjector( const QgsRasterProjector &projector );

    ~QgsRasterProjector();

    QgsRasterProjector & operator=( const QgsRasterProjector &projector );

    QgsRasterProjector *clone() const override;

    int bandCount() const override;

    QGis::DataType dataType( int bandNo ) const override;

    void setCRS( const QgsCoordinateReferenceSystem & theSrcCRS, const QgsCoordinateReferenceSystem & theDestCRS,
                 int srcDatumTransform = -1, int destDatumTransform = -1 );

    QgsCoordinateReferenceSystem srcCrs() const { return mSrcCRS; }

    QgsCoordinateReferenceSystem destCrs() const { return mDestCRS; }

    Q_DECL_DEPRECATED void setMaxSrcRes( double theMaxSrcXRes, double theMaxSrcYRes ) { Q_UNUSED( theMaxSrcXRes ); Q_UNUSED( theMaxSrcYRes ); }

    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 ) override;
    QgsRasterBlock *block2( int bandNo, const QgsRectangle & extent, int width, int height, QgsRasterBlockFeedback* feedback = nullptr ) override;

    bool destExtentSize( const QgsRectangle& theSrcExtent, int theSrcXSize, int theSrcYSize,
                         QgsRectangle& theDestExtent, int& theDestXSize, int& theDestYSize );

    static bool extentSize( const QgsCoordinateTransform* ct,
                            const QgsRectangle& theSrcExtent, int theSrcXSize, int theSrcYSize,
                            QgsRectangle& theDestExtent, int& theDestXSize, int& theDestYSize );

  private:

    QgsCoordinateReferenceSystem mSrcCRS;

    QgsCoordinateReferenceSystem mDestCRS;

    int mSrcDatumTransform;

    int mDestDatumTransform;

    Precision mPrecision;

};


class ProjectorData
{
  public:
    ProjectorData( const QgsRectangle &extent, int width, int height, QgsRasterInterface *input, const QgsCoordinateTransform *inverseCt, QgsRasterProjector::Precision precision );
    ~ProjectorData();

    bool srcRowCol( int theDestRow, int theDestCol, int *theSrcRow, int *theSrcCol );

    QgsRectangle srcExtent() const { return mSrcExtent; }
    int srcRows() const { return mSrcRows; }
    int srcCols() const { return mSrcCols; }

  private:
    void destPointOnCPMatrix( int theRow, int theCol, double *theX, double *theY );

    int matrixRow( int theDestRow );
    int matrixCol( int theDestCol );

    inline bool preciseSrcRowCol( int theDestRow, int theDestCol, int *theSrcRow, int *theSrcCol );

    inline bool approximateSrcRowCol( int theDestRow, int theDestCol, int *theSrcRow, int *theSrcCol );

    void insertRows( const QgsCoordinateTransform* ct );

    void insertCols( const QgsCoordinateTransform* ct );

    void calcCP( int theRow, int theCol, const QgsCoordinateTransform* ct );

    bool calcRow( int theRow, const QgsCoordinateTransform* ct );

    bool calcCol( int theCol, const QgsCoordinateTransform* ct );

    void calcSrcExtent();

    void calcSrcRowsCols();

    bool checkCols( const QgsCoordinateTransform* ct );

    bool checkRows( const QgsCoordinateTransform* ct );

    void calcHelper( int theMatrixRow, QgsPoint *thePoints );

    void nextHelper();

    QString cpToString();

    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<QgsPoint> > mCPMatrix;

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

    /* Warning: using QList is slow on access */
    QgsPoint *pHelperTop;

    /* Warning: using QList is slow on access */
    QgsPoint *pHelperBottom;

    int mHelperTopRow;

    int mCPCols;
    int mCPRows;

    double mSqrTolerance;

    double mMaxSrcXRes;
    double mMaxSrcYRes;

};


#endif