/usr/include/OTB-5.8/otbPointSetExtractROI.txx is in libotb-dev 5.8.0+dfsg-3.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 | /*=========================================================================
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 otbPointSetExtractROI_txx
#define otbPointSetExtractROI_txx
#include "otbPointSetExtractROI.h"
#include "itkMacro.h"
namespace otb
{
/**
*
*/
template <class TInputPointSet, class TOutputPointSet>
PointSetExtractROI<TInputPointSet, TOutputPointSet>
::PointSetExtractROI()
{
}
/**
*
*/
template <class TInputPointSet, class TOutputPointSet>
void
PointSetExtractROI<TInputPointSet, TOutputPointSet>
::PrintSelf(std::ostream& os, itk::Indent indent) const
{
Superclass::PrintSelf(os, indent);
}
/**
* This method causes the filter to generate its output.
*/
template <class TInputPointSet, class TOutputPointSet>
void
PointSetExtractROI<TInputPointSet, TOutputPointSet>
::GenerateData(void)
{
typedef typename TInputPointSet::PointsContainer InputPointsContainer;
typedef typename TInputPointSet::PointsContainerPointer InputPointsContainerPointer;
typedef typename TOutputPointSet::PointsContainerPointer OutputPointsContainerPointer;
typedef typename TInputPointSet::PointDataContainerPointer InputPointDataContainerPointer;
typedef typename TOutputPointSet::PointDataContainerPointer OutputPointDataContainerPointer;
InputPointSetPointer inputPointSet = this->GetInput();
OutputPointSetPointer outputPointSet = this->GetOutput();
if (!inputPointSet)
{
itkExceptionMacro(<< "Missing Input PointSet");
}
if (!outputPointSet)
{
itkExceptionMacro(<< "Missing Output PointSet");
}
outputPointSet->SetBufferedRegion(outputPointSet->GetRequestedRegion()); //TODO update outputRegion
InputPointsContainerPointer inPoints = inputPointSet->GetPoints();
InputPointDataContainerPointer inData = inputPointSet->GetPointData();
OutputPointsContainerPointer outPoints = outputPointSet->GetPoints();
OutputPointDataContainerPointer outData = outputPointSet->GetPointData();
typename InputPointsContainer::ConstIterator inputPoint = inPoints->Begin();
// Commented cause using iterator on the pointSetData crash
// Use a direct access to pointSet Data instead to avoid segfault
//typename InputPointDataContainer::ConstIterator inputData = inData->Begin();
//if (inData.IsNotNull())
// {
// inputData = inData->Begin();
// }
while (inputPoint != inPoints->End())
{
typename InputPointsContainer::Element point = inputPoint.Value();
if ((((point[0] >= m_StartX) && (point[0] < m_StartX + m_SizeX))
|| ((point[0] <= m_StartX) && (point[0] > m_StartX + m_SizeX))) //cover the case when size<0
&& (((point[1] >= m_StartY) && (point[1] < m_StartY + m_SizeY))
|| ((point[1] <= m_StartY) && (point[1] > m_StartY + m_SizeY))))
{
// Add the point
outPoints->push_back(point);
// Get & Add the data
typename InputPointSetType::PixelType data;
inputPointSet->GetPointData(inputPoint.Index(), &data);
outData->push_back(data/*inputData.Value()*/);
}
++inputPoint;
//++inputData;
}
}
} // end namespace otb
#endif
|