diff --git a/zUtil_Deformation/gt6DCreateProjDataFromGrainCluster.m b/zUtil_Deformation/gt6DCreateProjDataFromGrainCluster.m index 60ec8b01499aa8bac6f7ae4ea021ccacbdb36312..125a7269e4914ef05af1b692f6859be5705f3921 100644 --- a/zUtil_Deformation/gt6DCreateProjDataFromGrainCluster.m +++ b/zUtil_Deformation/gt6DCreateProjDataFromGrainCluster.m @@ -235,6 +235,12 @@ function [refor, estim_space_bbox_pix, estim_orient_bbox_rod] = gt6DCreateProjDa volume_verts = cat(1, volume_verts, ... bsxfun(@plus, grs(ii_g).fwd_sim.gv_verts, center_shift)); end + k = convhull(volume_verts); + volume_verts = volume_verts(k, :); + center_verts = max(volume_verts, [], 1) - min(volume_verts, [], 1); + volume_verts = bsxfun(@minus, volume_verts, center_verts); + volume_verts = bsxfun(@times, volume_verts, conf.ospace_oversize); + volume_verts = bsxfun(@plus, volume_verts, center_verts); proj_bl_masks = gt6DSpreadProjectVertices2Det(refor, volume_verts, fwdsim_sel_refl(~inconvenient_etas), p, conf.det_index); blobs = assign_masks(blobs, proj_bl_masks, sel_reflections, true);