ecat63h.c File Reference

Procedures for editing ECAT 6.3 header contents. More...

#include "libtpcimgio.h"

Go to the source code of this file.

Functions
int	ecat63CopyMainheader (ECAT63_mainheader *h1, ECAT63_mainheader *h2)
int	ecat63CopyScanheader (ECAT63_scanheader *h1, ECAT63_scanheader *h2)
int	ecat63EditMHeader (ECAT63_mainheader *h, char *field, char *value, int verbose)

Variables
int	ECAT63_TEST =0
char	ecat63errmsg [128]

Detailed Description

Procedures for editing ECAT 6.3 header contents.

Author

Vesa Oikonen

Definition in file ecat63h.c.

Function Documentation

ecat63CopyMainheader()

int ecat63CopyMainheader	(	ECAT63_mainheader *	h1,
		ECAT63_mainheader *	h2 )

Copy ECAT 6.3 main header from one struct into another.

Returns

Returns 0 if ok.

Parameters

h1	Pointer to source header
h2	Pointer to target header

Definition at line 16 of file ecat63h.c.

  {
  int i;
 
  if(h1==NULL || h2==NULL) return(1);
  memset(h2, 0, sizeof(ECAT63_mainheader));
 
  for(i=0; i<14; i++) h2->ecat_format[i]=h1->ecat_format[i];
  for(i=0; i<14; i++) h2->fill1[i]=h1->fill1[i];
  for(i=0; i<20; i++) h2->original_file_name[i]=h1->original_file_name[i];
  h2->sw_version=h1->sw_version;
  h2->data_type=h1->data_type;
  h2->system_type=h1->system_type;
  h2->file_type=h1->file_type;
  for(i=0; i<10; i++) h2->node_id[i]=h1->node_id[i];
  h2->scan_start_day=h1->scan_start_day;
  h2->scan_start_month=h1->scan_start_month;
  h2->scan_start_year=h1->scan_start_year;
  h2->scan_start_hour=h1->scan_start_hour;
  h2->scan_start_minute=h1->scan_start_minute;
  h2->scan_start_second=h1->scan_start_second;
  for(i=0; i<8; i++) h2->isotope_code[i]=h1->isotope_code[i];
  h2->isotope_halflife=h1->isotope_halflife;
  for(i=0; i<32; i++) h2->radiopharmaceutical[i]=h1->radiopharmaceutical[i];
  h2->gantry_tilt=h1->gantry_tilt;
  h2->gantry_rotation=h1->gantry_rotation;
  h2->bed_elevation=h1->bed_elevation;
  h2->rot_source_speed=h1->rot_source_speed;
  h2->wobble_speed=h1->wobble_speed;
  h2->transm_source_type=h1->transm_source_type;
  h2->axial_fov=h1->axial_fov;
  h2->transaxial_fov=h1->transaxial_fov;
  h2->transaxial_samp_mode=h1->transaxial_samp_mode;
  h2->coin_samp_mode=h1->coin_samp_mode;
  h2->axial_samp_mode=h1->axial_samp_mode;
  h2->calibration_factor=h1->calibration_factor;
  h2->calibration_units=h1->calibration_units;
  h2->compression_code=h1->compression_code;
  for(i=0; i<12; i++) h2->study_name[i]=h1->study_name[i];
  for(i=0; i<16; i++) h2->patient_id[i]=h1->patient_id[i];
  for(i=0; i<32; i++) h2->patient_name[i]=h1->patient_name[i];
  h2->patient_sex=h1->patient_sex;
  for(i=0; i<10; i++) h2->patient_age[i]=h1->patient_age[i];
  for(i=0; i<10; i++) h2->patient_height[i]=h1->patient_height[i];
  for(i=0; i<10; i++) h2->patient_weight[i]=h1->patient_weight[i];
  h2->patient_dexterity=h1->patient_dexterity;
  for(i=0; i<32; i++) h2->physician_name[i]=h1->physician_name[i];
  for(i=0; i<32; i++) h2->operator_name[i]=h1->operator_name[i];
  for(i=0; i<32; i++) h2->study_description[i]=h1->study_description[i];
  h2->acquisition_type=h1->acquisition_type;
  h2->bed_type=h1->bed_type;
  h2->septa_type=h1->septa_type;
  for(i=0; i<20; i++) h2->facility_name[i]=h1->facility_name[i];
  h2->num_planes=h1->num_planes;
  h2->num_frames=h1->num_frames;
  h2->num_gates=h1->num_gates;
  h2->num_bed_pos=h1->num_bed_pos;
  h2->init_bed_position=h1->init_bed_position;
  for(i=0; i<15; i++) h2->bed_offset[i]=h1->bed_offset[i];
  h2->plane_separation=h1->plane_separation;
  h2->lwr_sctr_thres=h1->lwr_sctr_thres;
  h2->lwr_true_thres=h1->lwr_true_thres;
  h2->upr_true_thres=h1->upr_true_thres;
  h2->collimator=h1->collimator;
  for(i=0; i<10; i++) h2->user_process_code[i]=h1->user_process_code[i];
  for(i=0; i<20; i++) h2->fill2[i]=h1->fill2[i];
  return(0);
}

Referenced by atnMake().

ecat63CopyScanheader()

int ecat63CopyScanheader	(	ECAT63_scanheader *	h1,
		ECAT63_scanheader *	h2 )

Copy ECAT 6.3 scan sub header from one struct into another.

Returns

Returns 0 if ok.

Parameters

h1	Pointer to source header
h2	Pointer to target header

Definition at line 94 of file ecat63h.c.

  {
  int i;
 
  if(h1==NULL || h2==NULL) return(1);
  memset(h2, 0, sizeof(ECAT63_scanheader));
 
  for(i=0; i<126; i++) h2->fill1[i]=h1->fill1[i];
  h2->data_type=h1->data_type;
  h2->dimension_1=h1->dimension_1;
  h2->dimension_2=h1->dimension_2;
  h2->smoothing=h1->smoothing;
  h2->processing_code=h1->processing_code;
  h2->sample_distance=h1->sample_distance;
  h2->isotope_halflife=h1->isotope_halflife;
  h2->frame_duration_sec=h1->frame_duration_sec;
  h2->gate_duration=h1->gate_duration;
  h2->r_wave_offset=h1->r_wave_offset;
  h2->scale_factor=h1->scale_factor;
  h2->scan_min=h1->scan_min;
  h2->scan_max=h1->scan_max;
  h2->prompts=h1->prompts;
  h2->delayed=h1->delayed;
  h2->multiples=h1->multiples;
  h2->net_trues=h1->net_trues;
  for(i=0; i<16; i++) h2->cor_singles[i]=h1->cor_singles[i];
  for(i=0; i<16; i++) h2->uncor_singles[i]=h1->uncor_singles[i];
  h2->tot_avg_cor=h1->tot_avg_cor;
  h2->tot_avg_uncor=h1->tot_avg_uncor;
  h2->total_coin_rate=h1->total_coin_rate;
  h2->frame_start_time=h1->frame_start_time;
  h2->frame_duration=h1->frame_duration;
  h2->loss_correction_fctr=h1->loss_correction_fctr;
  for(i=0; i<22; i++) h2->fill2[i]=h1->fill2[i];
  return(0);
}

ecat63EditMHeader()

int ecat63EditMHeader	(	ECAT63_mainheader *	h,
		char *	field,
		char *	value,
		int	verbose )

Edit ECAT 6.3 main header.

Returns

Returns 0, if ok, and 1 or 2, if field name or or value is invalid.

Parameters

h	Pointer to ECAT 6.3 mainheader structure
field	Field name to be changed
value	New value for the field
verbose	Verbose level; if <=0, then nothing is printed into stdout

Definition at line 140 of file ecat63h.c.

  {
//  int yy, mm, dd;
  short int si;
  float f;
 
  if(verbose>0) printf("ecat63EditMHeader('%s', '%s')\n", field, value);
  si=atoi(value); f=atof(value);
 
  if(strcmp(field, "ecat_format")==0 || strcmp(field, "magic_number")==0) {
    strlcpy(h->ecat_format, value, 14);
  } else if(strcmp(field, "fill1")==0) {
    strlcpy(h->fill1, value, 14);
  } else if(strcmp(field, "original_file_name")==0) {
    strlcpy(h->original_file_name, value, 20);
  } else if(strcmp(field, "sw_version")==0) {
    if(si<=0) return(2); else h->sw_version=si;
  } else if(strcmp(field, "data_type")==0) {
    if(si<0) return(2); else h->data_type=si;
  } else if(strcmp(field, "system_type")==0) {
    if(si<0) return(2); else h->system_type=si;
  } else if(strcmp(field, "file_type")==0) {
    if(si<0) return(2); else h->file_type=si;
  } else if(strcmp(field, "node_id")==0 || strcmp(field, "serial_number")==0) {
    strlcpy(h->node_id, value, 10);
  } else if(strcmp(field, "scan_start_day")==0) {
    if(si<0) return(2); else h->scan_start_day=si;
  } else if(strcmp(field, "scan_start_month")==0) {
    if(si<0) return(2); else h->scan_start_month=si;
  } else if(strcmp(field, "scan_start_year")==0) {
    if(si<0) return(2); else h->scan_start_year=si;
  } else if(strcmp(field, "scan_start_hour")==0) {
    if(si<0) return(2); else h->scan_start_hour=si;
  } else if(strcmp(field, "scan_start_minute")==0) {
    if(si<0) return(2); else h->scan_start_minute=si;
  } else if(strcmp(field, "scan_start_second")==0) {
    if(si<0) return(2); else h->scan_start_second=si;
  } else if(strcmp(field, "scan_start_time")==0) {
    int YYYY, MM, DD, hh, mm, ss, n;
    n=sscanf(value, "%d-%d-%d %d:%d:%d", &YYYY, &MM, &DD, &hh, &mm, &ss);
    if(n!=6 && n!=3) return(2);
    h->scan_start_year=YYYY;
    h->scan_start_month=MM;
    h->scan_start_day=DD;
    if(n==6) {
      h->scan_start_hour=hh;
      h->scan_start_minute=mm;
      h->scan_start_second=ss;
    }
  } else if(strcmp(field, "isotope_code")==0 || strcmp(field, "isotope_name")==0) {
    strlcpy(h->isotope_code, value, 8);
  } else if(strcmp(field, "isotope_halflife")==0) {
    if(f<=1.0E-3) return(2); else h->isotope_halflife=f;
  } else if(strcmp(field, "radiopharmaceutical")==0) {
    strlcpy(h->radiopharmaceutical, value, 32);
  } else if(strcmp(field, "gantry_tilt")==0) {
    h->gantry_tilt=f;
  } else if(strcmp(field, "gantry_rotation")==0) {
    h->gantry_rotation=f;
  } else if(strcmp(field, "bed_elevation")==0) {
    h->bed_elevation=f;
  } else if(strcmp(field, "rot_source_speed")==0) {
    h->rot_source_speed=si;
  } else if(strcmp(field, "wobble_speed")==0) {
    h->wobble_speed=si;
  } else if(strcmp(field, "transm_source_type")==0) {
    h->transm_source_type=si;
  } else if(strcmp(field, "axial_fov")==0) {
    h->axial_fov=f;
  } else if(strcmp(field, "transaxial_fov")==0) {
    h->transaxial_fov=f;
  } else if(strcmp(field, "transaxial_samp_mode")==0) {
    h->transaxial_samp_mode=si;
  } else if(strcmp(field, "coin_samp_mode")==0) {
    h->coin_samp_mode=si;
  } else if(strcmp(field, "axial_samp_mode")==0) {
    h->axial_samp_mode=si;
  } else if(strcmp(field, "calibration_factor")==0) {
    h->calibration_factor=f;
  } else if(strcmp(field, "calibration_units")==0) {
    h->calibration_units=si;
  } else if(strcmp(field, "compression_code")==0) {
    h->compression_code=si;
  } else if(strcmp(field, "study_name")==0) {
    strlcpy(h->study_name, value, 12);
  } else if(strcmp(field, "patient_id")==0) {
    strlcpy(h->patient_id, value, 16);
  } else if(strcmp(field, "patient_name")==0) {
    strlcpy(h->patient_name, value, 32);
  } else if(strcmp(field, "patient_sex")==0) {
    h->patient_sex=value[0];
  } else if(strcmp(field, "patient_age")==0) {
    strlcpy(h->patient_name, value, 10);
  } else if(strcmp(field, "patient_height")==0) {
    strlcpy(h->patient_height, value, 10);
  } else if(strcmp(field, "patient_weight")==0) {
    strlcpy(h->patient_weight, value, 10);
  } else if(strcmp(field, "patient_dexterity")==0) {
    h->patient_dexterity=value[0];
  } else if(strcmp(field, "physician_name")==0) {
    strlcpy(h->physician_name, value, 32);
  } else if(strcmp(field, "operator_name")==0) {
    strlcpy(h->operator_name, value, 32);
  } else if(strcmp(field, "study_description")==0) {
    strlcpy(h->study_description, value, 32);
  } else if(strcmp(field, "acquisition_type")==0) {
    h->acquisition_type=si;
  } else if(strcmp(field, "bed_type")==0) {
    h->bed_type=si;
  } else if(strcmp(field, "septa_type")==0) {
    h->septa_type=si;
  } else if(strcmp(field, "facility_name")==0) {
    strlcpy(h->facility_name, value, 20);
  } else if(strcmp(field, "num_planes")==0) {
    h->num_planes=si;
  } else if(strcmp(field, "num_frames")==0) {
    h->num_frames=si;
  } else if(strcmp(field, "num_gates")==0) {
    h->num_gates=si;
  } else if(strcmp(field, "num_bed_pos")==0) {
    h->num_bed_pos=si;
  } else if(strcmp(field, "init_bed_position")==0) {
    h->init_bed_position=f;
  } else if(strcmp(field, "bed_offset")==0) {
    sscanf(value, "%f %f %f %f %f %f %f %f %f %f %f %f %f %f %f",
      h->bed_offset+0, h->bed_offset+1, h->bed_offset+2,
      h->bed_offset+3, h->bed_offset+4, h->bed_offset+5,
      h->bed_offset+6, h->bed_offset+7, h->bed_offset+8,
      h->bed_offset+9, h->bed_offset+10, h->bed_offset+11,
      h->bed_offset+12, h->bed_offset+13, h->bed_offset+14
    );
  } else if(strcmp(field, "plane_separation")==0) {
    h->plane_separation=f;
  } else if(strcmp(field, "lwr_sctr_thres")==0) {
    h->lwr_sctr_thres=si;
  } else if(strcmp(field, "lwr_true_thres")==0) {
    h->lwr_true_thres=si;
  } else if(strcmp(field, "upr_true_thres")==0) {
    h->upr_true_thres=si;
  } else if(strcmp(field, "collimator")==0) {
    h->collimator=f;
  } else if(strcmp(field, "user_process_code")==0) {
    strlcpy(h->user_process_code, value, 10);
  } else if(strcasecmp(field, "FILL2")==0) {
    char *cptr;
    cptr=strtok(value, " \t,;\n\r");
    for(int i=0; i<20; i++) {
      if(cptr==NULL) break;
      h->fill2[i]=(short int)atoi(cptr);
      cptr=strtok(NULL, " \t,;\n\r");
    }
  } else
    return(1);
 
  return(0);
}

Variable Documentation

ECAT63_TEST

int ECAT63_TEST =0

Verbose prints from ECAT 6.3 functions

Definition at line 6 of file ecat63h.c.

Referenced by ecat63Create(), ecat63Matenter(), ecat63PrintAttnheader(), ecat63PrintImageheader(), ecat63PrintMainheader(), ecat63PrintNormheader(), ecat63PrintScanheader(), ecat63ReadAllToImg(), ecat63ReadAttnheader(), ecat63ReadAttnMatrix(), ecat63ReadImageMatrix(), ecat63ReadMainheader(), ecat63ReadMatdata(), ecat63ReadNormheader(), ecat63ReadScanheader(), ecat63ReadScanMatrix(), ecat63WriteAttn(), ecat63WriteAttnheader(), ecat63WriteImage(), ecat63WriteImageheader(), ecat63WriteImageMatrix(), ecat63WriteMainheader(), ecat63WriteMatdata(), ecat63WriteNorm(), ecat63WriteNormheader(), ecat63WriteScan(), ecat63WriteScanheader(), ecat63WriteScanMatrix(), ecat6PrintSubheader(), and imgWriteEcat63Frame().

ecat63errmsg

char ecat63errmsg[128]

Error message from ECAT 6.3 functions

Definition at line 7 of file ecat63h.c.

Referenced by ecat63AddImg(), ecat63ReadAllToImg(), ecat63ReadAttnMatrix(), ecat63ReadImageMatrix(), ecat63ReadPlaneToImg(), ecat63ReadScanMatrix(), ecat63WriteAllImg(), ecat63WriteImageMatrix(), ecat63WriteScanMatrix(), imgRead(), and imgWriteEcat63Frame().