You can also use mri_concat, like
mri_concat \
010130_vc6126/surf/lh.thickness.fsaverage.mgh \
zmaps/OCplusYC/lh_mean.mgh \
--paired-diff \
--o 010130_vc6126/surf/lh.thick_min_OCADmeanthick.mgh
Masking can be done by mri_mask, and thresholding with mri_binarize.
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
