src/core/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 you 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 <cmath>

//class QgsRectangle;
class QgsPoint;

//class CORE_EXPORT QgsRasterProjector
class QgsRasterProjector
{
//    Q_OBJECT
  public:
    QgsRasterProjector(
      QgsCoordinateReferenceSystem theSrcCRS,
      QgsCoordinateReferenceSystem theDestCRS,
      QgsRectangle theDestExtent,
      int theDestRows, int theDestCols,
      double theMaxSrcXRes, double theMaxSrcYRes,
      QgsRectangle theExtent
    );

    ~QgsRasterProjector();


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

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

    QgsPoint srcPoint( int theRow, int theCol );

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

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

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

    void insertRows();

    void insertCols();

    /* calculate single control point in current matrix */
    void calcCP( int theRow, int theCol );

    bool calcRow( int theRow );

    bool calcCol( int theCol );

    void calcSrcExtent();

    void calcSrcRowsCols();

    bool checkCols();

    bool checkRows();

    //void calcHelper ( int theMatrixRow, QList<QgsPoint> *thePoints );
    void calcHelper( int theMatrixRow, QgsPoint *thePoints );

    void nextHelper();

    QgsRectangle srcExtent() { return mSrcExtent; }

    int srcRows() { return mSrcRows; }
    int srcCols() { return mSrcCols; }
    void setSrcRows( int theRows ) { mSrcRows = theRows; mSrcXRes = mSrcExtent.height() / mSrcRows; }
    void setSrcCols( int theCols ) { mSrcCols = theCols; mSrcYRes = mSrcExtent.width() / mSrcCols; }

    QString cpToString();

  private:
    QgsCoordinateReferenceSystem mSrcCRS;

    QgsCoordinateReferenceSystem mDestCRS;

    QgsCoordinateTransform mCoordinateTransform;

    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;

    /* 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;

    bool mApproximate;
};

#endif