/******************************** LICENSE ********************************

 Copyright 2007 European Centre for Medium-Range Weather Forecasts (ECMWF)

 Licensed under the Apache License, Version 2.0 (the "License");
 you may not use this file except in compliance with the License.
 You may obtain a copy of the License at 

    http://www.apache.org/licenses/LICENSE-2.0

 Unless required by applicable law or agreed to in writing, software
 distributed under the License is distributed on an "AS IS" BASIS,
 WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
 See the License for the specific language governing permissions and
 limitations under the License.

 ******************************** LICENSE ********************************/

/*! \file RasterData.h
    \brief Definition of the Template class RasterData.
    
    Magics Team - ECMWF 2005
    
    Started: Tue 12-Apr-2005
    
    Changes:
    
*/

#ifndef RasterData_H
#define RasterData_H

#include <algorithm>
#include "magics.h"

class TeProjection;


namespace magics {

class Raster : public vector<double>
{
public:
	Raster() {}
	virtual ~Raster() {}
	double operator()(int row, int column) const
		{ return (*this)[row * columns_ + column]; }
	double min() const { if(this->empty()) return -1.; return *(std::min_element(begin(), end())); }
	double max() const { if(this->empty()) return -1.; return *(std::max_element(begin(), end())); }
        inline int   getRows()             const { return rows_; }
	inline int   getColumns()          const { return columns_; }	
	inline void setColumns(int columns)  { columns_ = columns; }
	inline void setRows(int rows)        { rows_ = rows; }

protected:
	 virtual void print(ostream& out) const { 
	 		out << "Raster[";
	 		out << size() << " points, ";
	 		out << min() << ", " << max() << "]"; 
	 }
	 int        rows_;
	 int        columns_;

	
private:
    //! Copy constructor - No copy allowed
	Raster(const Raster&);
    //! Overloaded << operator to copy - No copy allowed
	Raster& operator=(const Raster&);
    
// -- Friends
    //! Overloaded << operator to call print().
	friend ostream& operator<<(ostream& s,const Raster& p)
		{ p.print(s); return s; }

	
};


class RasterData  : public Raster {
	

public:
	RasterData() {}
	virtual ~RasterData() {}
	
	TeProjection& getProjection() { return *projection_; }
	
	inline double getXResolution()      const { return x_; }
	inline double getYResolution()      const { return y_; }

	inline const UserPoint& getLowerLeftCorner()  const { return lowerLeft_; }
	inline const UserPoint& getUpperRightCorner() const { return upperRight_; }
	
	inline void setXResolution(double x ) { x_ = x; }
	inline void setYResolution(double y ) { y_ = y; }
	
	inline void setUpperRightCorner(double lon, double lat) 
		{ upperRight_ = UserPoint(lon, lat); }
	inline void setLowerLeftCorner(double lon, double lat) 
		{ lowerLeft_ = UserPoint(lon, lat); }
	
	inline void setProjection(TeProjection* projection) { projection_ = projection; }
	

protected:
     //! Method to print string about this class on to a stream of type ostream (virtual).
	 virtual void print(ostream&) const {}
	 TeProjection* projection_;
	 double      x_;
	 double      y_;
	 UserPoint             lowerLeft_;
	 UserPoint             upperRight_;
	 	 
private:
    //! Copy constructor - No copy allowed
	RasterData(const RasterData&);
    //! Overloaded << operator to copy - No copy allowed
	RasterData& operator=(const RasterData&);
    
// -- Friends
    //! Overloaded << operator to call print().
	friend ostream& operator<<(ostream& s,const RasterData& p)
		{ p.print(s); return s; }

};

} // namespace magics


#endif
