ColsherFilter.h

Go to the documentation of this file.

//
//

/*
    Copyright (C) 2000 PARAPET partners
    Copyright (C) 2000- 2004, Hammersmith Imanet Ltd

    This file is part of STIR.

    SPDX-License-Identifier: Apache-2.0 AND License-ref-PARAPET-license

    See STIR/LICENSE.txt for details
*/

#ifndef __ColsherFilter_H__
#define __ColsherFilter_H__

#ifdef NRFFT
#  include "stir_experimental/Filter.h"
#else
#  include "stir/ArrayFilterUsingRealDFTWithPadding.h"
#  include "stir/TimedObject.h"
#endif

START_NAMESPACE_STIR
#ifdef NRFFT
template <typename elemT>
class Viewgram;
#endif

class ColsherFilter :
#ifdef NRFFT
    public Filter2D<float>
#else
    public ArrayFilterUsingRealDFTWithPadding<2, float>,
    public TimedObject
#endif
{

public:
#ifndef NRFFT

  ColsherFilter()
  {}
  explicit ColsherFilter(float theta_max,
                         float alpha_colsher_axial = 1.F,
                         float fc_colsher_axial = 0.5F,
                         float alpha_colsher_radial = 1.F,
                         float fc_colsher_radial = 0.5F,
                         const int stretch_factor_axial = 2,
                         const int stretch_factor_planar = 2);

  Succeeded set_up(int height, int width, float theta, float d_a, float d_b);
#else
  ColsherFilter(int height,
                int width,
                float gamma,
                float theta_max,
                float d_a,
                float d_b,
                float alpha_colsher_axial,
                float fc_colsher_axial,
                float alpha_colsher_radial,
                float fc_colsher_radial);
#endif

  virtual std::string parameter_info() const;

  ~ColsherFilter() override
  {}

private:
#ifdef NRFFT
  float gamma;
  /* d_a, d_b are used to convert to millimeter.*/
  float d_a;
  float d_b;
#endif
  float theta_max;
  float alpha_axial;
  float fc_axial;
  float alpha_planar;
  float fc_planar;
  int stretch_factor_axial;
  int stretch_factor_planar;
};

#ifdef NRFFT
void Filter_proj_Colsher(Viewgram<float>& view_i, Viewgram<float>& view_i1, ColsherFilter& CFilter, int PadS, int PadZ);
#endif
END_NAMESPACE_STIR

#endif //  __ColsherFilter_H__