1 #include <creaImageIOGimmickView.h>
2 #include <creaImageIOSystem.h>
3 #include "boost/filesystem.hpp"
5 #include <gdcmSerieHelper.h> // JPR
6 #include <gdcmFile.h> // JPR
7 #include <vtkGdcmReader.h> // JPR
10 namespace fs = boost::filesystem;
15 ///Class used to represent the actual state of the image selected and to perform comparisons on its values
19 ImageExtent(const std::string& x, const std::string& y, const std::string& z, const std::string& t)
21 sscanf(x.c_str(),"%d",&mExtent[0]);
22 sscanf(y.c_str(),"%d",&mExtent[1]);
23 sscanf(z.c_str(),"%d",&mExtent[2]);
24 sscanf(t.c_str(),"%d",&mExtent[3]);
25 if(x==""){mExtent[0]=1;}
26 if(y==""){mExtent[1]=1;}
27 if(z==""){mExtent[2]=1;}
28 if(t==""){mExtent[3]=1;}
30 if (mExtent[3]>1) mDim=4;
31 else if (mExtent[2]>1) mDim=3;
32 else if (mExtent[1]>1) mDim=2;
33 else if (mExtent[0]>1) mDim=1;
38 void Clear() { mExtent[0] = mExtent[1] = mExtent[2] = mExtent[3] = 1; }
40 ///Returns true if the two extents are compatible
41 bool IsCompatible( const ImageExtent& );
43 ///Adds the extent passed as a parameter to the current extent
44 void Add ( const ImageExtent& );
46 ///Returns the ieth position of the extent
47 int Get(int i) { return mExtent[i]; }
49 ///Returns the dimension of the current image
50 void SetDimension(int dim) { mDim=dim; }
52 ///Returns the dimension of the current image
53 int GetDimension() { return mDim; }
60 //======================================================================
62 //======================================================================
64 GimmickView::GimmickView(boost::shared_ptr<Gimmick> gimmick, int threads)
68 GimmickDebugMessage(1,"GimmickView::GimmickView"
70 // Anciently started the threads ...
71 // Threads now automatically start at first image request
74 //======================================================================
76 //======================================================================
78 GimmickView::~GimmickView()
80 GimmickDebugMessage(1,"GimmickView::~GimmickView"
83 //======================================================================
85 //======================================================================
86 /// Initializes the view :
87 /// Creates the TreeViews for all the TreeHandler of the Controller
89 void GimmickView::Initialize()
93 //======================================================================
95 //======================================================================
97 void GimmickView::Finalize()
101 //======================================================================
103 //======================================================================
104 /// Create the tree views
105 void GimmickView::CreateTreeViews()
107 GimmickMessage(2,"Creating the tree views"<<std::endl);
108 Gimmick::TreeHandlerMapType::const_iterator i;
109 for (i = mGimmick->GetTreeHandlerMap().begin();
110 i!= mGimmick->GetTreeHandlerMap().end();
113 this->CreateTreeView(i->second);
117 /// Create a tree view with a given name
118 void GimmickView::CreateSingleTreeView(std::string &i_name)
120 this->CreateTreeView(mGimmick->GetTreeHandlerMap()[i_name]);
124 //======================================================================
126 //======================================================================
127 /// Updates the TreeView of given name from level l to bottom
128 /// (calls the virtual method TreeView::Update())
129 void GimmickView::UpdateTreeViewLevel(const std::string& t, int l)
131 TreeViewMapType::iterator i;
132 i = GetTreeViewMap().find(t);
133 if ( i == GetTreeViewMap().end() )
135 GimmickError("INTERNAL ERROR : GimmickView::UpdateTreeView : '"
136 <<t<<"' is not in TreeViewMap");
138 i->second->UpdateLevel(l);
141 //======================================================================
142 /// Clears the status and begins a new selection process
143 void GimmickView::ResetExtent()
147 mImageExtent.reset();
152 //======================================================================
154 //======================================================================
155 bool ImageExtent::IsCompatible(const ImageExtent& ie)
157 bool compatible=true;
158 ImageExtent * extent= (ImageExtent*)&ie;
159 if((*extent).Get(0)!=Get(0)
160 || (*extent).Get(1)!=Get(1))
167 //======================================================================
169 //======================================================================
170 void ImageExtent::Add(const ImageExtent& ie)
172 ImageExtent * extent= (ImageExtent*)&ie;
173 mExtent[2]+=(*extent).Get(2);
180 //======================================================================
181 /// No selected image
182 bool GimmickView::NoValidateSelected ()
184 GimmickDebugMessage(2,"Validating selected"<<std::endl);
185 std::string mMessage;
186 mMessage="Cannot have 0 images selected!";
188 modifyValidationSignal(valid);
189 SetMessage(mMessage);
193 //======================================================================
194 ///Validates the dimension compliance of the images with the maximum and
195 ///minimum given, and between their sizes
196 bool GimmickView::ValidateSelected (tree::Node* sel, int min_dim, int max_dim)
198 GimmickDebugMessage(2,"Validating selected"<<std::endl);
199 std::string mMessage;
203 mMessage="Cannot have 0 images selected!";
208 boost::shared_ptr<ImageExtent> ie=boost::shared_ptr<ImageExtent>(new ImageExtent((*sel).GetAttribute("D0028_0010"),
209 (*sel).GetAttribute("D0028_0011"),
210 (*sel).GetAttribute("D0028_0012"),
215 if((mImageExtent->Get(min_dim-1)<2)||(mImageExtent->Get(max_dim)>1))
221 std::stringstream out;
222 out << mImageExtent->GetDimension() << "D image " << mImageExtent->Get(0) << "x"<< mImageExtent->Get(1) << "x"<< mImageExtent->Get(2) <<" selected";
223 mMessage = out.str();
224 mImageExtent->SetDimension(2);
230 if(mImageExtent->IsCompatible(*ie))
232 if(mImageExtent->GetDimension()==max_dim && mImageExtent->Get(max_dim)>2)
234 std::stringstream out;
235 out<<"Cannot add this image to selection : would result in a "<<mImageExtent->GetDimension()+1<<"D image!";
241 std::stringstream out;
242 out<<"Selecting "<<mImageExtent->GetDimension()<<"D images is not allowed !";
246 else if(min_dim==3 && (ie->Get(2)+mImageExtent->Get(2))<2)
248 std::stringstream out;
249 out << "Cannot build the selection as it would result in a ";
250 out << mImageExtent->GetDimension();
251 out << "D image, and the minimum is ";
259 mImageExtent->Add(*ie);
260 std::stringstream out;
261 out << mImageExtent->GetDimension() << "D image " << mImageExtent->Get(0) << "x"<< mImageExtent->Get(1) << "x"<< mImageExtent->Get(2) <<" selected";
262 mMessage = out.str();
267 mMessage="The selected images are not compatible.";
273 modifyValidationSignal(valid);
274 SetMessage(mMessage);
278 //======================================================================
279 void GimmickView::modifyValidationSignal(bool ivalid)
281 mValidationSignal(ivalid);
284 //======================================================================
285 ///Reads Images (Non Threaded)
286 void GimmickView::ReadImagesNotThreaded(std::vector<vtkImageData*>& s, std::vector<std::string> im, int dimension)
291 #define GIMMICK_NO_IMAGE_SELECTION 0
292 #define GIMMICK_2D_IMAGE_SELECTION 2
293 #define GIMMICK_3D_IMAGE_SELECTION 3
294 #define GIMMICK_4D_IMAGE_SELECTION 4
301 // Create the output data
304 // Only one image : give it
305 vtkImageData* out = vtkImageData::New();
306 GimmickDebugMessage(3, "State Check: Full Filename: "
309 out->ShallowCopy(mReader.GetImage(im.front()));
312 else if (im.size()>1) // Test inutile ? JPR
314 // Read the FIRST image to get info ... about all the other ones. ?!? // JPR
315 vtkImageData* first = mReader.GetImage( im.front());
320 std::vector<std::string>::iterator it;
321 for (it=im.begin(); it!=im.end(); ++it)
323 vtkImageData* out = vtkImageData::New();
324 out->ShallowCopy(mReader.GetImage(*it));
328 else if (dimension == 3)
330 // n*2D to 3D :: --> NON ! Ca va tres bien pour du 2D+T (le ZSpacing, on s'en tape)
331 // // ne convient pas pour du n*2D a lire comme un 3D
332 // // Il faut un gdcmSerieHelper !
333 vtkImageData* out = vtkImageData::New();
334 out->CopyStructure(first);
335 out->SetScalarType(first->GetScalarType());
337 //first->GetExtent(ext); // JPR
338 first->GetWholeExtent(ext); // renvoie egalement 0,0 en Z // JPR
341 std::cout <<"in GimmickView::ReadImagesNotThreaded GetWholeExtent ext =";
342 for (int jjj=0;jjj<6;jjj++)
343 std:cout << " [" << jjj << "]=" << ext[jjj];
344 std::cout << std::endl;
348 ext[5] = im.size()-1;
352 ext[5] = ext[5] * im.size()-1; // to deal with multiframes - JPR
356 // LG : TODO : Z Spacing ? // Voir + loin // JPR
358 out->AllocateScalars();
361 first->GetDimensions(dim);
362 unsigned long imsize = dim[0] * dim[1];
363 imsize = imsize * dim[2] ; // deal with multiframes // JPR
365 imsize = imsize * first->GetScalarSize();
368 std::vector<std::string>::iterator it;
369 for (it=im.begin(); it!=im.end(); ++it)
371 vtkImageData* cur = mReader.GetImage( (*it));
372 // void* src = cur->GetScalarPointer(0,0,0);
373 // void* dst = out->GetScalarPointer(0,0,slice);
374 memcpy(out->GetScalarPointer(0,0,slice), cur->GetScalarPointer(0,0,0), imsize);
381 else // ici, on traite le cas *vraiment* n*2D as 3D // JPR
383 typedef std::vector<GDCM_NAME_SPACE::File* > FileList;
385 //GDCM_NAME_SPACE::File *f = GDCM_NAME_SPACE::File::New();
386 GDCM_NAME_SPACE::SerieHelper *sh = GDCM_NAME_SPACE::SerieHelper::New();
387 std::vector<std::string>::iterator it;
388 for (it=im.begin(); it!=im.end(); ++it)
390 //std::cout << "avant Load ---------" << *it << std::endl;
392 ///\TODO liberer les GDCM_NAME_SPACE::File a la fin! // JPR
393 GDCM_NAME_SPACE::File *f = GDCM_NAME_SPACE::File::New();
398 fileVector.push_back(f);
400 std::cout << "File [" << *it << "] not GDCM Readable; Skipped!" << std::endl;
404 for (int i=0; i<fileVector.size();i++)
406 std::cout << "avant tri ---------" << fileVector[i]->GetFileName() << std::endl;
409 if (!sh->IsCoherent(&fileVector))
411 std::cout << "---------fileset is NOT coherent ?!?" << *it << std::endl;
416 sh->OrderFileList(&fileVector);
419 for (int i=0; i<fileVector.size();i++)
421 std::cout << "apres tri ---------" << fileVector[i]->GetFileName() << std::endl;
424 vtkGdcmReader *reader = vtkGdcmReader::New();
425 reader->SetCoherentFileList(&fileVector);
427 s.push_back(reader->GetOutput());
428 reader->GetOutput()->Print(std::cout);
433 //======================================================================
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 std::cout << "----------------------------------------n3D" << std::endl;
458 vtkImageData* out = vtkImageData::New();
459 out->ShallowCopy(mReader.GetImage(*it));
466 // n *2D + T in a vector :
467 std::cout << "---------------------------------------- n *2D + T" << std::endl;
468 std::vector<std::string>::iterator it;
469 for (it=im.begin(); it!=im.end(); ++it)
471 std::cout << "---------------------------------------- [" << *it << "]" << std::endl;
472 //vtkImageData* out = mReader.GetImage( (*it));
474 vtkImageData* out = vtkImageData::New();
475 out->ShallowCopy(mReader.GetImage(*it));
481 //======================================================================
483 //======================================================================
484 ///Requests the reading of an image
485 void GimmickView::RequestReading(tree::Node* n,
486 int prio, int selection_index, boost::shared_ptr<ImagePointerHolder> p)
493 ImageEventType t(n,selection_index);
495 mImageEventMap[n->GetAttribute("FullFileName")] = t;
496 mReader.Request(this,n->GetAttribute("FullFileName"),prio);
498 //======================================================================
500 //======================================================================
502 OnMultiThreadImageReaderEvent(const std::string& filename,
503 MultiThreadImageReaderUser::EventType e,
506 GimmickDebugMessage(7,
507 "MultiThreadImageReader event : "<<e<<std::endl);
510 if (filename.size()==0)
512 //What to do in this case?
514 GimmickDebugMessage(5,
515 "Pushing unknown image in queue"
517 mImageEventQueue.push_back(ImageEventType(image));*/
520 ImageEventTypeMap::iterator i;
521 //JCP 22-06-2009, test mImageEventMap.size() > 0
522 if(mImageEventMap.size()>0){
523 i = mImageEventMap.find(filename);
524 if (i!=mImageEventMap.end())
526 GimmickDebugMessage(5,
527 "Putting image of file '"<<filename<<"' on pointer"
529 ImageEventType ie(i->second);
531 ie.pointerHolder->Set(ie.image);
532 //mImageEventMap.erase(i);
538 std::string mess="ERROR: MultiThreadImageReader: Cannot read image in file ";
541 GimmickMessage(1,mess);
542 ImageEventTypeMap::iterator i;
543 i = mImageEventMap.find(filename);
544 if (i!=mImageEventMap.end())
546 ImageEventType ie(i->second);
548 ie.pointerHolder->Set(GetDefaultImage());
549 //mImageEventMap.erase(i);
553 else if (e==ImageUnloaded)
555 std::string mess="Unloaded image in file ";
558 GimmickMessage(1,mess);
559 ImageEventTypeMap::iterator i;
560 i = mImageEventMap.find(filename);
561 if (i!=mImageEventMap.end())
563 ImageEventType ie(i->second);
565 ie.pointerHolder->Set(GetDefaultImage());
566 //mImageEventMap.erase(i);
571 //====================================================================
573 //====================================================================
574 void GimmickView::ConnectValidationObserver(ValidationCallbackType callback)
576 mValidationSignal.connect(callback);
579 } // EO namespace creaImageIO