dispersion_polynomial.h

Go to the documentation of this file.

#ifndef DISPERSION_POLYNOMIAL_H
#define DISPERSION_POLYNOMIAL_H
#include "dispersion.h"
#include "unit.h"
#include "sub_state_vector_array.h"
#include "polynomial_eval.h"

namespace FullPhysics {
/****************************************************************/
class DispersionPolynomial: public SubStateVectorArray<Dispersion> {
public:
  DispersionPolynomial(const blitz::Array<double, 1>& Coeff, 
                       const blitz::Array<bool, 1>& Used_flag,
                       const Unit& Coeff_unit,
                       const std::string& Band_name,
                       int Number_pixel,
                       bool Is_one_based);
  DispersionPolynomial(const blitz::Array<double, 1>& Coeff, 
                       const blitz::Array<bool, 1>& Used_flag,
                       const std::string& Coeff_unit_name,
                       const std::string& Band_name,
                       int Number_pixel, 
                       bool Is_one_based);
  virtual ~DispersionPolynomial() {}

//-----------------------------------------------------------------------
//-----------------------------------------------------------------------

  double dispersion_offset() const { return coeff.value()(0); }

//-----------------------------------------------------------------------
//-----------------------------------------------------------------------

  double dispersion_spacing() const { return coeff.value()(1); }

//-----------------------------------------------------------------------
//-----------------------------------------------------------------------

  double dispersion_offset_uncertainty() const 
  { 
    if(sv_cov_sub.rows() < 1)
      return 0;
    double t = sv_cov_sub(0,0);
    return (t < 0 ? 0 : sqrt(t)); 
  }

//-----------------------------------------------------------------------
//-----------------------------------------------------------------------

  double dispersion_spacing_uncertainty() const 
  { 
    if(sv_cov_sub.rows() < 2)
      return 0;
    double t = sv_cov_sub(1,1);
    return (t < 0 ? 0 : sqrt(t)); 
  }

  virtual boost::shared_ptr<Dispersion> clone() const;

  virtual std::string sub_state_identifier() const { return "dispersion/" + band_name_; } 

  virtual std::string state_vector_name_i(int i) const;
  virtual SpectralDomain pixel_grid() const;
  virtual void print(std::ostream& Os) const;
private:
  void initialize();
  bool is_one_based;
  Unit coeff_unit;
  std::string band_name_;
  blitz::Array<double, 1> index_array;
  // Very similar to index_array, but always 1 based.
  blitz::Array<int, 1> spectral_index;
};
}
#endif