%This code takes the 3D dataset for the brain images, 
%and generates a mask. The data should be stored as ims
mag = abs(ims);
mmPerVox = [1.4 1.4 1.4];
imName = 'qt';
hdr = saveAnalyze(mag, imName, mmPerVox);

betLevel = 0.5;
brainMask = mrAnatExtractBrain(imName, mmPerVox, betLevel);


