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