otbStereoSensorModelToElevationMapFilter.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 __otbStereoSensorModelToElevationMapFilter_h
#define __otbStereoSensorModelToElevationMapFilter_h

#include "itkImageToImageFilter.h"
#include "itkInterpolateImageFunction.h"
#include "otbGenericRSTransform.h"
#include "itkNormalizedCorrelationImageToImageMetric.h"

namespace otb
{
template <class TInputImage, class TOutputHeight>
class ITK_EXPORT StereoSensorModelToElevationFilter
  : public itk::ImageToImageFilter<TInputImage, TOutputHeight>
{
public:
  typedef StereoSensorModelToElevationFilter                 Self;
  typedef itk::ImageToImageFilter<TInputImage, TOutputHeight> Superclass;
  typedef itk::SmartPointer<Self>                            Pointer;
  typedef itk::SmartPointer<const Self>                      ConstPointer;

  itkNewMacro(Self);

  itkTypeMacro(StereoSensorModelToElevationFilter, ImageToImageFilter);

  typedef TInputImage                                        InputImageType;
  typedef typename InputImageType::PixelType                 InputPixelType;
  typedef TOutputHeight                                      OutputImageType;
  typedef typename OutputImageType::RegionType               OutputRegionType;
  typedef typename OutputImageType::PixelType                HeightType;
  typedef typename OutputImageType::PointType                OutputPointType;

  typedef itk::InterpolateImageFunction<TInputImage, double> InterpolatorType;
  typedef typename InterpolatorType::Pointer                 InterpolatorPointerType;
  typedef GenericRSTransform<>                               GenericRSTransformType;
  typedef GenericRSTransform<double, 3, 3>                   GenericRSTransform3DType;
  typedef itk::NormalizedCorrelationImageToImageMetric
    <TInputImage, TInputImage>                               CorrelationMetricType;
  typedef typename InputImageType::SizeType                  RadiusType;

  void SetMasterInput(const TInputImage * image);

  void SetSlaveInput(const TInputImage * image);

  const TInputImage * GetMasterInput() const;

  const TInputImage * GetSlaveInput() const;

  TOutputHeight * GetCorrelationOutput();

  // Deprecated calls to elevation setters
  otbLegacyElevationMacro();

  itkLegacyMacro(itkSetMacro(UseDEM, bool));

  itkLegacyMacro(itkGetMacro(UseDEM, bool));

  itkLegacyMacro(itkBooleanMacro(UseDEM));

  itkSetMacro(LowerElevation, double);

  itkSetMacro(HigherElevation, double);

  itkSetMacro(ElevationStep, double);

  itkSetMacro(CorrelationThreshold, double);

  itkGetMacro(CorrelationThreshold, double);

  itkSetMacro(VarianceThreshold, double);

  itkGetMacro(VarianceThreshold, double);

  itkSetMacro(SubtractInitialElevation, bool);

  itkGetMacro(SubtractInitialElevation, bool);

  itkBooleanMacro(SubtractInitialElevation);

  void SetRadius(unsigned int radius)
  {
    m_Radius.Fill(radius);
    this->Modified();
  }

protected:
  StereoSensorModelToElevationFilter();

  virtual ~StereoSensorModelToElevationFilter();

  virtual void ThreadedGenerateData(const OutputRegionType& outputRegionForThread,
                                    int threadId);

  virtual void GenerateInputRequestedRegion(void);

  virtual void BeforeThreadedGenerateData();

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

  inline double Correlation(const std::vector<double>& master,
                            const std::vector<double>& slave) const;

  InterpolatorPointerType m_Interpolator;

  RadiusType m_Radius;

  double  m_LowerElevation;

  double  m_HigherElevation;

  double  m_ElevationStep;

  bool m_UseDEM;

  double m_CorrelationThreshold;

  double m_VarianceThreshold;

  bool m_SubtractInitialElevation;
};
} // End namespace otb

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

#endif