1 #include <creaImageIOGimmickView.h>
2 #include <creaImageIOSystem.h>
3 #include "boost/filesystem.hpp"
6 #include <gdcmGlobal.h>
7 #include <gdcmSerieHelper.h>
8 #include <vtkGdcmReader.h>
11 #if defined(USE_GDCM2)
12 #include <vtkGDCMImageReader.h>
16 namespace fs = boost::filesystem;
20 ///Class used to represent the actual state of the image selected and to perform comparisons on its values
24 ImageExtent(const std::string& x, const std::string& y, const std::string& z, const std::string& t)
26 sscanf(x.c_str(),"%d",&mExtent[0]);
27 sscanf(y.c_str(),"%d",&mExtent[1]);
28 sscanf(z.c_str(),"%d",&mExtent[2]);
29 sscanf(t.c_str(),"%d",&mExtent[3]);
30 if(x==""){mExtent[0]=1;}
31 if(y==""){mExtent[1]=1;}
32 if(z==""){mExtent[2]=1;}
33 if(t==""){mExtent[3]=1;}
35 if (mExtent[3]>1) mDim=4;
36 else if (mExtent[2]>1) mDim=3;
37 else if (mExtent[1]>1) mDim=2;
38 else if (mExtent[0]>1) mDim=1;
44 void Clear() { mExtent[0] = mExtent[1] = mExtent[2] = mExtent[3] = 1; }
46 ///Returns true if the two extents are compatible
47 bool IsCompatible( const ImageExtent& );
49 ///Adds the extent passed as a parameter to the current extent
50 void Add ( const ImageExtent& );
52 ///Returns the ieth position of the extent
53 int Get(int i) { return mExtent[i]; }
55 ///Returns the dimension of the current image
56 void SetDimension(int dim) { mDim=dim; }
58 ///Returns the dimension of the current image
59 int GetDimension() { return mDim; }
66 //======================================================================
68 //======================================================================
70 GimmickView::GimmickView(boost::shared_ptr<Gimmick> gimmick, int threads)
74 GimmickDebugMessage(1,"GimmickView::GimmickView"
76 // Anciently started the threads ...
77 // Threads now automatically start at first image request
81 //======================================================================
83 //======================================================================
85 GimmickView::~GimmickView()
87 GimmickDebugMessage(1,"GimmickView::~GimmickView"
90 //======================================================================
92 //======================================================================
93 /// Initializes the view :
94 /// Creates the TreeViews for all the TreeHandler of the Controller
96 void GimmickView::Initialize()
100 //======================================================================
102 //======================================================================
104 void GimmickView::Finalize()
108 //======================================================================
110 //======================================================================
111 /// Create the tree views
112 void GimmickView::CreateTreeViews()
114 GimmickMessage(2,"Creating the tree views"<<std::endl);
115 Gimmick::TreeHandlerMapType::const_iterator i;
116 for (i = mGimmick->GetTreeHandlerMap().begin();
117 i!= mGimmick->GetTreeHandlerMap().end();
120 this->CreateTreeView(i->second);
124 /// Create a tree view with a given name
125 void GimmickView::CreateSingleTreeView(std::string &i_name)
127 this->CreateTreeView(mGimmick->GetTreeHandlerMap()[i_name]);
132 //======================================================================
134 //======================================================================
135 /// Updates the TreeView of given name from level l to bottom
136 /// (calls the virtual method TreeView::Update())
137 void GimmickView::UpdateTreeViewLevel(const std::string& t, int l)
139 TreeViewMapType::iterator i;
140 i = GetTreeViewMap().find(t);
141 if ( i == GetTreeViewMap().end() )
143 GimmickError("INTERNAL ERROR : GimmickView::UpdateTreeView : '"
144 <<t<<"' is not in TreeViewMap");
146 i->second->UpdateLevel(l);
149 //======================================================================
150 /// Clears the status and begins a new selection process
151 void GimmickView::ResetExtent()
155 mImageExtent.reset();
160 //======================================================================
162 //======================================================================
163 bool ImageExtent::IsCompatible(const ImageExtent& ie)
165 bool compatible=true;
166 ImageExtent * extent= (ImageExtent*)&ie;
167 if((*extent).Get(0)!=Get(0)
168 || (*extent).Get(1)!=Get(1))
175 //======================================================================
177 //======================================================================
178 void ImageExtent::Add(const ImageExtent& ie)
180 ImageExtent * extent= (ImageExtent*)&ie;
181 mExtent[2]+=(*extent).Get(2);
188 //======================================================================
189 /// No selected image
190 bool GimmickView::NoValidateSelected ()
192 GimmickDebugMessage(2,"Validating selected"<<std::endl);
193 std::string mMessage;
194 mMessage="Cannot have 0 images selected!";
196 modifyValidationSignal(valid);
197 SetMessage(mMessage);
201 //======================================================================
202 ///Validates the dimension compliance of the images with the maximum and
203 ///minimum given, and between their sizes
204 bool GimmickView::ValidateSelected (tree::Node* sel, int min_dim, int max_dim)
206 GimmickDebugMessage(2,"Validating selected"<<std::endl);
207 std::string mMessage;
211 mMessage="Cannot have 0 images selected!";
216 boost::shared_ptr<ImageExtent> ie=boost::shared_ptr<ImageExtent>(new ImageExtent((*sel).GetAttribute("D0028_0010"),
217 (*sel).GetAttribute("D0028_0011"),
218 (*sel).GetAttribute("D0028_0012"),
223 if((mImageExtent->Get(min_dim-1)<2)||(mImageExtent->Get(max_dim)>1))
229 std::stringstream out;
230 out << mImageExtent->GetDimension() << "D image " << mImageExtent->Get(0) << "x"<< mImageExtent->Get(1) << "x"<< mImageExtent->Get(2) <<" selected";
231 mMessage = out.str();
232 mImageExtent->SetDimension(2);
238 if(mImageExtent->IsCompatible(*ie))
240 if(mImageExtent->GetDimension()==max_dim && mImageExtent->Get(max_dim)>2)
242 std::stringstream out;
243 out<<"Cannot add this image to selection : would result in a "<<mImageExtent->GetDimension()+1<<"D image!";
249 std::stringstream out;
250 out<<"Selecting "<<mImageExtent->GetDimension()<<"D images is not allowed !";
254 else if(min_dim==3 && (ie->Get(2)+mImageExtent->Get(2))<2)
256 std::stringstream out;
257 out << "Cannot build the selection as it would result in a ";
258 out << mImageExtent->GetDimension();
259 out << "D image, and the minimum is ";
267 mImageExtent->Add(*ie);
268 std::stringstream out;
269 out << mImageExtent->GetDimension() << "D image " << mImageExtent->Get(0) << "x"<< mImageExtent->Get(1) << "x"<< mImageExtent->Get(2) <<" selected";
270 mMessage = out.str();
275 mMessage="The selected images are not compatible.";
281 modifyValidationSignal(valid);
282 SetMessage(mMessage);
286 //======================================================================
287 void GimmickView::modifyValidationSignal(bool ivalid)
289 mValidationSignal(ivalid);
292 //======================================================================
293 ///Reads Images (Non Threaded)
294 void GimmickView::ReadImagesNotThreaded(std::vector<vtkImageData*>& s, std::vector<std::string> im, int dimension)
299 #define GIMMICK_NO_IMAGE_SELECTION 0
300 #define GIMMICK_2D_IMAGE_SELECTION 2
301 #define GIMMICK_3D_IMAGE_SELECTION 3
302 #define GIMMICK_4D_IMAGE_SELECTION 4
309 // Create the output data
312 vtkImageData * out=vtkImageData::New();
313 out->ShallowCopy(mReader.GetImage(im.front()));
316 else if (im.size()>1) // Test inutile ? JPR
318 vtkImageData* first = mReader.GetImage( im.front());
322 std::vector<std::string>::iterator it;
323 for (it=im.begin(); it!=im.end(); ++it)
325 vtkImageData* out = vtkImageData::New();
326 out->ShallowCopy(mReader.GetImage(*it));
333 vtkImageData* out = vtkImageData::New();
334 // out->CopyStructure(first);
335 out->SetScalarType(first->GetScalarType());
336 out->SetNumberOfScalarComponents(first->GetNumberOfScalarComponents());
338 //first->GetExtent(ext); // JPR
339 first->GetWholeExtent(ext); // renvoie egalement 0,0 en Z // JPR
343 ext[5] = im.size()-1;
347 ext[5] = ext[5] * im.size()-1; // to deal with multiframes - JPR
351 // LG : TODO : Z Spacing ?
354 first->GetDimensions(dim);
356 out->SetDimensions(dim[0], dim[1], im.size() );
357 out->AllocateScalars();
360 unsigned long imsize = dim[0] * dim[1];
361 imsize = imsize * dim[2] ; // deal with multiframes // JPR
365 // differents formats char , short, etc...
366 // differents components 1..3 ex. jpg ->RGB 3
367 imsize = imsize * first->GetScalarSize() * first->GetNumberOfScalarComponents();
370 // Order the file name vector
373 first->GetSpacing(spc);
374 spc[2]=OrderTheFileNameVector(im);
376 out->SetSpacing(spc);
379 std::vector<std::string>::iterator it;
381 for (it=im.begin(); it!=im.end(); ++it)
383 vtkImageData* cur = mReader.GetImage( (*it) );
384 memcpy(out->GetScalarPointer(0,0,slice), cur->GetScalarPointer(0,0,0), imsize);
394 //======================================================================
396 #if defined(USE_GDCM)
397 double GimmickView::OrderTheFileNameVector(std::vector<std::string> &im)
400 typedef std::vector<GDCM_NAME_SPACE::File* > FileList;
402 //GDCM_NAME_SPACE::File *f = GDCM_NAME_SPACE::File::New();
403 GDCM_NAME_SPACE::SerieHelper *sh = GDCM_NAME_SPACE::SerieHelper::New();
404 std::vector<std::string> lstAux;
405 std::vector<std::string>::iterator it;
406 for (it=im.begin(); it!=im.end(); ++it)
408 ///\TODO liberer les GDCM_NAME_SPACE::File a la fin! // JPR
409 GDCM_NAME_SPACE::File *f = GDCM_NAME_SPACE::File::New();
414 fileVector.push_back(f);
416 lstAux.push_back(*it);
421 if ((fileVector.size()>1) && (sh->IsCoherent( &fileVector )))
423 sh->OrderFileList(&fileVector);
424 spacing= sh->GetZSpacing();
427 for (i=0; i<fileVector.size(); i++)
429 im.push_back( (fileVector[i])->GetFileName() );
431 for (i=0; i<lstAux.size(); i++)
433 im.push_back( lstAux[i] );
436 std::sort( im.begin(), im.end() );
445 #if defined(USE_GDCM2)
447 double GimmickView::OrderTheFileNameVector(std::vector<std::string> &im)
452 //======================================================================
455 void GimmickView::ReadImagesNotThreadedInVector(std::vector<vtkImageData*>& s, std::vector<std::string> im, int dimension)
457 // Create the output data
460 // Only one image : give it
461 vtkImageData* out = vtkImageData::New();
462 GimmickDebugMessage(3, "State Check: Full Filename: "
465 out->ShallowCopy(mReader.GetImage(im.front()));
468 else if (im.size()>1) // Test inutile ? JPR
470 vtkImageData* first = mReader.GetImage( im.front());
474 std::vector<std::string>::iterator it;
475 for (it=im.begin(); it!=im.end(); ++it)
477 vtkImageData* out = vtkImageData::New();
478 out->ShallowCopy(mReader.GetImage(*it));
485 // n *2D + T in a vector :
487 std::vector<std::string>::iterator it;
488 for (it=im.begin(); it!=im.end(); ++it)
490 vtkImageData* out = mReader.GetImage( (*it));
496 //======================================================================
498 //======================================================================
499 ///Requests the reading of an image
500 void GimmickView::RequestReading(tree::Node* n,
501 int prio, int selection_index, boost::shared_ptr<ImagePointerHolder> p)
508 ImageEventType t(n,selection_index);
510 mImageEventMap[n->GetAttribute("FullFileName")] = t;
511 mReader.Request(this,n->GetAttribute("FullFileName"),prio);
513 //======================================================================
515 //======================================================================
517 OnMultiThreadImageReaderEvent(const std::string& filename,
518 MultiThreadImageReaderUser::EventType e,
521 GimmickDebugMessage(7,
522 "MultiThreadImageReader event : "<<e<<std::endl);
525 if (filename.size()==0)
527 //What to do in this case?
529 GimmickDebugMessage(5,
530 "Pushing unknown image in queue"
532 mImageEventQueue.push_back(ImageEventType(image));*/
535 ImageEventTypeMap::iterator i;
536 //JCP 22-06-2009, test mImageEventMap.size() > 0
537 if(mImageEventMap.size()>0){
538 i = mImageEventMap.find(filename);
539 if (i!=mImageEventMap.end())
541 GimmickDebugMessage(5,
542 "Putting image of file '"<<filename<<"' on pointer"
544 ImageEventType ie(i->second);
546 ie.pointerHolder->Set(ie.image);
547 //mImageEventMap.erase(i);
553 std::string mess="ERROR: MultiThreadImageReader: Cannot read image in file ";
556 GimmickMessage(1,mess);
557 ImageEventTypeMap::iterator i;
558 i = mImageEventMap.find(filename);
559 if (i!=mImageEventMap.end())
561 ImageEventType ie(i->second);
563 ie.pointerHolder->Set(GetDefaultImage());
564 //mImageEventMap.erase(i);
568 else if (e==ImageUnloaded)
570 std::string mess="Unloaded image in file ";
573 GimmickMessage(1,mess);
574 ImageEventTypeMap::iterator i;
575 i = mImageEventMap.find(filename);
576 if (i!=mImageEventMap.end())
578 ImageEventType ie(i->second);
580 ie.pointerHolder->Set(GetDefaultImage());
581 //mImageEventMap.erase(i);
586 //====================================================================
588 //====================================================================
589 void GimmickView::ConnectValidationObserver(ValidationCallbackType callback)
591 mValidationSignal.connect(callback);
594 } // EO namespace creaImageIO