otbVcaImageFilter.h

Go to the documentation of this file.

/*
 * Copyright (C) 1999-2011 Insight Software Consortium
 * 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 otbVCAImageFilter_h
#define otbVCAImageFilter_h
 
#include "otbImage.h"
#include "itkAbsImageFilter.h"
#include "otbDotProductImageFilter.h"
#include "otbProjectiveProjectionImageFilter.h"
#include "otbVectorImageToMatrixImageFilter.h"
#include "otbStreamingMinMaxImageFilter.h"
#include "otbStreamingStatisticsImageFilter.h"
#include "otbPCAImageFilter.h"
#include "otbVectorImageToAmplitudeImageFilter.h"
#include "otbConcatenateScalarValueImageFilter.h"
 
#include "itkMersenneTwisterRandomVariateGenerator.h"
#include "vnl/algo/vnl_svd.h"
 
namespace otb
{
 
template <class TVectorImage>
class ITK_EXPORT VCAImageFilter : public itk::ImageToImageFilter<TVectorImage, TVectorImage>
{
public:
  typedef VCAImageFilter Self;
  typedef itk::ImageToImageFilter<TVectorImage, TVectorImage> Superclass;
  typedef itk::SmartPointer<Self>       Pointer;
  typedef itk::SmartPointer<const Self> ConstPointer;
 
  typedef TVectorImage                                VectorImageType;
  typedef typename VectorImageType::IndexType         IndexType;
  typedef typename VectorImageType::SizeType          SizeType;
  typedef typename VectorImageType::RegionType        RegionType;
  typedef typename VectorImageType::PixelType         PixelType;
  typedef typename VectorImageType::InternalPixelType InternalPixelType;
  typedef InternalPixelType                           PrecisionType;
 
  typedef otb::Image<InternalPixelType> ImageType;
 
  typedef itk::AbsImageFilter<ImageType, ImageType> AbsImageFilterType;
  typedef otb::ProjectiveProjectionImageFilter<VectorImageType, VectorImageType, PrecisionType> ProjectiveProjectionImageFilterType;
  typedef otb::DotProductImageFilter<VectorImageType, ImageType>   DotProductImageFilterType;
  typedef otb::MatrixImageFilter<VectorImageType, VectorImageType> MatrixImageFilterType;
  typedef otb::VectorImageToMatrixImageFilter<VectorImageType> VectorImageToMatrixImageFilterType;
  typedef otb::StreamingMinMaxImageFilter<ImageType>           StreamingMinMaxImageFilterType;
  typedef otb::StreamingStatisticsVectorImageFilter<VectorImageType, PrecisionType> StreamingStatisticsVectorImageFilterType;
  typedef otb::StreamingStatisticsImageFilter<ImageType>         StreamingStatisticsImageFilterType;
  typedef itk::Statistics::MersenneTwisterRandomVariateGenerator RandomVariateGeneratorType;
  typedef otb::PCAImageFilter<VectorImageType, VectorImageType, otb::Transform::FORWARD> ForwardPCAImageFilterType;
  typedef otb::PCAImageFilter<VectorImageType, VectorImageType, otb::Transform::INVERSE> InversePCAImageFilterType;
  typedef otb::VectorImageToAmplitudeImageFilter<VectorImageType, ImageType>       VectorImageToAmplitudeImageFilterType;
  typedef otb::ConcatenateScalarValueImageFilter<VectorImageType, VectorImageType> ConcatenateScalarValueImageFilterType;
 
  // creation of SmartPointer
  itkNewMacro(Self);
  // runtime type information
  itkTypeMacro(VCAImageFilter, itk::ImageToImageFilter);
 
  itkGetMacro(NumberOfEndmembers, unsigned int);
  itkSetMacro(NumberOfEndmembers, unsigned int);
 
  void Update() override
  {
    this->GenerateData();
  }
 
  void EnlargeOutputRequestedRegion(itk::DataObject* itkNotUsed(output)) override
  {
    this->GetOutput()->SetRequestedRegion(this->GetOutput()->GetLargestPossibleRegion());
  }
 
protected:
  VCAImageFilter();
 
  ~VCAImageFilter() override;
 
  void PrintSelf(std::ostream& os, itk::Indent indent) const override;
 
  void GenerateOutputInformation() override;
 
  void GenerateData() override;
 
private:
  VCAImageFilter(const Self&) = delete;
  void operator=(const Self&) = delete;
 
  unsigned int m_NumberOfEndmembers;
};
 
} // end namesapce otb
 
#include "otbVcaImageFilter.hxx"
 
#endif /* otbVCAImageFilter_h */