enabling psf from parameters file for 3DTV reconstruction of grains

Signed-off-by: Wolfgang Ludwig <wolfgang.ludwig@esrf.fr>

5_reconstruction/gtAstraCreateProjections.m

-function proj = gtAstraCreateProjections(volume, proj_geom, proj_size, do_permute, proj_model)
+function proj = gtAstraCreateProjections(volume, proj_geom, proj_size, do_permute, proj_model, psf)
 % GTASTRACREATEPROJECTIONS  Calculate oblique angle projections of a volume
 %     proj = gtAstraCreateProjections(volume, proj_geom, proj_size, [do_permute], [proj_model])
 %     -----------------------------------------------------------------------------------------
@@ -26,10 +26,10 @@ if (~OK)
     throw(mexc)
 end
 
-if ~exist('do_permute','var') || isempty(do_permute)
+if (~exist('do_permute','var') || isempty(do_permute))
     do_permute = false;
 end
-if ~exist('proj_model','var') || isempty(proj_model)
+if (~exist('proj_model','var') || isempty(proj_model))
    proj_model = 'parallel3d_vec';
    disp(['using ' proj_model ' as projection model'])
 end
@@ -53,7 +53,13 @@ vol_geom = astra_create_vol_geom(yvol, xvol, zvol);
 proj_id = astra_create_sino3d_cuda(volume, proj_geom, vol_geom);
 proj    = astra_mex_data3d('get', proj_id);
 
-if do_permute
+if (exist('psf','var') && ~isempty(psf))
+    psf_obj = GtPSF();
+    psf_obj.set_psf_direct(proj.psf, [hsize, vsize]);
+    proj = psf_obj.apply_psf_direct(proj);
+end
+
+if (do_permute)
     proj = permute(proj, [3 1 2]);
     proj = reshape(proj, vsize, hsize, 1, nproj);
 end

5_reconstruction/gtAstraReconstructGrain.m

@@ -36,6 +36,11 @@ function varargout = gtAstraReconstructGrain(grainID, phaseID, parameters)
     gr.proj.stack = gr.proj.stack(:, gr.proj.selected, :);
     gr.proj.geom  = gr.proj.geom(gr.proj.selected, :);
 
+    if (isfield(parameters.rec.grains.options, 'psf') ...
+            && ~isempty(parameters.rec.grains.options.psf))
+        gr.proj.psf = parameters.rec.grains.options.psf;
+    end
+
     is_cone = isfield(parameters.labgeo, 'sourcepoint') ...
         && (~isempty(parameters.labgeo.sourcepoint));