1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
|
#include <stdio.h>
#include <stdlib.h>
typedef enum { ECAT6, ECAT7, Interfile } FileFormat;
typedef enum {
MAT_OK,
MAT_READ_ERROR,
MAT_WRITE_ERROR,
MAT_INVALID_DIRBLK,
MAT_ACS_FILE_NOT_FOUND,
MAT_INTERFILE_OPEN_ERR,
MAT_FILE_TYPE_NOT_MATCH,
MAT_READ_FROM_NILFPTR,
MAT_NOMHD_FILE_OBJECT,
MAT_NIL_SHPTR,
MAT_NIL_DATA_PTR,
MAT_MATRIX_NOT_FOUND,
MAT_UNKNOWN_FILE_TYPE,
MAT_ACS_CREATE_ERR,
MAT_BAD_ATTRIBUTE,
MAT_BAD_FILE_ACCESS_MODE,
MAT_INVALID_DIMENSION,
MAT_NO_SLICES_FOUND,
MAT_INVALID_DATA_TYPE,
MAT_INVALID_MBED_POSITION
} MatrixErrorCode;
typedef enum {
NoData, Sinogram, PetImage, AttenCor, Normalization,
PolarMap, ByteVolume, PetVolume, ByteProjection,
PetProjection, ByteImage, Short3dSinogram, Byte3dSinogram, Norm3d,
Float3dSinogram,InterfileImage, NumDataSetTypes
} DataSetType;
typedef enum {
UnknownMatDataType, ByteData, VAX_Ix2, VAX_Ix4,
VAX_Rx4, IeeeFloat, SunShort, SunLong, NumMatrixDataTypes, ColorData,
BitData
} MatrixDataType;
MatrixErrorCode matrix_errno;
char matrix_errtxt[132];
typedef struct XMAIN_HEAD {
char magic_number[14];
char original_file_name[32];
short sw_version;
short system_type;
short file_type;
char serial_number[10];
short align_0; /* 4 byte alignment purpose */
unsigned int scan_start_time;
char isotope_code[8];
float isotope_halflife;
char radiopharmaceutical[32];
float gantry_tilt;
float gantry_rotation;
float bed_elevation;
float intrinsic_tilt;
short wobble_speed;
short transm_source_type;
float distance_scanned;
float transaxial_fov;
short angular_compression;
short coin_samp_mode;
short axial_samp_mode;
short align_1;
float calibration_factor;
short calibration_units;
short calibration_units_label;
short compression_code;
char study_name[12];
char patient_id[16];
char patient_name[32];
char patient_sex[1];
char patient_dexterity[1];
float patient_age;
float patient_height;
float patient_weight;
int patient_birth_date;
char physician_name[32];
char operator_name[32];
char study_description[32];
short acquisition_type;
short patient_orientation;
char facility_name[20];
short num_planes;
short num_frames;
short num_gates;
short num_bed_pos;
float init_bed_position;
float bed_offset[15];
float plane_separation;
short lwr_sctr_thres;
short lwr_true_thres;
short upr_true_thres;
char user_process_code[10];
short acquisition_mode;
short align_2;
float bin_size;
float branching_fraction;
unsigned int dose_start_time;
float dosage;
float well_counter_factor;
char data_units[32];
short septa_state;
short align_3;
} Main_header;
typedef struct XIMAGE_SUB {
short data_type;
short num_dimensions;
short x_dimension;
short y_dimension;
short z_dimension;
short align_0;
float z_offset;
float x_offset;
float y_offset;
float recon_zoom;
float scale_factor;
short image_min;
short image_max;
float x_pixel_size;
float y_pixel_size;
float z_pixel_size;
unsigned int frame_duration;
unsigned int frame_start_time;
short filter_code;
short align_1;
float x_resolution;
float y_resolution;
float z_resolution;
float num_r_elements;
float num_angles;
float z_rotation_angle;
float decay_corr_fctr;
int processing_code;
unsigned int gate_duration;
int r_wave_offset;
int num_accepted_beats;
float filter_cutoff_frequency;
float filter_resolution;
float filter_ramp_slope;
short filter_order;
short align_2;
float filter_scatter_fraction;
float filter_scatter_slope;
char annotation[40];
float mt_1_1;
float mt_1_2;
float mt_1_3;
float mt_2_1;
float mt_2_2;
float mt_2_3;
float mt_3_1;
float mt_3_2;
float mt_3_3;
float rfilter_cutoff;
float rfilter_resolution;
short rfilter_code;
short rfilter_order;
float zfilter_cutoff;
float zfilter_resolution;
short zfilter_code;
short zfilter_order;
float mt_1_4;
float mt_2_4;
float mt_3_4;
short scatter_type;
short recon_type;
short recon_views;
short align_3;
} Image_subheader;
typedef struct matdirnode {
int matnum ;
int strtblk ;
int endblk ;
int matstat ;
struct matdirnode *next ;
} MatDirNode ;
typedef struct matdirlist {
int nmats ;
MatDirNode *first ;
MatDirNode *last ;
} MatDirList ;
typedef struct matrix_file {
char *fname ;
Main_header *mhptr ;
MatDirList *dirlist ;
FILE *fptr ;
int acs ;
FileFormat file_format;
char **interfile_header;
} MatrixFile;
typedef struct matrixdata {
int matnum ; /* matrix number */
MatrixFile *matfile ; /* pointer to parent */
DataSetType mat_type ; /* type of matrix? */
MatrixDataType data_type ; /* type of data */
void * shptr ; /* pointer to sub-header */
void * data_ptr ; /* pointer to data */
int data_size ; /* size of data in bytes */
int xdim; /* dimensions of data */
int ydim; /* y dimension */
int zdim; /* for volumes */
float scale_factor ; /* valid if data is int? */
float pixel_size; /* xdim data spacing (cm) */
float y_size; /* ydim data spacing (cm) */
float z_size; /* zdim data spacing (cm) */
float data_min; /* min value of data */
float data_max; /* max value of data */
float x_origin; /* x origin of data */
float y_origin; /* y origin of data */
float z_origin; /* z origin of data */
} MatrixData ;
MatrixFile *matrix_create(const char *fname, Main_header * proto_mhptr);
int matrix_write(MatrixFile *mptr, int matnum, MatrixData *data);
int matrix_close(MatrixFile *mptr);
int mat_numcod(int frame, int plane, int gate, int data, int bed);
|
