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;

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

    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;

    int mSrcDatumTransform = -1;

    int mDestDatumTransform = -1;

    Precision mPrecision = Approximate;

};


#ifndef SIP_RUN

class ProjectorData
{
  public:
    ProjectorData( const QgsRectangle &extent, int width, int height, QgsRasterInterface *input, const QgsCoordinateTransform &inverseCt, QgsRasterProjector::Precision precision );
    ~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();

    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