otbGCPsToRPCSensorModelImageFilter.h

Go to the documentation of this file.

/*
 * Copyright (C) 2005-2024 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 otbGCPsToRPCSensorModelImageFilter_h
#define otbGCPsToRPCSensorModelImageFilter_h
 
#include "itkUnaryFunctorImageFilter.h"
#include "itkCastImageFilter.h"
#include "otbDEMHandler.h"
 
namespace otb
{
 
template <class TImage>
class ITK_EXPORT GCPsToRPCSensorModelImageFilter : public itk::CastImageFilter<TImage, TImage>
{
public:
  typedef GCPsToRPCSensorModelImageFilter Self;
  typedef itk::CastImageFilter<TImage, TImage> Superclass;
  typedef itk::SmartPointer<Self>       Pointer;
  typedef itk::SmartPointer<const Self> ConstPointer;
 
  typedef itk::Point<double, 2>               Point2DType;
  typedef itk::Point<double, 3>               Point3DType;
  typedef std::pair<Point2DType, Point3DType> GCPType;
  typedef std::vector<GCPType> GCPsContainerType;
  typedef std::vector<double>  ErrorsContainerType;
 
  typedef itk::ContinuousIndex<> ContinuousIndexType;
  typedef itk::ContinuousIndex<double, 3> Continuous3DIndexType;
 
  typedef otb::DEMHandler                  DEMHandlerType;
 
  itkNewMacro(Self);
 
  itkTypeMacro(GCPsToRPCSensorModelImageFilter, CastImageFilter);
 
  itkStaticConstMacro(ImageDimension, unsigned int, TImage::ImageDimension);
 
  //  itkSetMacro(UseImageGCPs, bool);
  // itkGetMacro(UseImageGCPs, bool);
  //  itkBooleanMacro(UseImageGCPs);
  void SetUseImageGCPs(bool use);
  void LoadImageGCPs();
 
  itkSetMacro(UseDEM, bool);
  itkGetMacro(UseDEM, bool);
  itkBooleanMacro(UseDEM);
 
  itkSetMacro(MeanElevation, double);
  itkGetConstReferenceMacro(MeanElevation, double);
 
  itkGetConstReferenceMacro(RMSGroundError, double);
 
  ErrorsContainerType& GetErrorsContainer();
 
  itkGetConstReferenceMacro(MeanError, double);
 
  GCPsContainerType& GetGCPsContainer();
 
  void SetGCPsContainer(const GCPsContainerType& container);
 
  void AddGCP(const Point2DType& sensorPoint, const Point3DType& groundPoint);
 
  void AddGCP(const Point2DType& sensorPoint, const Point2DType& groundPoint);
 
  void RemoveGCP(unsigned int id);
 
  void ClearGCPs();
 
protected:
  GCPsToRPCSensorModelImageFilter();
 
  ~GCPsToRPCSensorModelImageFilter() override;
 
  void PrintSelf(std::ostream& os, itk::Indent indent) const override;
 
  void GenerateOutputInformation() override;
 
  void Modified() const override;
 
private:
  GCPsToRPCSensorModelImageFilter(const Self&) = delete;
  void operator=(const Self&) = delete;
 
  void ComputeErrors();
 
  bool m_UseImageGCPs;
 
  double m_RMSGroundError;
 
  ErrorsContainerType m_ErrorsContainer;
 
  double m_MeanError;
 
  bool m_UseDEM;
 
  double m_MeanElevation;
 
  GCPsContainerType m_GCPsContainer;
 
  ImageMetadata m_ImageMetadata;
 
  mutable bool m_ModelUpToDate;
 
}; // end of class
 
} // end of namespace otb
 
#ifndef OTB_MANUAL_INSTANTIATION
#include "otbGCPsToRPCSensorModelImageFilter.hxx"
#endif
 
#endif