X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;f=common%2FclitkCommon.txx;h=d7c91405dda8f45ff76138c341784f4d05b16b4c;hb=a9bddd3844a7b1fed036cbe9de2e1cc00b92598f;hp=f9eeb8c8bfd9280319f0ed0d38a47a393827b87b;hpb=0b7c9b1e1215634b02cbd38d4e4ba101d6111ba8;p=clitk.git

diff --git a/common/clitkCommon.txx b/common/clitkCommon.txx
index f9eeb8c..d7c9140 100644
--- a/common/clitkCommon.txx
+++ b/common/clitkCommon.txx
@@ -1,7 +1,7 @@
 /*=========================================================================
   Program:   vv                     http://www.creatis.insa-lyon.fr/rio/vv
 
-  Authors belong to: 
+  Authors belong to:
   - University of LYON              http://www.universite-lyon.fr/
   - Léon Bérard cancer center       http://oncora1.lyon.fnclcc.fr
   - CREATIS CNRS laboratory         http://www.creatis.insa-lyon.fr
@@ -22,13 +22,52 @@
    * @file   clitkCommon.txx
    * @author David Sarrut <david.sarrut@creatis.insa-lyon.fr>
    * @date   18 May 2006
-   * 
+   *
    -------------------------------------------------*/
 
+//-------------------------------------------------------
+// Utility functions for text file parsing (author: joel schaerer)
+
+template<class ElementType>
+ElementType parse_value(std::string str)
+{
+  std::istringstream parser(str);
+  ElementType value;
+  parser >> value;
+  assert(!parser.fail());
+  return value;
+}
+
+template<class ElementType>
+std::vector<ElementType> parse_string(std::string str,char delim)
+{
+  std::istringstream ss(str);
+  std::string token;
+  std::vector<ElementType> result;
+  while (std::getline(ss,token,delim)) {
+    result.push_back(parse_value<ElementType>(token));
+  }
+  return result;
+}
+
+template<class ElementType>
+std::vector<std::vector<ElementType> > parse_file(const char* filename,char delim)
+{
+  std::ifstream fs(filename);
+  std::string line;
+  std::vector<std::vector<ElementType> > result;
+  while (std::getline(fs,line)) {
+    if (line[0] != '#') //skip comments
+      result.push_back(parse_string<ElementType>(line,delim));
+  }
+  return result;
+}
+
 //--------------------------------------------------------------------
 // Convert float, double ... to string
 template<class T>
-std::string toString(const T & t) {
+std::string toString(const T & t)
+{
   std::ostringstream myStream;
   myStream << t << std::flush;
   return(myStream.str());
@@ -38,7 +77,8 @@ std::string toString(const T & t) {
 //--------------------------------------------------------------------
 // Convert float*, double* ... to string
 template<class T>
-std::string toStringVector(const T * t, const int n) {
+std::string toStringVector(const T * t, const int n)
+{
   std::ostringstream myStream;
   for(int i=0; i<n-1; i++)
     myStream << clitk::toString<T>(t[i]) << " ";
@@ -50,7 +90,8 @@ std::string toStringVector(const T * t, const int n) {
 //--------------------------------------------------------------------
 // Convert float*, double* ... to string
 template<class T>
-std::string toStringVector(const T & t, const int n) {
+std::string toStringVector(const T & t, const int n)
+{
   std::ostringstream myStream;
   for(int i=0; i<n-1; i++)
     myStream << t[i] << " ";
@@ -62,7 +103,8 @@ std::string toStringVector(const T & t, const int n) {
 //--------------------------------------------------------------------
 // Convert float*, double* ... to string
 template<class T>
-std::string toStringVector(const std::vector<T> & t) {
+std::string toStringVector(const std::vector<T> & t)
+{
   return toStringVector(&t[0], t.size());
 }
 //--------------------------------------------------------------------
@@ -70,7 +112,8 @@ std::string toStringVector(const std::vector<T> & t) {
 //--------------------------------------------------------------------
 // Convert a pixel type to another (downcast)
 template<class TPixelUp, class TPixelDown>
-TPixelDown PixelTypeDownCast(const TPixelUp & x) {
+TPixelDown PixelTypeDownCast(const TPixelUp & x)
+{
   return (TPixelDown)lrint(x);
 }
 //--------------------------------------------------------------------
@@ -79,7 +122,7 @@ TPixelDown PixelTypeDownCast(const TPixelUp & x) {
 template<class Type>
 struct vectorComparisonLowerThan: public std::binary_function<int, int, bool> {
   vectorComparisonLowerThan(const std::vector<Type> & v):vect(v) {};
-  bool operator()(int x, int y) { 
+  bool operator()(int x, int y) {
     return (vect[x] < vect[y]);
   }
   const std::vector<Type> & vect;
@@ -90,7 +133,7 @@ struct vectorComparisonLowerThan: public std::binary_function<int, int, bool> {
 template<class Type>
 struct vectorComparisonGreaterThan: public std::binary_function<int, int, bool> {
   vectorComparisonGreaterThan(const std::vector<Type> & v):vect(v) {};
-  bool operator()(int x, int y) { 
+  bool operator()(int x, int y) {
     return (vect[x] > vect[y]);
   }
   const std::vector<Type> & vect;
@@ -99,23 +142,25 @@ struct vectorComparisonGreaterThan: public std::binary_function<int, int, bool>
 
 //--------------------------------------------------------------------
 template<class Type>
-void GetSortedIndex(const std::vector<Type> & toSort, std::vector<int> & index, bool increasing) {
-  index.resize(toSort.size());  
+void GetSortedIndex(const std::vector<Type> & toSort, std::vector<int> & index, bool increasing)
+{
+  index.resize(toSort.size());
   for(unsigned int i=0; i<index.size(); i++) index[i] = i;
-  if (increasing) 
-    std::sort(index.begin(), 
-	      index.end(), 
-	      vectorComparisonLowerThan<double>(toSort));
-  else 
-    std::sort(index.begin(), 
-	      index.end(), 
-	      vectorComparisonGreaterThan<double>(toSort));
+  if (increasing)
+    std::sort(index.begin(),
+              index.end(),
+              vectorComparisonLowerThan<double>(toSort));
+  else
+    std::sort(index.begin(),
+              index.end(),
+              vectorComparisonGreaterThan<double>(toSort));
 }
 //--------------------------------------------------------------------
 
 //--------------------------------------------------------------------
 template<class TPixel>
-std::string GetTypeAsString() {
+std::string GetTypeAsString()
+{
   //  http://www.vtk.org/doc/release/3/html/vtkSetGet_8h-source.html
   // and
   // itkImageIOBase.cxx
@@ -137,17 +182,18 @@ std::string GetTypeAsString() {
 
 //--------------------------------------------------------------------
 template<class ImageType>
-void CloneImage(const typename ImageType::Pointer & input, typename ImageType::Pointer & output) {
+void CloneImage(const typename ImageType::Pointer & input, typename ImageType::Pointer & output)
+{
   output->SetRegions(input->GetLargestPossibleRegion());
   output->SetOrigin(input->GetOrigin());
   output->SetSpacing(input->GetSpacing());
   output->Allocate();
-  typedef itk::ImageRegionConstIterator<ImageType> ConstIteratorType; 
+  typedef itk::ImageRegionConstIterator<ImageType> ConstIteratorType;
   ConstIteratorType pi(input,input->GetLargestPossibleRegion());
-  pi.GoToBegin();  
-  typedef itk::ImageRegionIterator<ImageType> IteratorType; 
+  pi.GoToBegin();
+  typedef itk::ImageRegionIterator<ImageType> IteratorType;
   IteratorType po(output,input->GetLargestPossibleRegion());
-  po.GoToBegin(); 
+  po.GoToBegin();
   while (!pi.IsAtEnd()) {
     po.Set(pi.Get());
     ++pi;