OrientationSampling: made it less restrictive by default (to avoid problems with bad grains)

Signed-off-by: Nicola Vigano <nicola.vigano@esrf.fr>

zUtil_Deformation/GtOrientationSampling.m

@@ -310,7 +310,11 @@ classdef GtOrientationSampling < handle
     end
 
     methods (Access = protected)
-        function gvdm_tmp = guess_ODF_BB(self)
+        function gvdm_tmp = guess_ODF_BB(self, strict)
+            if (~exist('strict', 'var') || isempty(strict))
+                strict = false;
+            end
+
             omega_step = 180 / self.parameters.acq.nproj;
 
             bls = self.bl(self.selected);
@@ -456,8 +460,14 @@ classdef GtOrientationSampling < handle
             trasl_ex_points_n = sum(trasl_ex_points_n .* rws_ex_w, 2) ./ sqrt(sum(rws_ex_w .^ 2, 2));
             trasl_ex_points_n = r01s_ex_w + trasl_ex_points_n(:, [1 1 1]) .* rws_ex_w;
 
-            all_plane_normals = [trasl_ex_points_w - avg_R_vec_exp; trasl_ex_points_n - avg_R_vec_exp];
-%             all_plane_normals = trasl_ex_points_w - avg_R_vec_exp;
+            all_plane_normals = [ ...
+                trasl_ex_points_w - avg_R_vec_exp; ...
+                trasl_ex_points_n - avg_R_vec_exp ];
+
+            if (~strict)
+                sel_points = convhull(all_plane_normals);
+                all_plane_normals = all_plane_normals(sel_points);
+            end
 
             % Computing the circumscribing polyhedron in Orientation-space
             verts = gtMathsGetPolyhedronVerticesFromPlaneNormals(all_plane_normals);