otbVectorDataTransformFilter.hxx

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 otbVectorDataTransformFilter_hxx
#define otbVectorDataTransformFilter_hxx
 
#include "otbVectorDataTransformFilter.h"
#include "itkProgressReporter.h"
#include <itkContinuousIndex.h>
#include "otbStopwatch.h"
#include "vcl_legacy_aliases.h"
 
namespace otb
{
 
template <class TInputVectorData, class TOutputVectorData>
VectorDataTransformFilter<TInputVectorData, TOutputVectorData>::VectorDataTransformFilter()
{
  //  m_Transform = GenericTransformType::New();
  // with ITK v4 you can't instantiate a transform (virtual). We do NOT want to
  // use the otb::Transform or we loose the capability of using all the existing
  // itk transform, so we just keep it as NULL and you have to be careful not
  // to burn yourself when using it.
  m_Transform = nullptr;
}
 
template <class TInputVectorData, class TOutputVectorData>
typename VectorDataTransformFilter<TInputVectorData, TOutputVectorData>::PointType
VectorDataTransformFilter<TInputVectorData, TOutputVectorData>::ProcessPoint(PointType pointCoord) const
{
  itk::Point<double, 2> point;
  point = m_Transform->TransformPoint(pointCoord);
  return point;
}
 
 
template <class TInputVectorData, class TOutputVectorData>
typename VectorDataTransformFilter<TInputVectorData, TOutputVectorData>::LinePointerType
VectorDataTransformFilter<TInputVectorData, TOutputVectorData>::ProcessLine(LinePointerType line) const
{
  typedef typename LineType::VertexListType::ConstPointer VertexListConstPointerType;
  typedef typename LineType::VertexListConstIteratorType  VertexListConstIteratorType;
  VertexListConstPointerType                              vertexList = line->GetVertexList();
  VertexListConstIteratorType                             it         = vertexList->Begin();
  typename LineType::Pointer                              newLine    = LineType::New();
  while (it != vertexList->End())
  {
    itk::Point<double, 2>           point;
    itk::ContinuousIndex<double, 2> index;
    typename LineType::VertexType pointCoord = it.Value();
    point                                    = m_Transform->TransformPoint(pointCoord);
    index[0]                                 = point[0];
    index[1]                                 = point[1];
    if (!vnl_math_isnan(index[0]) && !vnl_math_isnan(index[1]))
      newLine->AddVertex(index);
    ++it;
  }
 
  return newLine;
}
 
template <class TInputVectorData, class TOutputVectorData>
typename VectorDataTransformFilter<TInputVectorData, TOutputVectorData>::PolygonPointerType
VectorDataTransformFilter<TInputVectorData, TOutputVectorData>::ProcessPolygon(PolygonPointerType polygon) const
{
  typedef typename PolygonType::VertexListType::ConstPointer VertexListConstPointerType;
  typedef typename PolygonType::VertexListConstIteratorType  VertexListConstIteratorType;
  VertexListConstPointerType                                 vertexList = polygon->GetVertexList();
  VertexListConstIteratorType                                it         = vertexList->Begin();
  typename PolygonType::Pointer                              newPolygon = PolygonType::New();
  while (it != vertexList->End())
  {
    itk::Point<double, 2>           point;
    itk::ContinuousIndex<double, 2> index;
    typename PolygonType::VertexType pointCoord = it.Value();
    point                                       = m_Transform->TransformPoint(pointCoord);
    index[0]                                    = point[0];
    index[1]                                    = point[1];
    if (!vnl_math_isnan(index[0]) && !vnl_math_isnan(index[1]))
      newPolygon->AddVertex(index);
    ++it;
  }
  return newPolygon;
}
 
template <class TInputVectorData, class TOutputVectorData>
typename VectorDataTransformFilter<TInputVectorData, TOutputVectorData>::PolygonListPointerType
VectorDataTransformFilter<TInputVectorData, TOutputVectorData>::ProcessPolygonList(PolygonListPointerType polygonList) const
{
 
  PolygonListPointerType newPolygonList = PolygonListType::New();
  for (typename PolygonListType::ConstIterator it = polygonList->Begin(); it != polygonList->End(); ++it)
  {
    newPolygonList->PushBack(this->ProcessPolygon(it.Get()));
  }
  return newPolygonList;
}
 
 
template <class TInputVectorData, class TOutputVectorData>
void VectorDataTransformFilter<TInputVectorData, TOutputVectorData>::GenerateData(void)
{
  Superclass::GenerateOutputInformation();
  this->AllocateOutputs();
 
  InputVectorDataPointer  inputPtr  = this->GetInput();
  OutputVectorDataPointer outputPtr = this->GetOutput();
 
  outputPtr->SetProjectionRef(inputPtr->GetProjectionRef());
 
  // // Create the output tree root
  // OutputDataNodePointerType newDataNode = OutputDataNodeType::New();
  // newDataNode->SetNodeType(inputPtr->GetRoot()->GetNodeType());
  // newDataNode->SetNodeId(inputPtr->GetRoot()->GetNodeId());
  // outputPtr->SetRoot(newDataNode);
 
  // Start recursive processing
  otb::Stopwatch chrono = otb::Stopwatch::StartNew();
  this->ProcessNode(inputPtr,inputPtr->GetRoot(),outputPtr,outputPtr->GetRoot());
  chrono.Stop();
  otbMsgDevMacro(<< "VectorDataTransformFilter: features processed in " << chrono.GetElapsedMilliseconds() << " ms.");
}
 
} // end namespace otb
 
#endif