qgsrasterblock.h

Go to the documentation of this file.

/***************************************************************************
    qgsrasterblock.h - Class representing a block of raster data
     --------------------------------------
    Date                 : Oct 9, 2012
    Copyright            : (C) 2012 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.                                   *
 *                                                                         *
 ***************************************************************************/

#ifndef QGSRASTERBLOCK_H
#define QGSRASTERBLOCK_H

#include "qgis_core.h"
#include "qgis_sip.h"
#include <limits>
#include <QImage>
#include "qgis.h"
#include "qgserror.h"
#include "qgslogger.h"
#include "qgsrasterrange.h"

class QgsRectangle;

class CORE_EXPORT QgsRasterBlock
{
  public:
    QgsRasterBlock();

    QgsRasterBlock( Qgis::DataType dataType, int width, int height );

    virtual ~QgsRasterBlock();

    bool reset( Qgis::DataType dataType, int width, int height );

    // TODO: consider if use isValid() at all, isEmpty() should be sufficient
    // and works also if block is valid but empty - difference between valid and empty?

    bool isValid() const { return mValid; }

    void setValid( bool valid ) { mValid = valid; }

    bool isEmpty() const;

    // Return data type size in bytes
    static int typeSize( int dataType )
    {
      // Modified and extended copy from GDAL
      switch ( dataType )
      {
        case Qgis::Byte:
          return 1;

        case Qgis::UInt16:
        case Qgis::Int16:
          return 2;

        case Qgis::UInt32:
        case Qgis::Int32:
        case Qgis::Float32:
        case Qgis::CInt16:
          return 4;

        case Qgis::Float64:
        case Qgis::CInt32:
        case Qgis::CFloat32:
          return 8;

        case Qgis::CFloat64:
          return 16;

        case Qgis::ARGB32:
        case Qgis::ARGB32_Premultiplied:
          return 4;

        default:
          return 0;
      }
    }

    // Data type in bytes
    int dataTypeSize() const
    {
      return typeSize( mDataType );
    }

    static bool typeIsNumeric( Qgis::DataType type );

    static bool typeIsColor( Qgis::DataType type );

    Qgis::DataType dataType() const { return mDataType; }

    static Qgis::DataType typeWithNoDataValue( Qgis::DataType dataType, double *noDataValue );

    bool hasNoDataValue() const { return mHasNoDataValue; }

    bool hasNoData() const;

    void setNoDataValue( double noDataValue );

    void resetNoDataValue();

    double noDataValue() const { return mNoDataValue; }

    static QByteArray valueBytes( Qgis::DataType dataType, double value );

    double value( int row, int column ) const;

    double value( qgssize index ) const;

    QRgb color( int row, int column ) const;

    QRgb color( qgssize index ) const;

    bool isNoData( int row, int column );

    bool isNoData( qgssize index );

    bool setValue( int row, int column, double value );

    bool setValue( qgssize index, double value );

    bool setColor( int row, int column, QRgb color );

    bool setColor( qgssize index, QRgb color );

    bool setIsNoData( int row, int column );

    bool setIsNoData( qgssize index );

    bool setIsNoData();

    bool setIsNoDataExcept( QRect exceptRect );

    void setIsData( int row, int column );

    void setIsData( qgssize index );

    QByteArray data() const;

    void setData( const QByteArray &data, int offset = 0 );

    char *bits( int row, int column ) SIP_SKIP;

    char *bits( qgssize index ) SIP_SKIP;

    char *bits() SIP_SKIP;

    static QString printValue( double value );

    static QString printValue( float value ) SIP_SKIP;

    bool convert( Qgis::DataType destDataType );

    QImage image() const;

    bool setImage( const QImage *image );

    inline static double readValue( void *data, Qgis::DataType type, qgssize index ) SIP_SKIP;

    inline static void writeValue( void *data, Qgis::DataType type, qgssize index, double value ) SIP_SKIP;

    void applyNoDataValues( const QgsRasterRangeList &rangeList );

    void applyScaleOffset( double scale, double offset );

    QgsError error() const { return mError; }

    void setError( const QgsError &error ) { mError = error;}

    QString toString() const;

    static QRect subRect( const QgsRectangle &extent, int width, int height, const QgsRectangle &subExtent );

    int width() const { return mWidth; }

    int height() const { return mHeight; }

  private:
    static QImage::Format imageFormat( Qgis::DataType dataType );
    static Qgis::DataType dataType( QImage::Format format );

    static bool isNoDataValue( double value, double noDataValue );

    bool isNoDataValue( double value ) const;

    bool createNoDataBitmap();

    static void *convert( void *srcData, Qgis::DataType srcDataType, Qgis::DataType destDataType, qgssize size );

    // Valid
    bool mValid = true;

    // Data type
    Qgis::DataType mDataType = Qgis::UnknownDataType;

    // Data type size in bytes, to make bits() fast
    int mTypeSize = 0;

    // Width
    int mWidth = 0;

    // Height
    int mHeight = 0;

    // Has no data value
    bool mHasNoDataValue = false;

    // No data value
    double mNoDataValue;

    static const QRgb NO_DATA_COLOR;

    // Data block for numerical data types, not used with image data types
    // QByteArray does not seem to be intended for large data blocks, does it?
    void *mData = nullptr;

    // Image for image data types, not used with numerical data types
    QImage *mImage = nullptr;

    // Bitmap of no data. One bit for each pixel. Bit is 1 if a pixels is no data.
    // Each row is represented by whole number of bytes (last bits may be unused)
    // to make processing rows easy.
    char *mNoDataBitmap = nullptr;

    // number of bytes in mNoDataBitmap row
    int mNoDataBitmapWidth = 0;

    // total size in bytes of mNoDataBitmap
    qgssize mNoDataBitmapSize = 0;

    // Error
    QgsError mError;
};

inline double QgsRasterBlock::readValue( void *data, Qgis::DataType type, qgssize index ) SIP_SKIP
{
  if ( !data )
  {
    return std::numeric_limits<double>::quiet_NaN();
  }

  switch ( type )
  {
    case Qgis::Byte:
      return static_cast< double >( ( static_cast< quint8 * >( data ) )[index] );
      break;
    case Qgis::UInt16:
      return static_cast< double >( ( static_cast< quint16 * >( data ) )[index] );
      break;
    case Qgis::Int16:
      return static_cast< double >( ( static_cast< qint16 * >( data ) )[index] );
      break;
    case Qgis::UInt32:
      return static_cast< double >( ( static_cast< quint32 * >( data ) )[index] );
      break;
    case Qgis::Int32:
      return static_cast< double >( ( static_cast< qint32 * >( data ) )[index] );
      break;
    case Qgis::Float32:
      return static_cast< double >( ( static_cast< float * >( data ) )[index] );
      break;
    case Qgis::Float64:
      return static_cast< double >( ( static_cast< double * >( data ) )[index] );
      break;
    default:
      QgsDebugMsg( QString( "Data type %1 is not supported" ).arg( type ) );
      break;
  }

  return std::numeric_limits<double>::quiet_NaN();
}

inline void QgsRasterBlock::writeValue( void *data, Qgis::DataType type, qgssize index, double value ) SIP_SKIP
{
  if ( !data ) return;

  switch ( type )
  {
    case Qgis::Byte:
      ( static_cast< quint8 * >( data ) )[index] = static_cast< quint8 >( value );
      break;
    case Qgis::UInt16:
      ( static_cast< quint16 * >( data ) )[index] = static_cast< quint16 >( value );
      break;
    case Qgis::Int16:
      ( static_cast< qint16 * >( data ) )[index] = static_cast< qint16 >( value );
      break;
    case Qgis::UInt32:
      ( static_cast< quint32 * >( data ) )[index] = static_cast< quint32 >( value );
      break;
    case Qgis::Int32:
      ( static_cast< qint32 * >( data ) )[index] = static_cast< qint32 >( value );
      break;
    case Qgis::Float32:
      ( static_cast< float * >( data ) )[index] = static_cast< float >( value );
      break;
    case Qgis::Float64:
      ( static_cast< double * >( data ) )[index] = value;
      break;
    default:
      QgsDebugMsg( QString( "Data type %1 is not supported" ).arg( type ) );
      break;
  }
}

inline double QgsRasterBlock::value( qgssize index ) const SIP_SKIP
{
  if ( !mData )
  {
    QgsDebugMsg( "Data block not allocated" );
    return std::numeric_limits<double>::quiet_NaN();
  }
  return readValue( mData, mDataType, index );
}

inline bool QgsRasterBlock::isNoDataValue( double value ) const SIP_SKIP
{
  return std::isnan( value ) || qgsDoubleNear( value, mNoDataValue );
}

#endif