I was able to get spm_mask to work with multiple input images and thresholds by changing line 81 from
msk=msk + (img>=thresh & finite(img));
to
msk=msk + (img>=thresh(i) & finite(img));
However, there may be something wrong in this too. The result is a mask that is much more conservative than what I thought would be the equivalent expression calculated with img_calc
In image calc I did
output = (i1>100 & i2>0.15 & i3>0.3) .* img_to_be_masked
Any ideas why this may be?
Stefano Marenco, MD
GCAP, Clinical Brain Disorders Branch,
NIMH
10 Center Drive, room 4S235
Bethesda, MD 20892
tel. (301) 435-8964
fax. (301) 480-7795