otbSatelliteRSR.h

Go to the documentation of this file.

/*=========================================================================

 Program:   ORFEO Toolbox
 Language:  C++
 Date:      $Date$
 Version:   $Revision$


 Copyright (c) Centre National d'Etudes Spatiales. All rights reserved.
 See OTBCopyright.txt for details.


 This software is distributed WITHOUT ANY WARRANTY; without even
 the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
 PURPOSE.  See the above copyright notices for more information.

 =========================================================================*/
#ifndef __otbSatelliteRSR_h
#define __otbSatelliteRSR_h

#include "itkDataObject.h"
#include <itkObjectFactory.h>
#include <vector>
#include <utility>
#include <fstream>
#include <limits>

#include "otbSpectralResponse.h"

namespace otb
{
template<class TPrecision = double, class TValuePrecision = double>
class SatelliteRSR: public itk::DataObject
{
public:
  typedef SatelliteRSR Self;
  typedef itk::DataObject Superclass;
  typedef itk::SmartPointer<Self> Pointer;
  typedef itk::SmartPointer<const Self> ConstPointer;

  itkNewMacro(Self)
; itkTypeMacro(SatelliteRSR, DataObject)
;

  itkGetConstMacro(NbBands, unsigned int)
; itkSetMacro(NbBands, unsigned int)
;

  itkSetMacro(SortBands, bool)
;

  typedef TPrecision PrecisionType;
  typedef TValuePrecision ValuePrecisionType;

  typedef SpectralResponse<TPrecision, TValuePrecision> SpectralResponseType;
  typedef typename SpectralResponseType::Pointer SpectralResponsePointerType;
  typedef std::vector<SpectralResponsePointerType> RSRVectorType;
  typedef typename SpectralResponseType::VectorPairType VectorPairType;

  virtual bool Clear();

  virtual int Size() const;

  void Load(const std::string & filename, ValuePrecisionType coefNormalization = 1.0);

  void Load(PrecisionType lambdaMin, PrecisionType lambdaMax, PrecisionType sampling,
            ValuePrecisionType coefNormalization = 1.0);

  struct sort_band
  {
    bool operator()(SpectralResponsePointerType a, SpectralResponsePointerType b)
    {
      PrecisionType aFirstNotNull;
      PrecisionType bFirstNotNull;

      typename VectorPairType::const_iterator it = a->GetResponse().begin();

      while ((*it).second == 0)
        {
        ++it;
        }
      aFirstNotNull = (*it).first;

      typename VectorPairType::const_iterator it2 = b->GetResponse().begin();

      while ((*it2).second == 0)
        {
        ++it2;
        }
      bFirstNotNull = (*it2).first;

      return aFirstNotNull < bFirstNotNull;
    }
  };

  inline ValuePrecisionType operator()(const PrecisionType & lambda, const unsigned int numBand);

  void PrintSelf(std::ostream& os, itk::Indent indent) const;

  RSRVectorType & GetRSR()
  {
    return m_RSR;
  }

protected:
  SatelliteRSR();
  //SatelliteRSR( const std::string & filename );
  virtual ~SatelliteRSR()
  {
  }
;

  bool m_SortBands;

private:
  SatelliteRSR(const Self&); //purposely not implemented
  void operator=(const Self&); //purposely not implemented

  RSRVectorType m_RSR;
  unsigned int m_NbBands;
};

}// end namespace otb


#ifndef OTB_MANUAL_INSTANTIATION
#include "otbSatelliteRSR.txx"
#endif

#endif