Freesurfer

freesurfer@nmr.mgh.harvard.edu

19 participants
30088 discussions

switching from v4.5.0-32bits to v5.0.0-64bits
by Daniel Ferreira Padilla 02 May '11

02 May '11

Dear experts, I have been processing my data on FreeSurfer v4.5.0, 32 bits (SUSE). Now we have changed our system to 64bits (UBUNTU) and we are planning to install FreeSurfer v5.0.0, 64 bits. Do I have to process my data again or can I carry on processing my new data on v5.0.0, 64 bits, and then mix all data (v4.5.0 and v5.0.0) for analysis in v5.0.0? Thank you very much in advance Daniel Ferreira

2 1

nii to label on a surface
by Michael Thomas Smith 02 May '11

02 May '11

Hello, I'm trying to get this ROI drawn on a brain surface. http://www.sal.mvm.ed.ac.uk/images/MT_roi.png To do this, I first ran a script (of the form)[1]: mri_cor2label --i ~/roi.nii --l ~/roi.label --id 1 Then opened the new roi.label file in tksurfer: :~$ tksurfer s01 lh orig The ROI seemed very small and in the wrong place: http://www.sal.mvm.ed.ac.uk/images/MT_scatter2.png http://www.sal.mvm.ed.ac.uk/images/MT_scatter3.png (the purple labelled areas are sulci contours) I'm guessing some form of registration went wrong. I tried to first convert it to .nii/.tgh but using the -rl option of mri_convert: :~$ mri_convert -rl ~/fs_subjects/s01/mri/T1.mgz -rt nearest ~/mainpref_mask_1_3.nii ~/test_mask.mgh This still looked ok, when I checked the ROI. I tried mri_cor2label again, but found the region was still in the wrong place. I figured I needed to give it more info about the surface that it's to be drawn on, so I tried: :~$ mri_cor2label --i ~/test_mask.mgh --l ~/test_mask.label --id 1 --surf s01 lh $Id: mri_cor2label.c,v 1.8.2.2 2009/02/20 16:37:50 greve Exp $ Loading mri /home/s0783410/test_mask.mgh Loading /home/s0783410/fs_subjects/s01/surf/lh.white ERROR: dim mismatch between surface (127260) and input (16777216) but that didn't seem to work. To see what might be going on, I tried out the whole thing on V1, which I got segmented by the v1_segment function I asked about in my last email. The output seemed reasonable: http://www.sal.mvm.ed.ac.uk/images/V1before.tiff and I was able to get it to a nice ROI in a nii file (thanks doug!): http://www.sal.mvm.ed.ac.uk/images/V1_nii.png and I got the nii file back to a surface using: :~$ mri_cor2label --i ~/imageA.nii --l ~/imageA.label --id 1 Then viewing it in: :~$ tksurfer s01 lh orig (and loading the imageA.label file). Which is a little odd: http://www.sal.mvm.ed.ac.uk/images/V1mess.tiff (although it seems to be lined up properly, it has a little extra region selected in the Lingual(?) Gyrus). I thought the difference might be because the V1 ROI came from the original T1 image used by freesurfer, while the ROI I'm trying to use is from a functional image. However, to check, I coregistered (resliced) the function ROI to the original anatomical image, but this didn't seem to make any different. Any suggestions or ideas about where I should be looking? This is I guess one of the most common tasks people use freesurfer for, so I guess I just need to look in the right place for help. Thanks for the help! Mike. [1] Actual output: s0783410@spatial2:~$ mri_cor2label --i ~/mainpref_mask_1_3.nii --l ~/mainpref_mask_1_3.label --id 1 $Id: mri_cor2label.c,v 1.8.2.2 2009/02/20 16:37:50 greve Exp $ Loading mri /home/s0783410/mainpref_mask_1_3.nii ------- Vox2RAS of input volume ----------- -1.719 0.000 0.000 110.003; 0.000 0.000 2.400 -30.001; 0.000 -1.719 0.000 110.003; 0.000 0.000 0.000 1.000; Scanning the volume Found 473 label voxels Writing label file /home/s0783410/mainpref_mask_1_3.label LabelWrite: saving to /home/s0783410/mainpref_mask_1_3.label Centroid: -0.91 -6.49 54.15 mri_cor2label completed SUCCESSFULLY -- The University of Edinburgh is a charitable body, registered in Scotland, with registration number SC005336.

1 0

Minimum n for glmfit maps?
by Adam Nitenson 02 May '11

02 May '11

Hi Freesurfers, Is there a minimum number of subjects that should be used in isxconcat groups and subsequent glmfit analysis? I recently generated activation maps from a group with only 2 subjects, but the resulting image (too large to attach) is very strage. The activation patterns are "whispy", almost looking like someone just doodled red and blue over parts of the brain, sometimes forming hollow circles. Is this just a consequence of only using 2 subjects or is there another problem I'm missing. I find it strange that both of the subjects' individual maps have rather robust (negative) activation)...much more than the group map would suggest. Another group of three subjects is also a bit strange (not nearly as bad as 2) and it doesn't seem until I reach a group of 4 or 5 that the patterns appear "normal". For the record here are the commands I used: isxconcat-sess -sf /cluster/roffman/users/Stable5_PerRun/Subject_Files/Apr2011/2_FB_DEP_TT_1T3T5T7T -d /cluster/roffman/users/Stable5_PerRun -analysis SIRP_Stable5 -c Cond1v3 -c Cond1vFix -c Cond3vFix -c Cond1v5 -c Cond1v7 -c Cond5vFix -c Cond7vFix -hemis -o /cluster/roffman/users/Stable5_PerRun/April_2011/2_FB_DEP_TT_1T3T5T7T mri_glmfit --y /cluster/roffman/users/Stable5_PerRun/April_2011/2_FB_DEP_TT_1T3T5T7T/SIRP_Stable5/Cond1v7/lh.ces.nii --fsgd /cluster/roffman/users/fsgd/April_2011/2_FB_DEP_TT_1T3T5T7T.fsgd --osgm --glmdir /cluster/roffman/users/Stable5_PerRun/April_2011/2_FB_DEP_TT_1T3T5T7T/SIRP_Stable5/osgm/LH1v7 --surf fsaverage lh --fwhm 4.6 Thanks! Adam Nitenson, B.S. Brain Genomics Lab Massachusetts General Hospital nitenson(a)nmr.mgh.harvard.edu Phone: 617-643-3215

2 1

inflated surface
by Tetiana Dadakova 30 Apr '11

30 Apr '11

Dear FreeSurfer list, When I look at inflated surface of my subjects (healthy males and females), almost all of them have a concavity in the posterior part of the brain (the arrow in attached image). It doesn't look correct to me. What can cause this problem? Thank you, Tanja.

3 4

distortion of flat map
by Yang Liu 30 Apr '11

30 Apr '11

Hi surfers, I am working on generating surface ROIs from the retinotopy surface. We cut the occipital lobe and draw paths on the flattened map to create labels. However, I found that there are distortions on the flat map that cause large errors when drawing these paths. I attached two pictures. The path_example2_dots.jpg are the dots that I selected to draw a path. Basically it is the method we use to demarcate the boundaries of visual areas. And the path_example2_line.jpg are the line that is supposed to connect the dots and to be the correct boundary. Notice the dots and the line do not overlap! so when I draw a line on the flat map, Freesurfer actually doesn't draw the line I want it to be. Has anyone met this problem before? Is there a better solution to create surface ROIs? Drawing the boundary on inflated surface? Yang

3 2

Compiling stats file for ROI label in QDec
by Krista Kelly 29 Apr '11

29 Apr '11

Using QDec, I've created an ROI label (lh.G_cuneus.label) and a .stats file (lh.G_cuneus.stats) that corresponds to this ROI for each participant in my group. Is there a quick way to compile cortical thickness data for this ROI from all participants? The manual suggests using aparcstats2table but I'm not sure how to select the stats for just this one ROI. Any help would be appreciated!

1 0

Run mri_ms_fitparms without motion correction
by Michael Waskom 29 Apr '11

29 Apr '11

Hi Bruce, et all, Is it possible to run mri_ms_fitparms without motion correction, if your input images are already in register? It seems like it would save quite a bit of processing time, but the usage help isn't all that useful (looks like there are a lot of hidden options buried in the source code). Thanks! Michael

2 1

tkmedit error
by Pom Sailasuta 29 Apr '11

29 Apr '11

Hi, I got an error with tkmedit command "freeglut (/home/pom/freesurfer/freesurfer/bin/tkmedit.bin): Unable to create direct context rendering for window 'Medit: ' This may hurt performance. Xlib: extension "XFree86-DRI" missing on display ":0.0". The only DRI related file I have is called /usr/X11R6/lib/modules/dri/fglrx_dri.so My freesurfer version is "freesurfer-Linux-rh9-stable-pub-v3.0.5". I am wondering if I need to update my linux, currently it is a very old OpenSuse 10.0 version? thank you -- Pom

1 0

HELP: recon-all error
by Yunjo Lee 28 Apr '11

28 Apr '11

Hi all, I sent an email a couple of days ago but didn't get any answer yet... I recently installed FreeSurfer in a virtual machine (Ubuntu). I installed it according to the freesurfer website and another source. http://mirko.windhoff.net/how_to/install_freesurfer_on_ubuntu When I ran recon-all, I kept getting the following error message (see below). Is there anything else that I have to install?? Thanks in advance, Yunjo *Here's the error message:* 0597: dt=0.056000, rms=0.417 (0.032%), neg=0, invalid=766 0598: dt=0.028000, rms=0.417 (-0.000%), neg=0, invalid=766 gcamFindOptimalTimeStep: Complete in 22111.002 ms 0599: dt=0.000000, rms=0.417 (0.000%), neg=0, invalid=766 GCAMwrite(transforms/talairach.m3z): gzip encountered error. writing output transformation to transforms/talairach.m3z... mri_ca_register: GCAMwrite(transforms/talairach.m3z) failed ERROR: mri_ca_register with non-zero status 252 but continuing despite the error #-------------------------------------- #@# CA Reg Inv Tue Apr 26 16:21:09 EDT 2011 /usr/local/freesurfer/subjects/YC7178/mri mri_ca_register -invert-and-save transforms/talairach.m3z Loading, Inverting, Saving, Exiting ... Reading transforms/talairach.m3z ERROR: cannot find or read transforms/talairach.m3z ERROR: mri_ca_register with non-zero status 0 Linux VMBronzeWing 2.6.35-22-generic #33-Ubuntu SMP Sun Sep 19 20:34:50 UTC 2010 i686 GNU/Linux recon-all -s YC7178 exited with ERRORS at Tue Apr 26 16:21:14 EDT 2011

2 1

Re: [Freesurfer] Freesurfer Digest, Vol 86, Issue 77
by Maryam Vaziri Pashkam 28 Apr '11

28 Apr '11

FOUND IT! The problem was with the startup.m file (nice catch dough). I had a change directory line (cd users/maryam/programs) in the end of the startup.m file that made matlab change directory to another directory that did not have the textfile necessary for the mcparams2extreg. I took that line out and now the preproc-sess works. Thanks for all the help. Maryam On Thu, Apr 28, 2011 at 5:32 PM, Douglas Greve <greve(a)nmr.mgh.harvard.edu>wrote: > That first one should have worked (or at least given you some other > error). That error indicates that it can't find mcp2x.m. But it's in the > same directory that you are running matlab from, right? I thought of > something else to check. Can you send me your ~/matlab/startup.m file? > > doug > > > On 4/28/11 3:51 PM, Maryam Vaziri Pashkam wrote: > > Hi, > > I tried " matlab -r mcp2x" and " matlab -r mcp2x.m" in both cased I get > an error that "undefined function or variable mcp2x. I tried cat mcp2x.m | > matlab, and the program runs and gives the exact same error as before when I > ran mcparams2extreg from within terminal (see bellow). It looks like an > error that occurs when you are in the wrong directory in matlab and you are > calling a file from another directory which does not make sense as the fmc > files are created in the correct path. > > Maryam > > >> >> >> >> >> >> >> >> ??? Error using ==> textread at 167 > File not found. > > >> ??? Undefined variable mc. > > >> ??? Undefined function or variable 'mc'. > > >> ??? Undefined function or variable 'ntrs'. > > >> >> ??? Undefined variable mc. > > >> >> ??? Undefined function or variable 'ntrs'. > > >> ??? Undefined function or variable 'x'. > > >> >> >> ??? Undefined function or variable 'x2'. > > >> >> >> >> >> >> >> >> >> >> ??? Reference to non-existent field 'vol'. > > >> >> ??? Reference to non-existent field 'vol'. > > Error in ==> MRIwrite at 70 > bmri.voldim = [size(mri.vol,1) size(mri.vol,2) size(mri.vol,3)]; > > > > > On Thu, Apr 28, 2011 at 8:52 AM, <greve(a)nmr.mgh.harvard.edu> wrote: > > I'm not sure where to proceed next. I can't create the error on my mac, so > it must be something about your version of matlab/macos. Let's go back to > trying the matlab -r command. Is this what you did: > > cd session/bold > mcparams2extreg -extreg junk.txt -mcfile RRR/fmc.mcdat -nkeep 6 -monly > mcp2x.m > matlab -r mcp2x > > It should not have generated the same error as when using mcp2x.m. > > Can you also try > cat mcp2x.m | matlab > > doug > > > > HI Doug, > > > > Did you get the file I sent you? > > Any new idea about how I can fix the problem? > > > > Maryam > > > > On Tue, Apr 26, 2011 at 5:25 AM, Douglas Greve > > <greve(a)nmr.mgh.harvard.edu>wrote: > > > >> > >> The file will be found in $FREESURFER_HOME/fsfast/bin/mcparams2extreg, > >> but > >> until we figure out how to run matlab in a way that does not crash, it > >> will > >> not do any good to edit this file. Actually, can you send me your > >> version of > >> mcparams2extreg? I ran your matlab file on my mac by catting it to > >> matlab, > >> but it failed with I tried to run it with -r (don't know why). > >> > >> doug > >> > >> On 4/25/11 4:41 PM, Maryam Vaziri Pashkam wrote: > >> > >> Here is the m file. > >> > >> I have two versions of Matlab: R2008a, R2009b. The one that freesurfer > >> runs when calling matlab codes in 2009b. > >> > >> Sorry I am not sure what I should copy to my home directory and edit. I > >> thought the line you asked me to change was in the mcparams2extreg > >> file > >> which I don't know where it is. Could you explain a bit more? > >> > >> Maryam > >> > >> > >> > >> On Mon, Apr 25, 2011 at 4:34 PM, Douglas Greve > >> <greve(a)nmr.mgh.harvard.edu>wrote: > >> > >>> Can you send me the matlab file? And what version of matlab are you > >>> running? On the mac, you'll probably have to copy it to your home > >>> directly, > >>> edit it, then copy it back. When you copy it back, you'll have to sudo > >>> to > >>> get administrator privileges. > >>> > >>> doug > >>> > >>> > >>> On 4/25/11 4:26 PM, Maryam Vaziri Pashkam wrote: > >>> > >>> There is an error: > >>> > >>> ??? Error: Unexpected MATLAB operator when I run matlab -r > >>> matlabfile.m > >>> I wanted to change the mcparams2extreg to see if it works but I don't > >>> know > >>> how to change it and where it is located. > >>> > >>> Maryam > >>> > >>> On Mon, Apr 25, 2011 at 3:22 PM, Douglas Greve > >>> <greve(a)nmr.mgh.harvard.edu > >>> > wrote: > >>> > >>>> > >>>> Yea, I vaguely remember this problem with 4.5 and the macs. Try > >>>> running > >>>> this from a terminal (not from matlab): > >>>> > >>>> matlab -r matlabfile.m > >>>> > >>>> where matlabfile.m is the file created with -monly. If this works you > >>>> can > >>>> change the call to matlab from mcparams2extreg from > >>>> > >>>> cat $MLF | $MATLAB -display iconic | tee -a $LF > >>>> > >>>> to > >>>> > >>>> $MATLAB -r $MLF -display iconic | tee -a $LF > >>>> > >>>> doug > >>>> > >>>> > >>>> > >>>> > >>>> > >>>> On 4/25/11 2:49 PM, Maryam Vaziri Pashkam wrote: > >>>> > >>>> Hi Dough, > >>>> > >>>> The generated matlab programs works from within matlab and generates > >>>> mcextreg.bhdr and mcextreg_000.bfloat and mcextreg_000.hdr. > >>>> > >>>> It seems that there is a problem calling matlab from terminal? any > >>>> idea > >>>> how I can fix this? I have two versions on Matlab on my computer. > >>>> Could this > >>>> be a source of the problem? I checked the program in both matlabs and > >>>> it > >>>> works when run from inside Matlab. > >>>> > >>>> Maryam > >>>> > >>>> > >>>> On Mon, Apr 25, 2011 at 4:32 AM, Douglas Greve < > >>>> greve(a)nmr.mgh.harvard.edu> wrote: > >>>> > >>>>> Try running mcparams2extreg with the -monly option to create a > >>>>> matlab > >>>>> file. Then run that matlab file from inside matlab. If that fails in > >>>>> the > >>>>> same way, send me the matlab file. > >>>>> > >>>>> doug > >>>>> > >>>>> > >>>>> On 4/25/11 4:02 AM, Maryam Vaziri Pashkam wrote: > >>>>> > >>>>> Hi Dough, > >>>>> > >>>>> I have attached the log file. It does not give any extra information. > >>>>> I > >>>>> ran the command you told me and gave me the exact same errors. It > >>>>> looks like > >>>>> there is a problem with that command and reading of the fmc.mcdat > >>>>> file. I > >>>>> have also attached the fmc.mcdat file from run 6. > >>>>> Any other idea? > >>>>> > >>>>> Thanks, > >>>>> Maryam > >>>>> > >>>>> On Mon, Apr 25, 2011 at 3:44 AM, < > >>>>> freesurfer-request(a)nmr.mgh.harvard.edu> wrote: > >>>>> > >>>>>> Send Freesurfer mailing list submissions to > >>>>>> freesurfer(a)nmr.mgh.harvard.edu > >>>>>> > >>>>>> To subscribe or unsubscribe via the World Wide Web, visit > >>>>>> https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer > >>>>>> or, via email, send a message with subject or body 'help' to > >>>>>> freesurfer-request(a)nmr.mgh.harvard.edu > >>>>>> > >>>>>> You can reach the person managing the list at > >>>>>> freesurfer-owner(a)nmr.mgh.harvard.edu > >>>>>> > >>>>>> When replying, please edit your Subject line so it is more specific > >>>>>> than "Re: Contents of Freesurfer digest..." > >>>>>> > >>>>>> > >>>>>> Today's Topics: > >>>>>> > >>>>>> 1. Fwd: error with motion correction (Maryam Vaziri Pashkam) > >>>>>> 2. Re: Fwd: error with motion correction (Douglas Greve) > >>>>>> > >>>>>> > >>>>>> > ---------------------------------------------------------------------- > >>>>>> > >>>>>> Message: 1 > >>>>>> Date: Sun, 24 Apr 2011 21:01:43 -0400 > >>>>>> From: Maryam Vaziri Pashkam <mvaziri.p(a)gmail.com> > >>>>>> Subject: [Freesurfer] Fwd: error with motion correction > >>>>>> To: freesurfer(a)nmr.mgh.harvard.edu > >>>>>> Message-ID: <BANLkTikkaJvqP5bxVD5e2moWQvwBrev9iA(a)mail.gmail.com> > >>>>>> Content-Type: text/plain; charset="iso-8859-1" > >>>>>> > >>>>>> Hi, > >>>>>> > >>>>>> I am running into trouble with motion correction on my Mac and > >>>>>> freesurfer > >>>>>> 4.5. > >>>>>> Originally I ran preproc-sess and got an error from mc-sess in the > >>>>>> part > >>>>>> (i > >>>>>> think) that runs the mcextreg. I ran mc-sess and the same error > >>>>>> happened (I > >>>>>> have copy pasted the log bellow). > >>>>>> It seems like the matlab program that runs mcparams2extreg cannot > >>>>>> read > >>>>>> the > >>>>>> fmc.mcdat file. when I did textread('fmc.mcdat') in matlab, I could > >>>>>> open the > >>>>>> file and it looked perfectly fine to me with the correct number of > >>>>>> TR's. > >>>>>> Any idea where the problem might be? > >>>>>> > >>>>>> Maryam > >>>>>> > >>>>>> preproc-sess -nosmooth -i f -sf sessid -df sessdir > >>>>>> -------------------------------------------------------------- > >>>>>> preproc-sess logfile is log/preproc-sess.log > >>>>>> -------------------------------------------------------------- > >>>>>> instem f > >>>>>> mc 1 f fmc > >>>>>> stc 0 fmc > >>>>>> sm 0 > >>>>>> mask 1 brain > >>>>>> 110419_xu_PDM1_02_PDM MC ----------------------------- > >>>>>> Sun Apr 24 20:42:29 EDT 2011 > >>>>>> /Users/maryam/mri-space/studies/PDM1/PDM1 > >>>>>> mc-sess -fstem f -fmcstem fmc -s 110419_xu_PDM1_02_PDM -d > >>>>>> /Users/maryam/mri-space/studies/PDM1/unpackdata -fsd bold -new > >>>>>> -update > >>>>>> -no-inorm > >>>>>> ---------------------------------------- > >>>>>> Logfile is > >>>>>> > >>>>>> > /Users/maryam/mri-space/studies/PDM1/PDM1/log/mc-110419_xu_PDM1_02_PDM-bold.log > >>>>>> ------------------------------------------- > >>>>>> > /Users/maryam/mri-space/studies/PDM1/unpackdata/110419_xu_PDM1_02_PDM > >>>>>> RunList: 005 006 007 008 <005%20006%20007%20008> 011 012 013 014 > >>>>>> Targ Stem is 005/f > >>>>>> Update not needed > >>>>>> -rw-r--r-- 1 maryam staff 129203552 Apr 24 20:12 005/fmc.nii > >>>>>> -rw-r--r-- 1 maryam staff 64601952 Apr 24 19:12 005/f.nii > >>>>>> ------- ************************** ----------- > >>>>>> ------- Motion Correcting Run 006 ----------- > >>>>>> ------- ************************** ----------- > >>>>>> sess = 110419_xu_PDM1_02_PDM > >>>>>> Sun Apr 24 20:42:30 EDT 2011 > >>>>>> mc-afni2 --i 006/f.nii --t 005/f.nii --frame 0 --o 006/fmc.nii > >>>>>> --mcdat > >>>>>> 006/fmc.mcdat > >>>>>> > >>>>>> > /Users/maryam/mri-space/studies/PDM1/unpackdata/110419_xu_PDM1_02_PDM/bold > >>>>>> mri_info --nframes --o /tmp/mriinfo.2758 005/f.nii > >>>>>> log file is 006/fmc.nii.mclog > >>>>>> mcdat file is 006/fmc.mcdat > >>>>>> tmp dir is 006/tmp.mc-afni2.2758 > >>>>>> #@# -------------------------------- > >>>>>> > >>>>>> > /Users/maryam/mri-space/studies/PDM1/unpackdata/110419_xu_PDM1_02_PDM/bold > >>>>>> mri_convert 005/f.nii 006/tmp.mc-afni2.2758/tempvol.nii.gz --nskip 0 > >>>>>> --ndrop > >>>>>> 205 -odt float > >>>>>> mri_convert 005/f.nii 006/tmp.mc-afni2.2758/tempvol.nii.gz --nskip 0 > >>>>>> --ndrop > >>>>>> 205 -odt float > >>>>>> nskip = 0 > >>>>>> ndrop = 205 > >>>>>> $Id: mri_convert.c,v 1.146.2.5 2009/04/08 18:40:23 nicks Exp $ > >>>>>> reading from 005/f.nii... > >>>>>> TR=2000.00, TE=0.00, TI=0.00, flip angle=0.00 > >>>>>> i_ras = (-1, -0.000171665, -3.15037e-05) > >>>>>> j_ras = (0.000174532, -0.983572, -0.180519) > >>>>>> k_ras = (-2.07852e-18, -0.180519, 0.983572) > >>>>>> changing data type from 4 to 3 (noscale = 0)... > >>>>>> Dropping last 205 frames > >>>>>> writing to 006/tmp.mc-afni2.2758/tempvol.nii.gz... > >>>>>> #@# -------------------------------- > >>>>>> > >>>>>> > /Users/maryam/mri-space/studies/PDM1/unpackdata/110419_xu_PDM1_02_PDM/bold > >>>>>> mri_convert 006/f.nii 006/tmp.mc-afni2.2758/invol.nii.gz -odt float > >>>>>> mri_convert 006/f.nii 006/tmp.mc-afni2.2758/invol.nii.gz -odt float > >>>>>> $Id: mri_convert.c,v 1.146.2.5 2009/04/08 18:40:23 nicks Exp $ > >>>>>> reading from 006/f.nii... > >>>>>> TR=2000.00, TE=0.00, TI=0.00, flip angle=0.00 > >>>>>> i_ras = (-1, -0.000171665, -3.15037e-05) > >>>>>> j_ras = (0.000174532, -0.983572, -0.180519) > >>>>>> k_ras = (-2.07852e-18, -0.180519, 0.983572) > >>>>>> changing data type from 4 to 3 (noscale = 0)... > >>>>>> writing to 006/tmp.mc-afni2.2758/invol.nii.gz... > >>>>>> #@# -------------------------------- > >>>>>> > >>>>>> > /Users/maryam/mri-space/studies/PDM1/unpackdata/110419_xu_PDM1_02_PDM/bold > >>>>>> 3dvolreg -verbose -dfile 006/fmc.mcdat -base > >>>>>> 006/tmp.mc-afni2.2758/tempvol.nii.gz -prefix > >>>>>> 006/tmp.mc-afni2.2758/outvol.nii.gz > >>>>>> 006/tmp.mc-afni2.2758/invol.nii.gz > >>>>>> ++ 3dvolreg: AFNI version=AFNI_2008_02_01_1144 (Jul 3 2008) > >>>>>> [32-bit] > >>>>>> ++ Authored by: RW Cox > >>>>>> *+ WARNING: If you are performing spatial transformations on an > >>>>>> oblique > >>>>>> dset, > >>>>>> such as 006/tmp.mc-afni2.2758/tempvol.nii.gz, > >>>>>> or viewing/combining it with volumes of differing obliquity, > >>>>>> you should consider running: > >>>>>> 3dWarp -deoblique > >>>>>> on this and other oblique datasets in the same session. > >>>>>> See 3dWarp -help for details. > >>>>>> ++ Oblique dataset:006/tmp.mc-afni2.2758/tempvol.nii.gz is 10.399981 > >>>>>> degrees > >>>>>> from plumb. > >>>>>> ++ Reading in base dataset 006/tmp.mc-afni2.2758/tempvol.nii.gz > >>>>>> ++ Oblique dataset:006/tmp.mc-afni2.2758/invol.nii.gz is 10.399981 > >>>>>> degrees > >>>>>> from plumb. > >>>>>> ++ Reading input dataset 006/tmp.mc-afni2.2758/invol.nii.gz > >>>>>> ++ Edging: x=4 y=4 z=2 > >>>>>> ++ Creating mask for -maxdisp > >>>>>> + Automask has 57667 voxels > >>>>>> + 10333 voxels left in -maxdisp mask after erosion > >>>>>> ++ Initializing alignment base > >>>>>> ++ Starting final pass on 206 sub-bricks: > >>>>>> > >>>>>> > 0..1..2..3..4..5..6..7..8..9..10..11..12..13..14..15..16..17..18..19..20..21..22..23..24..25..26..27..28..29..30..31..32..33..34..35..36..37..38..39..40..41..42..43..44..45..46..47..48..49..50..51..52..53..54..55..56..57..58..59..60..61..62..63..64..65..66..67..68..69..70..71..72..73..74..75..76..77..78..79..80..81..82..83..84..85..86..87..88..89..90..91..92..93..94..95..96..97..98..99..100..101..102..103..104..105..106..107..108..109..110..111..112..113..114..115..116..117..118..119..120..121..122..123..124..125..126..127..128..129..130..131..132..133..134..135..136..137..138..139..140..141..142..143..144..145..146..147..148..149..150..151..152..153..154..155..156..157..158..159..160..161..162..163..164..165..166..167..168..169..170..171..172..173..174..175..176..177..178..179..180..181..182..183..184..185..186..187..188..189..190..191..192..193..194..195..196..197..198..199..200..201..202..203..204..205.. > >>>>>> ++ CPU time for realignment=16.6 s [=0.0808 s/sub-brick] > >>>>>> ++ Min : roll=-0.008 pitch=-0.103 yaw=-0.096 dS=+0.255 dL=-0.019 > >>>>>> dP=-0.151 > >>>>>> ++ Mean: roll=+0.031 pitch=-0.011 yaw=-0.045 dS=+0.619 dL=+0.037 > >>>>>> dP=-0.055 > >>>>>> ++ Max : roll=+0.060 pitch=+0.112 yaw=-0.010 dS=+0.693 dL=+0.075 > >>>>>> dP=+0.015 > >>>>>> ++ Max displacements (mm) for each sub-brick: > >>>>>> 0.49 0.65 0.63 0.74 0.69 0.75 0.71 0.77 0.70 0.76 0.67 0.77 0.72 > >>>>>> 0.73 > >>>>>> 0.66 > >>>>>> 0.67 0.67 0.67 0.69 0.70 0.72 0.66 0.71 0.68 0.72 0.67 0.73 0.68 > >>>>>> 0.73 > >>>>>> 0.66 > >>>>>> 0.73 0.66 0.70 0.64 0.72 0.69 0.71 0.71 0.72 0.66 0.66 0.67 0.72 > >>>>>> 0.70 > >>>>>> 0.71 > >>>>>> 0.72 0.71 0.69 0.72 0.72 0.71 0.72 0.72 0.69 0.71 0.74 0.73 0.69 > >>>>>> 0.72 > >>>>>> 0.72 > >>>>>> 0.76 0.76 0.79 0.73 0.75 0.74 0.75 0.75 0.71 0.74 0.67 0.72 0.65 > >>>>>> 0.72 > >>>>>> 0.65 > >>>>>> 0.69 0.69 0.65 0.68 0.64 0.72 0.67 0.67 0.68 0.69 0.65 0.73 0.67 > >>>>>> 0.71 > >>>>>> 0.69 > >>>>>> 0.73 0.72 0.76 0.74 0.72 0.66 0.68 0.69 0.69 0.68 0.68 0.70 0.73 > >>>>>> 0.70 > >>>>>> 0.73 > >>>>>> 0.78 0.74 0.73 0.75 0.74 0.72 0.74 0.63 0.66 0.68 0.74 0.72 0.74 > >>>>>> 0.79 > >>>>>> 0.74 > >>>>>> 0.74 0.76 0.79 0.79 0.74 0.76 0.70 0.75 0.72 0.78 0.74 0.78 0.76 > >>>>>> 0.80 > >>>>>> 0.76 > >>>>>> 0.79 0.77 0.76 0.76 0.75 0.73 0.80 0.78 0.81 0.78 0.84 0.76 0.79 > >>>>>> 0.77 > >>>>>> 0.79 > >>>>>> 0.76 0.83 0.77 0.79 0.77 0.84 0.76 0.80 0.77 0.81 0.72 0.79 0.74 > >>>>>> 0.80 > >>>>>> 0.77 > >>>>>> 0.83 0.81 0.82 0.76 0.81 0.75 0.78 0.74 0.75 0.71 0.75 0.73 0.77 > >>>>>> 0.77 > >>>>>> 0.80 > >>>>>> 0.77 0.80 0.75 0.73 0.77 0.72 0.76 0.72 0.76 0.72 0.72 0.69 0.66 > >>>>>> 0.68 > >>>>>> 0.65 > >>>>>> 0.68 0.64 0.70 0.62 0.70 0.62 0.69 0.67 0.74 0.75 0.74 > >>>>>> ++ Max displacement in automask = 0.84 (mm) at sub-brick 145 > >>>>>> ++ Wrote dataset to disk in 006/tmp.mc-afni2.2758/outvol.nii.gz > >>>>>> #@# -------------------------------- > >>>>>> > >>>>>> > /Users/maryam/mri-space/studies/PDM1/unpackdata/110419_xu_PDM1_02_PDM/bold > >>>>>> mri_convert 006/tmp.mc-afni2.2758/outvol.nii.gz 006/fmc.nii > >>>>>> --in_like > >>>>>> 005/f.nii > >>>>>> mri_convert 006/tmp.mc-afni2.2758/outvol.nii.gz 006/fmc.nii > >>>>>> --in_like > >>>>>> 005/f.nii > >>>>>> reading info from 005/f.nii... > >>>>>> $Id: mri_convert.c,v 1.146.2.5 2009/04/08 18:40:23 nicks Exp $ > >>>>>> reading from 006/tmp.mc-afni2.2758/outvol.nii.gz... > >>>>>> TR=2000.00, TE=0.00, TI=0.00, flip angle=0.00 > >>>>>> i_ras = (-1, -0.000171665, -3.15037e-05) > >>>>>> j_ras = (0.000174532, -0.983572, -0.180519) > >>>>>> k_ras = (-2.07852e-18, -0.180519, 0.983572) > >>>>>> writing to 006/fmc.nii... > >>>>>> ... cleaning up .... > >>>>>> Sun Apr 24 20:43:59 EDT 2011 > >>>>>> mc-afni2 done > >>>>>> INFO: Making external regressor from mc params > >>>>>> mcparams2extreg -mcfile 006/fmc.mcdat -northog 6 -extreg > >>>>>> 006/mcextreg > >>>>>> Logfile is 006/mcparams2extreg.log > >>>>>> ----------- Matlab file -------------- > >>>>>> mcfile = '006/fmc.mcdat'; > >>>>>> extreg = '006/mcextreg'; > >>>>>> nkeep = []; > >>>>>> northog = [6]; > >>>>>> pctorthog = []; > >>>>>> monly = 0; > >>>>>> > >>>>>> mc = textread(mcfile); > >>>>>> mc = mc(:,2:7); > >>>>>> ntrs = size(mc,1); > >>>>>> if(ntrs < 6) > >>>>>> fprintf('ERROR: ntrs = %d < 6\n',ntrs); > >>>>>> if(~monly) > >>>>>> fprintf('Quiting matlab\n'); > >>>>>> quit; > >>>>>> fprintf('should not be here\n'); > >>>>>> end > >>>>>> end > >>>>>> > >>>>>> if(~isempty(nkeep)) > >>>>>> x = mc(:,1:nkeep); > >>>>>> elseif(~isempty(northog)) > >>>>>> [u s v] = svd(mc); > >>>>>> ds = diag(s); > >>>>>> pct = 100*cumsum(ds)/sum(ds); > >>>>>> x = u(:,1:northog); > >>>>>> nkeep = northog; > >>>>>> fprintf('INFO: northog = %d, pct = %g\n',northog,pct(nkeep)); > >>>>>> else > >>>>>> [u s v] = svd(mc); > >>>>>> ds = diag(s); > >>>>>> pct = 100*cumsum(ds)/sum(ds); > >>>>>> nkeep = min(find(pct > pctorthog)); > >>>>>> x = u(:,1:nkeep); > >>>>>> fprintf('INFO: pctorthog = %g, nkeep = %d\n',pctorthog,nkeep); > >>>>>> end > >>>>>> > >>>>>> x2 = zeros(1,1,nkeep,ntrs); > >>>>>> x2(1,1,:,:) = x'; %' > >>>>>> %fmri_svbvolume(x2,extreg); > >>>>>> > >>>>>> mri.vol = permute(x2,[1 3 2 4]); > >>>>>> mri.tr = 0; > >>>>>> mri.flip_angle = 0; > >>>>>> mri.te = 0; > >>>>>> mri.ti = 0; > >>>>>> mri.vox2ras0 = eye(4); > >>>>>> mri.xsize = 1; > >>>>>> mri.ysize = 1; > >>>>>> mri.zsize = 1; > >>>>>> mri.volres = [1 1 1]; > >>>>>> mri.volsize = [size(mri.vol,1) size(mri.vol,2) size(mri.vol,3)]; > >>>>>> fname = sprintf('%s.bhdr',extreg); > >>>>>> MRIwrite(mri,fname); > >>>>>> > >>>>>> quit; > >>>>>> > >>>>>> ----------------------------------- > >>>>>> ------------------------------------------ > >>>>>> ------- matlab output -------------------- > >>>>>> Warning: Unable to open display 'iconic'. You will not be able to > >>>>>> display > >>>>>> graphics on the screen. > >>>>>> > >>>>>> < M A T L A B (R) > > >>>>>> Copyright 1984-2009 The MathWorks, Inc. > >>>>>> Version 7.9.0.529 (R2009b) 64-bit (maci64) > >>>>>> August 12, 2009 > >>>>>> > >>>>>> To get started, type one of these: helpwin, helpdesk, or demo. > >>>>>> For product information, visit www.mathworks.com. > >>>>>> > >>>>>> >> >> >> >> >> >> >> >> ??? Error using ==> textread at 167 > >>>>>> File not found. > >>>>>> > >>>>>> >> ??? Undefined variable mc. > >>>>>> > >>>>>> >> ??? Undefined function or variable 'mc'. > >>>>>> > >>>>>> >> ??? Undefined function or variable 'ntrs'. > >>>>>> > >>>>>> >> >> ??? Undefined function or variable 'mc'. > >>>>>> > >>>>>> >> >> ??? Undefined function or variable 'ntrs'. > >>>>>> > >>>>>> >> ??? Undefined function or variable 'x'. > >>>>>> > >>>>>> >> >> >> ??? Undefined function or variable 'x2'. > >>>>>> > >>>>>> >> >> >> >> >> >> >> >> >> >> ??? Reference to non-existent field > >>>>>> 'vol'. > >>>>>> > >>>>>> >> >> ??? Reference to non-existent field 'vol'. > >>>>>> > >>>>>> Error in ==> MRIwrite at 70 > >>>>>> bmri.voldim = [size(mri.vol,1) size(mri.vol,2) size(mri.vol,3)]; > >>>>>> > >>>>>> >> >> ------------------------------------------ > >>>>>> ERROR: output not created > >>>>>> ERROR: mc-sess failed > >>>>>>

1 0

Jump to page:

2026

2025

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

2009

2008

2007

2006

2005

2004

Freesurfer