WiscPlan_v2/view_beamlet.m

blet1 = read_ryan_beamlets('F:\021_WiscPlan_data\DELETEME3\FRO_Beam1\scenario1\batch_dose.bin');
blet2 = read_ryan_beamlets('F:\021_WiscPlan_data\DELETEME3\FRO_Beam1\scenario2\batch_dose.bin');


idx = 420;

tabula = zeros([blet1{idx}.x_count, blet1{idx}.y_count, blet1{idx}.z_count]);
tabula(blet1{idx}.non_zero_indices)   = tabula(blet1{idx}.non_zero_indices)+blet1{idx}.non_zero_values;
orthoslice(tabula);

tabula = zeros([blet2{idx}.x_count, blet2{idx}.y_count, blet2{idx}.z_count]);
tabula(blet2{idx}.non_zero_indices)   = tabula(blet2{idx}.non_zero_indices)+blet2{idx}.non_zero_values;
orthoslice(tabula);