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

6D-Reconstruction/Twin: added possibility an even sampling of the grid and use...

6D-Reconstruction/Twin: added possibility an even sampling of the grid and use the previously computed orientation bboxes

Signed-off-by: default avatarNicola Vigano <nicola.vigano@esrf.fr>
parent dbe59b6d
No related branches found
No related tags found
No related merge requests found
...@@ -7,6 +7,9 @@ function gtReconstructGrainTwinCluster(grain_ids, phase_id, parameters, varargin ...@@ -7,6 +7,9 @@ function gtReconstructGrainTwinCluster(grain_ids, phase_id, parameters, varargin
end end
parameters.fsim.check_spot = false; parameters.fsim.check_spot = false;
conf = struct('recompute_osbb', true);
[conf, ~] = parse_pv_pairs(conf, varargin);
rec_opts = gtReconstruct6DGetParamenters(parameters); rec_opts = gtReconstruct6DGetParamenters(parameters);
if (isfield(parameters.rec, 'grains') ... if (isfield(parameters.rec, 'grains') ...
&& isfield(parameters.rec.grains, 'options') ... && isfield(parameters.rec.grains, 'options') ...
...@@ -29,10 +32,14 @@ function gtReconstructGrainTwinCluster(grain_ids, phase_id, parameters, varargin ...@@ -29,10 +32,14 @@ function gtReconstructGrainTwinCluster(grain_ids, phase_id, parameters, varargin
num_grains = numel(grs); num_grains = numel(grs);
for ii_g = 1:num_grains for ii_g = 1:num_grains
bb_gvdm = cat(1, grs(ii_g).bb_ors(:).R_vector)';
sampler(ii_g) = GtOrientationSampling(parameters, grs(ii_g)); %#ok<AGROW> sampler(ii_g) = GtOrientationSampling(parameters, grs(ii_g)); %#ok<AGROW>
sampler(ii_g).make_even_simple_grid('cubic', rec_opts.grid_edge, bb_gvdm, rec_opts.ospace_oversize)
if (conf.recompute_osbb)
sampler(ii_g).make_simple_grid_estim_ODF('cubic', rec_opts.grid_edge, false, rec_opts.ospace_oversize);
else
bb_gvdm = cat(1, grs(ii_g).bb_ors(:).R_vector)';
sampler(ii_g).make_even_simple_grid('cubic', rec_opts.grid_edge, bb_gvdm, rec_opts.ospace_oversize)
end
if (rec_opts.ospace_super_sampling > 1) if (rec_opts.ospace_super_sampling > 1)
sampler(ii_g).make_supersampling_simple_grid([1 2 3], rec_opts.ospace_super_sampling); sampler(ii_g).make_supersampling_simple_grid([1 2 3], rec_opts.ospace_super_sampling);
end end
......
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