2 # ---------------------------------------------------------------------
4 # Copyright (c) CREATIS (Centre de Recherche en Acquisition et Traitement de l'Image
6 # Authors : Eduardo Davila, Frederic Cervenansky, Claire Mouton
7 # Previous Authors : Laurent Guigues, Jean-Pierre Roux
8 # CreaTools website : www.creatis.insa-lyon.fr/site/fr/creatools_accueil
10 # This software is governed by the CeCILL-B license under French law and
11 # abiding by the rules of distribution of free software. You can use,
12 # modify and/ or redistribute the software under the terms of the CeCILL-B
13 # license as circulated by CEA, CNRS and INRIA at the following URL
14 # http://www.cecill.info/licences/Licence_CeCILL-B_V1-en.html
15 # or in the file LICENSE.txt.
17 # As a counterpart to the access to the source code and rights to copy,
18 # modify and redistribute granted by the license, users are provided only
19 # with a limited warranty and the software's author, the holder of the
20 # economic rights, and the successive licensors have only limited
23 # The fact that you are presently reading this means that you have had
24 # knowledge of the CeCILL-B license and that you accept its terms.
25 # ------------------------------------------------------------------------
29 #include <creaImageIOMultiThreadImageReader.h>
30 #include <creaImageIOImageReader.h>
32 #include <creaImageIOSystem.h>
34 #include <creaImageIODicomImageReader.h>
37 #include <creaImageIOGimmick.h>
44 //=====================================================================
45 void MultiThreadImageReaderUser::MultiThreadImageReaderSendEvent
46 ( const std::string& filename,
50 wxMutexLocker lock(mMultiThreadImageReaderUserMutex);
51 this->OnMultiThreadImageReaderEvent(filename,type,image);
53 //=====================================================================
58 void Delete() { printf("EED wxThreadEED2::Delete() \n"); }
59 int Run() { printf("EED wxThreadEED2::Run() \n"); return 0;}
60 void Pause() { printf("EED wxThreadEED2::Pause() \n"); }
61 void Create() { printf("EED wxThreadEED2::Create() \n"); }
62 bool IsAlive() { printf("EED wxThreadEED2::IsAlive() \n"); return false; }
63 bool TestDestroy() { printf("EED wxThreadEED2::TestDestroy() \n"); return false; }
64 int GetCurrentId() { printf("EED wxThreadEED2::GetCurrentId() \n"); return -999; }
67 //=====================================================================
69 // class ThreadedImageReader: public wxThread
70 class ThreadedImageReader: public wxThreadEED2
73 ThreadedImageReader(MultiThreadImageReader* tir) :
74 mMultiThreadImageReader(tir)
78 vtkImageData* Read(const std::string& filename);
81 void operator()(ThreadedImageReader* p)
86 friend struct deleter;
90 MultiThreadImageReader* mMultiThreadImageReader;
93 //=====================================================================
96 //=====================================================================
97 MultiThreadImageReader::MultiThreadImageReader(int number_of_threads)
98 : //mDoNotSignal(false),
101 mTotalMemMax(1000000)
103 // std::cout << "#### MultiThreadImageReader::MultiThreadImageReader("
104 // << " #threads= " << number_of_threads <<" )"<<std::endl;
107 // Create the threads
109 for (int i=0; i<number_of_threads; i++)
111 //ThreadedImageReader* t = new ThreadedImageReader(this);
112 boost::shared_ptr<ThreadedImageReader> t(new ThreadedImageReader(this), ThreadedImageReader::deleter());
113 mThreadedImageReaderList.push_back(t);
114 std::cout << " ===> Thread "<<i <<" successfully added"<< std::endl;
116 mNumberOfThreadedReadersRunning = 0;
118 mQueue.set(mComparator);
119 mQueue.set(mIndexer);
121 // no thread : alloc self reader
122 // if (number_of_threads==0)
124 mReader = new ImageReader();
127 //=====================================================================
130 //=====================================================================
131 bool MultiThreadImageReader::Start()
133 // std::cout << "#### MultiThreadImageReader::Start()"
135 if (mNumberOfThreadedReadersRunning > 0) return true;
137 ThreadedImageReaderListType::iterator i;
138 for (i =mThreadedImageReaderList.begin();
139 i!=mThreadedImageReaderList.end();
143 if ( (*i)->Run() != wxTHREAD_NO_ERROR )
145 std::cout << "ERROR starting a thread"<< std::endl;
148 std::cout << " ===> Thread "<<(*i)->GetCurrentId()
149 <<" successfully created"<< std::endl;
152 wxMutexLocker locker(GetMultiThreadImageReaderUserMutex());
153 // std::cout << "EO Start : #Threads running = "
154 // << mNumberOfThreadedReadersRunning<<std::endl;
157 //=====================================================================
159 //=====================================================================
160 void MultiThreadImageReader::Stop()
162 // std::cout << "#### MultiThreadImageReader::Stop()"
164 // std::cout << "Sending stop order to the threads..."<<std::endl;
168 ThreadedImageReaderListType::iterator i;
169 for (i =mThreadedImageReaderList.begin();
170 i!=mThreadedImageReaderList.end();
173 std::cout << " ===> Thread "<<(*i)->GetCurrentId()
174 <<" successfully stopped"<< std::endl;
182 mThreadedImageReaderList.clear();
183 // Wait a little to be sure that all threads have stopped
184 // A better way to do this ?
185 // wxMilliSleep(1000);
186 // New method : the threads generate a stop event when they have finished
187 // We wait until all threads have stopped
188 // std::cout << "Waiting for stop signals..."<<std::endl;
195 wxMutexLocker locker(GetMultiThreadImageReaderUserMutex());
196 // std::cout << "#Threads running = "
197 // << mNumberOfThreadedReadersRunning<<std::endl;
198 // Break if all readers have stopped
199 if (mNumberOfThreadedReadersRunning <= 0)
206 // std::cout << "All threads stopped : OK "<<std::endl;
214 //=====================================================================
217 void MultiThreadImageReader::CleanMImagesMap()
219 ImageMapType::iterator j;
220 for (j=mImages.begin(); j!=mImages.end(); ++j)
227 //=====================================================================
228 MultiThreadImageReader::~MultiThreadImageReader()
230 // std::cout << "#### MultiThreadImageReader::~MultiThreadImageReader()"
233 if (mReader) delete mReader;
234 mThreadedImageReaderList.clear();
239 ImageMapType::iterator it;
240 for (it=mImages.begin() ; it!=mImages.end(); it++)
242 printf("MultiThreadImageReader::~MultiThreadImageReader %s ", it.first->GetFilename() );
247 //=====================================================================
249 //=====================================================================
250 void MultiThreadImageReader::UpdateUnloadPriority(ImageToLoadPtr p,
253 // not in unload queue : ciao
254 if (p->UnloadIndex()<0) return;
255 int old_prio = p->GetPriority();
256 if (priority > old_prio)
258 p->SetPriority(priority);
259 mUnloadQueue.downsort(p->UnloadIndex());
261 else if ( old_prio > priority )
263 p->SetPriority(priority);
264 mUnloadQueue.upsort(p->UnloadIndex());
267 //=====================================================================
268 // function to read attributes for a file
269 void MultiThreadImageReader::getAttributes(const std::string filename,
270 std::map <std::string , std::string> &infos,std::vector<std::string> i_attr)
272 mReader->getAttributes(filename, infos, i_attr);
275 //=====================================================================
276 void MultiThreadImageReader::Request( MultiThreadImageReaderUser* user,
277 const std::string& filename,
280 wxMutexLocker lock(GetMultiThreadImageReaderUserMutex()); //mMutex);
282 if (mNumberOfThreadedReadersRunning==0)
283 // if (mThreadedImageReaderList.size()==0)
285 // no detached reader : use self reader
286 ImageToLoad itl(user,filename);
287 ImageMapType::iterator i = mImages.find(&itl);
288 if (i!=mImages.end())
290 ImageToLoadPtr pitl = const_cast<ImageToLoadPtr>(i->first);
292 if (pitl->GetImage() != 0)
296 UpdateUnloadPriority(pitl,priority);
297 SignalImageRead(pitl,false);
298 return; // pitl->GetImage();
301 ImageToLoadPtr pitl = new ImageToLoad(user,filename,0);
303 pitl->SetImage(mReader->ReadImage(filename));
304 UpdateUnloadPriority(pitl,priority);
305 SignalImageRead(pitl,true);
306 // return pitl->GetImage();
310 ImageToLoad itl(user,filename);
311 ImageMapType::iterator i = mImages.find(&itl);
312 if (i!=mImages.end())
315 if (i->first->GetImage() != 0)
317 // Already read : ok :signal the user
318 UpdateUnloadPriority(i->first,priority);
319 SignalImageRead(i->first,false);
322 /// Already requested : change the priority
323 ImageToLoadPtr pitl = const_cast<ImageToLoadPtr>(i->first);
324 pitl->SetPriority(priority);
326 if (pitl->Index()>=0)
329 mQueue.upsort(pitl->Index());
331 // Not read but not in queue = being read = ok
339 // Never requested before or unloaded
340 ImageToLoadPtr pitl = new ImageToLoad(user,filename,priority);
345 //=====================================================================
347 //=====================================================================
348 void MultiThreadImageReader::OnMultiThreadImageReaderEvent
349 (const std::string& filename,
350 MultiThreadImageReaderUser::EventType e,
353 if ((e==MultiThreadImageReaderUser::ImageLoaded) &&
354 (filename == mRequestedFilename))
356 mRequestedImage = image;
358 else if (e==MultiThreadImageReaderUser::ThreadedReaderStarted)
360 mNumberOfThreadedReadersRunning++;
361 // std::cout << "#TR=" << mNumberOfThreadedReadersRunning << std::endl;
363 else if (e==MultiThreadImageReaderUser::ThreadedReaderStopped)
366 mNumberOfThreadedReadersRunning--;
367 // std::cout << "#TR=" << mNumberOfThreadedReadersRunning << std::endl;
370 //=====================================================================
372 //=====================================================================
373 vtkImageData* MultiThreadImageReader::GetImage(const std::string& filename)
376 // std::cout << "** MultiThreadImageReader::GetImage('"<<filename<<"')"
381 // wxMutexLocker lock(GetMultiThreadImageReaderUserMutex()); //mMutex);
383 // std::cout << "** MultiThreadImageReader::GetImage('"<<filename
387 // if (mNumberOfThreadedReadersRunning==0)
388 // if (mThreadedImageReaderList.size()==0)
391 ImageToLoad itl(this,filename);
392 ImageMapType::iterator i = mImages.find(&itl);
393 if (i!=mImages.end())
395 ImageToLoadPtr pitl = const_cast<ImageToLoadPtr>(i->first);
397 if (pitl->GetImage() != NULL)
400 UpdateUnloadPriority( pitl, GetMaximalPriorityWithoutLocking()+1 );
401 return pitl->GetImage();
402 } // if pitl->GetImage()
405 ImageToLoadPtr pitl = new ImageToLoad(this,filename,0);
406 mImages[pitl] = NULL;
407 pitl->SetImage( mReader->ReadImage(filename) );
408 UpdateUnloadPriority( pitl, GetMaximalPriorityWithoutLocking()+1 );
409 return pitl->GetImage();
413 mRequestedFilename = filename;
415 ImageToLoad itl(this,filename);
416 ImageMapType::iterator i = mImages.find(&itl);
417 if (i!=mImages.end())
419 // Already inserted in queue
420 if (i->first->GetImage() != 0)
422 // Already read : ok : return it
423 return i->first->GetImage();
425 /// Already requested : change the priority
426 ImageToLoadPtr pitl = const_cast<ImageToLoadPtr>(i->first);
427 pitl->SetPriority( GetMaximalPriorityWithoutLocking() + 1 );
428 pitl->SetUser( this );
430 if (pitl->Index()>=0)
433 mQueue.upsort(pitl->Index());
435 // Not read but not in queue = being read = ok
438 pitl->SetUser( this );
444 // Never requested before or unloaded
445 ImageToLoadPtr pitl =
446 new ImageToLoad(this,filename,
447 GetMaximalPriorityWithoutLocking() + 1);
455 // std::cout << "Waiting..."<<std::endl;
458 // Waiting that it is read
462 // std::cout << n++ << std::endl;
466 // wxMutexLocker lock(mMutex);
467 wxMutexLocker lock(GetMultiThreadImageReaderUserMutex());
468 if (mRequestedImage!=0)
470 return mRequestedImage;
478 printf("EED MultiThreadImageReader::GetImage END\n");
481 //=====================================================================
483 //=====================================================================
484 void MultiThreadImageReader::SignalImageRead(ImageToLoadPtr p,
488 // std::cout << "MultiThreadImageReader::SignalImageRead" <<std::endl;
489 // std::cout << "this="<<this <<std::endl;
490 // std::cout << "user="<<p->GetUser() <<std::endl;
492 if ( p->GetUser() == this )
494 GetMultiThreadImageReaderUserMutex().Unlock();
497 p->GetUser()->MultiThreadImageReaderSendEvent
499 MultiThreadImageReaderUser::ImageLoaded,
503 AN ATTEMPT TO UNLOAD OLDEST IMAGE IF EXCEEDED A CERTAIN MEMORY QUOTA
507 GimmickMessage(5,"Image '"<<p->GetFilename()<<"' read"<<std::endl);
509 // wxMutexLocker lock(GetMultiThreadImageReaderUserMutex());
511 mUnloadQueue.insert(p);
514 //EED 2017-01-01 Migration VTK7
515 #if VTK_MAJOR_VERSION <= 5
516 p->GetImage()->UpdateInformation();
517 p->GetImage()->PropagateUpdateExtent();
518 long ImMem = p->GetImage()->GetEstimatedMemorySize();
522 p->GetImage()->GetExtent(ext);
523 dim[0] = ext[1]-ext[0]+1;
524 dim[1] = ext[3]-ext[2]+1;
525 dim[2] = ext[5]-ext[4]+1;
526 long ImMem = dim[0]*dim[1]*dim[2]*p->GetImage()->GetScalarSize();;
530 GimmickMessage(5,"==> Image in memory = "<<mUnloadQueue.size()<<std::endl);
531 GimmickMessage(5,"==> Total mem = "<<mTotalMem<<" Ko"<<std::endl);
535 while (mTotalMem > mTotalMemMax)
538 " ! Exceeded max of "
539 << mTotalMemMax << " Ko : unloading oldest image ... "
541 if ( mUnloadQueue.size() <= 1 )
544 " Only one image : cannot load AND unload it !!"
549 ImageToLoadPtr unload = mUnloadQueue.remove_top();
550 MultiThreadImageReaderUser* user = unload->GetUser();
553 if ((user!=0)&&(user!=this))
555 user->GetMultiThreadImageReaderUserMutex().Lock();
559 std::string filename = unload->GetFilename();
561 GimmickMessage(5,"'" << filename << "'" << std::endl);
563 //EED 2017-01-01 Migration VTK7
564 #if VTK_MAJOR_VERSION <= 5
565 mTotalMem -= unload->GetImage()->GetEstimatedMemorySize();
569 unload->GetImage()->GetExtent(ext);
570 dim[0] = ext[1]-ext[0]+1;
571 dim[1] = ext[3]-ext[2]+1;
572 dim[2] = ext[5]-ext[4]+1;
573 mTotalMem -= dim[0]*dim[1]*dim[2]*unload->GetImage()->GetScalarSize();
576 GimmickMessage(5," ==> Total mem = "<<mTotalMem<<" Ko "<<std::endl);
580 // std::cout << "unlock..."<<std::endl;
581 // user->GetMultiThreadImageReaderUserMutex().Unlock();
582 // std::cout << "event"<<std::endl;
583 user->MultiThreadImageReaderSendEvent
585 MultiThreadImageReaderUser::ImageUnloaded,
587 // std::cout << "event ok"<<std::endl;
590 if (unload->Index()>=0)
592 // GimmickMessage(5,"still in queue"<<std::endl);
594 unload->Index() = -1;
597 ImageMapType::iterator it = mImages.find(unload);
598 if (it!=mImages.end())
602 // std::cout << "delete..."<<std::endl;
604 // std::cout << "delete ok."<<std::endl;
608 //=====================================================================
610 //=====================================================================
611 int MultiThreadImageReader::GetMaximalPriority()
613 wxMutexLocker lock(GetMultiThreadImageReaderUserMutex()); //mMutex);
614 return GetMaximalPriorityWithoutLocking();
616 //=====================================================================
619 //=====================================================================
620 int MultiThreadImageReader::GetMaximalPriorityWithoutLocking()
625 max = mQueue.top()->GetPriority();
627 if (mUnloadQueue.size()>0)
629 int max2 = mUnloadQueue.top()->GetPriority();
630 if (max2>max) max=max2;
634 //=====================================================================
637 //=====================================================================
638 //=====================================================================
639 //=====================================================================
640 //=====================================================================
642 //=====================================================================
643 void* ThreadedImageReader::Entry()
645 // std::cout << "### Thread "<<GetCurrentId()<<"::Entry()"
648 mMultiThreadImageReader->MultiThreadImageReaderSendEvent
650 MultiThreadImageReaderUser::ThreadedReaderStarted,
653 // While was not deleted
654 while (!TestDestroy())
656 //std::cout << "### Thread "<<GetCurrentId()<<" still alive" << std::endl;
659 mMultiThreadImageReader->MultiThreadImageReaderEventLock();
662 if (mMultiThreadImageReader->mQueue.size()>0)
664 MultiThreadImageReader::ImageToLoadPtr i =
665 mMultiThreadImageReader->mQueue.remove_top();
667 mMultiThreadImageReader->MultiThreadImageReaderEventUnlock();
671 // std::cout << "### Thread "<<GetCurrentId()<<" : reading '"
672 // << i->GetFilename() << "'" << std::endl;
675 vtkImageData* im = Read(i->GetFilename());
677 // Store it in the map
678 mMultiThreadImageReader->MultiThreadImageReaderEventLock();
680 MultiThreadImageReader::ImageToLoad itl(0,i->GetFilename());
681 MultiThreadImageReader::ImageMapType::iterator it =
682 mMultiThreadImageReader->mImages.find(&itl);
683 MultiThreadImageReader::ImageToLoadPtr
684 pitl = const_cast<MultiThreadImageReader::ImageToLoadPtr>
687 mMultiThreadImageReader->SignalImageRead(pitl,true);//i->GetFilename());
688 mMultiThreadImageReader->MultiThreadImageReaderEventUnlock(); //mMutex.Unlock();
690 // std::cout << "### Thread "<<GetCurrentId()<<" : reading '"
691 // << i->GetFilename() << "' : DONE" << std::endl;
694 mMultiThreadImageReader->MultiThreadImageReaderEventUnlock();
696 // Wait a little to avoid blocking
700 // std::cout << "### Thread "<<GetCurrentId()<<" stopping"
705 //=====================================================================
707 //=====================================================================
708 void ThreadedImageReader::OnExit()
710 mMultiThreadImageReader->MultiThreadImageReaderSendEvent
712 MultiThreadImageReaderUser::ThreadedReaderStopped,
715 //=====================================================================
717 //=====================================================================
718 vtkImageData* ThreadedImageReader::Read(const std::string& filename)
720 return mReader.ReadImage(filename);
722 //=====================================================================
724 } // namespace creaImageIO