otbSinclairToReciprocalCovarianceMatrixFunctor.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 __otbSinclairToReciprocalCovarianceMatrixFunctor_h
#define __otbSinclairToReciprocalCovarianceMatrixFunctor_h

#include "vcl_complex.h"

namespace otb
{
namespace Functor
{
template <class TInput1, class TInput2, class TInput3, class TOutput>
class SinclairToReciprocalCovarianceMatrixFunctor
{
public:
  typedef typename std::complex <double>           ComplexType;
  typedef typename TOutput::ValueType              OutputValueType;
  inline TOutput operator ()(const TInput1& Shh, const TInput2& Shv, const TInput3& Svv)
  {
    TOutput result;

    result.SetSize(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);

    result[0] = static_cast<OutputValueType>( std::norm( S_hh ) );
    result[1] = static_cast<OutputValueType>( S_hh*vcl_conj(S_hv) );
    result[2] = static_cast<OutputValueType>( S_hh*vcl_conj(S_vv) );
    result[3] = static_cast<OutputValueType>( std::norm( S_hv ) );
    result[4] = static_cast<OutputValueType>( S_hv*vcl_conj(S_vv) );
    result[5] = static_cast<OutputValueType>( std::norm( S_vv ) );

    return (result);
  }

  unsigned int GetNumberOfComponentsPerPixel()
  {
    return NumberOfComponentsPerPixel;
  }

  SinclairToReciprocalCovarianceMatrixFunctor() {}

  virtual ~SinclairToReciprocalCovarianceMatrixFunctor() {}

protected:


private:
    itkStaticConstMacro(NumberOfComponentsPerPixel, unsigned int, 6 );

};

} // namespace Functor
} // namespace otb

#endif