Hi Rudolph,
we store things like thickness in 'volume' formats these days - the volumes
are nvertices x 1 x 1. Not sure why it didn't work, I would have thought it
would. Not sure. Akram: can you run mri_info on the input files and send
along the results?
Bruce
On Sat, 1 Nov 2008
[EMAIL PROTECTED] wrote:
mris_calc currently only works with flat curvature files (eg ?h.curv, ?h.sulc,
?h.area, ?h.thickness, etc) and not on volumes.
------Original Message------
From: Akram Bakkour
Sender: [EMAIL PROTECTED]
To: freesurfer@nmr.mgh.harvard.edu
Subject: [Freesurfer] mris_calc
Sent: Oct 31, 2008 20:13
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
_______________________________________________
Freesurfer mailing list
Freesurfer@nmr.mgh.harvard.edu
https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer