otbSinclairToReciprocalCircularCovarianceMatrixFunctor.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 __otbSinclairToReciprocalCircularCovarianceMatrixFunctor_h
#define __otbSinclairToReciprocalCircularCovarianceMatrixFunctor_h

#include "vcl_complex.h"
#include "otbSinclairToReciprocalCovarianceMatrixFunctor.h"

namespace otb
{
namespace Functor
{
template <class TInput1, class TInput2, class TInput3, class TOutput>
class SinclairToReciprocalCircularCovarianceMatrixFunctor
{
public:
  typedef double                                   RealType;
  typedef std::complex <RealType>                  ComplexType;
  typedef typename TOutput::ValueType              OutputValueType;
  typedef SinclairToReciprocalCovarianceMatrixFunctor<ComplexType, ComplexType, ComplexType, TOutput> SinclairToReciprocalCovarianceFunctorType;

  inline TOutput operator ()(const TInput1& Shh, const TInput2& Shv, const TInput3& Svv)
  {
    TOutput result;

    result.SetSize(m_NumberOfComponentsPerPixel);


    const ComplexType S_hh = static_cast<ComplexType>(Shh);
    const ComplexType S_hv = static_cast<ComplexType>(Shv);
    const ComplexType S_vv = static_cast<ComplexType>(Svv);

    const ComplexType coef(0.5);

    const ComplexType j2S_hv = S_hv * ComplexType(0.0, 2.0);

    const ComplexType Sll = coef * ( -S_hh-j2S_hv+S_vv );
    const ComplexType Slr = coef * ( -S_hh+-S_vv );
    const ComplexType Srr = coef * ( -S_hh+j2S_hv+S_vv );

    //const ComplexType conjSll = vcl_conj(Sll);
    //const ComplexType conjSlr = vcl_conj(Slr);
    //const ComplexType conjSrr = vcl_conj(Srr);

    SinclairToReciprocalCovarianceFunctorType funct;
    return ( funct(Sll, Slr, Srr ) );
  }

  unsigned int GetNumberOfComponentsPerPixel()
  {
    return m_NumberOfComponentsPerPixel;
  }

  SinclairToReciprocalCircularCovarianceMatrixFunctor() : m_NumberOfComponentsPerPixel(6) {}

  virtual ~SinclairToReciprocalCircularCovarianceMatrixFunctor() {}

protected:


private:
    unsigned int m_NumberOfComponentsPerPixel;
};

} // namespace Functor
} // namespace otb

#endif