SPECTUB_Weight3d.h

/*
   Copyright (c) 2013, Biomedical Image Group (GIB), Universitat de Barcelona, Barcelona, Spain.
   Copyright (c) 2013, University College London
   This file is part of STIR.
 
   SPDX-License-Identifier: Apache-2.0
 
   See STIR/LICENSE.txt for details
 
 \author Carles Falcon
 
*/
 
#ifndef _WEIGHT3D_SPECTUB_H
#define _WEIGHT3D_SPECTUB_H

 
namespace SPECTUB
{
 
void wm_calculation(const int kOS,
                    const SPECTUB::angle_type* const ang,
                    SPECTUB::voxel_type vox,
                    bin_type bin,
                    const SPECTUB::volume_type& vol,
                    const proj_type& prj,
                    const float* attmap,
                    const bool* msk_3d,
                    const bool* msk_2d,
                    const int maxszb,
                    const SPECTUB::discrf_type* const gaussdens,
                    const int* const NITEMS,
                    SPECTUB::wm_da_type& wm,
                    const SPECTUB::wmh_type& wmh,
                    const float* Rrad);
 
void wm_size_estimation(int kOS,
                        const SPECTUB::angle_type* const ang,
                        SPECTUB::voxel_type vox,
                        bin_type bin,
                        const SPECTUB::volume_type& vol,
                        const SPECTUB::proj_type& prj,
                        const bool* const msk_3d,
                        const bool* const msk_2d,
                        const int maxszb,
                        const SPECTUB::discrf_type* const gaussdens,
                        int* NITEMS,
                        const SPECTUB::wmh_type& wmh,
                        const float* Rrad);
 
//... geometric component ............................................
 
void calc_gauss(SPECTUB::discrf_type* gaussdens);
 
void calc_vxprj(SPECTUB::angle_type* ang);
 
void voxel_projection(SPECTUB::voxel_type* vox, float* eff, float lngcmd2, const SPECTUB::wmh_type& wmh);
 
void fill_psf_no(SPECTUB::psf2da_type* psf,
                 SPECTUB::psf1d_type* psf1d_h,
                 const SPECTUB::voxel_type& vox,
                 const angle_type* const ang,
                 float szdx,
                 const SPECTUB::wmh_type& wmh);
 
void fill_psf_2d(SPECTUB::psf2da_type* psf,
                 SPECTUB::psf1d_type* psf1d_h,
                 const SPECTUB::voxel_type& vox,
                 SPECTUB::discrf_type const* const gaussdens,
                 float szdx,
                 const SPECTUB::wmh_type& wmh);
 
void fill_psf_3d(SPECTUB::psf2da_type* psf,
                 SPECTUB::psf1d_type* psf1d_h,
                 SPECTUB::psf1d_type* psf1d_v,
                 const SPECTUB::voxel_type& vox,
                 SPECTUB::discrf_type const* const gaussdens,
                 float szdx,
                 float thdx,
                 float thcmd2,
                 const wmh_type& wmh);
 
void calc_psf_bin(float center_psf,
                  float binszcm,
                  SPECTUB::discrf_type const* const vxprj,
                  SPECTUB::psf1d_type* psf,
                  const SPECTUB::wmh_type& wmh);
 
//... attenuation...................................................
 
// not used
// void size_attpth_simple( psf2da_type *psf, voxel_type vox, volume_type vol, float *att, angle_type *ang );
 
// void size_attpth_full( psf2da_type *psf, voxel_type vox, volume_type vol, float *att, angle_type *ang );
 
void
calc_att_path(const bin_type& bin, const SPECTUB::voxel_type& vox, const SPECTUB::volume_type& vol, SPECTUB::attpth_type* attpth);
 
float calc_att(const SPECTUB::attpth_type* const attpth, const float* const attmap, int islc, const SPECTUB::wmh_type& wmh);
 
int comp_dist(float dx, float dy, float dz, float dlast);
 
void error_weight3d(int nerr, const std::string& txt); // error messages in weight3d_SPECT
 
} // namespace SPECTUB
 
#endif