Mise en place des classes otb::Image et otb::VectorImage (lecture des Meta...

Mise en place des classes otb::Image et otb::VectorImage (lecture des Meta données des images GDAL).

Testing/Code/IO/otbVectorImageTest.cxx

+/*=========================================================================
+
+  Programme :   OTB (ORFEO ToolBox)
+  Auteurs   :   CS - C.Ruffel
+  Language  :   C++
+  Date      :   09 mai 2006
+  Version   :   
+  Role      :   Test d'instanciation de la classe otb::VectorImage
+  		Lecture de meta donnes
+  $Id$
+
+=========================================================================*/
+#if defined(_MSC_VER)
+#pragma warning ( disable : 4786 )
+#endif
+
+
+
+#include "itkExceptionObject.h"
+#include <iostream>
+#include <fstream>
+#include <string>
+
+#include "otbVectorImage.h"
+#include "otbImageFileReader.h"
+
+
+int otbVectorImageTest(int argc, char* argv[])
+{
+  try
+  {
+        // Verify the number of parameters in the command line
+        const char * inputFilename  = argv[1];
+        const char * outputAsciiFilename  = argv[2];
+      //  unsigned int  GCPnum((unsigned int)::atoi(argv[2]));
+
+        typedef unsigned char  	                                InputPixelType;
+        const   unsigned int        	                        Dimension = 2;
+
+	std::ofstream file;
+	
+	file.open(outputAsciiFilename);
+
+        typedef otb::VectorImage< InputPixelType,  Dimension >        InputImageType;
+        
+        InputImageType::Pointer image = InputImageType::New();
+        
+        typedef otb::ImageFileReader< InputImageType >         ReaderType;
+        
+
+        ReaderType::Pointer reader = ReaderType::New();
+ 
+        reader->SetFileName( inputFilename  );
+	reader->Update();
+
+	image = reader->GetOutput();
+	
+	file << "------ IMAGE --------" << std::endl;
+	std::cout << "------ IMAGE --------" << std::endl;
+// Image n'est pas ecrit dans le fichier car des pointeurs sont affichs, donc la valeur dpend de l'instanciation 
+// -> pour pouvoir faire de tests de non regression sur la sortie ecrite dans le fichier ASCII
+	std::cout << image << std::endl;	
+	std::cout << "---------------------" << std::endl;
+	
+	file << "Spacing " << image->GetSpacing() << std::endl;
+	file << "Origin " << image->GetOrigin() << std::endl;
+	
+	file << "Projection REF " << image->GetProjectionRef() << std::endl;
+		 
+	file << "GCP Projection " << image->GetGCPProjection() << std::endl;
+	
+	unsigned int GCPCount = image->GetGCPCount();	
+	file << "GCP Count " << image->GetGCPCount() << std::endl; 
+	
+	for(unsigned int GCPnum = 0 ; GCPnum < GCPCount ; GCPnum++ )
+	{
+		file << "GCP[" << GCPnum << "] Id " << image->GetGCPId(GCPnum) << std::endl;
+		file << "GCP[" << GCPnum << "] Info " << image->GetGCPInfo(GCPnum) << std::endl;
+		file << "GCP[" << GCPnum << "] Row " << image->GetGCPRow(GCPnum) << std::endl;
+		file << "GCP[" << GCPnum << "] Col " << image->GetGCPCol(GCPnum) << std::endl;
+		file << "GCP[" << GCPnum << "] X " << image->GetGCPX(GCPnum) << std::endl;
+		file << "GCP[" << GCPnum << "] Y " << image->GetGCPY(GCPnum) << std::endl;
+		file << "GCP[" << GCPnum << "] Z " << image->GetGCPZ(GCPnum) << std::endl;
+		file << "----------------" << std::endl;
+	}
+	
+	InputImageType::VectorType tab = image->GetGeoTransform();
+
+	file << "Geo Transform " << std::endl;
+	for(int i = 0 ; i < tab.size() ; i++ )
+	{
+		file << " " <<i<<" -> "<<tab[i]<< std::endl;
+	}
+	tab.clear();
+	
+	tab = image->GetUpperLeftCorner();
+	file << "Corners " << std::endl;
+	for(int i = 0 ; i < tab.size() ; i++ )
+	{
+		file << " UL[" <<i<<"] -> "<<tab[i]<< std::endl;
+	}
+	tab.clear();
+		
+	tab = image->GetUpperRightCorner();
+	for(int i = 0 ; i < tab.size() ; i++ )
+	{
+		file << " UR[" <<i<<"] -> "<<tab[i]<< std::endl;
+	}
+	tab.clear();
+		
+	tab = image->GetLowerLeftCorner();
+	for(int i = 0 ; i < tab.size() ; i++ )
+	{
+		file << " LL[" <<i<<"] -> "<<tab[i]<< std::endl;
+	}
+	tab.clear();
+			
+	tab = image->GetLowerRightCorner();
+	for(int i = 0 ; i < tab.size() ; i++ )
+	{
+		file << " LR[" <<i<<"] -> "<<tab[i]<< std::endl;
+	}	
+	tab.clear();
+		
+	file.close();
+	
+  } 
+  catch( itk::ExceptionObject & err ) 
+  { 
+    std::cerr << "Exception OTB attrappee dans exception ITK !" << std::endl; 
+    std::cerr << err << std::endl; 
+    return EXIT_FAILURE;
+  } 
+  catch( ... )
+  {
+    std::cerr << "Exception OTB non attrappee !" << std::endl; 
+    return EXIT_FAILURE;
+  }
+  
+  return EXIT_SUCCESS;
+}
+