otbLogPolarTransform.hxx

Go to the documentation of this file.

/*
 * Copyright (C) 2005-2022 Centre National d'Etudes Spatiales (CNES)
 *
 * This file is part of Orfeo Toolbox
 *
 *     https://www.orfeo-toolbox.org/
 *
 * 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.
 */
 
#ifndef otbLogPolarTransform_hxx
#define otbLogPolarTransform_hxx
 
#include "otbLogPolarTransform.h"
#include "otbMacro.h"
#include "otbMath.h"
 
namespace otb
{
template <class TScalarType>
LogPolarTransform<TScalarType>::LogPolarTransform() : Superclass(4)
{
  m_Center[0] = 0.0;
  m_Center[1] = 0.0;
  m_Scale[0]  = 1.0;
  m_Scale[1]  = 1.0;
}
 
template <class TScalarType>
LogPolarTransform<TScalarType>::~LogPolarTransform()
{
}
 
template <class TScalarType>
void LogPolarTransform<TScalarType>::SetParameters(const ParametersType& parameters)
{
  m_Center[0] = parameters[0];
  m_Center[1] = parameters[1];
  m_Scale[0]  = parameters[2];
  m_Scale[1]  = parameters[3];
  otbMsgDebugMacro(<< "Call To SetParameters: Center=" << m_Center << ", Scale=" << m_Scale);
  this->m_Parameters = parameters;
  this->Modified();
}
 
template <class TScalarType>
typename LogPolarTransform<TScalarType>::ParametersType& LogPolarTransform<TScalarType>::GetParameters(void) const
{
  // Filling parameters vector
  this->m_Parameters[0] = m_Center[0];
  this->m_Parameters[1] = m_Center[1];
  this->m_Parameters[2] = m_Scale[0];
  this->m_Parameters[3] = m_Scale[1];
 
  return this->m_Parameters;
}
template <class TScalarType>
typename LogPolarTransform<TScalarType>::OutputPointType LogPolarTransform<TScalarType>::TransformPoint(const InputPointType& point) const
{
  OutputPointType result;
  double          theta  = point[0] * m_Scale[0] * CONST_PI_180;
  double          logRho = point[1] * m_Scale[1];
  result[0]              = m_Center[0];
  result[1]              = m_Center[1];
  result[0] += std::exp(logRho) * std::cos(theta);
  result[1] += std::exp(logRho) * std::sin(theta);
  return result;
}
 
template <class TScalarType>
typename LogPolarTransform<TScalarType>::OutputVectorType LogPolarTransform<TScalarType>::TransformVector(const InputVectorType& vector) const
{
  OutputVectorType result;
  double           theta  = vector[0] * m_Scale[0] * CONST_PI_180;
  double           logRho = vector[1] * m_Scale[1];
  result[0]               = 0.;
  result[1]               = 0.;
  result[0] += std::exp(logRho) * std::cos(theta);
  result[1] += std::exp(logRho) * std::sin(theta);
 
  return result;
}
template <class TScalarType>
typename LogPolarTransform<TScalarType>::OutputVnlVectorType LogPolarTransform<TScalarType>::TransformVector(const InputVnlVectorType& vector) const
{
  OutputVnlVectorType result;
  double              theta  = vector[0] * m_Scale[0] * CONST_PI_180;
  double              logRho = vector[1] * m_Scale[1];
  result[0]                  = 0.;
  result[1]                  = 0.;
  result[0] += std::exp(logRho) * std::cos(theta);
  result[1] += std::exp(logRho) * std::sin(theta);
 
  return result;
}
template <class TScalarType>
void LogPolarTransform<TScalarType>::PrintSelf(std::ostream& os, itk::Indent indent) const
{
  Superclass::PrintSelf(os, indent);
  os << indent << "Center: " << m_Center << std::endl;
  os << indent << "Scale: " << m_Scale << std::endl;
}
 
} // end namespace otb
#endif