diff --git a/Modules/DiffusionImaging/DiffusionCore/Algorithms/Reconstruction/itkMultiShellAdcAverageReconstructionImageFilter.cpp b/Modules/DiffusionImaging/DiffusionCore/Algorithms/Reconstruction/itkMultiShellAdcAverageReconstructionImageFilter.cpp
new file mode 100644
index 0000000000..d2644eda4f
--- /dev/null
+++ b/Modules/DiffusionImaging/DiffusionCore/Algorithms/Reconstruction/itkMultiShellAdcAverageReconstructionImageFilter.cpp
@@ -0,0 +1,95 @@
+/*===================================================================
+
+The Medical Imaging Interaction Toolkit (MITK)
+
+Copyright (c) German Cancer Research Center,
+Division of Medical and Biological Informatics.
+All rights reserved.
+
+This software is distributed WITHOUT ANY WARRANTY; without
+even the implied warranty of MERCHANTABILITY or FITNESS FOR
+A PARTICULAR PURPOSE.
+
+See LICENSE.txt or http://www.mitk.org for details.
+
+===================================================================*/
+/*=========================================================================
+
+Program:   Tensor ToolKit - TTK
+Module:    $URL: svn://scm.gforge.inria.fr/svn/ttk/trunk/Algorithms/itkElectrostaticRepulsionDiffusionGradientReductionFilter.txx $
+Language:  C++
+Date:      $Date: 2010-06-07 13:39:13 +0200 (Mo, 07 Jun 2010) $
+Version:   $Revision: 68 $
+
+Copyright (c) INRIA 2010. All rights reserved.
+See LICENSE.txt for details.
+
+This software is distributed WITHOUT ANY WARRANTY; without even
+the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
+PURPOSE.  See the above copyright notices for more information.
+
+=========================================================================*/
+#ifndef _itk_MultiShellAdcAverageReconstructionImageFilter_cpp_
+#define _itk_MultiShellAdcAverageReconstructionImageFilter_cpp_
+#endif
+
+#define _USE_MATH_DEFINES
+
+#include "itkMultiShellAdcAverageReconstructionImageFilter.h"
+#include <math.h>
+#include <time.h>
+#include <itkImageRegionIterator.h>
+#include <itkImageRegion.h>
+
+namespace itk
+{
+
+template <class TInputScalarType, class TOutputScalarType>
+MultiShellAdcAverageReconstructionImageFilter<TInputScalarType, TOutputScalarType>
+::MultiShellAdcAverageReconstructionImageFilter()
+{
+  this->SetNumberOfRequiredInputs( 1 );
+}
+
+template <class TInputScalarType, class TOutputScalarType>
+void
+MultiShellAdcAverageReconstructionImageFilter<TInputScalarType, TOutputScalarType>
+::ThreadedGenerateData(const OutputImageRegionType &outputRegionForThread, int /*threadId*/)
+{
+
+  // Get input gradient image pointer
+  typename InputImageType::Pointer inputImage = static_cast< InputImageType * >(ProcessObject::GetInput(0));
+  // ImageRegionIterator for the input image
+  ImageRegionConstIterator< InputImageType > iit(inputImage, outputRegionForThread);
+  iit.GoToBegin();
+
+  // Get output gradient image pointer
+  typename OutputImageType::Pointer outputImage = static_cast< OutputImageType * >(ProcessObject::GetOutput(0));
+  // ImageRegionIterator for the output image
+  ImageRegionIterator< OutputImageType > oit(outputImage, outputRegionForThread);
+  oit.GoToBegin();
+
+  const int numShells = m_BValueMap.size()-1;
+  BValueMap::iterator it = m_BValueMap.begin();
+  std::vector<double> adcVec = new Vector<double>(numShells);
+
+  // initialize output image
+  typename OutputImageType::Pointer outImage = OutputImageType::New();
+  outImage->SetSpacing( this->GetInput()->GetSpacing() );   // Set the image spacing
+  outImage->SetOrigin( this->GetInput()->GetOrigin() );     // Set the image origin
+  outImage->SetDirection( this->GetInput()->GetDirection() );  // Set the image direction
+  outImage->SetLargestPossibleRegion( this->GetInput()->GetLargestPossibleRegion());
+  outImage->SetBufferedRegion( this->GetInput()->GetLargestPossibleRegion() );
+  outImage->SetRequestedRegion( this->GetInput()->GetLargestPossibleRegion() );
+  outImage->SetVectorLength( vecLength ); // Set the vector length
+  outImage->Allocate();
+
+
+  this->SetNumberOfRequiredOutputs (1);
+  this->SetNthOutput (0, outImage);
+  MITK_INFO << "...done";
+}
+
+
+
+} // end of namespace
diff --git a/Modules/DiffusionImaging/DiffusionCore/Algorithms/Reconstruction/itkMultiShellAdcAverageReconstructionImageFilter.h b/Modules/DiffusionImaging/DiffusionCore/Algorithms/Reconstruction/itkMultiShellAdcAverageReconstructionImageFilter.h
new file mode 100644
index 0000000000..bbb6a7492d
--- /dev/null
+++ b/Modules/DiffusionImaging/DiffusionCore/Algorithms/Reconstruction/itkMultiShellAdcAverageReconstructionImageFilter.h
@@ -0,0 +1,116 @@
+/*===================================================================
+
+The Medical Imaging Interaction Toolkit (MITK)
+
+Copyright (c) German Cancer Research Center,
+Division of Medical and Biological Informatics.
+All rights reserved.
+
+This software is distributed WITHOUT ANY WARRANTY; without
+even the implied warranty of MERCHANTABILITY or FITNESS FOR
+A PARTICULAR PURPOSE.
+
+See LICENSE.txt or http://www.mitk.org for details.
+
+===================================================================*/
+/*=========================================================================
+
+Program:   Tensor ToolKit - TTK
+Module:    $URL: svn://scm.gforge.inria.fr/svn/ttk/trunk/Algorithms/itkElectrostaticRepulsionDiffusionGradientReductionFilter.h $
+Language:  C++
+Date:      $Date: 2010-06-07 13:39:13 +0200 (Mo, 07 Jun 2010) $
+Version:   $Revision: 68 $
+
+Copyright (c) INRIA 2010. All rights reserved.
+See LICENSE.txt for details.
+
+This software is distributed WITHOUT ANY WARRANTY; without even
+the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
+PURPOSE.  See the above copyright notices for more information.
+
+=========================================================================*/
+#ifndef _itk_MultiShellAdcAverageReconstructionImageFilter_h_
+#define _itk_MultiShellAdcAverageReconstructionImageFilter_h_
+
+#include <itkImageToImageFilter.h>
+#include <itkVectorImage.h>
+#include <itkPointShell.h>
+
+namespace itk
+{
+
+/**
+* \brief Select subset of the input vectors equally distributed over the sphere using an iterative electrostatic repulsion strategy.   */
+
+  template <class TInputScalarType, class TOutputScalarType>
+  class MultiShellAdcAverageReconstructionImageFilter
+    : public ImageToImageFilter<itk::VectorImage<TInputScalarType,3>, itk::VectorImage<TOutputScalarType,3> >
+  {
+
+  public:
+
+    typedef MultiShellAdcAverageReconstructionImageFilter Self;
+    typedef SmartPointer<Self>                      Pointer;
+    typedef SmartPointer<const Self>                ConstPointer;
+    typedef ImageToImageFilter< itk::VectorImage<TOutputScalarType,3>, itk::VectorImage<TOutputScalarType,3> >
+      Superclass;
+
+    /** Method for creation through the object factory. */
+    itkNewMacro(Self)
+
+    /** Runtime information support. */
+    itkTypeMacro(MultiShellAdcAverageReconstructionImageFilter, ImageToImageFilter)
+
+    typedef TInputScalarType                         InputScalarType;
+    typedef itk::VectorImage<InputScalarType,3>      InputImageType;
+    typedef typename InputImageType::PixelType       InputPixelType;
+
+    typedef TOutputScalarType                         OutputScalarType;
+    typedef itk::VectorImage<OutputScalarType,3>      OutputImageType;
+    typedef typename OutputImageType::PixelType       OutputPixelType;
+
+    typedef OutputScalarType                          BaselineScalarType;
+    typedef BaselineScalarType                        BaselinePixelType;
+    typedef typename itk::Image<BaselinePixelType,3>  BaselineImageType;
+
+    typedef vnl_vector_fixed< double, 3 >                               GradientDirectionType;
+    typedef itk::VectorContainer< unsigned int, GradientDirectionType > GradientDirectionContainerType;
+
+    typedef std::vector<unsigned int>           IndicesVector;
+    typedef std::map<double, IndicesVector>     BValueMap;
+
+    itkGetMacro(OriginalGradientDirections, GradientDirectionContainerType::Pointer)
+    itkSetMacro(OriginalGradientDirections, GradientDirectionContainerType::Pointer)
+
+    itkGetMacro(GradientDirections, GradientDirectionContainerType::Pointer)
+    itkSetMacro(GradientDirections, GradientDirectionContainerType::Pointer)
+
+    inline void SetOriginalBValueMap(BValueMap inp){m_OriginalBValueMap = inp;}
+
+  protected:
+    MultiShellAdcAverageReconstructionImageFilter();
+    ~MultiShellAdcAverageReconstructionImageFilter() {}
+
+    void ThreadedGenerateData(const OutputImageRegionType &outputRegionForThread, int threadId);
+
+    GradientDirectionContainerType::Pointer m_GradientDirections;   ///< container for the subsampled output gradient directions
+    GradientDirectionContainerType::Pointer m_OriginalGradientDirections;   ///< input gradient directions
+
+    IndicesVector m_UsedGradientIndices;
+    IndicesVector m_UnusedGradientIndices;
+    IndicesVector m_BaselineImageIndices;
+
+    BValueMap m_BValueMap;
+
+    };
+
+
+} // end of namespace
+
+
+#ifndef ITK_MANUAL_INSTANTIATION
+#include "itkMultiShellAdcAverageReconstructionImageFilter.cpp"
+#endif
+
+
+#endif
diff --git a/Modules/DiffusionImaging/DiffusionCore/files.cmake b/Modules/DiffusionImaging/DiffusionCore/files.cmake
index 07866a8d47..44d1bfa51a 100644
--- a/Modules/DiffusionImaging/DiffusionCore/files.cmake
+++ b/Modules/DiffusionImaging/DiffusionCore/files.cmake
@@ -1,116 +1,117 @@
 set(CPP_FILES
 
   # DicomImport
   DicomImport/mitkDicomDiffusionImageReader.cpp
   # DicomImport/mitkGroupDiffusionHeadersFilter.cpp
   DicomImport/mitkDicomDiffusionImageHeaderReader.cpp
   DicomImport/mitkGEDicomDiffusionImageHeaderReader.cpp
   DicomImport/mitkPhilipsDicomDiffusionImageHeaderReader.cpp
   DicomImport/mitkSiemensDicomDiffusionImageHeaderReader.cpp
   DicomImport/mitkSiemensMosaicDicomDiffusionImageHeaderReader.cpp
 
   # DataStructures
   IODataStructures/mitkDiffusionCoreObjectFactory.cpp
 
   # DataStructures -> DWI
   IODataStructures/DiffusionWeightedImages/mitkDiffusionImageHeaderInformation.cpp
   IODataStructures/DiffusionWeightedImages/mitkDiffusionImageSource.cpp
   IODataStructures/DiffusionWeightedImages/mitkNrrdDiffusionImageReader.cpp
   IODataStructures/DiffusionWeightedImages/mitkNrrdDiffusionImageWriter.cpp
   IODataStructures/DiffusionWeightedImages/mitkNrrdDiffusionImageIOFactory.cpp
   IODataStructures/DiffusionWeightedImages/mitkNrrdDiffusionImageWriterFactory.cpp
   IODataStructures/DiffusionWeightedImages/mitkDiffusionImageSerializer.cpp
 
   # DataStructures -> QBall
   IODataStructures/QBallImages/mitkQBallImageSource.cpp
   IODataStructures/QBallImages/mitkNrrdQBallImageReader.cpp
   IODataStructures/QBallImages/mitkNrrdQBallImageWriter.cpp
   IODataStructures/QBallImages/mitkNrrdQBallImageIOFactory.cpp
   IODataStructures/QBallImages/mitkNrrdQBallImageWriterFactory.cpp
   IODataStructures/QBallImages/mitkQBallImage.cpp
   IODataStructures/QBallImages/mitkQBallImageSerializer.cpp
 
   # DataStructures -> Tensor
   IODataStructures/TensorImages/mitkTensorImageSource.cpp
   IODataStructures/TensorImages/mitkNrrdTensorImageReader.cpp
   IODataStructures/TensorImages/mitkNrrdTensorImageWriter.cpp
   IODataStructures/TensorImages/mitkNrrdTensorImageIOFactory.cpp
   IODataStructures/TensorImages/mitkNrrdTensorImageWriterFactory.cpp
   IODataStructures/TensorImages/mitkTensorImage.cpp
   IODataStructures/TensorImages/mitkTensorImageSerializer.cpp
 
   # Rendering
   Rendering/vtkMaskedProgrammableGlyphFilter.cpp
   Rendering/mitkCompositeMapper.cpp
   Rendering/mitkVectorImageVtkGlyphMapper3D.cpp
   Rendering/vtkOdfSource.cxx
   Rendering/vtkThickPlane.cxx
   Rendering/mitkOdfNormalizationMethodProperty.cpp
   Rendering/mitkOdfScaleByProperty.cpp
   Rendering/mitkPlanarFigureMapper3D.cpp
 
   # Algorithms
   Algorithms/mitkPartialVolumeAnalysisHistogramCalculator.cpp
   Algorithms/mitkPartialVolumeAnalysisClusteringCalculator.cpp
 
   # Function Collection
   mitkDiffusionFunctionCollection.cpp
 )
 
 set(H_FILES
   # function Collection
   mitkDiffusionFunctionCollection.h
 
   # Rendering
   Rendering/mitkDiffusionImageMapper.h
   Rendering/mitkOdfVtkMapper2D.h
   Rendering/mitkPlanarFigureMapper3D.h
 
   # Reconstruction
   Algorithms/Reconstruction/itkDiffusionQballReconstructionImageFilter.h
   Algorithms/Reconstruction/mitkTeemDiffusionTensor3DReconstructionImageFilter.h
   Algorithms/Reconstruction/itkAnalyticalDiffusionQballReconstructionImageFilter.h
   Algorithms/Reconstruction/itkDiffusionMultiShellQballReconstructionImageFilter.h
   Algorithms/Reconstruction/itkPointShell.h
   Algorithms/Reconstruction/itkOrientationDistributionFunction.h
   Algorithms/Reconstruction/itkDiffusionIntravoxelIncoherentMotionReconstructionImageFilter.h
+  Algorithms/Reconstruction/itkMultiShellAdcAverageReconstructionImageFilter.h
 
   # IO Datastructures
   IODataStructures/DiffusionWeightedImages/mitkDiffusionImage.h
 
   # Algorithms
   Algorithms/itkDiffusionQballGeneralizedFaImageFilter.h
   Algorithms/itkDiffusionQballPrepareVisualizationImageFilter.h
   Algorithms/itkTensorDerivedMeasurementsFilter.h
   Algorithms/itkBrainMaskExtractionImageFilter.h
   Algorithms/itkB0ImageExtractionImageFilter.h
   Algorithms/itkB0ImageExtractionToSeparateImageFilter.h
   Algorithms/itkTensorImageToDiffusionImageFilter.h
   Algorithms/itkTensorToL2NormImageFilter.h
   Algorithms/itkGaussianInterpolateImageFunction.h
   Algorithms/mitkPartialVolumeAnalysisHistogramCalculator.h
   Algorithms/mitkPartialVolumeAnalysisClusteringCalculator.h
   Algorithms/itkDiffusionTensorPrincipalDirectionImageFilter.h
   Algorithms/itkCartesianToPolarVectorImageFilter.h
   Algorithms/itkPolarToCartesianVectorImageFilter.h
   Algorithms/itkDistanceMapFilter.h
   Algorithms/itkProjectionFilter.h
   Algorithms/itkResidualImageFilter.h
   Algorithms/itkExtractChannelFromRgbaImageFilter.h
   Algorithms/itkTensorReconstructionWithEigenvalueCorrectionFilter.h
   Algorithms/itkMergeDiffusionImagesFilter.h
   Algorithms/itkDwiPhantomGenerationFilter.h
   Algorithms/itkFiniteDiffOdfMaximaExtractionFilter.h
   Algorithms/itkMrtrixPeakImageConverter.h
   Algorithms/itkFslPeakImageConverter.h
   Algorithms/itkShCoefficientImageImporter.h
   Algorithms/itkOdfMaximaExtractionFilter.h
   Algorithms/itkResampleDwiImageFilter.h
 
 
 
 )
 
 set( TOOL_FILES
 )