Matlab Code for Reducing FA: Difference between revisions

From Psych 221 Image Systems Engineering
Jump to navigation Jump to search
Newer edit →
imported>Psych 204
Created page with 'Return to function alteredFAmasked = AlterFA(dt6, mask, densitymap, effectsize) % input: % dt6: ie dt.dt6 (can use dt=dtiLoadDt6('dt6.mat')) % mask: name of brain mask, ie …'
 
imported>Psych 204
No edit summary
Line 1: Line 1:
Return to  
Return to [[Christine McLeavey & Jessica Tsang]]
function alteredFAmasked = AlterFA(dt6, mask, densitymap, effectsize)
<br>
<br>function alteredFAmasked = AlterFA(dt6, mask, densitymap, effectsize)
<br>
<br>% input:
<br>%  dt6: ie dt.dt6 (can use dt=dtiLoadDt6('dt6.mat'))
<br>%  mask: name of brain mask, ie 'ab_mask.nii.gz'
<br>%  densitymap: image of fiber density (betw 0 and 1), i.e. using
<br>%      dtiComputeFiberDensityNOGUI
<br>%  effectsize: the factor by which to lower/raise FA (to be weighted by
<br>%      fiber density). With effectsize of 0.1, weighted effectsize =
<br>%      0.1*fiberdensity, and newFA = origFA*[1-(weighted effectsize)]
<br>%
<br>% returns:
<br>%  XxYxZ image with altered fa, in image space
<br>%
<br>% 2009.11.15 JT
<br>
<br>fa = dtiComputeFA(dt6);
<br>brainMask=readFileNifti(mask);
<br>brainMask.data=double(brainMask.data);
<br>% effectsize = .1;


% input:
<br>weightedESmap = effectsize*densitymap;
%  dt6: ie dt.dt6 (can use dt=dtiLoadDt6('dt6.mat'))
<br>alteredFA = fa.*(1-weightedESmap);
%  mask: name of brain mask, ie 'ab_mask.nii.gz'
<br>alteredFAmasked = alteredFA.*brainMask.data;
%  densitymap: image of fiber density (betw 0 and 1), i.e. using
<br>
%      dtiComputeFiberDensityNOGUI
<br>% fa = dtiComputeFA(dt6);
%  effectsize: the factor by which to lower/raise FA (to be weighted by
<br>% weightedESmap = effectsize*densitymap;
%      fiber density). With effectsize of 0.1, weighted effectsize =
<br>% alteredFA = fa.*(1-weightedESmap);
%      0.1*fiberdensity, and newFA = origFA*[1-(weighted effectsize)]
<br>
%
<br>return
% returns:
%  XxYxZ image with altered fa, in image space
%
% 2009.11.15 JT
 
fa = dtiComputeFA(dt6);
brainMask=readFileNifti(mask);
brainMask.data=double(brainMask.data);
% effectsize = .1;
 
weightedESmap = effectsize*densitymap;
alteredFA = fa.*(1-weightedESmap);
alteredFAmasked = alteredFA.*brainMask.data;
 
% fa = dtiComputeFA(dt6);
% weightedESmap = effectsize*densitymap;
% alteredFA = fa.*(1-weightedESmap);
 
return

Revision as of 20:36, 9 December 2009

Return to Christine McLeavey & Jessica Tsang

function alteredFAmasked = AlterFA(dt6, mask, densitymap, effectsize)

% input:
% dt6: ie dt.dt6 (can use dt=dtiLoadDt6('dt6.mat'))
% mask: name of brain mask, ie 'ab_mask.nii.gz'
% densitymap: image of fiber density (betw 0 and 1), i.e. using
% dtiComputeFiberDensityNOGUI
% effectsize: the factor by which to lower/raise FA (to be weighted by
% fiber density). With effectsize of 0.1, weighted effectsize =
% 0.1*fiberdensity, and newFA = origFA*[1-(weighted effectsize)]
%
% returns:
% XxYxZ image with altered fa, in image space
%
% 2009.11.15 JT

fa = dtiComputeFA(dt6);
brainMask=readFileNifti(mask);
brainMask.data=double(brainMask.data);
% effectsize = .1;


weightedESmap = effectsize*densitymap;
alteredFA = fa.*(1-weightedESmap);
alteredFAmasked = alteredFA.*brainMask.data;

% fa = dtiComputeFA(dt6);
% weightedESmap = effectsize*densitymap;
% alteredFA = fa.*(1-weightedESmap);

return

Retrieved from "http://vista.su.domains/psych221wiki/index.php?title=Matlab_Code_for_Reducing_FA&oldid=1553"