From 915bf2f85e1afdd1fdacfb7aa14e1f3b5bef269d Mon Sep 17 00:00:00 2001
From: schaerer <schaerer>
Date: Tue, 5 Oct 2010 13:26:37 +0000
Subject: [PATCH] allow the user to chose the axis along which to compute the
 mip

---
 vv/vvMaximumIntensityProjection.cxx | 5 ++++-
 vv/vvMaximumIntensityProjection.h   | 4 +++-
 2 files changed, 7 insertions(+), 2 deletions(-)

diff --git a/vv/vvMaximumIntensityProjection.cxx b/vv/vvMaximumIntensityProjection.cxx
index 7c2c3d7..4eb84f7 100644
--- a/vv/vvMaximumIntensityProjection.cxx
+++ b/vv/vvMaximumIntensityProjection.cxx
@@ -25,9 +25,12 @@
 #include "vvFromITK.h"
 #include "vvMaximumIntensityProjection.h"
 #include <QMessageBox>
+#include <QInputDialog>
 
 void vvMaximumIntensityProjection::Compute(vvSlicerManager * slicer_manager)
 {
+  mDimension = QInputDialog::getInteger(0, "MIP Axis","Choose the axis along which to perform the MIP",0,0,\
+                                     slicer_manager->GetImage()->GetNumberOfDimensions()-1,1);
 #define TRY_TYPE(TYPE)							\
   if (clitk::IsSameType<TYPE>(image->GetScalarTypeAsITKString())) { this->Update_WithPixelType<TYPE>(image); return; }
   std::string list = clitk::CreateListOfTypes<short>();
@@ -69,7 +72,7 @@ void vvMaximumIntensityProjection::Update_WithDimAndPixelType(vvImage::Pointer i
   typedef itk::Image<PixelType,Dim-1> OutputImageType;
   typedef itk::MaximumProjectionImageFilter<ImageType,OutputImageType> FilterType;
   typename FilterType::Pointer filter = FilterType::New();
-  filter->SetProjectionDimension(Dim-1);
+  filter->SetProjectionDimension(mDimension);
   typename ImageType::ConstPointer input = vvImageToITK<ImageType>(image);
   filter->SetInput(input);
   filter->Update();
diff --git a/vv/vvMaximumIntensityProjection.h b/vv/vvMaximumIntensityProjection.h
index 6da17f0..9d15164 100644
--- a/vv/vvMaximumIntensityProjection.h
+++ b/vv/vvMaximumIntensityProjection.h
@@ -23,7 +23,7 @@ class vvSlicerManager;
 class vvMaximumIntensityProjection
 {
 public:
-    vvMaximumIntensityProjection() {};
+    vvMaximumIntensityProjection(): mDimension(0) {};
     ///Computes the MIP image on the given vvImage
     void Compute(vvSlicerManager*);
     vvImage::Pointer GetOutput() {return mOutputImage;};
@@ -33,6 +33,8 @@ protected:
     template<class PixelType,int Dim> void Update_WithDimAndPixelType(vvImage::Pointer);
     template<class PixelType> void Update_WithPixelType(vvImage::Pointer);
     vvImage::Pointer mOutputImage;
+    //Dimension along which to compute the MIP
+    unsigned int mDimension;
 };
 
 #endif
-- 
2.47.1