X-Git-Url: http://git.kpe.io/?a=blobdiff_plain;f=include%2Fir.h;h=5a0b7291412afdd97914bfe0c1b1ca105525d8a7;hb=da11926c1039ae4c880915db508d69614ec7c79e;hp=7a3ab95ed581d5ce41e349179dedba2341d98c49;hpb=0ed0b45e2292346f54ea196825135ccc35536bfe;p=ctsim.git diff --git a/include/ir.h b/include/ir.h index 7a3ab95..5a0b729 100644 --- a/include/ir.h +++ b/include/ir.h @@ -2,8 +2,20 @@ ** This is part of the CTSim program ** Copyright (C) 1983-2000 Kevin Rosenberg ** -** $Id: ir.h,v 1.7 2000/05/03 19:51:41 kevin Exp $ +** $Id: ir.h,v 1.11 2000/05/05 02:37:31 kevin Exp $ ** $Log: ir.h,v $ +** Revision 1.11 2000/05/05 02:37:31 kevin +** renamed phmelm to pelm +** +** Revision 1.10 2000/05/04 18:16:34 kevin +** renamed filter definitions +** +** 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 ** @@ -96,38 +108,38 @@ struct image_st { float calctime; /* time to calculate voxels in seconds */ }; -struct phmelem_st { - int type; /* phmelem type (box, ellipse, etc) */ - double atten; /* X-ray attenuation coefficient */ - double cx,cy; /* center of phmelem */ - double u,v; /* size of phmelem */ - double rot; /* phmelem rotation angle (in radians) */ - double *x, *y; /* ptr to array of points in obj world coord */ - int pts; /* number of points in outline arrays */ - double xmin, xmax, ymin, ymax; /* phmelem limits */ - double radius; /* " */ - struct { /* transform matrices */ - GRFMTX_2D p_to_o; /* map from phantom to standard phmelem coords */ - GRFMTX_2D o_to_p; /* map from std phmelem coords to phantom coords */ +struct pelm_st { + int type; /* pelm type (box, ellipse, etc) */ + double atten; /* X-ray attenuation coefficient */ + double cx,cy; /* center of pelm */ + double u,v; /* size of pelm */ + double rot; /* pelm rotation angle (in radians) */ + double *x, *y; /* ptr to array of points in obj world coord */ + int pts; /* number of points in outline arrays */ + double xmin, xmax, ymin, ymax; /* pelm limits */ + double radius; /* " */ + struct { /* transform matrices */ + 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 */ } xform; - struct phmelem_st *next; /* pointer to next phmelem in phantom */ + struct pelm_st *next; /* pointer to next pelm in phantom */ }; -typedef struct phmelem_st PHMELEM; +typedef struct pelm_st PELM; -struct phm_st { /* Phantom structure */ - PHMELEM *objlist; /* phmelem list */ +struct phm_st { /* Phantom structure */ + PELM *pelm_list; /* pelm linked-list */ int type; - int n_pelem; /* number of phmelems in phantom */ - double xmin, xmax, ymin, ymax; /* extent of phmelems in phmelem coordinates */ - double radius; /* " " */ + int n_pelm; /* number of pelms in phantom */ + double xmin, xmax, ymin, ymax; /* extent of pelms in pelm coordinates */ + double radius; /* " " */ }; typedef struct image_st IMAGE; typedef struct phm_st PHANTOM; -#define P_PHMELEMS 0 /* PHANTOM made of phmelems */ -#define P_UNIT_PULSE 1 /* Special PHANTOM, not made of phmelems */ +#define P_PELMS 0 /* Phantom made of pelms */ +#define P_UNIT_PULSE 1 /* Special PHANTOM, not made of pelms */ #define P_FILTER 9 /* defined only by a type */ /*----------------------------------------------------------------------*/ @@ -150,7 +162,7 @@ struct detarray_st { #define DETECTOR_EQUIANGLE 2 #define DETECTOR_EQUILINEAR 3 -struct detect_st { +struct detector_st { int geometry; /* Geometry of detectory */ int ndet; /* Number of detectors in array */ int nview; /* Number of rotated views */ @@ -173,6 +185,7 @@ struct raysum_st { int fd; int file_mode; int header_size; + int geometry; struct detarray_st **view; /* Pointer to array of detarray_st pointers */ char remark[LENREMARK]; /* description of raysum data */ @@ -189,14 +202,14 @@ struct raysum_st { }; typedef struct detarray_st DETARRAY; -typedef struct detect_st DETECTOR; +typedef struct detector_st DETECTOR; typedef struct raysum_st RAYSUM; /*----------------------------------------------------------------------*/ /* USER SYMBOLS */ /*----------------------------------------------------------------------*/ -/* codes for phmelem types, passed to add_obj() */ +/* codes for pelm types, passed to phm_add_pelm() */ #define O_RECTANGLE 1 #define O_TRIANGLE 2 @@ -205,25 +218,23 @@ typedef struct raysum_st RAYSUM; #define O_SEGMENT 5 /* Codes for Coordinate Types */ -/* Defines coords for inside_obj() */ +/* Defines coords for pelm_is_point_inside() */ -#define OBJ_COORD -1 /* Normalized Phmelem Coordinates */ +#define PELM_COORD -1 /* Normalized Pelm Coordinates */ #define PHM_COORD -2 /* User's PHANTOM Coordinates */ /* Codes for Filter types */ -#define W_BANDLIMIT 1 /* filter types for genfilter() */ -#define W_SINC 2 -#define W_G_HAMMING 3 -#define W_COSINE 4 -#define W_TRIANGLE 5 - -#define W_A_BANDLIMIT 11 /* filters times abs() of function */ -#define W_A_SINC 12 -#define W_AG_HAMMING 13 -#define W_A_COSINE 14 - -#define W_SHEPP 21 +#define FILTER_BANDLIMIT 1 /* filter types for genfilter() */ +#define FILTER_SINC 2 +#define FILTER_G_HAMMING 3 +#define FILTER_COSINE 4 +#define FILTER_TRIANGLE 5 +#define FILTER_ABS_BANDLIMIT 11 /* filters times abs() of function */ +#define FILTER_ABS_SINC 12 +#define FILTER_ABS_G_HAMMING 13 +#define FILTER_ABS_COSINE 14 +#define FILTER_SHEPP 21 /* function domains */ @@ -251,7 +262,7 @@ typedef struct raysum_st RAYSUM; /* Constants for sizing PHANTOM */ #define PERCENT_PHM_SIZE_INCR 0.0 /* Fractional increase in PHANTOM */ - /* limits compared to phmelem size */ + /* limits compared to pelm size */ #define N_EXTRA_DETECTORS 4 /* Number of extra detectors */ /* widths when calculating detlen */ @@ -307,17 +318,17 @@ typedef struct raysum_st RAYSUM; #define O_BPROJ_DIFF2_STR "diff2" #define O_BPROJ_IDIFF2_STR "idiff2" -#define RS_TRACE_ROW_TITLE 1 -#define RS_TRACE_ROW_TITLE2 2 -#define RS_TRACE_ROW_PHANT_ID 4 -#define RS_TRACE_ROW_CHROMATIC 7 -#define RS_TRACE_ROW_SCATTER 8 -#define RS_TRACE_ROW_PHOT_STAT 9 -#define RS_TRACE_ROW_NDET 12 -#define RS_TRACE_ROW_NVIEW 13 -#define RS_TRACE_ROW_SAMPLES 14 -#define RS_TRACE_ROW_CURR_VIEW 17 -#define RS_TRACE_ROW_ATTEN 18 +#define RAYSUM_TRACE_ROW_TITLE 1 +#define RAYSUM_TRACE_ROW_TITLE2 2 +#define RAYSUM_TRACE_ROW_PHANT_ID 4 +#define RAYSUM_TRACE_ROW_CHROMATIC 7 +#define RAYSUM_TRACE_ROW_SCATTER 8 +#define RAYSUM_TRACE_ROW_PHOT_STAT 9 +#define RAYSUM_TRACE_ROW_NDET 12 +#define RAYSUM_TRACE_ROW_NVIEW 13 +#define RAYSUM_TRACE_ROW_SAMPLES 14 +#define RAYSUM_TRACE_ROW_CURR_VIEW 17 +#define RAYSUM_TRACE_ROW_ATTEN 18 /*----------------------------------------------------------------------*/ @@ -390,31 +401,32 @@ int main(const int argc, char * const argv[]); /* 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 int func_type); double convolve_both(const double f1[], const double f2[], const double dx, const int n, const int np); -void rotate2d(double x[], double y[], int pts, double angle); -void xlat2d(double x[], double y[], int pts, double xoffset, double yoffset); -void scale2d(double x[], double y[], int pts, double xfact, double yfact); -int circle_pts(double theta); + /* dialogs.c */ -int phm_add_pelem_kb(PHANTOM *phm); +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(int filt_type, double bw, double xmin, double xmax, int n, double param, int 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 abscos_int(double u, double w); +double integral_abscos(double u, double w); + /* image.c */ IMAGE *image_create(const char *fname, const int nx, const int ny); int image_clear(IMAGE *im); int image_save(IMAGE *im); IMAGE *image_load(const char *fname); void image_filter_response(IMAGE *im, int domain, double bw, int filt_type, double filt_param, int opt_trace); + /* options.c */ int opt_set_trace(const char *optarg, const char *program); const char *name_of_phantom(const int phmid); @@ -427,45 +439,49 @@ int opt_set_filter_domain(const char *optarg, const char *program); const char *name_of_filter_domain(const int domain); int opt_set_backproj(const char *optarg, const char *program); const char *name_of_backproj(const int backproj); + /* phm.c */ PHANTOM *phm_create(const int phmid); PHANTOM *phm_create_from_file(const char *fname); PHANTOM *phm_init(void); -int phm_add_pelem_file(PHANTOM *phm, const char *fname); -void phm_add_pelem (PHANTOM *phm, const int type, const double cx, const double cy, +int phm_add_pelm_file(PHANTOM *phm, const char *fname); +void phm_add_pelm (PHANTOM *phm, const int type, const double cx, const double cy, const double u, const double v, const double rot, const double atten); -int pelem_make_points(PHMELEM *obj); -void pelem_make_xform (PHMELEM *obj); -PHMELEM *pelem_alloc(void); +int pelm_make_points(PELM *obj); +void pelm_make_xform (PELM *obj); +PELM *pelm_alloc(void); void calc_arc(double x[], double y[], const int pts, const double xcent, const double ycent, const double r, const double start, const double stop); void calc_ellipse(double x[], double y[], const int pts, const double u, const double v); +int circle_pts(double theta); void phm_print(PHANTOM *phm); #if HAVE_INTERACTIVE_GRAPHICS void phm_show(const PHANTOM *phm); void phm_draw(const PHANTOM *phm); #endif + /* phm2image.c */ void phm_to_image(const PHANTOM *phm, IMAGE *im, const int col_start, const int col_count, const int nsample, const int trace); -int inside_obj(const PHMELEM *obj, double x, double y, const int coord_type); +int pelm_is_point_inside(const PELM *obj, double x, double y, const int coord_type); + /* phmstd.c */ -void herm_head(PHANTOM *phm); -void row_head(PHANTOM *phm); -void row_bord_head(PHANTOM *phm); -/* ray.c */ -DETECTOR *detect_create(const PHANTOM *phm, int ndet, int nview, int nsample, const double rot_anglen); -void detect_free(DETECTOR *det); -double calc_rsum(const PHANTOM *phm, const double x1, const double y1, const double x2, const double y2); -double calc_objsum(const PHMELEM *obj, const double x1, const double y1, const double x2, const double y2); -int clipobj(const PHMELEM *obj, double *x1, double *y1, double *x2, double *y2); +void phm_std_herman (PHANTOM *phm); +void phm_std_rowland (PHANTOM *phm); +void phm_std_rowland_bordered (PHANTOM *phm); + /* raycollect.c */ -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); -void rs_trace_showprm (const char *label, const char *fmt, int row, int color, ...); +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 (const PELM *obj, const double x1, const double y1, const double x2, const double y2); +int pelm_clip_line (const PELM *obj, 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); +void detector_free(DETECTOR *det); + /* rayio.c */ RAYSUM *raysum_create(const char *fname, const int nview, const int ndet); RAYSUM *raysum_create_from_det(const char *fname, const DETECTOR *det); @@ -484,5 +500,4 @@ 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); - #endif