/*****************************************************************************
-** This is part of the CTSim program
-** Copyright (C) 1983-2000 Kevin Rosenberg
-**
-** $Id: ir.h,v 1.22 2000/06/07 07:43:39 kevin Exp $
-** $Log: ir.h,v $
-** Revision 1.22 2000/06/07 07:43:39 kevin
-** *** empty log message ***
-**
-** Revision 1.21 2000/06/07 03:49:54 kevin
-** *** empty log message ***
-**
-** Revision 1.20 2000/06/07 02:30:27 kevin
-** Added C++ image files
-**
-** Revision 1.19 2000/06/07 00:59:38 kevin
-** added imagefiles
-**
-** Revision 1.18 2000/06/05 01:32:45 kevin
-** Added C++ compatibility
-**
-** Revision 1.17 2000/06/03 06:29:08 kevin
-** *** empty log message ***
-**
-** Revision 1.16 2000/05/24 22:48:17 kevin
-** First functional version of SDF library for X-window
-**
-** Revision 1.15 2000/05/16 04:33:17 kevin
-** Updated documentation
-**
-** Revision 1.14 2000/05/11 14:07:00 kevin
-** Added support for Windows NT
-**
-** Revision 1.13 2000/05/08 20:00:48 kevin
-** ANSI C changes
-**
-** Revision 1.12 2000/05/07 12:46:19 kevin
-** made c++ compatible
-**
-** Revision 1.11 2000/05/05 02:37:31 kevin
-** renamed phmelm to pelm
+** FILE IDENTIFICATION
**
-** Revision 1.10 2000/05/04 18:16:34 kevin
-** renamed filter definitions
+** Name: ir.h
+** Purpose: Master Image Reconstruction Header
+** Programmer: Kevin Rosenberg
+** Date Started: July 1, 1984
**
-** Revision 1.9 2000/05/04 04:29:18 kevin
-** *** empty log message ***
-**
-** Revision 1.8 2000/05/04 04:25:55 kevin
-** Renamed phantom and phantom-element functions/variables
-**
-** Revision 1.7 2000/05/03 19:51:41 kevin
-** function renaming for phantoms and phantom elements
-**
-** Revision 1.6 2000/05/03 08:49:49 kevin
-** Code cleanup
-**
-** Revision 1.5 2000/05/02 20:00:25 kevin
-** *** empty log message ***
-**
-** Revision 1.4 2000/05/02 15:31:39 kevin
-** code cleaning
-**
-** Revision 1.3 2000/04/29 23:24:29 kevin
-** *** empty log message ***
-**
-** Revision 1.2 2000/04/28 14:14:16 kevin
-** *** empty log message ***
-**
-** Revision 1.1.1.1 2000/04/28 13:02:43 kevin
-** Initial CVS import for first public release
+** This is part of the CTSim program
+** Copyright (C) 1983-2000 Kevin Rosenberg
**
+** $Id: ir.h,v 1.26 2000/06/13 16:20:31 kevin Exp $
**
**
** This program is free software; you can redistribute it and/or modify
** Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
******************************************************************************/
-/* ir.h Header File for Image Reconstruction System
- * Programmer: Kevin Rosenberg
- * Date Started: 7-1-84
- */
-
#ifndef IR_H
#define IR_H
-struct histo_st {
- int *b; /* Histogram array (# of elements in each bin) */
- int nbin; /* Number of histogram bins */
- double xmin, xmax, xinc; /* Limits of histogram boundaries */
-};
-
-typedef struct histo_st HISTOGRAM;
-
-/*---------------------------------------------------------------------------*/
-
static const int POINTS_PER_CIRCLE=90;
#define MAXREMARK 99
SEGMENT
} PElmType;
-struct pelm_st {
+
+class PhmElement
+{
+ public:
PElmType type; /* pelm type (box, ellipse, etc) */
double atten; /* X-ray attenuation coefficient */
double cx,cy; /* center of pelm */
double radius; /* " */
GRFMTX_2D p_to_o; /* map from phantom to standard pelm coords */
GRFMTX_2D o_to_p; /* map from std pelm coords to phantom coords */
- struct pelm_st *next; /* pointer to next pelm in phantom */
+ class PhmElement *next; /* pointer to next pelm in phantom */
};
-typedef struct pelm_st PELM;
+typedef class PhmElement PELM;
typedef enum {
P_PELMS, /* Phantom made of Pelms */
P_FILTER /* defined only by this type */
} PhmType;
-struct phm_st { /* Phantom structure */
+
+/* Phantom class */
+class Phantom
+{
+public:
PELM *pelm_list; /* pelm linked-list */
PhmType type;
int n_pelm; /* number of pelms in phantom */
double xmin, xmax, ymin, ymax; /* extent of pelms in pelm coordinates */
double radius; /* " " */
};
-typedef struct phm_st PHANTOM;
+typedef class Phantom PHANTOM;
/*----------------------------------------------------------------------*/
typedef float DETECT_TYPE;
-struct detarray_st {
+struct DetectorArray
+{
DETECT_TYPE *detval; /* Pointer to array of values recorded by detector */
int ndet; /* Number of detectors in array */
double view_angle; /* View angle in radians */
};
-typedef struct detarray_st DETARRAY;
+typedef struct DetectorArray DETARRAY;
typedef enum {
DETECTOR_EQUILINEAR
} ScannerGeometry;
-struct detector_st {
+
+class Detector
+{
+ public:
ScannerGeometry geometry; /* Geometry of detectory */
int ndet; /* Number of detectors in array */
int nview; /* Number of rotated views */
double angle; /* Starting angle */
} init;
};
-typedef struct detector_st DETECTOR;
+typedef class Detector DETECTOR;
-struct raysum_st {
+
+class Projections
+{
+ public:
int fd; /* Disk file descriptor */
int file_mode; /* Current file mode (read or write) */
int header_size; /* Size of disk file header */
int geometry; /* Geometry of scanner */
- struct detarray_st **view; /* Pointer to array of detarray_st pointers */
+ struct DetectorArray **view; /* Pointer to array of detarray_st pointers */
char remark[MAXREMARK+1]; /* description of raysum data */
double calctime; /* time required to calculate raysums */
double det_inc; /* increment between detectors */
double phmlen; /* Length of PHANTOM edge (phm is square) */
};
-typedef struct raysum_st RAYSUM;
+typedef class Projections RAYSUM;
+
/*----------------------------------------------------------------------*/
/* USER SYMBOLS */
const static int RAYSUM_TRACE_ROW_ATTEN=18;
+
/*************************************************************************
* FUNCTION DECLARATIONS
************************************************************************/
-#ifdef __cplusplus
-extern "C" {
-#endif /* __cplusplus */
-
-#ifdef MPI_CT
-#define MPI_MAX_PROCESS 128
-struct mpi_ct_st
-{
- int my_rank;
- int nproc;
- int base_local_work_units;
- int remainder_work_units;
- int local_work_units[MPI_MAX_PROCESS];
- int start_work_unit[MPI_MAX_PROCESS];
- MPI_Comm comm;
-};
-
-extern struct mpi_ct_st mpi_ct;
-void mpi_ct_calc_work_units(const unsigned int global_work_units);
-#endif
-
-
-/* bspline.c */
-int bspline(int samples, int zoom_factor, int spline_order, double input[], double output[]);
-
-/* convolve.c */
-double convolve(const double f1[], const double f2[], const double dx, const int n, const int np, const FunctionSymmetry func_type);
-double convolve_both(const double f1[], const double f2[], const double dx, const int n, const int np);
-
-/* dialogs.c */
-int phm_add_pelm_kb(PHANTOM *phm);
-PHANTOM *phm_select(void);
-int interpolation_select(void);
-int filter_select(double *filter_param);
-
-/* filter.c */
-double *filter_generate(const FilterType filt_type, double bw, double xmin, double xmax, int n, double param, const DomainType domain, int numint);
-double filter_spatial_response_calc(int filt_type, double x, double bw, double param, int n);
-double filter_spatial_response_analytic(int filt_type, double x, double bw, double param);
-double filter_frequency_response(int filt_type, double u, double bw, double param);
-double sinc(double x, double mult);
+/* convolve.cpp */
+double convolve (const double f1[], const double f2[], const double dx, const int n, const int np, const FunctionSymmetry func_type);
+
+/* dialogs.cpp */
+int phm_add_pelm_kb (PHANTOM *phm);
+PHANTOM *phm_select (void);
+int interpolation_select (void);
+int filter_select (double *filter_param);
+
+/* filter.cpp */
+double *filter_generate (const FilterType filt_type, double bw, double xmin, double xmax, int n, double param, const DomainType domain, int numint);
+double filter_spatial_response_calc (int filt_type, double x, double bw, double param, int n);
+double filter_spatial_response_analytic (int filt_type, double x, double bw, double param);
+double filter_frequency_response (int filt_type, double u, double bw, double param);
+double sinc (double x, double mult);
double integral_abscos(double u, double w);
-/* options.c */
+/* options.cpp */
int opt_set_trace(const char *optarg);
const char *name_of_phantom(const int phmid);
int opt_set_phantom(const char *optarg);
BackprojType opt_set_backproj(const char *optarg);
const char *name_of_backproj(const BackprojType backproj);
-/* phm.c */
+/* phm.cpp */
PHANTOM *phm_create(const int phmid);
PHANTOM *phm_create_from_file(const char *fname);
PHANTOM *phm_init(void);
void phm_draw(const PHANTOM *phm);
#endif
-/* phmstd.c */
+/* phmstd.cpp */
void phm_std_herman (PHANTOM *phm);
void phm_std_rowland (PHANTOM *phm);
void phm_std_rowland_bordered (PHANTOM *phm);
-/* raycollect.c */
+/* raycollect.cpp */
int raysum_collect(RAYSUM *rs, const DETECTOR *det, const PHANTOM *phm, const int start_view, const int trace, const int unit_pulse);
void rayview(const PHANTOM *phm, DETARRAY *darray, const DETECTOR *det, const double xd1, const double yd1, const double xd2, const double yd2, const double xs1, const double ys1, const double xs2, const double ys2, const int unit_pulse);
double phm_ray_attenuation (const PHANTOM *phm, const double x1, const double y1, const double x2, const double y2);
double pelm_ray_attenuation (PELM *pelm, const double x1, const double y1, const double x2, const double y2);
-int pelm_clip_line (const PELM *pelm, double *x1, double *y1, double *x2, double *y2);
+int pelm_clip_line (const PELM *pelm, double& x1, double& y1, double& x2, double& y2);
void raysum_trace_show_param (const char *label, const char *fmt, int row, int color, ...);
-/* scanner.c */
-DETECTOR *detector_create(const PHANTOM *phm, int geometry, int ndet, int nview, int nsample, const double rot_anglen);
+/* scanner.cpp */
+DETECTOR *detector_create(const PHANTOM *phm, const ScannerGeometry geometry, int ndet, int nview, int nsample, const double rot_anglen);
void detector_free(DETECTOR *det);
-/* rayio.c */
+/* rayio.cpp */
RAYSUM *raysum_create(const char *fname, const int nview, const int ndet);
RAYSUM *raysum_create_from_det(const char *fname, const DETECTOR *det);
RAYSUM *raysum_open(const char *filename);
int raysum_write_header(RAYSUM *rs);
int raysum_read(RAYSUM *rs);
int raysum_write(RAYSUM *rs);
+void raysum_print(const RAYSUM *rs);
+void raysum_print_info(const RAYSUM *rs);
DETARRAY *detarray_alloc(const int n);
void detarray_free(DETARRAY *darray);
int detarray_read(RAYSUM *rs, DETARRAY *darray, const int view_num);
int detarray_write(RAYSUM *rs, const DETARRAY *darray, const int view_num);
-int raysum_print(const RAYSUM *rs);
-
-#ifdef __cplusplus
-}
/* From phm2image.cpp */
void phm_to_imagefile (const PHANTOM *phm, ImageFile& im, const int col_start, const int col_count, const int nsample, const int trace);
const double view_angle, const int interp_type);
void backproj_term_id2 (const RAYSUM *rs, ImageFile& im);
-
-#endif /* __cplusplus */
-
-
#endif
projects / ctsim.git / blobdiff