OK, the roi analysis assumes that the annot exists for each subject. I see you have provided a full path instead of just an annotation name. Your annotation is in fsaverage space, you need to convert it to each subject's space using mri_surf2surf (you'll need the --sval-annot flag).
doug
On 07/12/2013 09:45 PM, Joseph Dien wrote:
I followed the tutorial:
mkanalysis-sess -fsd bold -analysis CPA.sm05.lh -surface fsaverage lh -fwhm 5 -event-related -paradigm CPA1.par -nconditions 20 -spmhrf 0 -TR 2 -refeventdur 20.265 -polyfit 2 -per-run -force -nuisreg nuisreg2.dat 7 -tpexclude tpexclude.dat
On Jul 12, 2013, at 7:34 PM, Douglas N Greve <greve@nmr.mgh.harvard.edu mailto:greve@nmr.mgh.harvard.edu> wrote:
Hi Joe,
On 07/12/2013 06:51 PM, Joseph Dien wrote:
I'm using Freesurfer 5.3.0
(well, actually I ran the original analyses with 5.2.0 and then am running this ROI follow-up with 5.2.0, could that be a problem? From the release notes, it didn't look like it would be an issue)
No, that is fine.
Anyway, I've run the following set of commands:
funcroi-config -roi /Volumes/Data/CP1/fsfast/secondLevel/CPA.sm05.lh/AR-v-CS/AR-v-CS_HC.wls/HC/roi.AR-v-CS_HC.HC.th2.neg.cluster1.roicfg
-annot /Volumes/Data/CP1/fsfast/secondLevel/CPA.sm05.lh/AR-v-CS/AR-v-CS_HC.wls/HC/cache.th20.neg.sig.ocn.annot
cluster-001 -analysis CPA.sm05.lh
funcroi-sess -roi /Volumes/Data/CP1/fsfast/secondLevel/CPA.sm05.lh/AR-v-CS/AR-v-CS_HC.wls/HC/roi.AR-v-CS_HC.HC.th2.neg.cluster1.roicfg
-s CPA003 -d .
I get the following output.
When I checked the results with:
tksurfer fsaverage lh inflated -overlay roi.AR-v-CS_HC.HC.th2.neg.cluster1.nii.gz -fthresh .5
sure enough, it didn't look right. Instead of the cluster, I got speckles all over the surface.
I used the Matlab command:
[vertex_coords, faces] = read_surf('lh.white');
and found sure enough that the size of the vertex array was only 143827 whereas mri_info indicates that the ocn file has 163842
When you made the analysis (mkanalysis-sess), do you specify the subject as "self" or "fsaverage"? doug
Any suggestions what I might be doing wrong?
Thanks!
Joe
/Volumes/Data/CP1/fsfast/CPA003 Fri Jul 12 18:23:28 EDT 2013 mri_annotation2label --subject CPA003 --hemi lh --annotation /Volumes/Data/CP1/fsfast/secondLevel/CPA.sm05.lh/AR-v-CS/AR-v-CS_HC.wls/HC/cache.th20.neg.sig.ocn.annot
--outdir /Volumes/Data/CP1/fsfast/CPA003/bold/CPA.sm05.lh/tmp.funcroi-sess.roi.AR-v-CS_HC.HC.th2.neg.cluster1/label subject = CPA003 annotation = /Volumes/Data/CP1/fsfast/secondLevel/CPA.sm05.lh/AR-v-CS/AR-v-CS_HC.wls/HC/cache.th20.neg.sig.ocn.annot hemi = lh outdir = /Volumes/Data/CP1/fsfast/CPA003/bold/CPA.sm05.lh/tmp.funcroi-sess.roi.AR-v-CS_HC.HC.th2.neg.cluster1/label surface = white
Reading surface /Volumes/Data/CP1/freesurfer/CPA003/surf/lh.white Loading annotations from /Volumes/Data/CP1/fsfast/secondLevel/CPA.sm05.lh/AR-v-CS/AR-v-CS_HC.wls/HC/cache.th20.neg.sig.ocn.annot MRISreadAnnotationIntoArray: vertex index out of range: 143827 i=00000000, in_array_size=143827 annot file: /Volumes/Data/CP1/fsfast/secondLevel/CPA.sm05.lh/AR-v-CS/AR-v-CS_HC.wls/HC/cache.th20.neg.sig.ocn.annot MRISreadAnnotationIntoArray: vertex index out of range: 143828 i=00000000, in_array_size=143827 annot file: /Volumes/Data/CP1/fsfast/secondLevel/CPA.sm05.lh/AR-v-CS/AR-v-CS_HC.wls/HC/cache.th20.neg.sig.ocn.annot MRISreadAnnotationIntoArray: vertex index out of range: 143829 i=00000000, in_array_size=143827 annot file: /Volumes/Data/CP1/fsfast/secondLevel/CPA.sm05.lh/AR-v-CS/AR-v-CS_HC.wls/HC/cache.th20.neg.sig.ocn.annot MRISreadAnnotationIntoArray: vertex index out of range: 143830 i=00000000, in_array_size=143827 annot file: /Volumes/Data/CP1/fsfast/secondLevel/CPA.sm05.lh/AR-v-CS/AR-v-CS_HC.wls/HC/cache.th20.neg.sig.ocn.annot MRISreadAnnotationIntoArray: vertex index out of range: 143831 i=00000000, in_array_size=143827 annot file: /Volumes/Data/CP1/fsfast/secondLevel/CPA.sm05.lh/AR-v-CS/AR-v-CS_HC.wls/HC/cache.th20.neg.sig.ocn.annot MRISreadAnnotationIntoArray: vertex index out of range: 143832 i=00000000, in_array_size=143827 annot file: /Volumes/Data/CP1/fsfast/secondLevel/CPA.sm05.lh/AR-v-CS/AR-v-CS_HC.wls/HC/cache.th20.neg.sig.ocn.annot MRISreadAnnotationIntoArray: vertex index out of range: 143833 i=00000000, in_array_size=143827
…
MRISreadAnnotationIntoArray: vertex index out of range: 163838 i=00000000, in_array_size=143827 annot file: /Volumes/Data/CP1/fsfast/secondLevel/CPA.sm05.lh/AR-v-CS/AR-v-CS_HC.wls/HC/cache.th20.neg.sig.ocn.annot MRISreadAnnotationIntoArray: vertex index out of range: 163839 i=00000000, in_array_size=143827 annot file: /Volumes/Data/CP1/fsfast/secondLevel/CPA.sm05.lh/AR-v-CS/AR-v-CS_HC.wls/HC/cache.th20.neg.sig.ocn.annot MRISreadAnnotationIntoArray: vertex index out of range: 163840 i=00000000, in_array_size=143827 annot file: /Volumes/Data/CP1/fsfast/secondLevel/CPA.sm05.lh/AR-v-CS/AR-v-CS_HC.wls/HC/cache.th20.neg.sig.ocn.annot MRISreadAnnotationIntoArray: vertex index out of range: 163841 i=00000000, in_array_size=143827 annot file: /Volumes/Data/CP1/fsfast/secondLevel/CPA.sm05.lh/AR-v-CS/AR-v-CS_HC.wls/HC/cache.th20.neg.sig.ocn.annot reading colortable from annotation file... colortable with 3 entries read (originally none) Seg base 0 max index = 2 1 952 /Volumes/Data/CP1/fsfast/CPA003/bold/CPA.sm05.lh/tmp.funcroi-sess.roi.AR-v-CS_HC.HC.th2.neg.cluster1/label/lh.cluster-001.label 2 532 /Volumes/Data/CP1/fsfast/CPA003/bold/CPA.sm05.lh/tmp.funcroi-sess.roi.AR-v-CS_HC.HC.th2.neg.cluster1/label/lh.cluster2.label mri_label2label --hemi lh --srclabel /Volumes/Data/CP1/fsfast/CPA003/bold/CPA.sm05.lh/tmp.funcroi-sess.roi.AR-v-CS_HC.HC.th2.neg.cluster1/label/lh.cluster-001.label
--srcsubject CPA003 --trglabel /Volumes/Data/CP1/fsfast/CPA003/bold/CPA.sm05.lh/tmp.funcroi-sess.roi.AR-v-CS_HC.HC.th2.neg.cluster1/target.label
--trgsubject fsaverage --regmethod surface
srclabel = /Volumes/Data/CP1/fsfast/CPA003/bold/CPA.sm05.lh/tmp.funcroi-sess.roi.AR-v-CS_HC.HC.th2.neg.cluster1/label/lh.cluster-001.label srcsubject = CPA003 trgsubject = fsaverage trglabel = /Volumes/Data/CP1/fsfast/CPA003/bold/CPA.sm05.lh/tmp.funcroi-sess.roi.AR-v-CS_HC.HC.th2.neg.cluster1/target.label regmethod = surface
srchemi = lh trghemi = lh trgsurface = white srcsurfreg = sphere.reg trgsurfreg = sphere.reg usehash = 1 Use ProjAbs = 0, 0 Use ProjFrac = 0, 0 DoPaint 0
SUBJECTS_DIR /Volumes/Data/CP1/freesurfer FREESURFER_HOME /Applications/freesurfer Loading source label. Found 952 points in source label. Starting surface-based mapping Reading source registration /Volumes/Data/CP1/freesurfer/CPA003/surf/lh.sphere.reg Rescaling ... original radius = 100 Reading target surface /Volumes/Data/CP1/freesurfer/fsaverage/surf/lh.white Reading target registration /Volumes/Data/CP1/freesurfer/fsaverage/surf/lh.sphere.reg Rescaling ... original radius = 100 Building target registration hash (res=16). Building source registration hash (res=16). INFO: found 952 nlabel points Performing mapping from target back to the source label 163842 Number of reverse mapping hits = 415 Checking for and removing duplicates Writing label file /Volumes/Data/CP1/fsfast/CPA003/bold/CPA.sm05.lh/tmp.funcroi-sess.roi.AR-v-CS_HC.HC.th2.neg.cluster1/target.label
1367 mri_label2label: Done
mri_label2label --hemi lh --srclabel /Volumes/Data/CP1/fsfast/CPA003/bold/CPA.sm05.lh/tmp.funcroi-sess.roi.AR-v-CS_HC.HC.th2.neg.cluster1/target.label
--outmask /Volumes/Data/CP1/fsfast/CPA003/bold/CPA.sm05.lh/roi.AR-v-CS_HC.HC.th2.neg.cluster1.nii.gz
--regmethod surface --trglabel /Volumes/Data/CP1/fsfast/CPA003/bold/CPA.sm05.lh/tmp.funcroi-sess.roi.AR-v-CS_HC.HC.th2.neg.cluster1/junk.label
--s fsaverage
srclabel = /Volumes/Data/CP1/fsfast/CPA003/bold/CPA.sm05.lh/tmp.funcroi-sess.roi.AR-v-CS_HC.HC.th2.neg.cluster1/target.label srcsubject = fsaverage trgsubject = fsaverage trglabel = /Volumes/Data/CP1/fsfast/CPA003/bold/CPA.sm05.lh/tmp.funcroi-sess.roi.AR-v-CS_HC.HC.th2.neg.cluster1/junk.label regmethod = surface
srchemi = lh trghemi = lh trgsurface = white srcsurfreg = sphere.reg trgsurfreg = sphere.reg usehash = 1 Use ProjAbs = 0, 0 Use ProjFrac = 0, 0 DoPaint 0
SUBJECTS_DIR /Volumes/Data/CP1/freesurfer FREESURFER_HOME /Applications/freesurfer Loading source label. Found 1311 points in source label. Starting surface-based mapping Reading source registration /Volumes/Data/CP1/freesurfer/fsaverage/surf/lh.sphere.reg Rescaling ... original radius = 100 Reading target surface /Volumes/Data/CP1/freesurfer/fsaverage/surf/lh.white Reading target registration /Volumes/Data/CP1/freesurfer/fsaverage/surf/lh.sphere.reg Rescaling ... original radius = 100 Building target registration hash (res=16). Building source registration hash (res=16). INFO: found 1311 nlabel points Performing mapping from target back to the source label 163842 Number of reverse mapping hits = 0 Checking for and removing duplicates Creating output /Volumes/Data/CP1/fsfast/CPA003/bold/CPA.sm05.lh/roi.AR-v-CS_HC.HC.th2.neg.cluster1.nii.gz Writing label file /Volumes/Data/CP1/fsfast/CPA003/bold/CPA.sm05.lh/tmp.funcroi-sess.roi.AR-v-CS_HC.HC.th2.neg.cluster1/junk.label
1311 mri_label2label: Done
Ended at Fri Jul 12 18:23:48 EDT 2013 funcroi-sess completed
Joseph Dien, Senior Research Scientist University of Maryland
E-mail: jdien07@mac.com mailto:jdien07@mac.com mailto:jdien07@mac.com Phone: 202-297-8117 http://joedien.com//
Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
-- Douglas N. Greve, Ph.D. MGH-NMR Center greve@nmr.mgh.harvard.edu mailto:greve@nmr.mgh.harvard.edu Phone Number: 617-724-2358 Fax: 617-726-7422
Bugs: surfer.nmr.mgh.harvard.edu/fswiki/BugReporting FileDrop: https://gate.nmr.mgh.harvard.edu/filedrop2 www.nmr.mgh.harvard.edu/facility/filedrop/index.html Outgoing: ftp://surfer.nmr.mgh.harvard.edu/transfer/outgoing/flat/greve/
Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
The information in this e-mail is intended only for the person to whom it is addressed. If you believe this e-mail was sent to you in error and the e-mail contains patient information, please contact the Partners Compliance HelpLine at http://www.partners.org/complianceline . If the e-mail was sent to you in error but does not contain patient information, please contact the sender and properly dispose of the e-mail.
Joseph Dien, Senior Research Scientist University of Maryland
E-mail: jdien07@mac.com mailto:jdien07@mac.com Phone: 202-297-8117 http://joedien.com//