X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;f=src2%2FcreaImageIOGimmickView.cpp;h=30ca0d5d094b9bd430e6e1824da392bb64057bf7;hb=a13b9d55cc9b8ca9b886e99a73946dd8eccefe65;hp=3e958f51b6b53655ed1234faabc17e7d9e5c9f31;hpb=f1caae70392f96f30dd90eb073be6be10fbf1caa;p=creaImageIO.git

diff --git a/src2/creaImageIOGimmickView.cpp b/src2/creaImageIOGimmickView.cpp
index 3e958f5..30ca0d5 100644
--- a/src2/creaImageIOGimmickView.cpp
+++ b/src2/creaImageIOGimmickView.cpp
@@ -1,15 +1,75 @@
 #include <creaImageIOGimmickView.h>
 #include <creaImageIOSystem.h>
+#include "boost/filesystem.hpp"
+
+#include <gdcmSerieHelper.h>    // JPR
+#include <gdcmFile.h>           // JPR
+#include <vtkGdcmReader.h>      // JPR
+#include <gdcmDebug.h>
+
+namespace fs = boost::filesystem;
 
 namespace creaImageIO
 {
+
+	///Class used to represent the actual state of the image selected and to perform comparisons on its values
+  class ImageExtent
+  {
+  public:
+    ImageExtent(const std::string& x, const std::string& y, const std::string& z, const std::string& t)
+	{
+	   sscanf(x.c_str(),"%d",&mExtent[0]);
+	   sscanf(y.c_str(),"%d",&mExtent[1]);
+	   sscanf(z.c_str(),"%d",&mExtent[2]);
+	   sscanf(t.c_str(),"%d",&mExtent[3]);
+	   if(x==""){mExtent[0]=1;}
+	   if(y==""){mExtent[1]=1;}
+	   if(z==""){mExtent[2]=1;}
+	   if(t==""){mExtent[3]=1;}
+
+	   if (mExtent[3]>1) mDim=4;
+	   else if (mExtent[2]>1) mDim=3;
+	   else if (mExtent[1]>1) mDim=2;
+	   else if (mExtent[0]>1) mDim=1;
+	   else mDim=0;
+	}
+
+	///Clears the extent
+    void Clear() { mExtent[0] = mExtent[1] = mExtent[2] = mExtent[3] = 1; }
+
+	///Returns true if the two extents are compatible
+    bool IsCompatible( const ImageExtent& );
+
+	///Adds the extent passed as a parameter to the current extent
+    void Add ( const ImageExtent& );
+		
+	///Returns the ieth position of the extent
+    int Get(int i) { return mExtent[i]; }
+    
+	///Returns the dimension of the current image
+    void SetDimension(int dim) { mDim=dim; }
+
+	///Returns the dimension of the current image
+    int GetDimension() { return mDim; }
+
+  private:
+    int mExtent[4];
+    int mDim;
+  };
+
+  //======================================================================
+
   //======================================================================
   // CTor
-  GimmickView::GimmickView(Gimmick* gimmick)
-    : mGimmick(gimmick)
+  GimmickView::GimmickView(boost::shared_ptr<Gimmick> gimmick, int threads)
+    : mGimmick(gimmick),
+	  mReader(threads)
   {
     GimmickDebugMessage(1,"GimmickView::GimmickView"
 			<<std::endl);
+	// Anciently started the threads ...
+    // Threads now automatically start at first image request
+    //mReader.Start();
   }
   //======================================================================
 
@@ -28,9 +88,10 @@ namespace creaImageIO
   /// 
   void GimmickView::Initialize()
   {
+	mReaderStarted=false;
   }
   //======================================================================
-  
+
   //======================================================================
   /// Finalize 
   void GimmickView::Finalize()
@@ -52,6 +113,14 @@ namespace creaImageIO
 	this->CreateTreeView(i->second);
       }
   }
+
+  /// Create a tree view with a given name
+  void GimmickView::CreateSingleTreeView(std::string &i_name)
+  {
+	this->CreateTreeView(mGimmick->GetTreeHandlerMap()[i_name]);
+     
+  }
+
   //======================================================================
 
   //======================================================================
@@ -66,10 +135,450 @@ namespace creaImageIO
 	GimmickError("INTERNAL ERROR : GimmickView::UpdateTreeView : '"
 		     <<t<<"' is not in TreeViewMap");
       }
-    i->second->UpdateLevel(l);    
+    i->second->UpdateLevel(l);
+  }
+
+  //======================================================================
+  /// Clears the status and begins a new selection process
+  void GimmickView::ResetExtent()
+  {
+	  if(mImageExtent!=0)
+	  {
+		  mImageExtent.reset();
+	  }
+	  valid=true;
+  }
+
+  //======================================================================
+
+  //======================================================================
+  bool ImageExtent::IsCompatible(const ImageExtent& ie)
+  {
+	  bool compatible=true;
+	  ImageExtent * extent= (ImageExtent*)&ie;
+	  if((*extent).Get(0)!=Get(0)
+		 || (*extent).Get(1)!=Get(1))
+	  {
+		  compatible=false;
+	  }
+	  return compatible;
+  }
+
+  //======================================================================
+
+  //======================================================================
+  void ImageExtent::Add(const ImageExtent& ie)
+  {
+	  ImageExtent * extent= (ImageExtent*)&ie;
+	  mExtent[2]+=(*extent).Get(2);
+	  if(mExtent[2]>1)
+	  {
+	  SetDimension(3);
+	  }
+  }
+
+  //======================================================================
+  /// No selected image
+  bool GimmickView::NoValidateSelected ()
+  {
+	GimmickDebugMessage(2,"Validating selected"<<std::endl);
+	std::string mMessage;
+	mMessage="Cannot have 0 images selected!";
+	valid=false;
+	modifyValidationSignal(valid);
+	SetMessage(mMessage);
+	return valid;
+  }
+
+  //======================================================================
+  ///Validates the dimension compliance of the images with the maximum and 
+  ///minimum given, and between their sizes
+  bool GimmickView::ValidateSelected (tree::Node* sel, int min_dim, int max_dim)
+  {
+	GimmickDebugMessage(2,"Validating selected"<<std::endl);
+	std::string mMessage;
+	
+	if(sel==0)
+	{
+		mMessage="Cannot have 0 images selected!";
+		valid=false;
+	}
+	else
+	{
+		boost::shared_ptr<ImageExtent> ie=boost::shared_ptr<ImageExtent>(new ImageExtent((*sel).GetAttribute("D0028_0010"),
+			                         			(*sel).GetAttribute("D0028_0011"),
+									(*sel).GetAttribute("D0028_0012"), 
+									""));
+	if(mImageExtent==0)
+	{
+		mImageExtent=ie;
+		if((mImageExtent->Get(min_dim-1)<2)||(mImageExtent->Get(max_dim)>1))
+		{
+			valid=false;
+		}
+		else
+		{
+			std::stringstream out;
+			out << mImageExtent->GetDimension() << "D image " << mImageExtent->Get(0) << "x"<< mImageExtent->Get(1) << "x"<< mImageExtent->Get(2) <<" selected";
+	        	mMessage = out.str();
+			mImageExtent->SetDimension(2);
+			valid=true;
+		}
+	}
+	else
+	{
+		if(mImageExtent->IsCompatible(*ie))
+		{
+			if(mImageExtent->GetDimension()==max_dim && mImageExtent->Get(max_dim)>2)
+			{
+				std::stringstream out;
+				out<<"Cannot add this image to selection : would result in a "<<mImageExtent->GetDimension()+1<<"D image!";
+				mMessage=out.str();
+				valid=false;
+			}
+			else if(max_dim<3)
+			{
+				std::stringstream out;
+				out<<"Selecting "<<mImageExtent->GetDimension()<<"D images is not allowed !";
+				mMessage=out.str();
+				valid=false;
+			}
+			else if(min_dim==3 && (ie->Get(2)+mImageExtent->Get(2))<2)
+			{
+				std::stringstream out;
+				out << "Cannot build the selection as it would result in a ";
+				out << mImageExtent->GetDimension();
+				out << "D image, and the minimum is ";
+				out << min_dim;
+				out << "D!";
+				mMessage=out.str();
+				valid=false;
+			}
+			else
+			{
+            			mImageExtent->Add(*ie);
+				std::stringstream out;
+				out << mImageExtent->GetDimension() << "D image " << mImageExtent->Get(0) << "x"<< mImageExtent->Get(1) << "x"<< mImageExtent->Get(2) <<" selected";
+	        		mMessage = out.str();
+			}
+		}
+		else
+		{
+			mMessage="The selected images are not compatible.";
+			valid=false;
+		}
+	  }
+	}
+
+	modifyValidationSignal(valid);
+	SetMessage(mMessage);
+	return valid;
+  }
+
+  //======================================================================
+  void GimmickView::modifyValidationSignal(bool ivalid)
+  {
+	  mValidationSignal(ivalid);
+  }
+
+   //======================================================================
+  ///Reads Images (Non Threaded)
+void GimmickView::ReadImagesNotThreaded(std::vector<vtkImageData*>& s, std::vector<std::string> im, int dimension)
+{
+
+/* remember!
+
+#define GIMMICK_NO_IMAGE_SELECTION 0
+#define GIMMICK_2D_IMAGE_SELECTION 2
+#define GIMMICK_3D_IMAGE_SELECTION 3
+#define GIMMICK_4D_IMAGE_SELECTION 4
+
+#define NATIVE 0
+#define _2D    2
+#define _3D    3
+
+*/
+	// Create the output data
+	if (im.size()==1) 
+	{
+		// O0.
+	}
+    else if (im.size()>1) // Test inutile ? JPR
+	{
+		// Read the FIRST image to get info ... about all the other ones. ?!? // JPR
+		vtkImageData* first = mReader.GetImage( im.front());
+
+		if (dimension == 2)
+		{
+		 // n3D
+		    std::vector<std::string>::iterator it;
+			for (it=im.begin(); it!=im.end(); ++it) 
+			{
+				vtkImageData* out = vtkImageData::New();
+				out->ShallowCopy(mReader.GetImage(*it));
+				s.push_back(out);
+			}
+		}	  
+		else if (dimension == 3)
+		{
+			// n*2D to 3D  :: --> NON ! Ca va tres bien pour du 2D+T (le ZSpacing, on s'en tape)
+			//			// ne convient pas pour du n*2D a  lire comme un 3D
+			//			// Il faut un gdcmSerieHelper !
+			vtkImageData* out = vtkImageData::New();
+			out->CopyStructure(first);	
+			out->SetScalarType(first->GetScalarType());
+			out->SetNumberOfScalarComponents(first->GetNumberOfScalarComponents());
+			int ext[6];
+			//first->GetExtent(ext);     // JPR
+			first->GetWholeExtent(ext);  // renvoie egalement 0,0 en Z // JPR
+						     // Normal :
+
+			if(ext[5] == 0)
+			{
+				ext[5] = im.size()-1;
+			}
+			else
+			{
+				ext[5] = ext[5] * im.size()-1; // to deal with multiframes - JPR
+			}
+			out->SetExtent(ext);
+
+			// LG : TODO : Z Spacing  ?  // Voir + loin // JPR
+
+			int dim[3];
+			first->GetDimensions(dim);
+
+			out->SetDimensions(dim[0], dim[1], im.size() );
+			out->AllocateScalars();
+			out->Update();
+
+			unsigned long imsize = dim[0] * dim[1];
+			imsize = imsize * dim[2] ;  // deal with multiframes // JPR
+
+
+//EED 03-11-2009
+			// differents formats char , short, etc...
+			// differents components 1..3  ex. jpg ->RGB 3
+			imsize = imsize * first->GetScalarSize() * first->GetNumberOfScalarComponents();
+
+
+			// Order the file name vector
+
+			double spc[3];
+			out->GetSpacing(spc);
+			spc[2]=OrderTheFileNameVector(im);	
+			out->SetSpacing(spc);
+			int slice = 0;
+			std::vector<std::string>::iterator it;
+
+			for (it=im.begin(); it!=im.end(); ++it) 
+			{
+				vtkImageData* cur = mReader.GetImage( (*it) );
+				memcpy(out->GetScalarPointer(0,0,slice), cur->GetScalarPointer(0,0,0), imsize);
+				slice++;
+			}	
+			s.push_back(out);
+
+		}  // dimension == 3
+
+	}  // size >1
+
+}
+  //======================================================================
+
+
+double GimmickView::OrderTheFileNameVector(std::vector<std::string> &im)
+{
+	double spacing=1;
+	typedef std::vector<GDCM_NAME_SPACE::File* > FileList;
+	FileList fileVector;
+	//GDCM_NAME_SPACE::File *f = GDCM_NAME_SPACE::File::New();
+	GDCM_NAME_SPACE::SerieHelper *sh = GDCM_NAME_SPACE::SerieHelper::New();
+	std::vector<std::string> lstAux;
+	std::vector<std::string>::iterator it;
+	for (it=im.begin(); it!=im.end(); ++it)
+	{
+		///\TODO liberer les GDCM_NAME_SPACE::File a la fin!  // JPR
+		GDCM_NAME_SPACE::File *f = GDCM_NAME_SPACE::File::New();
+		f->SetFileName(*it);
+		f->Load();
+		if (f->IsReadable())
+		{
+			fileVector.push_back(f);
+		} else {
+			lstAux.push_back(*it);
+		}
+	} // for
+			
+	
+	if ((fileVector.size()>1) && (sh->IsCoherent( &fileVector )))
+	{
+			sh->OrderFileList(&fileVector);
+			spacing= sh->GetZSpacing();
+			im.clear();
+			int i;
+			for (i=0; i<fileVector.size(); i++)
+			{
+				im.push_back( (fileVector[i])->GetFileName() );
+			}
+			for (i=0; i<lstAux.size(); i++)
+			{
+				im.push_back( lstAux[i] );
+			}
+	}else {
+		std::sort( im.begin(), im.end() );
+	}
+			
+   return spacing;
+}
+
+
+//======================================================================
+
+
+void GimmickView::ReadImagesNotThreadedInVector(std::vector<vtkImageData*>& s, std::vector<std::string> im, int dimension)
+{
+	// Create the output data
+	if (im.size()==1)
+	{
+		// Only one image : give it
+		vtkImageData* out = vtkImageData::New();
+		GimmickDebugMessage(3, "State Check: Full Filename: "
+						<<im.front()
+						<<std::endl);
+		out->ShallowCopy(mReader.GetImage(im.front()));
+		s.push_back( out );
+	}
+    	else if (im.size()>1) // Test inutile ? JPR
+	{
+		vtkImageData* first = mReader.GetImage( im.front());
+		if (dimension == 2)
+		{
+		 // n3D
+		    std::vector<std::string>::iterator it;
+			for (it=im.begin(); it!=im.end(); ++it)
+			{
+			std::cout << "----------------------------------------n3D" << std::endl;
+				vtkImageData* out = vtkImageData::New();
+				out->ShallowCopy(mReader.GetImage(*it));
+				s.push_back(out);
+			}
+		}
+		else
+		{
+			// n2D to 3D // NO!
+			// n *2D + T in a vector :
+			std::cout << "---------------------------------------- n *2D + T" << std::endl;			
+			std::vector<std::string>::iterator it;
+			for (it=im.begin(); it!=im.end(); ++it) 
+			{
+			std::cout << "---------------------------------------- [" << *it << "]" << std::endl;
+				//vtkImageData* out = mReader.GetImage( (*it));
+				//s.push_back(out);
+				vtkImageData* out = vtkImageData::New();
+				out->ShallowCopy(mReader.GetImage(*it));
+				s.push_back(out);
+			}
+		}
+	}
+}
+  //======================================================================
+
+  //======================================================================
+  ///Requests the reading of an image
+  void GimmickView::RequestReading(tree::Node* n, 
+	  int prio, int selection_index, boost::shared_ptr<ImagePointerHolder> p)
+  {
+     if(!mReaderStarted)
+     {
+        mReader.Start();
+        mReaderStarted=true;
+      }
+      ImageEventType t(n,selection_index);
+      t.pointerHolder = p;
+      mImageEventMap[n->GetAttribute("FullFileName")] = t;
+      mReader.Request(this,n->GetAttribute("FullFileName"),prio);
   }
   //======================================================================
+
+  //======================================================================
+  void GimmickView::
+  OnMultiThreadImageReaderEvent(const std::string& filename,
+				MultiThreadImageReaderUser::EventType e,
+				vtkImageData* image)
+  {
+    GimmickDebugMessage(7,
+			"MultiThreadImageReader event : "<<e<<std::endl);
+	if (e==ImageLoaded)
+	{
+		if (filename.size()==0)
+		{
+		  //What to do in this case?
+		  /*
+			GimmickDebugMessage(5,
+					"Pushing unknown image in queue"
+					<<std::endl);
+			mImageEventQueue.push_back(ImageEventType(image));*/
+			return;
+		}
+		ImageEventTypeMap::iterator i;
+//JCP 22-06-2009, test mImageEventMap.size() > 0
+		if(mImageEventMap.size()>0){
+			i = mImageEventMap.find(filename);
+			if (i!=mImageEventMap.end())
+			{
+				GimmickDebugMessage(5,
+						"Putting image of file '"<<filename<<"' on pointer"
+						<<std::endl);
+				ImageEventType ie(i->second);
+				ie.image = image;
+				ie.pointerHolder->Set(ie.image);
+				//mImageEventMap.erase(i);
+			}
+		}    
+	}
+	else if (e==Error)
+	{
+		std::string mess="ERROR: MultiThreadImageReader: Cannot read image in file ";
+		mess+=filename;
+		mess+="\n";
+		GimmickMessage(1,mess);
+		ImageEventTypeMap::iterator i;
+		i = mImageEventMap.find(filename);
+		if (i!=mImageEventMap.end())
+		{
+		ImageEventType ie(i->second);
+		ie.image = image;
+		ie.pointerHolder->Set(GetDefaultImage());
+		//mImageEventMap.erase(i);
+		}
+	}
+
+	else if (e==ImageUnloaded)
+	{
+		std::string mess="Unloaded image in file ";
+		mess+=filename;
+		mess+="\n";
+		GimmickMessage(1,mess);
+				ImageEventTypeMap::iterator i;
+		i = mImageEventMap.find(filename);
+		if (i!=mImageEventMap.end())
+		{
+		ImageEventType ie(i->second);
+		ie.image = image;
+		ie.pointerHolder->Set(GetDefaultImage());
+		//mImageEventMap.erase(i);
+		}
+	}
+  }
+
+  //====================================================================
+
+  //====================================================================
+  void GimmickView::ConnectValidationObserver(ValidationCallbackType callback)
+  {
+    mValidationSignal.connect(callback);
+  }	  
   
 } // EO namespace creaImageIO
 
-