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