Beauchamp:MapIco: Difference between revisions

From OpenWetWare
Jump to navigationJump to search
No edit summary
No edit summary
Line 3: Line 3:
==Old and New==
==Old and New==
In March 2017, the AFNI team discovered a bug in the MapIcosahedron program. Here is a comparison of group average results with the old and new (fixed) versions.
In March 2017, the AFNI team discovered a bug in the MapIcosahedron program. Here is a comparison of group average results with the old and new (fixed) versions.
If you have old surface averages and would like to fix them, you can use a script like this:
  set ss = (MS MT MU MW MX MY MZ NA NB NC ND NE NF NG NH NU NV NW OA )
  foreach s ( $ss )
  set sd = /Volumes/data/BCM/{$s}/fs/SUMA
  cd $sd 
  pwd
  mkdir oldstd
  mv *141* *60* oldstd/
  set sid = fs
  foreach size (60 141)
    foreach hem ( lh rh )
      MapIcosahedron -spec ${sid}_$hem.spec -ld $size  \
            -dset_map $hem.thickness.gii.dset            \
            -dset_map $hem.curv.gii.dset                \
            -dset_map $hem.sulc.gii.dset                \
            -prefix std.$size.
    end
    inspec -LRmerge std.$size.${sid}_lh.spec std.$size.${sid}_rh.spec \
          -prefix std.$size.${sid}_both.spec
end
end
==Anatomical Data==
==Anatomical Data==
New
New

Revision as of 08:20, 30 March 2017

Return to Notes on Cortical Surface Averages

Old and New

In March 2017, the AFNI team discovered a bug in the MapIcosahedron program. Here is a comparison of group average results with the old and new (fixed) versions. If you have old surface averages and would like to fix them, you can use a script like this:

 set ss = (MS MT MU MW MX MY MZ NA NB NC ND NE NF NG NH NU NV NW OA )
 foreach s ( $ss )
 set sd = /Volumes/data/BCM/{$s}/fs/SUMA
 cd $sd  
 pwd 
 mkdir oldstd
 mv *141* *60* oldstd/
 set sid = fs
 foreach size (60 141)
   foreach hem ( lh rh )
      MapIcosahedron -spec ${sid}_$hem.spec -ld $size   \
           -dset_map $hem.thickness.gii.dset            \
           -dset_map $hem.curv.gii.dset                 \
           -dset_map $hem.sulc.gii.dset                 \
           -prefix std.$size.
   end
   inspec -LRmerge std.$size.${sid}_lh.spec std.$size.${sid}_rh.spec \
          -prefix std.$size.${sid}_both.spec 
end
end

Anatomical Data

New

     

Old

      

Functional Data

Functional data of the Mouth vs. Eye contrast mapped to the average surface. It is hard to quantify the difference, but one way is to compare the maximum statistical values; with better intersubject alignment, these would presumably increase.

  1. OLD: Max t-value for Mouth: 7.9 Eye:8.8 Eye vs Mouth: 5.8
  2. NEW: Max t-value for Mouth: 9.3 Eye:10.1 Eye vs Mouth: 5.6


New

 

Old