qMRinfo('vfa_t1'); % Describe the model
Model = vfa_t1;
The current
Model
is an instance ofvfa_t1
class.
You can manually download the sample data for vfa_t1
by clicking here.
dataDir = downloadData(Model,pwd);
Protocol is set according to the example data
% |- vfa_t1 object needs 1 protocol field(s) to be assigned: % |- VFAData
FlipAngle = [3.0000; 20.0000];
% FlipAngle is a vector of [2X1]
TR = [0.0150; 0.0150];
% TR is a vector of [2X1]
Model.Prot.VFAData.Mat = [ FlipAngle TR];
%%
% VFAData.nii.gz contains [128 128 1 2] data.
data.VFAData=double(load_nii_data('vfa_t1_data/VFAData.nii.gz'));
% B1map.nii.gz contains [128 128] data.
data.B1map=double(load_nii_data('vfa_t1_data/B1map.nii.gz'));
% Mask.nii.gz contains [128 128] data.
data.Mask=double(load_nii_data('vfa_t1_data/Mask.nii.gz'));
This section will fit data.
FitResults = FitData(data,Model,0);
- Output map will be displayed.
- If available, a graph will be displayed to show fitting in a voxel.
qMRshowOutput(FitResults,data,Model);
- qMR maps are saved in NIFTI and in a structure
FitResults.mat
that can be loaded in qMRLab graphical user interface.- Model object stores all the options and protocol
- These objects can be easily shared or be used for simulation.
FitResultsSave_nii(FitResults, 'vfa_t1_data/VFAData.nii.gz');