nii_reslice_target.m

function [outhdr, outimg] = nii_reslice_target(inhdr, inimg, tarhdr, interp) 
%Reslice input image to match dimensions of target image (either to disk or memory)
%  inhdr: image to reslice- either filename of NIfTI header or loaded NIfTI header structure
%  inimg: (optonal) NIfTI image data (only if inhdr is a structure)
%  tarhdr: image to match- either filename of NIfTI header or loaded NIfTI header structure
%  interp: (optional) if 0 nearest neighbor interpolation, else trilinear interpolation (default)
% Outputs: if not specified, resliced image saved to disk, else returns resliced header and image
%Chris Rorden (2014) see John Ashburner and Ged Ridgway's reorient.m 
% http://opensource.org/licenses/BSD-2-Clause
%Examples
%  nii_reslice_target('T1_LM1001.nii','','jhu.nii',0);
% Next example shows how to use this without writing to disk
%  inhdr = spm_vol('T1_LM1001.nii'); %load header
%  inimg = spm_read_vols(inhdr); %load volume
%  tarhdr = spm_vol('jhu.nii'); %load header
%  [outhdr,outimg] = nii_reslice_target(inhdr,inimg, tarhdr); %resize in memory
%  spm_write_vol(outhdr,outimg); %save resized image
% Next example: specify target dimensions without explicitly providing a full header
%  tarhdr.mat = [-1 0 0 79; 0 1 0 -113; 0 0 1 -51; 0 0 0 1];
%  tarhdr.dim = [157 189 136];
%  nii_reslice_target('qCBV.nii','',tarhdr);
% Next example: rescale for maximum resolution (in voxels)
%  fnm = 'single_subj_T1.nii';
%  outres = 256; %rescale so maximum dimension has 256 voxels
%  tarhdr = spm_vol(fnm);
%  mx = max(tarhdr.dim(:));
%  scale = outres/mx;
%  tarhdr.dim = floor(tarhdr.dim * scale);
%  tarhdr.mat(1:3,1:3) = tarhdr.mat(1:3,1:3)/scale;
%  nii_reslice_target(fnm,'',tarhdr);

if ~exist('inhdr','var')
    inhdr = spm_select(1,'image','Select source image that will be resliced');
end
if ~exist('tarhdr','var')
    tarhdr = spm_select(1,'image','Select target image (source will be resliced to match target)');
end
if ~exist('interp','var')
    interp = 1;%linear
end
if ~isnumeric(interp), interp = interp + 0; end; %change false/true to 0/1
if ~isstruct(tarhdr)
    tarhdr = spm_vol(tarhdr); %load target header
end
imgdim = tarhdr.dim(1:3);
if ~isstruct(inhdr)
    inhdr = spm_vol(inhdr); %load input header
    %if numel(inhdr) > 1
    %    fprintf('%s warning: only reslicing first volume of 4D image %s\n', mfilename, inhdr(1).fname);
    %end
    %inhdr = inhdr(1); %if 4D, only process first volume
	inimg = spm_read_vols(inhdr); %load input image
end
nvol = size(inimg,4);
inhdr = inhdr(1);
outhdr            = inhdr;
[pth,nam,ext] = fileparts(outhdr.fname);
outhdr.fname      = fullfile(pth,['r' nam ext]);
outhdr.dim(1:3)   = imgdim(1:3);
%outhdr.dim(4) = nvol;
outhdr.mat        = tarhdr.mat;
if isequal(inhdr.mat, outhdr.mat) && isequal(outhdr.dim(1:3), inhdr.dim(1:3))
    fprintf('%s no need to reslice: input image already aligned to target image\n', mfilename);
    outimg = inimg;
else %if reslicing is required
    outimg = zeros([outhdr.dim(1:3), nvol]);
    for vol = 1:nvol 
        for i = 1:imgdim(3)
            M = inv(spm_matrix([0 0 -i])*inv(outhdr.mat)*inhdr.mat);
            outimg(:,:,i,vol) = spm_slice_vol(inimg(:,:,:,vol), M, imgdim(1:2), interp); % (linear interp)
        end
    end
end
if nargout < 2
    for vol=1:nvol
        outhdr.n(1)=vol;
        spm_write_vol(outhdr,outimg(:, :, :, vol));
    end;
    %outhdr = spm_create_vol(outhdr); %save header to disk
    %spm_write_vol(outhdr,outimg); %save image to disk
end
%end nii_reslice_target()