[CCPPETMR/Hackathon-STIR] Reading HDF5 accurately (#2)

[pwadhwa351]

—
You are receiving this because you are subscribed to this thread.
Reply to this email directly, view it on GitHub, or mute the thread.

[Kris Thielemans]

assuming that GE always stores segments as 0, -1, +1, ...
schematically.

std::vector<int> get_ax_pos_offset(const ProjDataInfo&)
{
  //std::vector<int> segment_sequence_in_stream(2*max_segment_num+1);
  std::vector<int> offset(min_segment_num,  max_segment_num);
  //segment_sequence_in_stream[0] = 0;
  int last_offset = get_num_axial_poss(segment);
  offset[0] = 0;
  for (int segment_num = 1; segment_num<=max_segment_num; ++segment_num)
  {
    //segment_sequence_in_stream[2*segment_num-1] = -segment_num;
    //segment_sequence_in_stream[2*segment_num] = segment_num;
    offset[-segment_num] = last_offset;
    last_offset += get_num_axial_poss(-segment_num);
    offset[segment_num] = last_offset;
    last_offset += get_num_axial_poss(segment_num); 
    }
  return offset;
}

[pwadhwa351]

Closed #2.