1 #include <creaImageIOGimmickView.h>
2 #include <creaImageIOSystem.h>
3 #include "boost/filesystem.hpp"
4 #include <gdcmGlobal.h>
5 #include <gdcmSerieHelper.h>
6 #include <vtkGdcmReader.h>
7 namespace fs = boost::filesystem;
11 ///Class used to represent the actual state of the image selected and to perform comparisons on its values
15 ImageExtent(const std::string& x, const std::string& y, const std::string& z, const std::string& t)
17 sscanf(x.c_str(),"%d",&mExtent[0]);
18 sscanf(y.c_str(),"%d",&mExtent[1]);
19 sscanf(z.c_str(),"%d",&mExtent[2]);
20 sscanf(t.c_str(),"%d",&mExtent[3]);
21 if(x==""){mExtent[0]=1;}
22 if(y==""){mExtent[1]=1;}
23 if(z==""){mExtent[2]=1;}
24 if(t==""){mExtent[3]=1;}
26 if (mExtent[3]>1) mDim=4;
27 else if (mExtent[2]>1) mDim=3;
28 else if (mExtent[1]>1) mDim=2;
29 else if (mExtent[0]>1) mDim=1;
35 void Clear() { mExtent[0] = mExtent[1] = mExtent[2] = mExtent[3] = 1; }
37 ///Returns true if the two extents are compatible
38 bool IsCompatible( const ImageExtent& );
40 ///Adds the extent passed as a parameter to the current extent
41 void Add ( const ImageExtent& );
43 ///Returns the ieth position of the extent
44 int Get(int i) { return mExtent[i]; }
46 ///Returns the dimension of the current image
47 void SetDimension(int dim) { mDim=dim; }
49 ///Returns the dimension of the current image
50 int GetDimension() { return mDim; }
57 //======================================================================
59 //======================================================================
61 GimmickView::GimmickView(boost::shared_ptr<Gimmick> gimmick, int threads)
65 GimmickDebugMessage(1,"GimmickView::GimmickView"
67 // Anciently started the threads ...
68 // Threads now automatically start at first image request
72 //======================================================================
74 //======================================================================
76 GimmickView::~GimmickView()
78 GimmickDebugMessage(1,"GimmickView::~GimmickView"
81 //======================================================================
83 //======================================================================
84 /// Initializes the view :
85 /// Creates the TreeViews for all the TreeHandler of the Controller
87 void GimmickView::Initialize()
91 //======================================================================
93 //======================================================================
95 void GimmickView::Finalize()
99 //======================================================================
101 //======================================================================
102 /// Create the tree views
103 void GimmickView::CreateTreeViews()
105 GimmickMessage(2,"Creating the tree views"<<std::endl);
106 Gimmick::TreeHandlerMapType::const_iterator i;
107 for (i = mGimmick->GetTreeHandlerMap().begin();
108 i!= mGimmick->GetTreeHandlerMap().end();
111 this->CreateTreeView(i->second);
115 /// Create a tree view with a given name
116 void GimmickView::CreateSingleTreeView(std::string &i_name)
118 this->CreateTreeView(mGimmick->GetTreeHandlerMap()[i_name]);
123 //======================================================================
125 //======================================================================
126 /// Updates the TreeView of given name from level l to bottom
127 /// (calls the virtual method TreeView::Update())
128 void GimmickView::UpdateTreeViewLevel(const std::string& t, int l)
130 TreeViewMapType::iterator i;
131 i = GetTreeViewMap().find(t);
132 if ( i == GetTreeViewMap().end() )
134 GimmickError("INTERNAL ERROR : GimmickView::UpdateTreeView : '"
135 <<t<<"' is not in TreeViewMap");
137 i->second->UpdateLevel(l);
140 //======================================================================
141 /// Clears the status and begins a new selection process
142 void GimmickView::ResetExtent()
146 mImageExtent.reset();
151 //======================================================================
153 //======================================================================
154 bool ImageExtent::IsCompatible(const ImageExtent& ie)
156 bool compatible=true;
157 ImageExtent * extent= (ImageExtent*)&ie;
158 if((*extent).Get(0)!=Get(0)
159 || (*extent).Get(1)!=Get(1))
166 //======================================================================
168 //======================================================================
169 void ImageExtent::Add(const ImageExtent& ie)
171 ImageExtent * extent= (ImageExtent*)&ie;
172 mExtent[2]+=(*extent).Get(2);
179 //======================================================================
180 /// No selected image
181 bool GimmickView::NoValidateSelected ()
183 GimmickDebugMessage(2,"Validating selected"<<std::endl);
184 std::string mMessage;
185 mMessage="Cannot have 0 images selected!";
187 modifyValidationSignal(valid);
188 SetMessage(mMessage);
192 //======================================================================
193 ///Validates the dimension compliance of the images with the maximum and
194 ///minimum given, and between their sizes
195 bool GimmickView::ValidateSelected (tree::Node* sel, int min_dim, int max_dim)
197 GimmickDebugMessage(2,"Validating selected"<<std::endl);
198 std::string mMessage;
202 mMessage="Cannot have 0 images selected!";
207 boost::shared_ptr<ImageExtent> ie=boost::shared_ptr<ImageExtent>(new ImageExtent((*sel).GetAttribute("D0028_0010"),
208 (*sel).GetAttribute("D0028_0011"),
209 (*sel).GetAttribute("D0028_0012"),
214 if((mImageExtent->Get(min_dim-1)<2)||(mImageExtent->Get(max_dim)>1))
220 std::stringstream out;
221 out << mImageExtent->GetDimension() << "D image " << mImageExtent->Get(0) << "x"<< mImageExtent->Get(1) << "x"<< mImageExtent->Get(2) <<" selected";
222 mMessage = out.str();
223 mImageExtent->SetDimension(2);
229 if(mImageExtent->IsCompatible(*ie))
231 if(mImageExtent->GetDimension()==max_dim && mImageExtent->Get(max_dim)>2)
233 std::stringstream out;
234 out<<"Cannot add this image to selection : would result in a "<<mImageExtent->GetDimension()+1<<"D image!";
240 std::stringstream out;
241 out<<"Selecting "<<mImageExtent->GetDimension()<<"D images is not allowed !";
245 else if(min_dim==3 && (ie->Get(2)+mImageExtent->Get(2))<2)
247 std::stringstream out;
248 out << "Cannot build the selection as it would result in a ";
249 out << mImageExtent->GetDimension();
250 out << "D image, and the minimum is ";
258 mImageExtent->Add(*ie);
259 std::stringstream out;
260 out << mImageExtent->GetDimension() << "D image " << mImageExtent->Get(0) << "x"<< mImageExtent->Get(1) << "x"<< mImageExtent->Get(2) <<" selected";
261 mMessage = out.str();
266 mMessage="The selected images are not compatible.";
272 modifyValidationSignal(valid);
273 SetMessage(mMessage);
277 //======================================================================
278 void GimmickView::modifyValidationSignal(bool ivalid)
280 mValidationSignal(ivalid);
283 //======================================================================
284 ///Reads Images (Non Threaded)
285 void GimmickView::ReadImagesNotThreaded(std::vector<vtkImageData*>& s, std::vector<std::string> im, int dimension)
290 #define GIMMICK_NO_IMAGE_SELECTION 0
291 #define GIMMICK_2D_IMAGE_SELECTION 2
292 #define GIMMICK_3D_IMAGE_SELECTION 3
293 #define GIMMICK_4D_IMAGE_SELECTION 4
300 // Create the output data
303 vtkImageData * out=vtkImageData::New();
304 out->ShallowCopy(mReader.GetImage(im.front()));
307 else if (im.size()>1) // Test inutile ? JPR
309 vtkImageData* first = mReader.GetImage( im.front());
313 std::vector<std::string>::iterator it;
314 for (it=im.begin(); it!=im.end(); ++it)
316 vtkImageData* out = vtkImageData::New();
317 out->ShallowCopy(mReader.GetImage(*it));
324 vtkImageData* out = vtkImageData::New();
325 // out->CopyStructure(first);
326 out->SetScalarType(first->GetScalarType());
327 out->SetNumberOfScalarComponents(first->GetNumberOfScalarComponents());
329 //first->GetExtent(ext); // JPR
330 first->GetWholeExtent(ext); // renvoie egalement 0,0 en Z // JPR
334 ext[5] = im.size()-1;
338 ext[5] = ext[5] * im.size()-1; // to deal with multiframes - JPR
342 // LG : TODO : Z Spacing ?
345 first->GetDimensions(dim);
347 out->SetDimensions(dim[0], dim[1], im.size() );
348 out->AllocateScalars();
351 unsigned long imsize = dim[0] * dim[1];
352 imsize = imsize * dim[2] ; // deal with multiframes // JPR
356 // differents formats char , short, etc...
357 // differents components 1..3 ex. jpg ->RGB 3
358 imsize = imsize * first->GetScalarSize() * first->GetNumberOfScalarComponents();
361 // Order the file name vector
364 first->GetSpacing(spc);
365 spc[2]=OrderTheFileNameVector(im);
366 out->SetSpacing(spc);
368 std::vector<std::string>::iterator it;
370 for (it=im.begin(); it!=im.end(); ++it)
372 vtkImageData* cur = mReader.GetImage( (*it) );
373 memcpy(out->GetScalarPointer(0,0,slice), cur->GetScalarPointer(0,0,0), imsize);
383 //======================================================================
386 double GimmickView::OrderTheFileNameVector(std::vector<std::string> &im)
389 typedef std::vector<GDCM_NAME_SPACE::File* > FileList;
391 //GDCM_NAME_SPACE::File *f = GDCM_NAME_SPACE::File::New();
392 GDCM_NAME_SPACE::SerieHelper *sh = GDCM_NAME_SPACE::SerieHelper::New();
393 std::vector<std::string> lstAux;
394 std::vector<std::string>::iterator it;
395 for (it=im.begin(); it!=im.end(); ++it)
397 ///\TODO liberer les GDCM_NAME_SPACE::File a la fin! // JPR
398 GDCM_NAME_SPACE::File *f = GDCM_NAME_SPACE::File::New();
403 fileVector.push_back(f);
405 lstAux.push_back(*it);
410 if ((fileVector.size()>1) && (sh->IsCoherent( &fileVector )))
412 sh->OrderFileList(&fileVector);
413 spacing= sh->GetZSpacing();
416 for (i=0; i<fileVector.size(); i++)
418 im.push_back( (fileVector[i])->GetFileName() );
420 for (i=0; i<lstAux.size(); i++)
422 im.push_back( lstAux[i] );
425 std::sort( im.begin(), im.end() );
432 //======================================================================
435 void GimmickView::ReadImagesNotThreadedInVector(std::vector<vtkImageData*>& s, std::vector<std::string> im, int dimension)
437 // Create the output data
440 // Only one image : give it
441 vtkImageData* out = vtkImageData::New();
442 GimmickDebugMessage(3, "State Check: Full Filename: "
445 out->ShallowCopy(mReader.GetImage(im.front()));
448 else if (im.size()>1) // Test inutile ? JPR
450 vtkImageData* first = mReader.GetImage( im.front());
454 std::vector<std::string>::iterator it;
455 for (it=im.begin(); it!=im.end(); ++it)
457 vtkImageData* out = vtkImageData::New();
458 out->ShallowCopy(mReader.GetImage(*it));
465 // n *2D + T in a vector :
467 std::vector<std::string>::iterator it;
468 for (it=im.begin(); it!=im.end(); ++it)
470 vtkImageData* out = mReader.GetImage( (*it));
476 //======================================================================
478 //======================================================================
479 ///Requests the reading of an image
480 void GimmickView::RequestReading(tree::Node* n,
481 int prio, int selection_index, boost::shared_ptr<ImagePointerHolder> p)
488 ImageEventType t(n,selection_index);
490 mImageEventMap[n->GetAttribute("FullFileName")] = t;
491 mReader.Request(this,n->GetAttribute("FullFileName"),prio);
493 //======================================================================
495 //======================================================================
497 OnMultiThreadImageReaderEvent(const std::string& filename,
498 MultiThreadImageReaderUser::EventType e,
501 GimmickDebugMessage(7,
502 "MultiThreadImageReader event : "<<e<<std::endl);
505 if (filename.size()==0)
507 //What to do in this case?
509 GimmickDebugMessage(5,
510 "Pushing unknown image in queue"
512 mImageEventQueue.push_back(ImageEventType(image));*/
515 ImageEventTypeMap::iterator i;
516 //JCP 22-06-2009, test mImageEventMap.size() > 0
517 if(mImageEventMap.size()>0){
518 i = mImageEventMap.find(filename);
519 if (i!=mImageEventMap.end())
521 GimmickDebugMessage(5,
522 "Putting image of file '"<<filename<<"' on pointer"
524 ImageEventType ie(i->second);
526 ie.pointerHolder->Set(ie.image);
527 //mImageEventMap.erase(i);
533 std::string mess="ERROR: MultiThreadImageReader: Cannot read image in file ";
536 GimmickMessage(1,mess);
537 ImageEventTypeMap::iterator i;
538 i = mImageEventMap.find(filename);
539 if (i!=mImageEventMap.end())
541 ImageEventType ie(i->second);
543 ie.pointerHolder->Set(GetDefaultImage());
544 //mImageEventMap.erase(i);
548 else if (e==ImageUnloaded)
550 std::string mess="Unloaded image in file ";
553 GimmickMessage(1,mess);
554 ImageEventTypeMap::iterator i;
555 i = mImageEventMap.find(filename);
556 if (i!=mImageEventMap.end())
558 ImageEventType ie(i->second);
560 ie.pointerHolder->Set(GetDefaultImage());
561 //mImageEventMap.erase(i);
566 //====================================================================
568 //====================================================================
569 void GimmickView::ConnectValidationObserver(ValidationCallbackType callback)
571 mValidationSignal.connect(callback);
574 } // EO namespace creaImageIO