From bd85849d11fcd4ef01da1966e9c86abde5485a75 Mon Sep 17 00:00:00 2001
From: Patrick Imbo <patrick.imbo@c-s.fr>
Date: Mon, 20 Mar 2006 14:57:06 +0000
Subject: [PATCH] HuImage invariants
---
.../FeatureExtraction/otbFlusserImage.cxx | 87 +++++++++++++++++++
1 file changed, 87 insertions(+)
create mode 100644 Testing/Code/FeatureExtraction/otbFlusserImage.cxx
diff --git a/Testing/Code/FeatureExtraction/otbFlusserImage.cxx b/Testing/Code/FeatureExtraction/otbFlusserImage.cxx
new file mode 100644
index 0000000000..1bc1b73874
--- /dev/null
+++ b/Testing/Code/FeatureExtraction/otbFlusserImage.cxx
@@ -0,0 +1,87 @@
+/*=========================================================================
+
+ Programme : OTB (ORFEO ToolBox)
+ Auteurs : CS - P.Imbo
+ Language : C++
+ Date : 20 mars 2006
+ Version :
+ Role :
+ $Id$
+
+=========================================================================*/
+#if defined(_MSC_VER)
+#pragma warning ( disable : 4786 )
+#endif
+
+#include "itkExceptionObject.h"
+#include "itkImage.h"
+
+#include "otbImageFileReader.h"
+#include "otbFlusserImageFunction.h"
+
+int otbFlusserImage( int argc, char ** argv )
+{
+ try
+ {
+ const char * inputFilename = argv[1];
+ unsigned int Number = 1;
+
+ typedef unsigned char InputPixelType;
+ const unsigned int Dimension = 2;
+
+ typedef itk::Image< InputPixelType, Dimension > InputImageType;
+ typedef otb::ImageFileReader< InputImageType > ReaderType;
+ typedef std::complex<float> ComplexType;
+ typedef float RealType;
+ typedef otb::FlusserImageFunction<InputImageType,float,float> FunctionType;
+
+ InputImageType::RegionType region;
+ InputImageType::SizeType size;
+ InputImageType::IndexType start;
+
+ start.Fill( 0 );
+ size[0] = 50;
+ size[1] = 50;
+
+ ReaderType::Pointer reader = ReaderType::New();
+ FunctionType::Pointer function =FunctionType::New();
+
+ reader->SetFileName( inputFilename );
+
+ InputImageType::Pointer image = reader->GetOutput();
+
+ region.SetIndex( start );
+ region.SetSize( size );
+
+ image->SetRegions(region);
+ image->Update();
+ function->SetInputImage( image );
+
+ InputImageType::IndexType index;
+ index[0]=10;
+ index[1]=10;
+
+ RealType Result;
+
+ for (Number = 1 ;Number<12;Number++)
+ {
+ function->SetNumber(Number);
+ Result = function->EvaluateAtIndex( index );
+ std::cout << "Flusser("<<Number<<") = "<< Result <<std::endl;
+ }
+
+ }
+ catch( itk::ExceptionObject & err )
+ {
+ std::cout << "Exception itk::ExceptionObject levee !" << std::endl;
+ std::cout << err << std::endl;
+ return EXIT_FAILURE;
+ }
+ catch( ... )
+ {
+ std::cout << "Exception levee inconnue !" << std::endl;
+ return EXIT_FAILURE;
+ }
+ return EXIT_SUCCESS;
+}
+
--
GitLab