recon_array_functions.h

Go to the documentation of this file.

//
//
/*
    Copyright (C) 2000 PARAPET partners
    Copyright (C) 2000- 2008, 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 __stir_recon_array_functions_h_
#define __stir_recon_array_functions_h_

#include "stir/common.h"

START_NAMESPACE_STIR

// TODO template this all in data type of Viewgram et al

template <typename elemT>
class SegmentByView;
template <typename elemT>
class SegmentBySinogram;
template <typename elemT>
class Viewgram;
template <typename elemT>
class RelatedViewgrams;
template <int num_dimensions, typename elemT>
class DiscretisedDensity;

#if 0
void multiply_and_add(DiscretisedDensity<3,float> &image_res, const DiscretisedDensity<3,float> &image_scaled, float scalar);
#endif

float neg_trunc(float x);

// MJ 17/05/2000 disabled
//  gives 1 if negative, 0 otherwise
// float neg_indicate(float x){return (x<=0.0)?1.0:0.0;}

#if 0

void threshold_min_to_small_positive_value_and_truncate_rim(DiscretisedDensity<3,float>& input_image, const int rim_truncation_image=0);

void divide_and_truncate(DiscretisedDensity<3,float>& numerator, 
                         const DiscretisedDensity<3,float>& denominator,
                         const int rim_truncation,
                         int & count);
#endif

void divide_and_truncate(Viewgram<float>& numerator,
                         const Viewgram<float>& denominator,
                         const int rim_truncation_sino,
                         int& count,
                         int& count2,
                         double* f = NULL);

void divide_and_truncate(RelatedViewgrams<float>& numerator,
                         const RelatedViewgrams<float>& denominator,
                         const int rim_truncation_sino,
                         int& count,
                         int& count2,
                         double* f = NULL);

void truncate_rim(DiscretisedDensity<3, float>& image_input,
                  const int rim_truncation_image,
                  const bool strictly_less_than_radius = true);

void truncate_rim(SegmentByView<float>& seg, const int rim_truncation_sino);

void truncate_rim(Viewgram<float>& viewgram, const int rim_truncation_sino);

void truncate_end_planes(DiscretisedDensity<3, float>& input_image, int input_num_planes = 1);

void divide_array(SegmentByView<float>& numerator, const SegmentByView<float>& denominator);

void divide_array(DiscretisedDensity<3, float>& numerator, const DiscretisedDensity<3, float>& denominator);

// MJ 03/01/2000  Trying to adhoc parallelize a loglikelihood computation

void accumulate_loglikelihood(Viewgram<float>& projection_data,
                              const Viewgram<float>& estimated_projections,
                              const int rim_truncation_sino,
                              double* accum);

END_NAMESPACE_STIR
#endif // __recon_array_functions_h_