FlusserMomentsImageFunctionExample.cxxΒΆ
Example usage:
./FlusserMomentsImageFunctionExample Input/ROISpot5.png 2
Example source code (FlusserMomentsImageFunctionExample.cxx):
#include "itkMacro.h"
#include "otbImage.h"
#include "otbImageFileReader.h"
// This example illustrates the use of the \doxygen{otb}{FlusserMomentsImageFunction}.
//
// The first step required to use this filter is to include its header file.
#include "otbFlusserMomentsImageFunction.h"
int main(int argc, char* argv[])
{
if (argc != 3)
{
std::cerr << "Usage: " << argv[0] << " inputImageFile ";
std::cerr << " neighborhood_radius" << std::endl;
return EXIT_FAILURE;
}
const char* inputFilename = argv[1];
const unsigned int radius = atoi(argv[2]);
using InputPixelType = unsigned char;
const unsigned int Dimension = 2;
using InputImageType = otb::Image<InputPixelType, Dimension>;
using ReaderType = otb::ImageFileReader<InputImageType>;
ReaderType::Pointer reader = ReaderType::New();
reader->SetFileName(inputFilename);
// The \doxygen{otb}{FlusserMomentsImageFunction} is templated over the
// input image type and the output (real) type value, so we start by
// defining:
using FlusserType = otb::FlusserMomentsImageFunction<InputImageType>;
using MomentType = FlusserType::OutputType;
FlusserType::Pointer fmFunction = FlusserType::New();
// We can choose the region and the pixel of the image which will
// used as coordinate origin
// for the moment computation
InputImageType::RegionType region;
InputImageType::SizeType size;
InputImageType::IndexType start;
start[0] = 0;
start[1] = 0;
size[0] = 50;
size[1] = 50;
reader->Update();
InputImageType::Pointer image = reader->GetOutput();
region.SetIndex(start);
region.SetSize(size);
image->SetRegions(region);
image->Update();
InputImageType::IndexType center;
center[0] = start[0] + size[0] / 2;
center[1] = start[1] + size[1] / 2;
// Next, we plug the input image into the complex moment function
// and we set its parameters.
fmFunction->SetInputImage(image);
fmFunction->SetNeighborhoodRadius(radius);
// In order to get the value of the moment, we call the
// \code{EvaluateAtIndex} method.
MomentType Result = fmFunction->EvaluateAtIndex(center);
for (unsigned int j = 0; j < 11; ++j)
{
std::cout << "The moment of order " << j + 1 << " is equal to " << Result[j] << std::endl;
}
// \relatedClasses
// \begin{itemize}
// \item \doxygen{otb}{FlusserPathFunction}
// \end{itemize}
return EXIT_SUCCESS;
}