I just modified fscalc.fsl to take surfaces. This is a front-end to
FSL's fslmaths command. You give it the exact same arguments, but the
files can be of any format. If the underlying data are on a surface,
then add --surf subject hemi, and it will reshape the surface to fit
in nifti (as long as the largest prime factor of the number of
verticies is less than 2^15, but that was probably obvious).
doug
On Fri, 31 Oct 2008, Akram Bakkour wrote:
Hello,
I have two thickness files I want to subtract one from the other, so I did:
on thigpen
source /usr/local/freesurfer/nmr-std-env
cd /space/elvin_005/users/akram/AGING_MASSIVE
setenv SUBJECTS_DIR $PWD
mris_surface_stats -mask fsaverage/label/lh.cort_thick.label \
-surf_name $SUBJECTS_DIR/fsaverage/surf/lh.white \
-src_type w \
-out_name $SUBJECTS_DIR/zmaps/OCplusYC/lh_std.mgh \
-mean $SUBJECTS_DIR/zmaps/OCplusYC/lh_mean.mgh \
OAS1_????_MR1/surf/lh.thickness.fsaverage.mgh
#OAS1_????_MR1/surf/lh.thickness.fsaverage.mgh were generated using -qcache
mris_calc mris_calc -o 010130_vc6126/surf/lh.thick_min_OCADmeanthick.mgh \
010130_vc6126/surf/lh.thickness.fsaverage.mgh sub zmaps/OCplusYC/lh_mean.mgh
\
and I get this weird error:
"mris_calc: curvature file '010130_vc6126/surf/lh.thickness.fsaverage.mgh'
has wrong magic number."
Any ideas what this means, or if there is any way around this?
I suppose I could do this in MATLAB, but I was curious about this magic
number business.
Thanks!
Happy Halloween!
A
--
Douglas N. Greve, Ph.D.
MGH-NMR Center
[EMAIL PROTECTED]
Phone Number: 617-724-2358
Fax: 617-726-7422
In order to help us help you, please follow the steps in:
surfer.nmr.mgh.harvard.edu/fswiki/BugReporting
_______________________________________________
Freesurfer mailing list
Freesurfer@nmr.mgh.harvard.edu
https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
