gtSetupReconstruction: added possibility to only reconstruct a set of grains

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

5_reconstruction/gtSetupReconstruction.m

-function gtSetupReconstruction(phaseID)
+function gtSetupReconstruction(phaseID, first, last)
 % GTSETUPRECONSTRUCTION  Simple helper function to launch reconstruction and segmentation of
 %                        absorption volume - as well as reconstruction and segmentation of grain volumes
 %
@@ -108,7 +108,8 @@ if ~exist(absorption_mask,'file')
             abs = load(absorption_volume, 'abs');
             abs = abs.abs;
         end
-        disp('Please create file volume_mask.mat from absorption volume (right click in the volume to select a seed point and then press ''Segment'')');
+        disp([ 'Please create file volume_mask.mat from absorption volume ' ...
+            '(right click in the volume to select a seed point and then press ''Segment'')' ]);
         GtGuiThresholdVolume(abs);
     end
 else
@@ -123,9 +124,12 @@ else
     end
 end
 
-
-first = 1;
-last = sample.phases{phaseID}.getNumberOfGrains();
+if (~exist('first', 'var') || isempty(first))
+    first = 1;
+end
+if (~exist('last', 'var') || isempty(last))
+    last = sample.phases{phaseID}.getNumberOfGrains();
+end
 
 check = inputwdefault('Launch reconstruction? [y/n]', 'y');
 rec_on_OAR = false;