qMRinfo('mt_sat'); % Describe the model
Model = mt_sat;
The current
Model
is an instance ofmt_sat
class.
You can manually download the sample data for mt_sat
by clicking here.
dataDir = downloadData(Model,pwd);
Protocol is set according to the example data
% |- mt_sat object needs 3 protocol field(s) to be assigned: % |- MTw % |- T1w % |- PDw
FlipAngle = 6;
TR = 0.028;
Model.Prot.MTw.Mat = [ FlipAngle TR ];
%%
FlipAngle = 20;
TR = 0.018;
Model.Prot.T1w.Mat = [ FlipAngle TR];
%%
FlipAngle = 6;
TR = 0.028;
Model.Prot.PDw.Mat = [ FlipAngle TR];
%%
% MTw.nii.gz contains [128 128 96] data.
data.MTw=double(load_nii_data('mt_sat_data/MTw.nii.gz'));
% T1w.nii.gz contains [128 128 96] data.
data.T1w=double(load_nii_data('mt_sat_data/T1w.nii.gz'));
% PDw.nii.gz contains [128 128 96] data.
data.PDw=double(load_nii_data('mt_sat_data/PDw.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, 'mt_sat_data/MTw.nii.gz');