1 /*=========================================================================
5 Author : Pierre Seroul (pierre.seroul@gmail.com)
8 Léon Bérard cancer center http://oncora1.lyon.fnclcc.fr
9 CREATIS-LRMN http://www.creatis.insa-lyon.fr
11 This program is free software: you can redistribute it and/or modify
12 it under the terms of the GNU General Public License as published by
13 the Free Software Foundation, version 3 of the License.
15 This program is distributed in the hope that it will be useful,
16 but WITHOUT ANY WARRANTY; without even the implied warranty of
17 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
18 GNU General Public License for more details.
20 You should have received a copy of the GNU General Public License
21 along with this program. If not, see <http://www.gnu.org/licenses/>.
23 =========================================================================*/
24 #include "vvSlicerManager.h"
28 #include "vvSlicerManagerCommand.h"
29 #include "vvInteractorStyleNavigator.h"
30 #include "vvLandmarks.h"
31 #include "vvImageReader.h"
32 #include "vvImageReader.h"
34 #include "vvImageMapToWLColors.h"
36 #include "vtkImageActor.h"
37 #include "vtkImageData.h"
38 #include "vtkRenderWindow.h"
39 #include "vtkRendererCollection.h"
40 #include "vtkRenderWindowInteractor.h"
41 #include "vtkImageMapToWindowLevelColors.h"
42 #include "vtkWindowLevelLookupTable.h"
43 #include "vtkColorTransferFunction.h"
44 #include "vtkImageClip.h"
45 #include <vtkLODActor.h>
46 #include <vtkPointData.h>
48 #include <vtksys/SystemTools.hxx>
49 //----------------------------------------------------------------------------
51 vvSlicerManager::vvSlicerManager(int numberOfSlicers)
60 mType = UNDEFINEDIMAGETYPE;
74 mOverlayReader = NULL;
79 for ( int i = 0; i < numberOfSlicers; i++)
81 vvSlicer *slicer = vvSlicer::New();
82 mSlicers.push_back(slicer);
86 vvSlicerManager::~vvSlicerManager()
88 for ( unsigned int i = 0; i < mSlicers.size(); i++)
90 if (mSlicers[i] != NULL)
91 mSlicers[i]->Delete();
103 delete mOverlayReader;
107 delete mFusionReader;
113 void vvSlicerManager::AddContour(vvMesh::Pointer contour,bool propagate)
115 for ( unsigned int i = 0; i < mSlicers.size(); i++)
117 mSlicers[i]->AddContour(contour,propagate);
121 void vvSlicerManager::ToggleContourSuperposition()
123 for ( unsigned int i = 0; i < mSlicers.size(); i++)
124 mSlicers[i]->ToggleContourSuperposition();
127 bool vvSlicerManager::SetImage(std::string filename,LoadedImageType type)
129 mFileName = filename;
132 mReader = new vvImageReader;
133 std::vector<std::string> filenames;
134 filenames.push_back(filename);
135 mReader->SetInputFilenames(filenames);
136 mReader->Update(type);
137 if (mReader->GetLastError().size() == 0)
139 mImage=mReader->GetOutput();
140 for ( unsigned int i = 0; i < mSlicers.size(); i++)
142 mSlicers[i]->SetFileName(vtksys::SystemTools::GetFilenameWithoutLastExtension(filename));
143 mSlicers[i]->SetImage(mReader->GetOutput());
148 mLastError = mReader->GetLastError();
154 void vvSlicerManager::SetImage(vvImage::Pointer image)
157 for (unsigned int i = 0; i < mSlicers.size();i++)
159 mSlicers[i]->SetImage(image);
163 bool vvSlicerManager::SetImages(std::vector<std::string> filenames,LoadedImageType type)
166 std::string fileWithoutExtension = vtksys::SystemTools::GetFilenameWithoutExtension(filenames[0]);
168 fileWithoutExtension += "_dicom";
169 else if (type == MERGED)
170 fileWithoutExtension += "_merged";
171 else if (type == MERGEDWITHTIME)
172 fileWithoutExtension += "_merged_wt";
174 mFileName = fileWithoutExtension + vtksys::SystemTools::GetFilenameExtension(filenames[0]);
176 mReader = new vvImageReader;
177 mReader->SetInputFilenames(filenames);
178 mReader->Update(type);
181 if (mReader->GetLastError().size() == 0)
183 mImage=mReader->GetOutput();
184 for ( unsigned int i = 0; i < mSlicers.size(); i++)
186 mSlicers[i]->SetFileName(fileWithoutExtension);
187 mSlicers[i]->SetImage(mReader->GetOutput());
192 mLastError = mReader->GetLastError();
198 bool vvSlicerManager::SetOverlay(std::string filename,int dim, std::string component)
200 mOverlayName = filename;
201 mOverlayComponent = component;
202 if (dim > mImage->GetNumberOfDimensions())
204 mLastError = " Overlay dimension cannot be greater then reference image!";
207 if (mOverlayReader == NULL)
208 mOverlayReader = new vvImageReader;
209 std::vector<std::string> filenames;
210 filenames.push_back(filename);
211 mOverlayReader->SetInputFilenames(filenames);
212 mOverlayReader->Update(mImage->GetNumberOfDimensions(),component.c_str(),mType);
213 if (mOverlayReader->GetLastError().size() == 0)
215 for ( unsigned int i = 0; i < mSlicers.size(); i++)
217 mSlicers[i]->SetOverlay(mOverlayReader->GetOutput());
222 mLastError = mOverlayReader->GetLastError();
228 bool vvSlicerManager::SetFusion(std::string filename,int dim, std::string component)
230 mFusionName = filename;
231 mFusionComponent = component;
232 if (dim > mImage->GetNumberOfDimensions())
234 mLastError = " Overlay dimension cannot be greater then reference image!";
237 if (mFusionReader == NULL)
238 mFusionReader = new vvImageReader;
239 std::vector<std::string> filenames;
240 filenames.push_back(filename);
241 mFusionReader->SetInputFilenames(filenames);
242 mFusionReader->Update(mImage->GetNumberOfDimensions(),component.c_str(),mType);
243 if (mFusionReader->GetLastError().size() == 0)
245 for ( unsigned int i = 0; i < mSlicers.size(); i++)
247 mSlicers[i]->SetFusion(mFusionReader->GetOutput());
252 mLastError = mFusionReader->GetLastError();
255 double *fusRange = mFusionReader->GetOutput()->GetVTKImages()[0]->GetScalarRange();
256 mFusionLevel = (fusRange[0]+fusRange[1])/2;
257 mFusionWindow = fusRange[1]-fusRange[0];
261 bool vvSlicerManager::SetVF(std::string filename)
263 if (mVectorReader == NULL)
264 mVectorReader = new vvImageReader;
265 mVectorReader->SetInputFilename(filename);
266 mVectorReader->Update(VECTORFIELD);
267 if (mVectorReader->GetLastError().size() != 0)
269 mLastError = mVectorReader->GetLastError();
273 return SetVF(mVectorReader->GetOutput(),filename);
276 bool vvSlicerManager::SetVF(vvImage::Pointer vf,std::string filename)
278 if (vf->GetNumberOfDimensions() > mImage->GetNumberOfDimensions())
280 mLastError = " Vector field dimension cannot be greater then reference image!";
285 for ( unsigned int i = 0; i < mSlicers.size(); i++)
287 mSlicers[i]->SetVF(vf);
292 void vvSlicerManager::SetExtractedImage(std::string filename,vvImage::Pointer image, int slice)
294 mFileName = filename;
295 mImage = vvImage::New();
296 if (image->GetNumberOfDimensions() == 4)
298 mImage->AddImage(image->GetVTKImages()[slice]);
299 for ( unsigned int i = 0; i < mSlicers.size(); i++)
301 mSlicers[i]->SetFileName(vtksys::SystemTools::GetFilenameWithoutLastExtension(filename));
302 mSlicers[i]->SetImage(mImage);
307 vtkImageClip* clipper = vtkImageClip::New();
309 image->GetVTKImages()[0]->GetWholeExtent(extent);
310 clipper->SetInput(image->GetVTKImages()[0]);
311 clipper->SetOutputWholeExtent(extent[0],extent[1],extent[2],extent[3],slice,slice);
313 mImage->AddImage(clipper->GetOutput());
314 for ( unsigned int i = 0; i < mSlicers.size(); i++)
316 mSlicers[i]->SetFileName(vtksys::SystemTools::GetFilenameWithoutLastExtension(filename));
317 mSlicers[i]->SetImage(mImage);
323 vvSlicer* vvSlicerManager::GetSlicer(int i)
328 void vvSlicerManager::UpdateSlicer(int num, bool state)
330 if (mSlicers[num]->GetImage())
331 mSlicers[num]->SetDisplayMode(state);
334 void vvSlicerManager::SetSlicerWindow(int i, vtkRenderWindow* RW)
336 mSlicers[i]->SetRenderWindow(i,RW);
339 void vvSlicerManager::SetInteractorStyleNavigator(int i, vtkInteractorStyle* style)
341 vvSlicerManagerCommand *smc = vvSlicerManagerCommand::New();
343 smc->SetSlicerNumber(i);
344 mSlicers[i]->GetRenderWindow()->GetInteractor()->SetInteractorStyle(style);
346 mSlicers[i]->GetRenderWindow()->GetInteractor()->
347 GetInteractorStyle()->AddObserver(vtkCommand::KeyPressEvent, smc);
348 mSlicers[i]->GetRenderWindow()->GetInteractor()->
349 GetInteractorStyle()->AddObserver(vtkCommand::WindowLevelEvent, smc);
350 mSlicers[i]->GetRenderWindow()->GetInteractor()->
351 GetInteractorStyle()->AddObserver(vtkCommand::EndWindowLevelEvent, smc);
352 mSlicers[i]->GetRenderWindow()->GetInteractor()->
353 GetInteractorStyle()->AddObserver(vtkCommand::StartWindowLevelEvent, smc);
354 mSlicers[i]->GetRenderWindow()->GetInteractor()->
355 GetInteractorStyle()->AddObserver(vtkCommand::PickEvent, smc);
356 mSlicers[i]->GetRenderWindow()->GetInteractor()->
357 GetInteractorStyle()->AddObserver(vtkCommand::StartPickEvent, smc);
358 mSlicers[i]->GetRenderWindow()->GetInteractor()->
359 GetInteractorStyle()->AddObserver(vtkCommand::LeaveEvent, smc);
360 mSlicers[i]->GetRenderWindow()->GetInteractor()->
361 GetInteractorStyle()->AddObserver(vtkCommand::UserEvent, smc);
362 mSlicers[i]->GetRenderWindow()->GetInteractor()->
363 GetInteractorStyle()->AddObserver(vtkCommand::MouseWheelForwardEvent, smc);
364 mSlicers[i]->GetRenderWindow()->GetInteractor()->
365 GetInteractorStyle()->AddObserver(vtkCommand::MouseWheelBackwardEvent, smc);
369 void vvSlicerManager::SetTSlice(int slice)
373 else if (slice > mSlicers[0]->GetTMax())
374 slice = mSlicers[0]->GetTMax();
376 mLandmarks->SetTime(slice);
377 for ( unsigned int i = 0; i < mSlicers.size(); i++)
379 mSlicers[i]->SetTSlice(slice);
380 if (mSlicers[i]->GetImageActor()->GetVisibility())
385 void vvSlicerManager::SetNextTSlice(int originating_slicer)
387 int t = mSlicers[0]->GetTSlice();
389 if (t > mSlicers[0]->GetTMax())
391 emit UpdateTSlice(originating_slicer,t);
394 void vvSlicerManager::SetPreviousTSlice(int originating_slicer)
396 int t = mSlicers[0]->GetTSlice();
399 t = mSlicers[0]->GetTMax();
400 emit UpdateTSlice(originating_slicer,t);
403 void vvSlicerManager::ToggleInterpolation()
405 bool interpolate=!(mSlicers[0]->GetImageActor()->GetInterpolate());
406 for ( unsigned int i = 0; i < mSlicers.size(); i++)
408 mSlicers[i]->GetImageActor()->SetInterpolate(interpolate);
413 void vvSlicerManager::SetTSliceInSlicer(int tslice, int slicer)
417 else if (tslice > mSlicers[slicer]->GetTMax())
418 tslice = mSlicers[slicer]->GetTMax();
420 mLandmarks->SetTime(tslice);
421 mSlicers[slicer]->SetTSlice(tslice);
422 if (mSlicers[slicer]->GetImageActor()->GetVisibility())
423 UpdateTSlice(slicer);
426 void vvSlicerManager::SetColorWindow(double s)
428 for ( unsigned int i = 0; i < mSlicers.size(); i++)
430 mSlicers[i]->SetColorWindow(s);
434 void vvSlicerManager::SetColorLevel(double s)
436 for ( unsigned int i = 0; i < mSlicers.size(); i++)
438 mSlicers[i]->SetColorLevel(s);
442 void vvSlicerManager::SetCursorVisibility(int s)
444 for ( unsigned int i = 0; i < mSlicers.size(); i++)
446 mSlicers[i]->SetCursorVisibility(s);
450 void vvSlicerManager::SetOpacity(int i, double factor)
452 mSlicers[i]->SetOpacity(1/factor);
455 void vvSlicerManager::UpdateViews(int current,int slicer)
457 double x = (mSlicers[slicer]->GetCurrentPosition()[0] - mSlicers[slicer]->GetInput()->GetOrigin()[0])
458 /mSlicers[slicer]->GetInput()->GetSpacing()[0];
459 double y = (mSlicers[slicer]->GetCurrentPosition()[1] - mSlicers[slicer]->GetInput()->GetOrigin()[1])
460 /mSlicers[slicer]->GetInput()->GetSpacing()[1];
461 double z = (mSlicers[slicer]->GetCurrentPosition()[2] - mSlicers[slicer]->GetInput()->GetOrigin()[2])
462 /mSlicers[slicer]->GetInput()->GetSpacing()[2];
464 if (x >= mSlicers[slicer]->GetInput()->GetWholeExtent()[0] &&
465 x <= mSlicers[slicer]->GetInput()->GetWholeExtent()[1] &&
466 y >= mSlicers[slicer]->GetInput()->GetWholeExtent()[2] &&
467 y <= mSlicers[slicer]->GetInput()->GetWholeExtent()[3] &&
468 z >= mSlicers[slicer]->GetInput()->GetWholeExtent()[4] &&
469 z <= mSlicers[slicer]->GetInput()->GetWholeExtent()[5])
471 mSlicers[slicer]->UpdateCursorPosition();
472 mSlicers[slicer]->SetCursorColor(10,212,255);
474 switch (mSlicers[slicer]->GetSliceOrientation())
476 case vtkImageViewer2::SLICE_ORIENTATION_XY:
477 if (mSlicers[slicer]->GetSlice() == (int)floor(z))
478 mSlicers[slicer]->Render();
480 mSlicers[slicer]->SetSlice((int)floor(z));
483 case vtkImageViewer2::SLICE_ORIENTATION_XZ:
484 if (mSlicers[slicer]->GetSlice() == (int)floor(y))
485 mSlicers[slicer]->Render();
487 mSlicers[slicer]->SetSlice((int)floor(y));
490 case vtkImageViewer2::SLICE_ORIENTATION_YZ:
491 if (mSlicers[slicer]->GetSlice() == (int)floor(x))
492 mSlicers[slicer]->Render();
494 mSlicers[slicer]->SetSlice((int)floor(x));
498 for ( unsigned int i = 0; i < mSlicers.size(); i++)
500 if (i != (unsigned int)slicer && mSlicers[i]->GetImageActor()->GetVisibility()
501 && mSlicers[i]->GetRenderWindow()->GetSize()[0] > 2
502 && mSlicers[i]->GetRenderWindow()->GetSize()[1] > 2)
504 mSlicers[i]->SetCurrentPosition(mSlicers[slicer]->GetCurrentPosition()[0],
505 mSlicers[slicer]->GetCurrentPosition()[1],
506 mSlicers[slicer]->GetCurrentPosition()[2],
507 mSlicers[slicer]->GetTSlice());
508 mSlicers[i]->UpdateCursorPosition();
509 if (current) //do not display corner annotation if image is the one picked
511 mSlicers[i]->SetCurrentPosition(-VTK_DOUBLE_MAX,-VTK_DOUBLE_MAX,
512 -VTK_DOUBLE_MAX, mSlicers[slicer]->GetTSlice());
513 mSlicers[i]->SetCursorColor(255,10,212);
517 mSlicers[i]->SetCursorColor(150,10,282);
519 switch (mSlicers[i]->GetSliceOrientation())
521 case vtkImageViewer2::SLICE_ORIENTATION_XY:
522 if (mSlicers[i]->GetSlice() == (int)floor(z))
523 mSlicers[i]->Render();
525 mSlicers[i]->SetSlice((int)floor(z));
528 case vtkImageViewer2::SLICE_ORIENTATION_XZ:
529 if (mSlicers[i]->GetSlice() == (int)floor(y))
530 mSlicers[i]->Render();
532 mSlicers[i]->SetSlice((int)floor(y));
535 case vtkImageViewer2::SLICE_ORIENTATION_YZ:
536 if (mSlicers[i]->GetSlice() == (int)floor(x))
537 mSlicers[i]->Render();
539 mSlicers[i]->SetSlice((int)floor(x));
549 void vvSlicerManager::UpdateLinked(int slicer)
551 double x = (mSlicers[slicer]->GetCurrentPosition()[0] - mSlicers[slicer]->GetInput()->GetOrigin()[0])
552 /mSlicers[slicer]->GetInput()->GetSpacing()[0];
553 double y = (mSlicers[slicer]->GetCurrentPosition()[1] - mSlicers[slicer]->GetInput()->GetOrigin()[1])
554 /mSlicers[slicer]->GetInput()->GetSpacing()[1];
555 double z = (mSlicers[slicer]->GetCurrentPosition()[2] - mSlicers[slicer]->GetInput()->GetOrigin()[2])
556 /mSlicers[slicer]->GetInput()->GetSpacing()[2];
558 if (x >= mSlicers[slicer]->GetInput()->GetWholeExtent()[0] &&
559 x <= mSlicers[slicer]->GetInput()->GetWholeExtent()[1] &&
560 y >= mSlicers[slicer]->GetInput()->GetWholeExtent()[2] &&
561 y <= mSlicers[slicer]->GetInput()->GetWholeExtent()[3] &&
562 z >= mSlicers[slicer]->GetInput()->GetWholeExtent()[4] &&
563 z <= mSlicers[slicer]->GetInput()->GetWholeExtent()[5])
565 for (std::list<std::string>::const_iterator i = mLinkedId.begin(); i != mLinkedId.end(); i++)
567 emit UpdateLinkManager(*i, slicer,mSlicers[slicer]->GetCurrentPosition()[0],
568 mSlicers[slicer]->GetCurrentPosition()[1],
569 mSlicers[slicer]->GetCurrentPosition()[2],mSlicers[slicer]->GetTSlice());
573 double vvSlicerManager::GetColorWindow()
576 return mSlicers[0]->GetColorWindow();
580 double vvSlicerManager::GetColorLevel()
583 return mSlicers[0]->GetColorLevel();
587 void vvSlicerManager::Render()
589 for ( unsigned int i = 0; i < mSlicers.size(); i++)
591 mSlicers[i]->Render();
595 void vvSlicerManager::GenerateDefaultLookupTable()
598 SetColorMap(mColorMap);
601 void vvSlicerManager::Reload()
603 mReader->Update(mType);
604 mImage=mReader->GetOutput();
605 for ( unsigned int i = 0; i < mSlicers.size(); i++)
607 mSlicers[i]->SetImage(mImage);
611 void vvSlicerManager::ReloadFusion()
613 mFusionReader->Update();
614 for ( unsigned int i = 0; i < mSlicers.size(); i++)
616 mSlicers[i]->SetFusion(mFusionReader->GetOutput());
617 mSlicers[i]->Render();
621 void vvSlicerManager::ReloadOverlay()
623 mOverlayReader->Update();
624 for ( unsigned int i = 0; i < mSlicers.size(); i++)
626 mSlicers[i]->SetOverlay(mOverlayReader->GetOutput());
627 mSlicers[i]->Render();
631 void vvSlicerManager::ReloadVF()
633 mVectorReader->Update(VECTORFIELD); //deletes the old images through the VF::Init() function
634 mVF=mVectorReader->GetOutput();
635 for ( unsigned int i = 0; i < mSlicers.size(); i++)
637 mSlicers[i]->SetVF(mVF);
638 mSlicers[i]->Render();
642 void vvSlicerManager::RemoveActor(const std::string& actor_type, int overlay_index)
644 for (unsigned int i = 0; i < mSlicers.size();i++)
646 mSlicers[i]->RemoveActor(actor_type,overlay_index);
648 if (actor_type=="vector")
652 delete mVectorReader;
658 void vvSlicerManager::RemoveActors()
660 ///This method leaks a few objects. See RemoveActor for what a correct implementation would look like
661 for ( unsigned int i = 0; i < mSlicers.size(); i++)
663 mSlicers[i]->SetDisplayMode(0);
664 mSlicers[i]->GetRenderer()->RemoveActor(mSlicers[i]->GetImageActor());
668 void vvSlicerManager::UpdateInfoOnCursorPosition(int slicer)
670 // int view = mSlicers[slicer]->GetSliceOrientation();
671 // int slice = mSlicers[slicer]->GetSlice();
672 double x = mSlicers[slicer]->GetCursorPosition()[0];
673 double y = mSlicers[slicer]->GetCursorPosition()[1];
674 double z = mSlicers[slicer]->GetCursorPosition()[2];
675 double X = (x - mSlicers[slicer]->GetInput()->GetOrigin()[0])/
676 mSlicers[slicer]->GetInput()->GetSpacing()[0];
677 double Y = (y - mSlicers[slicer]->GetInput()->GetOrigin()[1])/
678 mSlicers[slicer]->GetInput()->GetSpacing()[1];
679 double Z = (z - mSlicers[slicer]->GetInput()->GetOrigin()[2])/
680 mSlicers[slicer]->GetInput()->GetSpacing()[2];
681 double value = -VTK_DOUBLE_MAX;
683 double xVec=0, yVec=0, zVec=0, valueVec=0;
686 double valueOver=0, valueFus=0;
687 if (X >= mSlicers[slicer]->GetInput()->GetWholeExtent()[0] &&
688 X <= mSlicers[slicer]->GetInput()->GetWholeExtent()[1] &&
689 Y >= mSlicers[slicer]->GetInput()->GetWholeExtent()[2] &&
690 Y <= mSlicers[slicer]->GetInput()->GetWholeExtent()[3] &&
691 Z >= mSlicers[slicer]->GetInput()->GetWholeExtent()[4] &&
692 Z <= mSlicers[slicer]->GetInput()->GetWholeExtent()[5])
694 value = mSlicers[slicer]->GetInput()->GetScalarComponentAsDouble(
698 if (mSlicers[slicer]->GetVFActor() && mSlicers[slicer]->GetVFActor()->GetVisibility())
701 unsigned int currentTime = mSlicers[slicer]->GetTSlice();
702 vtkImageData *vf = NULL;
704 if (mSlicers[slicer]->GetVF()->GetVTKImages().size() > currentTime)
705 vf = mSlicers[slicer]->GetVF()->GetVTKImages()[currentTime];
707 vf = mSlicers[slicer]->GetVF()->GetVTKImages()[0];
711 double Xvf = (x - vf->GetOrigin()[0])/ vf->GetSpacing()[0];
712 double Yvf = (y - vf->GetOrigin()[1])/ vf->GetSpacing()[1];
713 double Zvf = (z - vf->GetOrigin()[2])/ vf->GetSpacing()[2];
714 xVec = vf->GetScalarComponentAsDouble( (int)floor(Xvf), (int)floor(Yvf), (int)floor(Zvf),0);
715 yVec = vf->GetScalarComponentAsDouble( (int)floor(Xvf), (int)floor(Yvf), (int)floor(Zvf),1);
716 zVec = vf->GetScalarComponentAsDouble( (int)floor(Xvf), (int)floor(Yvf), (int)floor(Zvf),2);
717 valueVec = sqrt(xVec*xVec + yVec*yVec + zVec*zVec);
720 if (mSlicers[slicer]->GetOverlayActor() && mSlicers[slicer]->GetOverlayActor()->GetVisibility())
723 double Xover = (x - mSlicers[slicer]->GetOverlay()->GetOrigin()[0])
724 /mSlicers[slicer]->GetOverlay()->GetSpacing()[0];
725 double Yover = (y - mSlicers[slicer]->GetOverlay()->GetOrigin()[1])
726 /mSlicers[slicer]->GetOverlay()->GetSpacing()[1];
727 double Zover = (z - mSlicers[slicer]->GetOverlay()->GetOrigin()[2])
728 /mSlicers[slicer]->GetOverlay()->GetSpacing()[2];
729 if (Xover >= mSlicers[slicer]->GetOverlayMapper()->GetInput()->GetWholeExtent()[0] &&
730 Xover <= mSlicers[slicer]->GetOverlayMapper()->GetInput()->GetWholeExtent()[1] &&
731 Yover >= mSlicers[slicer]->GetOverlayMapper()->GetInput()->GetWholeExtent()[2] &&
732 Yover <= mSlicers[slicer]->GetOverlayMapper()->GetInput()->GetWholeExtent()[3] &&
733 Zover >= mSlicers[slicer]->GetOverlayMapper()->GetInput()->GetWholeExtent()[4] &&
734 Zover <= mSlicers[slicer]->GetOverlayMapper()->GetInput()->GetWholeExtent()[5])
736 valueOver = static_cast<vtkImageData*>(mSlicers[slicer]->GetOverlayMapper()->GetInput())->
737 GetScalarComponentAsDouble(
740 (int)floor(Zover),0);
743 if (mSlicers[slicer]->GetFusionActor() && mSlicers[slicer]->GetFusionActor()->GetVisibility())
746 double Xfus = (x - mSlicers[slicer]->GetFusion()->GetOrigin()[0])
747 /mSlicers[slicer]->GetFusion()->GetSpacing()[0];
748 double Yfus = (y - mSlicers[slicer]->GetFusion()->GetOrigin()[1])
749 /mSlicers[slicer]->GetFusion()->GetSpacing()[1];
750 double Zfus = (z - mSlicers[slicer]->GetFusion()->GetOrigin()[2])
751 /mSlicers[slicer]->GetFusion()->GetSpacing()[2];
752 if (Xfus >= mSlicers[slicer]->GetFusionMapper()->GetInput()->GetWholeExtent()[0] &&
753 Xfus <= mSlicers[slicer]->GetFusionMapper()->GetInput()->GetWholeExtent()[1] &&
754 Yfus >= mSlicers[slicer]->GetFusionMapper()->GetInput()->GetWholeExtent()[2] &&
755 Yfus <= mSlicers[slicer]->GetFusionMapper()->GetInput()->GetWholeExtent()[3] &&
756 Zfus >= mSlicers[slicer]->GetFusionMapper()->GetInput()->GetWholeExtent()[4] &&
757 Zfus <= mSlicers[slicer]->GetFusionMapper()->GetInput()->GetWholeExtent()[5])
759 valueFus = static_cast<vtkImageData*>(mSlicers[slicer]->GetFusionMapper()->GetInput())->
760 GetScalarComponentAsDouble(
766 emit UpdatePosition(mSlicers[slicer]->GetCursorVisibility(),
768 emit UpdateVector(displayVec,xVec, yVec, zVec, valueVec);
769 emit UpdateOverlay(displayOver,valueOver,value);
770 emit UpdateFusion(displayFus,valueFus);
771 for (unsigned int i = 0; i < mSlicers.size(); i++)
773 if (mSlicers[i]->GetImageActor()->GetVisibility() == 1)
774 emit UpdateWindows(i,mSlicers[i]->GetSliceOrientation(),mSlicers[i]->GetSlice());
776 emit UpdateWindows(i,-1,-1);
781 void vvSlicerManager::Activated()
783 emit currentImageChanged(mId);
786 void vvSlicerManager::UpdateWindowLevel()
788 emit WindowLevelChanged(mSlicers[0]->GetColorWindow(),mSlicers[0]->GetColorLevel(),mPreset,mColorMap);
791 void vvSlicerManager::UpdateSlice(int slicer)
793 emit UpdateSlice(slicer, mSlicers[slicer]->GetSlice());
796 void vvSlicerManager::UpdateTSlice(int slicer)
798 emit UpdateTSlice(slicer,mSlicers[0]->GetTSlice());
801 void vvSlicerManager::UpdateSliceRange(int slicer)
803 emit UpdateSliceRange(slicer,
804 mSlicers[slicer]->GetSliceRange()[0], mSlicers[slicer]->GetSliceRange()[1],
805 0,mSlicers[slicer]->GetTMax());
808 void vvSlicerManager::SetPreset(int preset)
810 //vtkLookupTable* LUT = static_cast<vtkLookupTable*>(mSlicers[0]->GetWindowLevel()->GetLookupTable());
811 double window = mSlicers[0]->GetColorWindow();
812 double level = mSlicers[0]->GetColorLevel();
814 std::string component_type=mImage->GetScalarTypeAsString();
818 if (component_type == "unsigned_char")
823 else if (component_type == "short")
831 mImage->GetScalarRange(range);
832 window = range[1] - range[0];
833 level = (range[1] + range[0])* 0.5;
864 this->SetColorWindow(window);
865 this->SetColorLevel(level);
873 void vvSlicerManager::SetLocalColorWindowing(const int slicer)
876 this->mSlicers[slicer]->GetExtremasAroundMousePointer(min, max);
877 this->SetColorWindow(max-min);
878 this->SetColorLevel(0.5*(min+max));
882 void vvSlicerManager::SetColorMap()
884 SetColorMap(mColorMap);
887 void vvSlicerManager::SetColorMap(int colormap)
890 range[0] = mSlicers[0]->GetInput()->GetScalarRange()[0];
891 range[1] = mSlicers[0]->GetInput()->GetScalarRange()[1];
893 double window = mSlicers[0]->GetWindowLevel()->GetWindow();
894 double level = mSlicers[0]->GetWindowLevel()->GetLevel();
896 vtkLookupTable* LUT = static_cast<vtkLookupTable*>(mSlicers[0]->GetWindowLevel()->GetLookupTable());
906 LUT = vtkLookupTable::New();
907 LUT->SetValueRange(0,1);
908 LUT->SetSaturationRange(1,1);
909 LUT->SetHueRange(0,0.18);
913 LUT = vtkLookupTable::New();
914 LUT->SetValueRange(0,1);
915 LUT->SetSaturationRange(1,1);
916 LUT->SetHueRange(0.4,0.80);
920 LUT = vtkLookupTable::New();
921 LUT->SetValueRange(0,1);
922 LUT->SetSaturationRange(1,1);
923 LUT->SetHueRange(0,1);
927 LUT = vtkLookupTable::New();
928 LUT->SetValueRange(0.,1);
929 LUT->SetSaturationRange(1,1);
930 LUT->SetHueRange(1,0.1);
931 //LUT->SetRampToLinear();
936 LUT->SetTableRange(level-fabs(window)/4,level+fabs(window)/4);
939 vtkLookupTable* fusLUT = NULL;
940 if (mSlicers[0]->GetFusion())
942 fusLUT = vtkLookupTable::New();
944 fusRange[0] = mFusionLevel - mFusionWindow/2;
945 fusRange[1] = mFusionLevel + mFusionWindow/2;
946 fusLUT->SetTableRange(fusRange[0],fusRange[1]);
947 fusLUT->SetValueRange(1,1);
948 fusLUT->SetSaturationRange(1,1);
949 if (mFusionColorMap == 1)
950 fusLUT->SetHueRange(0,0.18);
951 else if (mFusionColorMap == 2)
952 fusLUT->SetHueRange(0.4,0.80);
953 else if (mFusionColorMap == 3)
954 fusLUT->SetHueRange(0,1);
956 if (mFusionColorMap == 0)
959 for ( unsigned int i = 0; i < mSlicers.size(); i++) {
960 if (mSlicers[i]->GetOverlay() && mSlicers[i]->GetOverlayActor()->GetVisibility()) {
961 vtkLookupTable* supLUT = vtkLookupTable::New();
962 supLUT->SetTableRange(range[0],range[1]);
963 supLUT->SetValueRange(1,1);
964 supLUT->SetSaturationRange(1,1);
965 supLUT->SetHueRange(double(mOverlayColor)/360,double(mOverlayColor)/360);
967 vtkLookupTable* invLUT = vtkLookupTable::New();
968 invLUT->SetTableRange(range[0],range[1]);
969 invLUT->SetValueRange(1,1);
970 invLUT->SetSaturationRange(1,1);
971 invLUT->SetHueRange(double((mOverlayColor+180)%360)/360,double((mOverlayColor+180)%360)/360);
973 dynamic_cast<vvImageMapToWLColors*>(mSlicers[i]->GetWindowLevel())
974 ->SetWindowLevelMode(true);
975 mSlicers[i]->GetWindowLevel()->SetLookupTable(supLUT);
976 mSlicers[i]->GetOverlayMapper()->SetLookupTable(invLUT);
980 else if (mSlicers[i]->GetOverlay())
982 //dynamic_cast<vvImageMapToWLColors*>(mSlicers[i]->GetWindowLevel())
983 //->SetWindowLevelMode(false);
984 mSlicers[i]->GetWindowLevel()->SetLookupTable(LUT);
988 mSlicers[i]->GetWindowLevel()->SetLookupTable(LUT);
990 if (mSlicers[i]->GetFusion() && mSlicers[i]->GetFusionActor()->GetVisibility())
992 mSlicers[i]->GetFusionActor()->SetOpacity(double(mFusionOpacity)/100);
993 mSlicers[i]->GetFusionMapper()->SetLookupTable(fusLUT);
999 mColorMap = colormap;
1002 vvLandmarks* vvSlicerManager::GetLandmarks()
1004 if (mLandmarks == NULL)
1006 mLandmarks = new vvLandmarks(mSlicers[0]->GetTMax()+1);
1007 for (unsigned int i = 0; i < mSlicers.size(); i++)
1008 mSlicers[i]->SetLandmarks(mLandmarks);
1013 void vvSlicerManager::AddLandmark(float x,float y,float z,float t)
1015 double x_index = (x - mSlicers[0]->GetInput()->GetOrigin()[0])/mSlicers[0]->GetInput()->GetSpacing()[0];
1016 double y_index = (y - mSlicers[0]->GetInput()->GetOrigin()[1])/mSlicers[0]->GetInput()->GetSpacing()[1];
1017 double z_index = (z - mSlicers[0]->GetInput()->GetOrigin()[2])/mSlicers[0]->GetInput()->GetSpacing()[2];
1018 double value = mSlicers[0]->GetInput()->GetScalarComponentAsDouble(
1022 this->GetLandmarks()->AddLandmark(x,y,z,t,value);
1023 emit LandmarkAdded();