view_beamlet.m 509 B

123456789101112131415
  1. % data = read_ryan_beamlets('F:\021_WiscPlan_data\Hippocampus_HD_2\batch_dose.bin')
  2. idx = 3000;
  3. tabula = zeros([data{idx}.x_count, data{idx}.y_count, data{idx}.z_count]);
  4. tabula(data{idx}.non_zero_indices) = tabula(data{idx}.non_zero_indices)+data{idx}.non_zero_values;
  5. tabula(data{idx+1}.non_zero_indices) = tabula(data{idx+1}.non_zero_indices)+data{idx+1}.non_zero_values;
  6. tabula(data{idx+2}.non_zero_indices) = tabula(data{idx+2}.non_zero_indices)+data{idx+2}.non_zero_values;
  7. orthoslice(tabula);