NonseparableConvolutionUsingRealDFTImageFilter.h

Go to the documentation of this file.

//
//
/*
    Copyright (C) 2007- 2007, Hammersmith Imanet
    This file is part of STIR.

    SPDX-License-Identifier: Apache-2.0
    See STIR/LICENSE.txt for details
*/

#ifndef __stir_NonseparableConvolutionUsingRealDFTImageFilter_H__
#define __stir_NonseparableConvolutionUsingRealDFTImageFilter_H__

#include "stir/NonseparableConvolutionUsingRealDFTImageFilter.h"
#include "stir/ArrayFilterUsingRealDFTWithPadding.h"
#include "stir/RegisteredParsingObject.h"
#include "stir/DataProcessor.h"
#include "stir/DiscretisedDensity.h"
#include "stir/VectorWithOffset.h"

START_NAMESPACE_STIR

// TODO!! remove define
// currently fixed at 3 because I didn't really have a good idea for the parsing
// keywords in n dimensions.
//
#define num_dimensions 3

template <typename elemT>
class NonseparableConvolutionUsingRealDFTImageFilter
    : public RegisteredParsingObject<NonseparableConvolutionUsingRealDFTImageFilter<elemT>,
                                     DataProcessor<DiscretisedDensity<3, elemT>>,
                                     DataProcessor<DiscretisedDensity<3, elemT>>>
{
private:
  typedef RegisteredParsingObject<NonseparableConvolutionUsingRealDFTImageFilter<elemT>,
                                  DataProcessor<DiscretisedDensity<3, elemT>>,
                                  DataProcessor<DiscretisedDensity<3, elemT>>>
      base_type;

public:
  static const char* const registered_name;

  NonseparableConvolutionUsingRealDFTImageFilter();


  NonseparableConvolutionUsingRealDFTImageFilter(const Array<num_dimensions, elemT>& filter_coefficients);

private:
  std::string _kernel_filename;
  shared_ptr<DiscretisedDensity<num_dimensions, elemT>> _kernel_sptr;
  shared_ptr<ArrayFilterUsingRealDFTWithPadding<num_dimensions, elemT>> _array_filter_sptr; // ChT::float
  Array<num_dimensions, elemT> _filter_coefficients;

  void set_defaults() override;
  void initialise_keymap() override;
  bool post_processing() override;

  Succeeded virtual_set_up(const DiscretisedDensity<num_dimensions, elemT>& image) override;
  void virtual_apply(DiscretisedDensity<num_dimensions, elemT>& out_density,
                     const DiscretisedDensity<num_dimensions, elemT>& in_density) const override;
  void virtual_apply(DiscretisedDensity<num_dimensions, elemT>& density) const override;
};

#undef num_dimensions

END_NAMESPACE_STIR

#endif