SeparableGaussianArrayFilter.h

Go to the documentation of this file.

/*
    Copyright (C) 2000 - 2009-06-22, Hammersmith Imanet Ltd
    Copyright (C) 2011-07-01 - 2011, Kris Thielemans
    Copyright (C) 2018, 2019, 2020, UCL
    See STIR/LICENSE.txt for details
*/

#ifndef __stir_SeparableGaussianArrayFilter_H__
#define __stir_SeparableGaussianArrayFilter_H__

#include "stir/DiscretisedDensity.h"
#include "stir/SeparableArrayFunctionObject.h"
#include "stir/Array.h"
#include "stir/BasicCoordinate.h"

#include <vector>

START_NAMESPACE_STIR

template <int num_dimensions, typename elemT>
class SeparableGaussianArrayFilter : public SeparableArrayFunctionObject<num_dimensions, elemT>
{
public:
  SeparableGaussianArrayFilter();


  SeparableGaussianArrayFilter(const BasicCoordinate<num_dimensions, float>& fwhm,
                               const BasicCoordinate<num_dimensions, int>& max_kernel_sizes,
                               bool normalise = true);

  SeparableGaussianArrayFilter(const float fwhm, const float max_kernel_sizes, bool normalise = true);

private:
  void construct_filter(bool normalise = true);

  void calculate_coefficients(VectorWithOffset<elemT>& filter_coefficients,
                              const int max_kernel_sizes,
                              const float fwhm,
                              bool normalise);

  BasicCoordinate<num_dimensions, float> fwhms;
  BasicCoordinate<num_dimensions, int> max_kernel_sizes;
};

END_NAMESPACE_STIR

#endif