continued the fusion sequence visualization mode

[clitk.git] / vv / vvSlicerManager.cxx

diff --git a/vv/vvSlicerManager.cxx b/vv/vvSlicerManager.cxx
index 6f5c55abff3cbb809394638bbaf937c901f021c3..01821d9baca187961efc7556492c3b39f7b1278a 100644 (file)
--- a/vv/vvSlicerManager.cxx
+++ b/vv/vvSlicerManager.cxx
@@ -40,6 +40,7 @@
 #include <vtkCamera.h>\r
 \r
 #include <qfileinfo.h>\r
+#include <QMessageBox>\r
 //----------------------------------------------------------------------------\r
 vvSlicerManager::vvSlicerManager(int numberOfSlicers)\r
 {\r
@@ -62,10 +63,12 @@ vvSlicerManager::vvSlicerManager(int numberOfSlicers)
   mFusionLevel = 1000;\r
   mFusionShowLegend = true;\r
   \r
+  mFusionSequenceInvolvementCode = -1;\r
+  mFusionSequenceIndexLinkedManager = -1;\r
   mFusionSequenceFrameIndex = -1;\r
-  mFusionSequenceSpatialSyncFlag = false;\r
   mFusionSequenceNbFrames = 0;\r
-  mFusionSequenceIndexLinkedManager = -1;\r
+  mFusionSequenceSpatialSyncFlag = false;\r
+  mFusionSequenceTemporalSyncFlag = false;\r
 \r
   mLandmarks = NULL;\r
   mLinkedId.resize(0);\r
@@ -307,8 +310,11 @@ bool vvSlicerManager::SetFusion(std::string filename,int dim, std::string compon
 //----------------------------------------------------------------------------\r
 \r
 //----------------------------------------------------------------------------\r
-bool vvSlicerManager::SetFusionSequence(std::vector<std::string> filenames,int dim, std::string component, vvImageReader::LoadedImageType type)\r
+//this function is called by vvMainWindow::AddFusionSequence for the primary sequence (CT), while the given files constitute the secondary sequence.\r
+bool vvSlicerManager::SetFusionSequence(std::vector<std::string> filenames, int dim, std::string component, vvImageReader::LoadedImageType type)\r
 {\r
+       mFusionSequenceInvolvementCode = 0;\r
+\r
        mFusionName = filenames[0];\r
        mFusionComponent = component;\r
 \r
@@ -326,7 +332,7 @@ bool vvSlicerManager::SetFusionSequence(std::vector<std::string> filenames,int d
 \r
        if (mFusionSequenceReader->GetLastError().size() == 0) {\r
                for ( unsigned int i = 0; i < mSlicers.size(); i++) {\r
-                       mSlicers[i]->SetFusion(mFusionSequenceReader->GetOutput(), true);\r
+                       mSlicers[i]->SetFusion(mFusionSequenceReader->GetOutput(), mFusionSequenceInvolvementCode);\r
                }\r
        } else {\r
                mLastError = mFusionSequenceReader->GetLastError();\r
@@ -411,8 +417,9 @@ vvSlicer* vvSlicerManager::GetSlicer(int i)
 //----------------------------------------------------------------------------\r
 void vvSlicerManager::UpdateSlicer(int num, bool state)\r
 {\r
-  if (mSlicers[num]->GetImage())\r
+  if (mSlicers[num]->GetImage()) {\r
     mSlicers[num]->SetDisplayMode(state);\r
+  }\r
 }\r
 //----------------------------------------------------------------------------\r
 \r
@@ -777,8 +784,9 @@ void vvSlicerManager::UpdateLinked(int slicer)
       z >= mSlicers[slicer]->GetInput()->GetWholeExtent()[4]-0.5 &&\r
       z <= mSlicers[slicer]->GetInput()->GetWholeExtent()[5]+0.5) {\r
     for (std::list<std::string>::const_iterator i = mLinkedId.begin(); i != mLinkedId.end(); i++) {\r
-               if (mFusionSequenceIndexLinkedManager>0) {\r
+               if (this->IsInvolvedInFusionSequence()) {\r
                        //this SlicerManager is involved in fusionSequence => do not synchronize the times\r
+      //TODO: check is something more specific should be done ...\r
                        emit UpdateLinkManager(*i, slicer, p[0], p[1], p[2], -1);\r
                }\r
                else {\r
@@ -935,8 +943,8 @@ void vvSlicerManager::Reload()
     mSlicers[i]->SetImage(mImage);\r
   }\r
   \r
-  //check if this image is involved in a fusion sequence, then the main transform matrix should be updated.\r
-  if (mFusionSequenceReader.IsNotNull()) {\r
+  //if this image is the primary sequence of a fusion sequence, then the main transform matrix should be updated.\r
+  if (this->IsMainSequenceOfFusionSequence()) {\r
          SetFusionSequenceMainTransformMatrix( mImage->GetTransform()[0]->GetMatrix() );\r
   }\r
 }\r
@@ -956,12 +964,13 @@ void vvSlicerManager::ReloadFusion()
 //----------------------------------------------------------------------------\r
 \r
 //----------------------------------------------------------------------------\r
+//the secondary sequence is being reloaded.\r
 void vvSlicerManager::ReloadFusionSequence()\r
 {\r
   mFusionSequenceReader->Update(mImage->GetNumberOfDimensions(),mFusionComponent.c_str(),vvImageReader::MERGEDWITHTIME);\r
 \r
   for ( unsigned int i = 0; i < mSlicers.size(); i++) {\r
-    mSlicers[i]->SetFusion(mFusionSequenceReader->GetOutput(), true);\r
+    mSlicers[i]->SetFusion(mFusionSequenceReader->GetOutput(), 1);\r
     mSlicers[i]->Render();\r
   }\r
 \r
@@ -972,7 +981,6 @@ void vvSlicerManager::ReloadFusionSequence()
   }\r
 \r
   //Update the list of initial transforms\r
-  //Warning, the main transform will not be updated on reload.........\r
   mFusionSequenceListInitialTransformMatrices.clear();\r
   for (unsigned i=0 ; i<mFusionSequenceNbFrames ; i++) {\r
          this->AddFusionSequenceInitialTransformMatrices( mFusionSequenceReader->GetOutput()->GetTransform()[i]->GetMatrix() );\r
@@ -1113,7 +1121,9 @@ void vvSlicerManager::UpdateInfoOnCursorPosition(int slicer)
                double Zover = (z - fusion->GetOrigin()[2]) / fusion->GetSpacing()[2];\r
                valueFus = this->GetScalarComponentAsDouble(fusion, Xover, Yover, Zover);\r
        }\r
-       else if (mFusionSequenceIndexLinkedManager>=0) {\r
+       else if (this->IsInvolvedInFusionSequence()) { \r
+    //if the cursor moves over the 'independent' version of the secondary sequence\r
+    //do not update the panel, just keep it as it is.\r
                displayFus = 1;\r
                valueFus = std::numeric_limits<double>::quiet_NaN();\r
        }\r
@@ -1171,9 +1181,9 @@ void vvSlicerManager::UpdateSlice(int slicer)
 void vvSlicerManager::UpdateTSlice(int slicer)\r
 {\r
   int slice = mSlicers[slicer]->GetSlice();\r
+\r
   int tslice = mSlicers[slicer]->GetMaxCurrentTSlice();\r
-  \r
-  if (mFusionSequenceIndexLinkedManager>=0) tslice = mSlicers[slicer]->GetTSlice();\r
+  if (this->IsInvolvedInFusionSequence()) tslice = mSlicers[slicer]->GetTSlice();\r
 \r
   if (mPreviousSlice[slicer] == slice) {\r
     if (mPreviousTSlice[slicer] == tslice) {\r
@@ -1184,8 +1194,6 @@ void vvSlicerManager::UpdateTSlice(int slicer)
   mPreviousSlice[slicer] = slice;\r
   mPreviousTSlice[slicer] = tslice;\r
 \r
-  if (mFusionSequenceIndexLinkedManager>=0) return;\r
-\r
   emit UpdateTSlice(slicer, tslice);\r
 }\r
 //----------------------------------------------------------------------------\r
@@ -1392,18 +1400,22 @@ void vvSlicerManager::SetColorMap(int colormap)
     LUT->Build();\r
   }\r
   vtkWindowLevelLookupTable* fusLUT = NULL;\r
+\r
+  //FUSION / FUSION SEQUENCE\r
   if (mSlicers[0]->GetFusion()) { // && mFusionColorMap >= 0) {\r
-       fusLUT = vtkWindowLevelLookupTable::New();\r
+    fusLUT = vtkWindowLevelLookupTable::New();\r
     double fusRange [2];\r
     fusRange[0] = mFusionLevel - mFusionWindow/2;\r
     fusRange[1] = mFusionLevel + mFusionWindow/2;\r
 \r
-       //check whether it is actually a fusionSequence or a fusion, before invoking mFusionReader...\r
-       double* frange;\r
-       if (mFusionReader.IsNull()) frange = mFusionSequenceReader->GetOutput()->GetVTKImages()[0]->GetScalarRange();\r
-    else                        frange = mFusionReader->GetOutput()->GetVTKImages()[0]->GetScalarRange();\r
+    //check whether it is actually a fusionSequence or a fusion, before invoking mFusionReader...\r
+    double* frange;\r
+    if (this->IsInvolvedInFusionSequence()) \r
+      frange = mFusionSequenceReader->GetOutput()->GetVTKImages()[0]->GetScalarRange();\r
+    else\r
+      frange = mFusionReader->GetOutput()->GetVTKImages()[0]->GetScalarRange();\r
 \r
-       fusLUT->SetTableRange(frange);\r
+    fusLUT->SetTableRange(frange);\r
     fusLUT->SetValueRange(1,1);\r
     fusLUT->SetSaturationRange(1,1);\r
     fusLUT->SetAlphaRange(1, 1);\r
@@ -1425,10 +1437,10 @@ void vvSlicerManager::SetColorMap(int colormap)
       fusLUT->SetValueRange(0,1);\r
       fusLUT->SetSaturationRange(0,0);\r
     }\r
-    \r
+\r
     fusLUT->ForceBuild();\r
     double v[4];\r
-    \r
+\r
     // set color table transparency\r
     //double alpha_range=(double)mFusionThresOpacity/10;\r
     double range_end = fusRange[0] + (double)mFusionThresOpacity*(fusRange[1] - fusRange[0])/100;\r
@@ -1444,7 +1456,7 @@ void vvSlicerManager::SetColorMap(int colormap)
     }\r
   }\r
   for ( unsigned int i = 0; i < mSlicers.size(); i++) {\r
-    \r
+\r
     if (mSlicers[i]->GetOverlay()) {\r
       vtkLookupTable* supLUT = vtkLookupTable::New();\r
       supLUT->SetTableRange(range[0],range[1]);\r
@@ -1467,7 +1479,7 @@ void vvSlicerManager::SetColorMap(int colormap)
     } else {\r
       mSlicers[i]->GetWindowLevel()->SetLookupTable(LUT);\r
     }\r
-    \r
+\r
     if (mSlicers[i]->GetFusion()) {\r
       mSlicers[i]->ShowFusionLegend(mFusionShowLegend);\r
       mSlicers[i]->GetFusionMapper()->SetLookupTable(fusLUT);\r