Skip to content
Snippets Groups Projects
Commit 0a1f8e6c authored by Nicola Vigano's avatar Nicola Vigano
Browse files

OrientationSampling: fixed bug related to the switch from included to selected

parent 4df1fc09
No related branches found
No related tags found
No related merge requests found
......@@ -639,7 +639,8 @@ classdef GtOrientationSampling < handle
ab = [grs(:).allblobs];
w_tab = [ab(:).omega] / gtGetOmegaStepDeg(self.parameters, self.detector_index);
ref_int_ws = self.ref_gr.allblobs.detector(self.detector_index).uvw(self.ondet(self.included), 3);
sel = self.ondet(self.included(self.selected));
ref_int_ws = self.ref_gr.allblobs.detector(self.detector_index).uvw(sel, 3);
w_tab = gtGrainAnglesTabularFix360deg(w_tab, ref_int_ws, self.parameters);
......@@ -651,7 +652,7 @@ classdef GtOrientationSampling < handle
function [delta_uvs, proj_uv_lims] = get_uv_deviations(self)
grs = [self.lattice(1).gr{:}];
num_blobs = numel(self.selected);
num_blobs = numel(find(self.selected));
num_oris = numel(grs);
uv_tab = zeros(num_blobs, 2, num_oris);
for ii_o = 1:num_oris
......@@ -1076,7 +1077,7 @@ classdef GtOrientationSampling < handle
for y_ii = 1:size(self.lattice(o_ii).gr, 2)
for x_ii = 1:size(self.lattice(o_ii).gr, 1)
ab = self.lattice(o_ii).gr{x_ii, y_ii, z_ii}.allblobs;
omegas = ab.omega(self.selected) / omega_step;
omegas = ab.omega / omega_step;
within_lims = find(omegas >= limits(:, 1) & omegas <= limits(:, 2));
pos_in_blobs = round(omegas(within_lims)) - limits(within_lims, 1) + 1;
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment