/****************************************************************************

 Copyright (C) 2010-2012 the Office National des Forets (ONF), France
                     and the Association de Recherche Technologie et Sciences (ARTS), Ecole Nationale Suprieure d'Arts et Metiers (ENSAM), Cluny, France.
                     All rights reserved.

 Contact : alexandre.piboule@onf.fr

 Developers : Alexandre PIBOULE (ONF)

 This file is part of PluginShared library 2.0.

 PluginShared is free library: you can redistribute it and/or modify
 it under the terms of the GNU Lesser General Public License as published by
 the Free Software Foundation, either version 3 of the License, or
 (at your option) any later version.

 PluginShared is distributed in the hope that it will be useful,
 but WITHOUT ANY WARRANTY; without even the implied warranty of
 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 GNU General Public License for more details.

 You should have received a copy of the GNU Lesser General Public License
 along with PluginShared.  If not, see <http://www.gnu.org/licenses/lgpl.html>.

*****************************************************************************/

#ifndef LVOX3_GRID3D_H
#define LVOX3_GRID3D_H

#include "tools/3dgrid/abstract/lvox3_abstractgrid3d.h"

template< typename DataT > class LVOX3_StandardGrid3DDrawManager;

/*!
 * \class CT_Grid3D
 * \ingroup PluginShared_Items
 * \brief <b>ItemDrawable managing a 3D grid</b>
 *
 * It represents grid in 3D, with templated values in cells.
 *
 */
template< typename DataT>
class LVOX3_Grid3D : public LVOX3_AbstractGrid3D
{
    CT_TYPE_TEMPLATED_IMPL_MACRO(LVOX3_Grid3D, DataT, LVOX3_AbstractGrid3D, 3D grid LVOX)
public:

    /*!
     *  \brief How to manage central cell when using neighboursValues method
     *
     */
    enum CenterMode
    {
        CM_KeepCenter = 1,      /*!< Gives neighbours and central value */
        CM_NAasCenter = 2,      /*!< Gives neighbours and NA for central value */
        CM_DropCenter = 3       /*!< Gives neighbours only */
    };

    /**
      * \brief Empty Contructor vide
      */
    LVOX3_Grid3D();

    /*!
     * \brief Contructor with integer column and row coordinates
     *
     * \param xmin Minimum X coordinate (bottom left corner)
     * \param ymin Minimum Y coordinate (bottom left corner)
     * \param zmin Minimum Z coordinate (bottom left corner)
     * \param dimx Number of colums
     * \param dimy Number of rows
     * \param dimz Number of zlevels
     * \param resolution X Size of a
     * \param resolution Y Size of a cell
     * \param resolution Z Size of a cell
     * \param na Value used to code NA
     * \param initValue Initialisation value for grid cells
     */

    LVOX3_Grid3D(
              double xmin,
              double ymin,
              double zmin,
              size_t dimx,
              size_t dimy,
              size_t dimz,
              double resolutionx,
              double resolutiony,
              double resolutionz,
              DataT na,
              DataT initValue);

    /*!
     * \brief Constructor with min and max (X,Y) coordinates
     *
     * \param xmin Minimum X coordinate (bottom left corner)
     * \param ymin Minimum Y coordinate (bottom left corner)
     * \param zmin Minimum Z coordinate (bottom left corner)
     * \param xmax Maximum X coordinate (upper right corner)
     * \param ymax Maximum Y coordinate (upper right corner)
     * \param zmax Maximum Z coordinate (upper right corner)
     * \param resolution X Size of a
     * \param resolution Y Size of a cell
     * \param resolution Z Size of a cell
     * \param na Value used to code NA
     * \param initValue Initialisation value for grid cells
     * \param coordConstructor Not used, only to ensure constructor different signatures
     */

    LVOX3_Grid3D(
              double xmin,
              double ymin,
              double zmin,
              double xmax,
              double ymax,
              double zmax,
              double resolutionx,
              double resolutiony,
              double resolutionz,
              DataT na,
              DataT initValue);

    /*!
     * \brief Factory with min and max (X,Y,Z) coordinates
     *
     * \param xmin Minimum X coordinate (bottom left corner)
     * \param ymin Minimum Y coordinate (bottom left corner)
     * \param zmin Minimum Z coordinate (bottom left corner)
     * \param xmax Maximum X coordinate (upper right corner)
     * \param ymax Maximum Y coordinate (upper right corner)
     * \param zmax Maximum Z coordinate (upper right corner)
     * \param resolution X Size of a
     * \param resolution Y Size of a cell
     * \param resolution Z Size of a cell
     * \param na Value used to code NA
     * \param initValue Initialisation value for grid cells
     * \param coordConstructor Not used, only to ensure constructor different signatures
     */

    static LVOX3_Grid3D<DataT>* createGrid3DFromXYZCoords(
                                                          double xmin,
                                                          double ymin,
                                                          double zmin,
                                                          double xmax,
                                                          double ymax,
                                                          double zmax,
                                                          double resolutionx,
                                                          double resolutiony,
                                                          double resolutionz,
                                                          DataT na,
                                                          DataT initValue,
                                                          bool extends = true);
    /*!
     * \brief Destructor
     */
    virtual ~LVOX3_Grid3D() override;

    /*!
     * \brief Initialize all grid cells values with val
     *
     * Used by Constructors.
     *
     * \param val
     */
    void initGridWithValue(const DataT val);

    /*!
     * \brief Generic setter taking a double as input value
     * \param index Index of cell to modify
     * \param value Double value to cast in effective type
     */
    virtual void setValueAtIndexFromDouble(const size_t &index, const double &value) override;

    /*!
     * \brief Compute min and max values
     */
    void computeMinMax() override;

    /**
      * \brief Gives the NA value
      * \return Valeur NA value
      */
    inline DataT NA() const {return _NAdata;}

    /*!
     * \brief Gives the max value of the grid
     * \return Max value
     */
    inline DataT dataMax() const { return _dataMax; }

    /*!
     * \brief Gives the min value of the grid
     * \return Min value
     */
    inline DataT dataMin() const { return _dataMin; }


    // CT_ITemplatedData3DArray
    size_t xArraySize() const { return xdim(); }

    // CT_ITemplatedData3DArray
    size_t yArraySize() const { return ydim(); }

    // CT_ITemplatedData3DArray
    size_t zArraySize() const { return zdim(); }


    /*!
     * \brief Set value at specified index
     * \param index Index
     * \param value Value
     * \return True if the value has actually been set
     */
    bool setValueAtIndex(const size_t index, const DataT value);


    /**
      * \brief Set the value at row liny, column colx and z level levz to value
      * \param colx column, first one is 0
      * \param liny row, first one is 0
      * \param levz z level, first one is 0
      * \param value Value
      * \return True if the value has actually been set
      */
    bool setValue(const size_t colx, const size_t liny, const size_t levz, DataT value);

    /**
      * \brief Set the value for the cell containing (x,y,z) coordinate to value
      * \param x X coordinate
      * \param y Y coordinate
      * \param z Z coordinate
      * \param value Value
      * \return True if the value has actually been set
      */
    bool setValueAtXYZ(const double x, const double y, const double z, const DataT value);


    /*!
     * \brief ives the value at specified index
     * \param index Index
     * \return Value
     */
    DataT valueAtIndex(const size_t index) const;

    /*!
     * \brief Generic [0;1]  (or -1 for NA) value accessor for use as CT_AbstractGrid3D
     * \param index index in the grid
     * \return A double value between 0 (min value) and 1 (max value), or -1 for NA
     */
    virtual double ratioValueAtIndex(const size_t index) const override;


    /*!
     * \brief Return a double value for any type (or NAN for NA)
     * \param index index in the grid
     * \return A double value
     */
    virtual double valueAtIndexAsDouble(const size_t index) const override;


    /*!
     * \brief Return a QString value for any type (or -1 for NA)
     * \param index index in the grid
     * \return A QString représenting value
     */
    virtual QString valueAtIndexAsString(const size_t index) const override;

    /*!
     * \brief return na value as a string
     *
     */
    virtual QString NAAsString() const override;

    virtual bool isNA(size_t index) const override;


    /**
      * \brief Gives the value at row liny and column colx and z level levz
      * \param colx column, first one is 0
      * \param liny row, first one is 0
      * \param levz z level, first one is 0
      * \return Value at row liny and column clox
      */
    DataT value(const size_t colx, const size_t liny, const size_t levz) const;

    // CT_ITemplatedData3DArray
    //DataT dataFromArray(const size_t &x, const size_t &y, const size_t &z) const;

    // CT_ITemplatedData3DArray
    //DataT dataFromArray(const size_t &index) const;

    /**
      * \brief Gives the value at (x,y,z) coordinate
      * \param x X coordinate
      * \param y Y coordinate
      * \param z Z coordinate
      * \return Value at (x,y,z)
      */
    DataT valueAtXYZ(const double x, const double y, const double z) const;


    /**
      * \brief Set the value for the cell at specified index, if value is superior to actual one
      * \param index Index
      * \param value Value
      * \return True if the value has actually been set
      */
    bool setMaxValueAtIndex(const size_t index, const DataT value);

    /**
      * \brief Set the value for the cell containing (x,y,z) coordinate to value, if value is superior to actual one
      * \param x X coordinate
      * \param y Y coordinate
      * \param z Z coordinate
      * \param value Value
      * \return True if the value has actually been set
      */
    bool setMaxValueAtXYZ(const double x, const double y, const double z, const DataT value);

    /**
      * \brief Set the value for the cell at specified index, if value is inferior to actual one
      * \param index Index
      * \param value Value
      * \return True if the value has actually been set
      */
    bool setMinValueAtIndex(const size_t index, const DataT value);

    /**
      * \brief Set the value for the cell containing (x,y,z) coordinate to value, if value is inferior to actual one
      * \param x X coordinate
      * \param y Y coordinate
      * \param z Z coordinate
      * \param value Value
      * \return True if the value has actually been set
      */
    bool setMinValueAtXYZ(const double x, const double y, const double z, const DataT value);

    /**
      * \brief Add value to the cell at the specified index
      * \param index Index
      * \param value Value
      * \return True if the value has actually been set
      */
    bool addValueAtIndex(const size_t index, const DataT value);

    /**
      * \brief Add value to the cell containing (x,y,z) coordinate
      * \param x X coordinate
      * \param y Y coordinate
      * \param z Z coordinate
      * \param value Value
      * \return True if the value has actually been updated
      */
    bool addValueAtXYZ(const double x, const double y, const double z, const DataT value);

    /**
      * \brief Gives neighbours values
      * \param colx column, first one is 0
      * \param liny row, first one is 0
      * \param levz z level, first one is 0
      * \param distance : Distance for neighbours search in cells
      * \param keepNAs : Should the NA values be kept ?
      * \param centermode : How should be treated the central cell ?
      * \return List of neighbours values
      */
    QList<DataT> neighboursValues(const size_t colx, const size_t liny, const size_t levz, const size_t distance, const bool keepNAs = false, const CenterMode centermode = LVOX3_Grid3D::CM_DropCenter) const;

    //Méthodes pour la sommation des distances
    /*!
     * \brief setSommaDelta Set Delta shot summation value
     * \param index index to assign to
     * \param value Value to assign
     */
    void setSommaDelta(size_t index, float value);

    /*!
     * \brief setSommaFree Set Free shot summation value
     * \param index index to assign to
     * \param value Value to assign
     */
    void setSommaFree(size_t index,float value);

    /*!
     * \brief setSommaEffectiveFree
     * \param index index to assign to
     * \param value Value to assign
     */
    void setSommaEffectiveFree(size_t index, float value);

    /*!
     * \brief setSommaEffectiveDelta
     * \param index index to assign to
     * \param value Value to assign
     */
    void setSommaEffectiveDelta(size_t index,float value);

    /*!
     * \brief setSommaEffectiveHits
     * \param index index to assign to
     * \param value Value to assign
     */
    void setSommaEffectiveHits(size_t index,float value);

    /*!
     * \brief setSommaDeltaSquare
     * \param index index to assign to
     * \param value Value to assign
     */
    void setSommaDeltaSquare(size_t index, float value);


    /*!
     * \brief addSommaDelta
     * \param index index to assign to
     * \param value Value to add to the sum of the delta distances
     */
    void addSommaDelta(size_t index, float value);

    /*!
     * \brief addSommaFree
     * \param index index to assign to
     * \param value Value to add to the sum of the free distances
     */
    void addSommaFree(size_t index,float value);

    /*!
     * \brief addSommaEffectiveFree
     * \param index index to assign to
     * \param value Value to add to the sum of the effective free distances
     */
    void addSommaEffectiveFree(size_t index, float value);

    /*!
     * \brief addSommaEffectiveDelta
     * \param index index to assign to
     * \param value Value to add to the sum of the effective delta distances
     */
    void addSommaEffectiveDelta(size_t index,float value);

    /*!
     * \brief addSommaEffectiveHits
     * \param index index to assign to
     * \param value Value to add to the sum of the effective hits distances
     */
    void addSommaEffectiveHits(size_t index, float value);

    /*!
     * \brief addSommaDeltaSquare
     * \param index index to assign to
     * \param value Value to add to the sum of the delta square distances
     */
    void addSommaDeltaSquare(size_t index,float value);


    /*!
     * \brief getSommaDelta
     * \param index index to get form
     * \return Sum of the delta distances
     */
    float getSommaDelta(size_t index)const;

    /*!
     * \brief getSommaFree
     * \param index index to get form
     * \return Sum of the free distances
     */
    float getSommaFree(size_t index)const;

    /*!
     * \brief getSommaEffectiveFree
     * \param index index to get form
     * \return Sum of the effective free distances
     */
    float getSommaEffectiveFree(size_t index)const;

    /*!
     * \brief getSommaEffectiveDelta
     * \param index index to get form
     * \return Sum of the effective delta distances
     */
    float getSommaEffectiveDelta(size_t index)const;

    /*!
     * \brief getSommaEffectiveHits
     * \param index index to get form
     * \return Sum of all the hit effective distances
     */
    float getSommaEffectiveHits(size_t index)const;

    /*!
     * \brief getSommaDeltaSquare
     * \param index index to get form
     * \return Sum of all the delta square distances
     */
    float getSommaDeltaSquare(size_t index)const;

    /*!
     * \brief setLambda1 Set the attenuation coefficient divided by the volume of a voxel
     * \param vegetation Vegetation element to divide by the volume
     */
    void setLambda1(float vegetation);

    /*!
     * \brief getLambda1 Get the attenuation coefficient
     * \return The attenuation coefficient
     */
    float getLambda1()const;

    /*!
     * \brief doSommation Call the other summation functions
     * \param index index to assign to
     * \param delta Distance a ray would parcour inside the voxel
     * \param free Real distance the ray parcour inside the voxel
     */
    void doSommation(size_t index, float delta, float free);

    /** Methods for colors management
     */
    inline bool colorsDefined() const {return !_colorMap.isEmpty();}

    inline void clearColors() {_colorMap.clear();}

    inline void setDefaultColor(const QColor &color) {_defaultColor = color;}

    inline void addColorForValue(DataT value, const QColor &color) {_colorMap.insert(value, color);}

    inline QColor getColorForValue(DataT value) const {return _colorMap.value(value, _defaultColor);}

    CT_ITEM_COPY_IMP(LVOX3_Grid3D<DataT>)


    inline size_t xdim() const {return LVOX3_AbstractGrid3D::xdim();}
    inline size_t ydim() const {return LVOX3_AbstractGrid3D::ydim();}
    inline size_t zdim() const {return LVOX3_AbstractGrid3D::zdim();}
    inline double minX() const {return LVOX3_AbstractGrid3D::minX();}
    inline double minY() const {return LVOX3_AbstractGrid3D::minY();}
    inline double minZ() const {return LVOX3_AbstractGrid3D::minZ();}
    inline double xresolution() const {return LVOX3_AbstractGrid3D::xresolution();}
    inline double yresolution() const {return LVOX3_AbstractGrid3D::yresolution();}
    inline double zresolution() const {return LVOX3_AbstractGrid3D::zresolution();}

    std::vector<float> _sommaDelta;
    std::vector<float> _sommaFree;
    //std::vector<double> _sommaI;
    std::vector<float> _sommaEffectiveFree;
    std::vector<float> _sommaEffectiveDelta;
    std::vector<float> _sommaEffectiveHits;
    std::vector<float> _sommaDeltaSquare;

    double _lambda1;

    DataT       _NAdata;            /*!< Valeur codant NA */

    DataT       _dataMax;           /*!< valeur maximale du grid*/
    DataT       _dataMin;           /*!< valeur minimale du grid*/
    std::vector<DataT> _data;       /*!< Tableau contenant les donnees pour chaque case de la grille*/

    QMap<DataT, QColor> _colorMap;
    QColor              _defaultColor;

    CT_DEFAULT_IA_BEGIN(LVOX3_Grid3D<DataT>)    
    CT_DEFAULT_IA_V2(LVOX3_Grid3D<DataT>, CT_AbstractCategory::staticInitDataXDimension(), &LVOX3_Grid3D<DataT>::xdim, QObject::tr("X dimension"))
    CT_DEFAULT_IA_V2(LVOX3_Grid3D<DataT>, CT_AbstractCategory::staticInitDataYDimension(), &LVOX3_Grid3D<DataT>::ydim, QObject::tr("Y dimension"))
    CT_DEFAULT_IA_V2(LVOX3_Grid3D<DataT>, CT_AbstractCategory::staticInitDataZDimension(), &LVOX3_Grid3D<DataT>::zdim, QObject::tr("Z dimension"))
    CT_DEFAULT_IA_V2(LVOX3_Grid3D<DataT>, CT_AbstractCategory::staticInitDataX(), &LVOX3_Grid3D<DataT>::minX, QObject::tr("X min"))
    CT_DEFAULT_IA_V2(LVOX3_Grid3D<DataT>, CT_AbstractCategory::staticInitDataY(), &LVOX3_Grid3D<DataT>::minY, QObject::tr("Y min"))
    CT_DEFAULT_IA_V2(LVOX3_Grid3D<DataT>, CT_AbstractCategory::staticInitDataZ(), &LVOX3_Grid3D<DataT>::minZ, QObject::tr("Z min"))
    CT_DEFAULT_IA_V2(LVOX3_Grid3D<DataT>, CT_AbstractCategory::staticInitDataXResolution(), &LVOX3_Grid3D<DataT>::xresolution, QObject::tr("X Resolution"))
    CT_DEFAULT_IA_V2(LVOX3_Grid3D<DataT>, CT_AbstractCategory::staticInitDataYResolution(), &LVOX3_Grid3D<DataT>::yresolution, QObject::tr("Y Resolution"))
    CT_DEFAULT_IA_V2(LVOX3_Grid3D<DataT>, CT_AbstractCategory::staticInitDataZResolution(), &LVOX3_Grid3D<DataT>::zresolution, QObject::tr("Z Resolution"))
    CT_DEFAULT_IA_V2(LVOX3_Grid3D<DataT>, CT_AbstractCategory::staticInitDataNa(), &LVOX3_Grid3D<DataT>::NA, QObject::tr("NA"))
    CT_DEFAULT_IA_V2(LVOX3_Grid3D<DataT>, CT_AbstractCategory::staticInitDataValue(), &LVOX3_Grid3D<DataT>::dataMin, QObject::tr("Min Value"))
    CT_DEFAULT_IA_V2(LVOX3_Grid3D<DataT>, CT_AbstractCategory::staticInitDataValue(), &LVOX3_Grid3D<DataT>::dataMax, QObject::tr("Max Value"))
    CT_DEFAULT_IA_END(LVOX3_Grid3D<DataT>)

    const static LVOX3_StandardGrid3DDrawManager<DataT> ABSGRID3D_DRAW_MANAGER;

};

// Spécialisations

template<>
inline bool LVOX3_Grid3D<bool>::NA() const {return false;}

template<>
void LVOX3_Grid3D<bool>::computeMinMax();

template<>
double LVOX3_Grid3D<bool>::ratioValueAtIndex(const size_t index) const;

template<>
double LVOX3_Grid3D<bool>::valueAtIndexAsDouble(const size_t index) const;

template<>
QString LVOX3_Grid3D<bool>::valueAtIndexAsString(const size_t index) const;

template<>
bool LVOX3_Grid3D<bool>::setMaxValueAtIndex(const size_t index, const bool value);

template<>
bool LVOX3_Grid3D<bool>::setMinValueAtIndex(const size_t index, const bool value);

template<>
bool LVOX3_Grid3D<bool>::addValueAtIndex(const size_t index, const bool value);

template<>
QList<bool> LVOX3_Grid3D<bool>::neighboursValues(const size_t colx, const size_t liny, const size_t levz, const size_t distance, const bool keepNAs, const CenterMode centermode) const;

template<>
QString LVOX3_Grid3D<unsigned long>::valueAtIndexAsString(const size_t index) const;

template<>
QString LVOX3_Grid3D<unsigned long>::NAAsString() const;

template<>
QString LVOX3_Grid3D<long int>::valueAtIndexAsString(const size_t index) const;

template<>
QString LVOX3_Grid3D<long int>::NAAsString() const;

// fin des spécialisations

// Includes the template implementations
#include "lvox3_grid3d.hpp"

#endif // LVOX3_GRID3D_H