Function calculateSos

Computes combined image for all coils calculated according to Tran-Gia et. al. (2014)

Contents

Input

Output

Copyrights

(C) All rights reserved.

The code may be used free of charge for non-commercial and educational purposes, the only requirement is that this text is preserved within the derivative work. For any other purpose you must contact the authors for permission. This code may not be redistributed without written permission from the authors.

ABOUT: This software implements basic functionalities of the FIR-MAP algorithm

IMPORTANT: If you use this software you should cite the following in any resulting publication: [1] Michal Staniszewski and Uwe Klose. Improvements of Fast Model-based Acceleration of Parameter Look-Locker T1 Mapping

function M_sos = calculateSos(meanPhase, cons_model_coil)

    nc = size(cons_model_coil,4); % number of coils
    np = size(cons_model_coil,1); % number of projections
    nr = size(cons_model_coil,2); % number of image rows
    O_nc=zeros(np,nr,nr);

    % iterate for all coils
    for c=1:nc
        % build imspace for each coil
        fprintf('Calculate SoS for coil#%d\n',c)
        M_Re=cons_model_coil(:,:,:,c);
        temp_meanPhase = meanPhase(c,:,:);
        temp_prod = M_Re.*exp(-1i*temp_meanPhase);
        M_Re = real(temp_prod); % (6) Tran-Gia et. al. (2014)
        % combine all coils
        O_nc=O_nc+sign(M_Re).*abs(M_Re).^2; % (8) Tran-Gia et. al. (2014)
    end
    M_sos=sign(O_nc).*sqrt(abs(O_nc)); % (7) Tran-Gia et. al. (2014)

end