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 vvSlicerManagerCallback *smc = vvSlicerManagerCallback::New();
343 mSlicers[i]->GetRenderWindow()->GetInteractor()->SetInteractorStyle(style);
345 mSlicers[i]->GetRenderWindow()->GetInteractor()->
346 GetInteractorStyle()->AddObserver(vtkCommand::KeyPressEvent, smc);
347 mSlicers[i]->GetRenderWindow()->GetInteractor()->
348 GetInteractorStyle()->AddObserver(vtkCommand::WindowLevelEvent, smc);
349 mSlicers[i]->GetRenderWindow()->GetInteractor()->
350 GetInteractorStyle()->AddObserver(vtkCommand::EndWindowLevelEvent, smc);
351 mSlicers[i]->GetRenderWindow()->GetInteractor()->
352 GetInteractorStyle()->AddObserver(vtkCommand::StartWindowLevelEvent, smc);
353 mSlicers[i]->GetRenderWindow()->GetInteractor()->
354 GetInteractorStyle()->AddObserver(vtkCommand::PickEvent, smc);
355 mSlicers[i]->GetRenderWindow()->GetInteractor()->
356 GetInteractorStyle()->AddObserver(vtkCommand::StartPickEvent, smc);
357 mSlicers[i]->GetRenderWindow()->GetInteractor()->
358 GetInteractorStyle()->AddObserver(vtkCommand::LeaveEvent, smc);
359 mSlicers[i]->GetRenderWindow()->GetInteractor()->
360 GetInteractorStyle()->AddObserver(vtkCommand::UserEvent, smc);
361 mSlicers[i]->GetRenderWindow()->GetInteractor()->
362 GetInteractorStyle()->AddObserver(vtkCommand::MouseWheelForwardEvent, smc);
363 mSlicers[i]->GetRenderWindow()->GetInteractor()->
364 GetInteractorStyle()->AddObserver(vtkCommand::MouseWheelBackwardEvent, smc);
368 void vvSlicerManager::SetTSlice(int slice)
372 else if (slice > mSlicers[0]->GetTMax())
373 slice = mSlicers[0]->GetTMax();
375 mLandmarks->SetTime(slice);
376 for ( unsigned int i = 0; i < mSlicers.size(); i++)
378 mSlicers[i]->SetTSlice(slice);
379 if (mSlicers[i]->GetImageActor()->GetVisibility())
384 void vvSlicerManager::SetNextTSlice(int originating_slicer)
386 int t = mSlicers[0]->GetTSlice();
388 if (t > mSlicers[0]->GetTMax())
390 emit UpdateTSlice(originating_slicer,t);
393 void vvSlicerManager::SetPreviousTSlice(int originating_slicer)
395 int t = mSlicers[0]->GetTSlice();
398 t = mSlicers[0]->GetTMax();
399 emit UpdateTSlice(originating_slicer,t);
402 void vvSlicerManager::ToggleInterpolation()
404 bool interpolate=!(mSlicers[0]->GetImageActor()->GetInterpolate());
405 for ( unsigned int i = 0; i < mSlicers.size(); i++)
407 mSlicers[i]->GetImageActor()->SetInterpolate(interpolate);
412 void vvSlicerManager::SetTSliceInSlicer(int tslice, int slicer)
416 else if (tslice > mSlicers[slicer]->GetTMax())
417 tslice = mSlicers[slicer]->GetTMax();
419 mLandmarks->SetTime(tslice);
420 mSlicers[slicer]->SetTSlice(tslice);
421 if (mSlicers[slicer]->GetImageActor()->GetVisibility())
422 UpdateTSlice(slicer);
425 void vvSlicerManager::SetColorWindow(double s)
427 for ( unsigned int i = 0; i < mSlicers.size(); i++)
429 mSlicers[i]->SetColorWindow(s);
433 void vvSlicerManager::SetColorLevel(double s)
435 for ( unsigned int i = 0; i < mSlicers.size(); i++)
437 mSlicers[i]->SetColorLevel(s);
441 void vvSlicerManager::SetCursorVisibility(int s)
443 for ( unsigned int i = 0; i < mSlicers.size(); i++)
445 mSlicers[i]->SetCursorVisibility(s);
449 void vvSlicerManager::SetOpacity(int i, double factor)
451 mSlicers[i]->SetOpacity(1/factor);
454 void vvSlicerManager::UpdateViews(int current,int slicer)
456 double x = (mSlicers[slicer]->GetCurrentPosition()[0] - mSlicers[slicer]->GetInput()->GetOrigin()[0])
457 /mSlicers[slicer]->GetInput()->GetSpacing()[0];
458 double y = (mSlicers[slicer]->GetCurrentPosition()[1] - mSlicers[slicer]->GetInput()->GetOrigin()[1])
459 /mSlicers[slicer]->GetInput()->GetSpacing()[1];
460 double z = (mSlicers[slicer]->GetCurrentPosition()[2] - mSlicers[slicer]->GetInput()->GetOrigin()[2])
461 /mSlicers[slicer]->GetInput()->GetSpacing()[2];
463 if (x >= mSlicers[slicer]->GetInput()->GetWholeExtent()[0] &&
464 x <= mSlicers[slicer]->GetInput()->GetWholeExtent()[1] &&
465 y >= mSlicers[slicer]->GetInput()->GetWholeExtent()[2] &&
466 y <= mSlicers[slicer]->GetInput()->GetWholeExtent()[3] &&
467 z >= mSlicers[slicer]->GetInput()->GetWholeExtent()[4] &&
468 z <= mSlicers[slicer]->GetInput()->GetWholeExtent()[5])
470 mSlicers[slicer]->UpdateCursorPosition();
471 mSlicers[slicer]->SetCursorColor(10,212,255);
473 switch (mSlicers[slicer]->GetSliceOrientation())
475 case vtkImageViewer2::SLICE_ORIENTATION_XY:
476 if (mSlicers[slicer]->GetSlice() == (int)floor(z))
477 mSlicers[slicer]->Render();
479 mSlicers[slicer]->SetSlice((int)floor(z));
482 case vtkImageViewer2::SLICE_ORIENTATION_XZ:
483 if (mSlicers[slicer]->GetSlice() == (int)floor(y))
484 mSlicers[slicer]->Render();
486 mSlicers[slicer]->SetSlice((int)floor(y));
489 case vtkImageViewer2::SLICE_ORIENTATION_YZ:
490 if (mSlicers[slicer]->GetSlice() == (int)floor(x))
491 mSlicers[slicer]->Render();
493 mSlicers[slicer]->SetSlice((int)floor(x));
497 for ( unsigned int i = 0; i < mSlicers.size(); i++)
499 if (i != (unsigned int)slicer && mSlicers[i]->GetImageActor()->GetVisibility()
500 && mSlicers[i]->GetRenderWindow()->GetSize()[0] > 2
501 && mSlicers[i]->GetRenderWindow()->GetSize()[1] > 2)
503 mSlicers[i]->SetCurrentPosition(mSlicers[slicer]->GetCurrentPosition()[0],
504 mSlicers[slicer]->GetCurrentPosition()[1],
505 mSlicers[slicer]->GetCurrentPosition()[2],
506 mSlicers[slicer]->GetTSlice());
507 mSlicers[i]->UpdateCursorPosition();
508 if (current) //do not display corner annotation if image is the one picked
510 mSlicers[i]->SetCurrentPosition(-VTK_DOUBLE_MAX,-VTK_DOUBLE_MAX,
511 -VTK_DOUBLE_MAX, mSlicers[slicer]->GetTSlice());
512 mSlicers[i]->SetCursorColor(255,10,212);
516 mSlicers[i]->SetCursorColor(150,10,282);
518 switch (mSlicers[i]->GetSliceOrientation())
520 case vtkImageViewer2::SLICE_ORIENTATION_XY:
521 if (mSlicers[i]->GetSlice() == (int)floor(z))
522 mSlicers[i]->Render();
524 mSlicers[i]->SetSlice((int)floor(z));
527 case vtkImageViewer2::SLICE_ORIENTATION_XZ:
528 if (mSlicers[i]->GetSlice() == (int)floor(y))
529 mSlicers[i]->Render();
531 mSlicers[i]->SetSlice((int)floor(y));
534 case vtkImageViewer2::SLICE_ORIENTATION_YZ:
535 if (mSlicers[i]->GetSlice() == (int)floor(x))
536 mSlicers[i]->Render();
538 mSlicers[i]->SetSlice((int)floor(x));
548 void vvSlicerManager::UpdateLinked(int slicer)
550 double x = (mSlicers[slicer]->GetCurrentPosition()[0] - mSlicers[slicer]->GetInput()->GetOrigin()[0])
551 /mSlicers[slicer]->GetInput()->GetSpacing()[0];
552 double y = (mSlicers[slicer]->GetCurrentPosition()[1] - mSlicers[slicer]->GetInput()->GetOrigin()[1])
553 /mSlicers[slicer]->GetInput()->GetSpacing()[1];
554 double z = (mSlicers[slicer]->GetCurrentPosition()[2] - mSlicers[slicer]->GetInput()->GetOrigin()[2])
555 /mSlicers[slicer]->GetInput()->GetSpacing()[2];
557 if (x >= mSlicers[slicer]->GetInput()->GetWholeExtent()[0] &&
558 x <= mSlicers[slicer]->GetInput()->GetWholeExtent()[1] &&
559 y >= mSlicers[slicer]->GetInput()->GetWholeExtent()[2] &&
560 y <= mSlicers[slicer]->GetInput()->GetWholeExtent()[3] &&
561 z >= mSlicers[slicer]->GetInput()->GetWholeExtent()[4] &&
562 z <= mSlicers[slicer]->GetInput()->GetWholeExtent()[5])
564 for (std::list<std::string>::const_iterator i = mLinkedId.begin(); i != mLinkedId.end(); i++)
566 emit UpdateLinkManager(*i, slicer,mSlicers[slicer]->GetCurrentPosition()[0],
567 mSlicers[slicer]->GetCurrentPosition()[1],
568 mSlicers[slicer]->GetCurrentPosition()[2],mSlicers[slicer]->GetTSlice());
572 double vvSlicerManager::GetColorWindow()
575 return mSlicers[0]->GetColorWindow();
579 double vvSlicerManager::GetColorLevel()
582 return mSlicers[0]->GetColorLevel();
586 void vvSlicerManager::Render()
588 for ( unsigned int i = 0; i < mSlicers.size(); i++)
590 mSlicers[i]->Render();
594 void vvSlicerManager::GenerateDefaultLookupTable()
597 SetColorMap(mColorMap);
600 void vvSlicerManager::Reload()
602 mReader->Update(mType);
603 mImage=mReader->GetOutput();
604 for ( unsigned int i = 0; i < mSlicers.size(); i++)
606 mSlicers[i]->SetImage(mImage);
610 void vvSlicerManager::ReloadFusion()
612 mFusionReader->Update();
613 for ( unsigned int i = 0; i < mSlicers.size(); i++)
615 mSlicers[i]->SetFusion(mFusionReader->GetOutput());
616 mSlicers[i]->Render();
620 void vvSlicerManager::ReloadOverlay()
622 mOverlayReader->Update();
623 for ( unsigned int i = 0; i < mSlicers.size(); i++)
625 mSlicers[i]->SetOverlay(mOverlayReader->GetOutput());
626 mSlicers[i]->Render();
630 void vvSlicerManager::ReloadVF()
632 mVectorReader->Update(VECTORFIELD); //deletes the old images through the VF::Init() function
633 mVF=mVectorReader->GetOutput();
634 for ( unsigned int i = 0; i < mSlicers.size(); i++)
636 mSlicers[i]->SetVF(mVF);
637 mSlicers[i]->Render();
641 void vvSlicerManager::RemoveActor(const std::string& actor_type, int overlay_index)
643 for (unsigned int i = 0; i < mSlicers.size();i++)
645 mSlicers[i]->RemoveActor(actor_type,overlay_index);
647 if (actor_type=="vector")
651 delete mVectorReader;
657 void vvSlicerManager::RemoveActors()
659 ///This method leaks a few objects. See RemoveActor for what a correct implementation would look like
660 for ( unsigned int i = 0; i < mSlicers.size(); i++)
662 mSlicers[i]->SetDisplayMode(0);
663 mSlicers[i]->GetRenderer()->RemoveActor(mSlicers[i]->GetImageActor());
667 void vvSlicerManager::UpdateInfoOnCursorPosition(int slicer)
669 // int view = mSlicers[slicer]->GetSliceOrientation();
670 // int slice = mSlicers[slicer]->GetSlice();
671 double x = mSlicers[slicer]->GetCursorPosition()[0];
672 double y = mSlicers[slicer]->GetCursorPosition()[1];
673 double z = mSlicers[slicer]->GetCursorPosition()[2];
674 double X = (x - mSlicers[slicer]->GetInput()->GetOrigin()[0])/
675 mSlicers[slicer]->GetInput()->GetSpacing()[0];
676 double Y = (y - mSlicers[slicer]->GetInput()->GetOrigin()[1])/
677 mSlicers[slicer]->GetInput()->GetSpacing()[1];
678 double Z = (z - mSlicers[slicer]->GetInput()->GetOrigin()[2])/
679 mSlicers[slicer]->GetInput()->GetSpacing()[2];
680 double value = -VTK_DOUBLE_MAX;
682 double xVec=0, yVec=0, zVec=0, valueVec=0;
685 double valueOver=0, valueFus=0;
686 if (X >= mSlicers[slicer]->GetInput()->GetWholeExtent()[0] &&
687 X <= mSlicers[slicer]->GetInput()->GetWholeExtent()[1] &&
688 Y >= mSlicers[slicer]->GetInput()->GetWholeExtent()[2] &&
689 Y <= mSlicers[slicer]->GetInput()->GetWholeExtent()[3] &&
690 Z >= mSlicers[slicer]->GetInput()->GetWholeExtent()[4] &&
691 Z <= mSlicers[slicer]->GetInput()->GetWholeExtent()[5])
693 value = mSlicers[slicer]->GetInput()->GetScalarComponentAsDouble(
697 if (mSlicers[slicer]->GetVFActor() && mSlicers[slicer]->GetVFActor()->GetVisibility())
700 unsigned int currentTime = mSlicers[slicer]->GetTSlice();
701 vtkImageData *vf = NULL;
703 if (mSlicers[slicer]->GetVF()->GetVTKImages().size() > currentTime)
704 vf = mSlicers[slicer]->GetVF()->GetVTKImages()[currentTime];
706 vf = mSlicers[slicer]->GetVF()->GetVTKImages()[0];
710 double Xvf = (x - vf->GetOrigin()[0])/ vf->GetSpacing()[0];
711 double Yvf = (y - vf->GetOrigin()[1])/ vf->GetSpacing()[1];
712 double Zvf = (z - vf->GetOrigin()[2])/ vf->GetSpacing()[2];
713 xVec = vf->GetScalarComponentAsDouble( (int)floor(Xvf), (int)floor(Yvf), (int)floor(Zvf),0);
714 yVec = vf->GetScalarComponentAsDouble( (int)floor(Xvf), (int)floor(Yvf), (int)floor(Zvf),1);
715 zVec = vf->GetScalarComponentAsDouble( (int)floor(Xvf), (int)floor(Yvf), (int)floor(Zvf),2);
716 valueVec = sqrt(xVec*xVec + yVec*yVec + zVec*zVec);
719 if (mSlicers[slicer]->GetOverlayActor() && mSlicers[slicer]->GetOverlayActor()->GetVisibility())
722 double Xover = (x - mSlicers[slicer]->GetOverlay()->GetOrigin()[0])
723 /mSlicers[slicer]->GetOverlay()->GetSpacing()[0];
724 double Yover = (y - mSlicers[slicer]->GetOverlay()->GetOrigin()[1])
725 /mSlicers[slicer]->GetOverlay()->GetSpacing()[1];
726 double Zover = (z - mSlicers[slicer]->GetOverlay()->GetOrigin()[2])
727 /mSlicers[slicer]->GetOverlay()->GetSpacing()[2];
728 if (Xover >= mSlicers[slicer]->GetOverlayMapper()->GetInput()->GetWholeExtent()[0] &&
729 Xover <= mSlicers[slicer]->GetOverlayMapper()->GetInput()->GetWholeExtent()[1] &&
730 Yover >= mSlicers[slicer]->GetOverlayMapper()->GetInput()->GetWholeExtent()[2] &&
731 Yover <= mSlicers[slicer]->GetOverlayMapper()->GetInput()->GetWholeExtent()[3] &&
732 Zover >= mSlicers[slicer]->GetOverlayMapper()->GetInput()->GetWholeExtent()[4] &&
733 Zover <= mSlicers[slicer]->GetOverlayMapper()->GetInput()->GetWholeExtent()[5])
735 valueOver = static_cast<vtkImageData*>(mSlicers[slicer]->GetOverlayMapper()->GetInput())->
736 GetScalarComponentAsDouble(
739 (int)floor(Zover),0);
742 if (mSlicers[slicer]->GetFusionActor() && mSlicers[slicer]->GetFusionActor()->GetVisibility())
745 double Xfus = (x - mSlicers[slicer]->GetFusion()->GetOrigin()[0])
746 /mSlicers[slicer]->GetFusion()->GetSpacing()[0];
747 double Yfus = (y - mSlicers[slicer]->GetFusion()->GetOrigin()[1])
748 /mSlicers[slicer]->GetFusion()->GetSpacing()[1];
749 double Zfus = (z - mSlicers[slicer]->GetFusion()->GetOrigin()[2])
750 /mSlicers[slicer]->GetFusion()->GetSpacing()[2];
751 if (Xfus >= mSlicers[slicer]->GetFusionMapper()->GetInput()->GetWholeExtent()[0] &&
752 Xfus <= mSlicers[slicer]->GetFusionMapper()->GetInput()->GetWholeExtent()[1] &&
753 Yfus >= mSlicers[slicer]->GetFusionMapper()->GetInput()->GetWholeExtent()[2] &&
754 Yfus <= mSlicers[slicer]->GetFusionMapper()->GetInput()->GetWholeExtent()[3] &&
755 Zfus >= mSlicers[slicer]->GetFusionMapper()->GetInput()->GetWholeExtent()[4] &&
756 Zfus <= mSlicers[slicer]->GetFusionMapper()->GetInput()->GetWholeExtent()[5])
758 valueFus = static_cast<vtkImageData*>(mSlicers[slicer]->GetFusionMapper()->GetInput())->
759 GetScalarComponentAsDouble(
765 emit UpdatePosition(mSlicers[slicer]->GetCursorVisibility(),
767 emit UpdateVector(displayVec,xVec, yVec, zVec, valueVec);
768 emit UpdateOverlay(displayOver,valueOver,value);
769 emit UpdateFusion(displayFus,valueFus);
770 for (unsigned int i = 0; i < mSlicers.size(); i++)
772 if (mSlicers[i]->GetImageActor()->GetVisibility() == 1)
773 emit UpdateWindows(i,mSlicers[i]->GetSliceOrientation(),mSlicers[i]->GetSlice());
775 emit UpdateWindows(i,-1,-1);
780 void vvSlicerManager::Activated()
782 emit currentImageChanged(mId);
785 void vvSlicerManager::UpdateWindowLevel()
787 emit WindowLevelChanged(mSlicers[0]->GetColorWindow(),mSlicers[0]->GetColorLevel(),mPreset,mColorMap);
790 void vvSlicerManager::UpdateSlice(int slicer)
792 emit UpdateSlice(slicer, mSlicers[slicer]->GetSlice());
795 void vvSlicerManager::UpdateTSlice(int slicer)
797 emit UpdateTSlice(slicer,mSlicers[0]->GetTSlice());
800 void vvSlicerManager::UpdateSliceRange(int slicer)
802 emit UpdateSliceRange(slicer,
803 mSlicers[slicer]->GetSliceRange()[0], mSlicers[slicer]->GetSliceRange()[1],
804 0,mSlicers[slicer]->GetTMax());
807 void vvSlicerManager::SetPreset(int preset)
810 mImage->GetScalarRange(range);
811 //vtkLookupTable* LUT = static_cast<vtkLookupTable*>(mSlicers[0]->GetWindowLevel()->GetLookupTable());
812 double window = mSlicers[0]->GetColorWindow();
813 double level = mSlicers[0]->GetColorLevel();
815 std::string component_type=mImage->GetScalarTypeAsString();
819 if (component_type == "unsigned_char")
824 else if (component_type == "short")
831 window = range[1] - range[0];
832 level = (range[1] + range[0])* 0.5;
863 this->SetColorWindow(window);
864 this->SetColorLevel(level);
871 void vvSlicerManager::SetColorMap()
873 SetColorMap(mColorMap);
876 void vvSlicerManager::SetColorMap(int colormap)
879 range[0] = mSlicers[0]->GetInput()->GetScalarRange()[0];
880 range[1] = mSlicers[0]->GetInput()->GetScalarRange()[1];
882 double window = mSlicers[0]->GetWindowLevel()->GetWindow();
883 double level = mSlicers[0]->GetWindowLevel()->GetLevel();
885 vtkLookupTable* LUT = static_cast<vtkLookupTable*>(mSlicers[0]->GetWindowLevel()->GetLookupTable());
895 LUT = vtkLookupTable::New();
896 LUT->SetValueRange(0,1);
897 LUT->SetSaturationRange(1,1);
898 LUT->SetHueRange(0,0.18);
902 LUT = vtkLookupTable::New();
903 LUT->SetValueRange(0,1);
904 LUT->SetSaturationRange(1,1);
905 LUT->SetHueRange(0.4,0.80);
909 LUT = vtkLookupTable::New();
910 LUT->SetValueRange(0,1);
911 LUT->SetSaturationRange(1,1);
912 LUT->SetHueRange(0,1);
916 LUT = vtkLookupTable::New();
917 LUT->SetValueRange(0.,1);
918 LUT->SetSaturationRange(1,1);
919 LUT->SetHueRange(1,0.1);
920 //LUT->SetRampToLinear();
925 LUT->SetTableRange(level-fabs(window)/4,level+fabs(window)/4);
928 vtkLookupTable* fusLUT = NULL;
929 if (mSlicers[0]->GetFusion())
931 fusLUT = vtkLookupTable::New();
933 fusRange[0] = mFusionLevel - mFusionWindow/2;
934 fusRange[1] = mFusionLevel + mFusionWindow/2;
935 fusLUT->SetTableRange(fusRange[0],fusRange[1]);
936 fusLUT->SetValueRange(1,1);
937 fusLUT->SetSaturationRange(1,1);
938 if (mFusionColorMap == 1)
939 fusLUT->SetHueRange(0,0.18);
940 else if (mFusionColorMap == 2)
941 fusLUT->SetHueRange(0.4,0.80);
942 else if (mFusionColorMap == 3)
943 fusLUT->SetHueRange(0,1);
945 if (mFusionColorMap == 0)
948 for ( unsigned int i = 0; i < mSlicers.size(); i++) {
949 if (mSlicers[i]->GetOverlay() && mSlicers[i]->GetOverlayActor()->GetVisibility()) {
950 vtkLookupTable* supLUT = vtkLookupTable::New();
951 supLUT->SetTableRange(range[0],range[1]);
952 supLUT->SetValueRange(1,1);
953 supLUT->SetSaturationRange(1,1);
954 supLUT->SetHueRange(double(mOverlayColor)/360,double(mOverlayColor)/360);
956 vtkLookupTable* invLUT = vtkLookupTable::New();
957 invLUT->SetTableRange(range[0],range[1]);
958 invLUT->SetValueRange(1,1);
959 invLUT->SetSaturationRange(1,1);
960 invLUT->SetHueRange(double((mOverlayColor+180)%360)/360,double((mOverlayColor+180)%360)/360);
962 dynamic_cast<vvImageMapToWLColors*>(mSlicers[i]->GetWindowLevel())
963 ->SetWindowLevelMode(true);
964 mSlicers[i]->GetWindowLevel()->SetLookupTable(supLUT);
965 mSlicers[i]->GetOverlayMapper()->SetLookupTable(invLUT);
969 else if (mSlicers[i]->GetOverlay())
971 //dynamic_cast<vvImageMapToWLColors*>(mSlicers[i]->GetWindowLevel())
972 //->SetWindowLevelMode(false);
973 mSlicers[i]->GetWindowLevel()->SetLookupTable(LUT);
977 mSlicers[i]->GetWindowLevel()->SetLookupTable(LUT);
979 if (mSlicers[i]->GetFusion() && mSlicers[i]->GetFusionActor()->GetVisibility())
981 mSlicers[i]->GetFusionActor()->SetOpacity(double(mFusionOpacity)/100);
982 mSlicers[i]->GetFusionMapper()->SetLookupTable(fusLUT);
988 mColorMap = colormap;
991 vvLandmarks* vvSlicerManager::GetLandmarks()
993 if (mLandmarks == NULL)
995 mLandmarks = new vvLandmarks(mSlicers[0]->GetTMax()+1);
996 for (unsigned int i = 0; i < mSlicers.size(); i++)
997 mSlicers[i]->SetLandmarks(mLandmarks);
1002 void vvSlicerManager::AddLandmark(float x,float y,float z,float t)
1004 double x_index = (x - mSlicers[0]->GetInput()->GetOrigin()[0])/mSlicers[0]->GetInput()->GetSpacing()[0];
1005 double y_index = (y - mSlicers[0]->GetInput()->GetOrigin()[1])/mSlicers[0]->GetInput()->GetSpacing()[1];
1006 double z_index = (z - mSlicers[0]->GetInput()->GetOrigin()[2])/mSlicers[0]->GetInput()->GetSpacing()[2];
1007 double value = mSlicers[0]->GetInput()->GetScalarComponentAsDouble(
1011 this->GetLandmarks()->AddLandmark(x,y,z,t,value);
1012 emit LandmarkAdded();