SeparableCartesianMetzImageFilter.h

Go to the documentation of this file.

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

    SPDX-License-Identifier: Apache-2.0

    See STIR/LICENSE.txt for details
*/

#ifndef __stir_SeparableCartesianMetzImageFilter_H__
#define __stir_SeparableCartesianMetzImageFilter_H__

#include "stir/SeparableMetzArrayFilter.h"
#include "stir/RegisteredParsingObject.h"
#include "stir/DataProcessor.h"
#include "stir/DiscretisedDensity.h"

START_NAMESPACE_STIR

// TODO!! remove define

#define num_dimensions 3

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

public:
  static const char* const registered_name;

  SeparableCartesianMetzImageFilter();

  // Construct metz filter given parameters
  // SeparableCartesianMetzImageFilter(const double fwhm_x,const double fwhm_y, const double fwhm_z,const int metz_power_x,const
  // int metz_power_y, const int metz_power_z);

  // Succeeded consistency_check( const DiscretisedDensity<num_dimensions,elemT>& image) const;

  VectorWithOffset<float> get_metz_fwhms() const;
  VectorWithOffset<float> get_metz_powers() const;

  VectorWithOffset<int> get_max_kernel_sizes() const;

private:
  VectorWithOffset<float> fwhms;
  VectorWithOffset<float> metz_powers;
  VectorWithOffset<int> max_kernel_sizes;

  SeparableMetzArrayFilter<num_dimensions, elemT> metz_filter;

  void set_defaults() override;
  void initialise_keymap() override;

  Succeeded virtual_set_up(const DiscretisedDensity<num_dimensions, elemT>& image) override;
  // new
  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