Beauchamp:IfCortModExists

From OpenWetWare

(Difference between revisions)
Jump to: navigation, search
(New page: == What if I have scanned a subject for whom a surface model already exists? == == Automation == All of the steps below are contained in the scripts file /Volumes/data9/surface/scripts/@...)
(Automation)
Line 11: Line 11:
subjname is the subject's name in the format lastname_firstname.
subjname is the subject's name in the format lastname_firstname.
original_surface_ec is the two letter experiment code which the surface was originally constructed from
original_surface_ec is the two letter experiment code which the surface was originally constructed from
-
(this must be determined by looking in the surfaces directory)
+
(this must be determined by looking in the surfaces directory).
 +
Note that an anatavg must exist in the afni directory. If it doesn't create with
 +
  3dcopy 3DSAG_T1.nii ${ec}anatavg
== Aligning to Experiment ==
== Aligning to Experiment ==

Revision as of 09:46, 16 March 2009

What if I have scanned a subject for whom a surface model already exists?

Automation

All of the steps below are contained in the scripts file /Volumes/data9/surface/scripts/@alreadymade

To run this command type,

 /Volumes/data9/surfaces/scripts/@alreadymade current_ec subjname original_surface_ec

Where current_ec is the current two-letter experiment code subjname is the subject's name in the format lastname_firstname. original_surface_ec is the two letter experiment code which the surface was originally constructed from (this must be determined by looking in the surfaces directory). Note that an anatavg must exist in the afni directory. If it doesn't create with

 3dcopy 3DSAG_T1.nii ${ec}anatavg

Aligning to Experiment

The scripts assumes that all of the steps performed the script file @prep_dir have been performed (see Preparation for Creating Cortical Surface Models page) e.g.

 /Volumes/data9/surfaces/scripts/@recon BD doe_jane

First, check to make sure that the already created surface has the necessary components created by the @finish script

 FSread_annot -input ../label/lh.aparc.a2005s.annot -roi_1D ./lh.aparc.a2005s.1D.roi
 FSread_annot -input ../label/rh.aparc.annot -roi_1D ./rh.aparc.1D.roi
 FSread_annot -input ../label/lh.aparc.annot -roi_1D ./lh.aparc.1D.roi
 mris_convert -c ../surf/lh.sulc ../surf/lh.smoothwm ./lh.sulc.asc
 mris_convert -c ../surf/rh.sulc ../surf/rh.smoothwm ./rh.sulc.asc
 set ec = TA203
 SurfSmooth -spec {$ec}_both.spec -surf_A lh.smoothwm.asc -met NN_geom -surf_out lh.smoothwm.SS500.ply -Niter 500 -match_area 0.01
 SurfSmooth -spec {$ec}_both.spec -surf_A rh.smoothwm.asc -met NN_geom -surf_out rh.smoothwm.SS500.ply -Niter 500 -match_area 0.01
 cp /Volumes/data9/surfaces/scripts/both.spec ./both.spec

Next, go to the surfaces directory, create a copy of the new anatomy, and align it the new data with the old surface model.

 3dcopy  /Volumes/data1/UT/CZ/afni/CZanatr1+orig ./CZanatr1
 @SUMA_AlignToExperiment -exp_anat CZanatr1+orig -surf_anat TA203_SurfVol+orig -prefix CZ_SurfVol_Alnd_Exp 
 mv CZ_SurfVol_Alnd_Exp+orig.* /Volumes/data1/UT/CZ/afni/

Next, create an @EC file that points to the appropriate directory:

 cat >> @CZ
 /Volumes/data9/surfaces/scripts/@ec  CZ /Volumes/data9/surfaces/TandonLang/TA203/SUMA
Personal tools