Hello, Freesurfer group, I registered the example_func.nii.gz to brain.mgz in tkregister2 using FS/anat2exf.register.dat, and --fslregout fsl.feat/reg/example_func2highres.mat. I have a couple of questions.
1) could anybody confirm that using--fslregout does apply .mat file for registration based on the following information on the screen? surfname set to orig target volume PT01/mri/brain.mgz movable volume Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/example_func.nii.gz reg file Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/reg/freesurfer/anat2exf.register.dat LoadVol 1 ZeroCRAS 0 $Id: tkregister2.c,v 1.86.2.2 2008/03/10 14:02:35 nicks Exp $ Diagnostic Level -1 ---- Input registration matrix -------- 0.995 -0.028 0.093 4.032; -0.091 0.069 0.994 1.194; -0.035 -0.997 0.066 0.264; 0.000 0.000 0.000 1.000; float2int = 0 --------------------------------------- INFO: loading target PT01/mri/brain.mgz INFO: changing target type to float Ttarg: -------------------- -1.000 0.000 0.000 128.000; 0.000 0.000 1.000 -128.000; 0.000 -1.000 0.000 128.000; 0.000 0.000 0.000 1.000; INFO: loading movable Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/example_func.nii.gz Tmov: -------------------- -3.438 0.000 0.000 110.000; 0.000 0.000 3.000 -60.000; 0.000 -3.438 0.000 110.000; 0.000 0.000 0.000 1.000; mkheaderreg = 0, float2int = 0 ---- Input registration matrix -------- 0.995 -0.028 0.093 4.032; -0.091 0.069 0.994 1.194; -0.035 -0.997 0.066 0.264; 0.000 0.000 0.000 1.000; subject = PT01 Reading lh surface orig Done reading surface Reading rh surface orig Done reading surface Zoom Factor = 2 FOV = 256 Opening window PT01 NVIDIA: could not open the device file /dev/nvidiactl (Permission denied). NVIDIA: Direct rendering failed; attempting indirect rendering. Setting scale Tagging mov volume tkregister2: interface: /net/dysthymia/freesurfer/lib/tcl/tkregister2.tcl tkregister.tcl: startup done tkregister.tcl: default macro interface (to change: macro,mini,micro) tkregister.tcl: or: setenv tkregisterinterface {macro,mini,micro} %
2) The results are quite the same (and good) in the tkregister2 with a slightly difference as shown below. Could anyone please tell me what these numbers mean?
with .bat file Screen: 256 257 (113,113,113), inorm = 1, mov_frame = 0 Anat: (128 126 127) ( 0.0 -1.0 2.0) 25.0000 32.7 Func: ( 31 32 21) ( 4.2 3.1 1.4) 756.0000 112.5
with .mat file Screen: 257 257 (175,175,175), inorm = 1, mov_frame = 0 Anat: (128 126 127) ( 0.0 -1.0 2.0) 25.0000 28.3 Func: ( 31 32 21) ( 4.2 3.1 1.4) 756.0000 106.2
Thanks,
Xin
-----Original Message----- From: Douglas N Greve [mailto:greve@nmr.mgh.harvard.edu] Sent: Mon 1/12/2009 4:02 PM To: Wang, Xin Cc: analysis-bugs@nmr.mgh.harvard.edu Subject: Re: [Freesurfer] reg-feat2anat failed
reg-feat2anat calls fslregister, so they should produce the same result. There is no need for an "intermediate" volume since you are registering from the subject's EPI to the subject's structure. Does the registration not look good?
Wang, Xin wrote:
Thank you, Dr. Greve. It seems work ok. what is the difference between reg-feat2anat and fslregister? although both commands are working for me, I don't think either of them serve my propose. Freesurfer does not use initial_highres to improve the reg of example_func to highres, but FSL/FEAT uses this inter-mediate step. I try to use registration.mat files of FSL/FEAT step-wise registrations to improve the reg. I think flag --fsl allow me to use a fsl.mat file as reg file. How to make it work in the way I need?
Thank you again,
Xin
*From:* Douglas N Greve [mailto:greve@nmr.mgh.harvard.edu] *Sent:* Mon 1/12/2009 3:29 PM *To:* Wang, Xin *Cc:* analysis-bugs@nmr.mgh.harvard.edu *Subject:* Re: [Freesurfer] reg-feat2anat failed
when you say everything looked ok, do you mean that the orientation was correct? If so, try running this:
fslregister --s subject --mov example_func.nii.gz --reg newreg.dat
then
tkregister2 --mov example_func.nii.gz --reg newreg.dat --surfs
Wang, Xin wrote:
Sorry, Dr. Greve. I copied addtional parts after I killed the tkmedit. The following is on the screen when I saw the image.
[]$ tkmedit -f Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/example_func.nii.gz Reading /net/dysthymia/freesurfer/lib/tcl/tkm_common.tcl Reading /net/dysthymia/freesurfer/lib/tcl/tkm_wrappers.tcl Reading /net/dysthymia/freesurfer/lib/tcl/fsgdfPlot.tcl Reading /net/dysthymia/freesurfer/lib/tcl/tkUtils.tcl %
*From:* Douglas N Greve [mailto:greve@nmr.mgh.harvard.edu] *Sent:* Mon 1/12/2009 3:10 PM *To:* Wang, Xin *Cc:* analysis-bugs@nmr.mgh.harvard.edu *Subject:* Re: [Freesurfer] reg-feat2anat failed
I don't understand. If it segfaulted how did you look at the images?
Wang, Xin wrote:
Thank you, Dr. Greve. This time tkmedit read the image, and the image looks ok to me. But following information appeared.
[]$ tkmedit -f Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/example_func.nii.gz Reading /net/dysthymia/freesurfer/lib/tcl/tkm_common.tcl Reading /net/dysthymia/freesurfer/lib/tcl/tkm_wrappers.tcl Reading /net/dysthymia/freesurfer/lib/tcl/fsgdfPlot.tcl Reading /net/dysthymia/freesurfer/lib/tcl/tkUtils.tcl % Killed
=================================== ERROR: A segfault has occurred. This is not your fault, : but is most likely an unrecoverable error and has : made the program unstable. : : Please send the contents of the file .xdebug_tkmedit : that should be in this directory to analysis-bugs@nmr.mgh.harvard.edu mailto:analysis-bugs@nmr.mgh.harvard.edu : : Now exiting... : []$ cat .xdebug_tkmedit tkmedit started: Mon Jan 12 15:14:52 2009 tkmedit.bin -f Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/example_func.nii.gz $Id: tkmedit.c,v 1.320.2.5 2008/03/10 14:02:35 nicks Exp $ $Name: stable4 $ Set user home dir to /net/data4/CCMB08/PTSD/Thickness Set subject home dir to Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/example_func.nii.gz Using interface file /net/dysthymia/freesurfer/lib/tcl/tkmedit.tcl Tcl: Successfully parsed tkmedit.tcl Tcl: Successfully parsed tkm_functional.tcl Segfault Entering main loop xDebug stack (length: 2) 01: main() 01: Entering main loop 00: 00:
*From:* Douglas N Greve [mailto:greve@nmr.mgh.harvard.edu] *Sent:* Mon 1/12/2009 2:55 PM *To:* Wang, Xin *Cc:* analysis-bugs@nmr.mgh.harvard.edu *Subject:* Re: [Freesurfer] reg-feat2anat failed
what is your tkmedit cmd?
It should be
tkmedit -f example_func.nii.gz
Wang, Xin wrote:
Hello, Dr. Greve, tkmedit can not open the example_func.nii.gzI attached the header of example_func.nii.gz (after mri_convert), and the contents of .xdebug_tkmedit file
[]$ fslhd
Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/example_func.nii.gz
filename Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/example_func.nii.gz sizeof_hdr 348 data_type FLOAT32 dim0 3 dim1 64 dim2 64 dim3 40 dim4 1 dim5 1 dim6 1 dim7 1 vox_units mm time_units s datatype 16 nbyper 4 bitpix 32 pixdim0 0.0000000000 pixdim1 3.4375000000 pixdim2 3.4375000000 pixdim3 3.0000000000 pixdim4 1.0000000000 pixdim5 1.0000000000 pixdim6 1.0000000000 pixdim7 1.0000000000 vox_offset 352 cal_max 0.0000 cal_min 0.0000 scl_slope 0.000000 scl_inter 0.000000 phase_dim 0 freq_dim 0 slice_dim 0 slice_name Unknown slice_code 0 slice_start 0 slice_end 0 slice_duration 0.000000 time_offset 0.000000 intent Unknown intent_code 0 intent_name intent_p1 0.000000 intent_p2 0.000000 intent_p3 0.000000 qform_name Scanner Anat qform_code 1 qto_xyz:1 -3.437500 0.000000 -0.000000 220.000000 qto_xyz:2 0.000000 3.437500 -0.000000 0.000000 qto_xyz:3 0.000000 0.000000 3.000000 0.000000 qto_xyz:4 0.000000 0.000000 0.000000 1.000000 qform_xorient Right-to-Left qform_yorient Posterior-to-Anterior qform_zorient Inferior-to-Superior sform_name Scanner Anat sform_code 1 sto_xyz:1 -3.437500 0.000000 0.000000 220.000000 sto_xyz:2 0.000000 3.437500 0.000000 0.000000 sto_xyz:3 0.000000 0.000000 3.000000 0.000000 sto_xyz:4 0.000000 0.000000 0.000000 1.000000 sform_xorient Right-to-Left sform_yorient Posterior-to-Anterior sform_zorient Inferior-to-Superior file_type NIFTI-1+ file_code 1 descrip FreeSurfer Jun 14 2008 aux_file .xdebug_tkmedit tkmedit started: Mon Jan 12 15:05:46 2009 tkmedit.bin Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/example_func.nii.gz $Id: tkmedit.c,v 1.320.2.5 2008/03/10 14:02:35 nicks Exp $ $Name: stable4 $ Set user home dir to /net/data4/CCMB08/PTSD/Thickness ERROR: Parsing arguments Subject and image not declared Fatal error in parsing command line args. Segfault Copying argv[1] xDebug stack (length: 3) 02: ParseCmdLineArgs( argc=2, argv=tkmedit.bin ) 02: Copying argv[1] 01: main() 01: Parsing command line arguments 00: 00:
*From:* Douglas N Greve [mailto:greve@nmr.mgh.harvard.edu] *Sent:* Mon 1/12/2009 2:38 PM *To:* Wang, Xin *Cc:* analysis-bugs@nmr.mgh.harvard.edu *Subject:* Re: [Freesurfer] reg-feat2anat failed
If the orientation of the example func in tkmedit is wrong, then the reg-feat2anat is destined to fail. tkmedit (usually) reads the nifti header correctly. How do you know the header is right?
Wang, Xin wrote:
My problem about tkregister2 is that the flag -fsl example_func2highres.mat does not work properly. The tkregister2 catastrophically failed to register two brains. I submit the bugreport of tkregister2 (please see the COMMAND
3). The
two early steps (mri_convert and reg-feat2anat) are included for
your
references.
Thanks,
Xin
FREESURFER_HOME: /net/dysthymia/freesurfer
Build stamp: freesurfer-Linux-centos4_x86_64-stable-pub-v4.0.5
RedHat release: Red Hat Enterprise Linux Client release 5.2
(Tikanga)
Kernel info: Linux 2.6.18-92.1.1.el5 x86_64
Subject name: PT01
COMMAND 1: []$ mri_convert Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/example_func.nii.gz Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/example_func.nii.gz mri_convert Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/example_func.nii.gz Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/example_func.nii.gz $Id: mri_convert.c,v 1.146.2.2 2008/03/02 02:53:20 nicks Exp $ reading from Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/example_func.nii.gz... TR=1000.00, TE=0.00, TI=0.00, flip angle=0.00 i_ras = (-1, 0, 0) j_ras = (0, 1, 0) k_ras = (0, 0, 1) writing to Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/example_func.nii.gz...
COMMAND 2: []$ reg-feat2anat --feat
Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat
--subject PT01 FeatDir is Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat template volume is Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/example_func.nii.gz log file is
Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/reg/freesurfer/reg-feat2anat.log
Mon Jan 12 14:08:56 CST 2009 /net/data4/CCMB08/PTSD/Thickness /net/dysthymia/freesurfer/bin/reg-feat2anat --feat Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat --subject PT01 Linux schizo.psych.med.umich.edu 2.6.18-92.1.1.el5 #1 SMP Thu
May 22
09:01:47 EDT 2008 x86_64 x86_64 x86_64 GNU/Linux example_func ras_good_flag 1
Initializing exf2anat with header tkregister2_cmdl --targ /net/data4/CCMB08/PTSD/Thickness/PT01/mri/brainmask.mgz --mov Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/example_func.nii.gz
--reg
Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/reg/freesurfer/anat2exf.register.dat.init
--fslregout
Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/reg/freesurfer/exf2anat.init.fsl.mat
--regheader --noedit target volume
/net/data4/CCMB08/PTSD/Thickness/PT01/mri/brainmask.mgz
movable volume Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/example_func.nii.gz reg file
Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/reg/freesurfer/anat2exf.register.dat.init
LoadVol 0 ZeroCRAS 0 $Id: tkregister2.c,v 1.86.2.2 2008/03/10 14:02:35 nicks Exp $ Diagnostic Level -1 INFO: loading target /net/data4/CCMB08/PTSD/Thickness/PT01/mri/brainmask.mgz Ttarg: -------------------- -1.000 0.000 0.000 128.000; 0.000 0.000 1.000 -128.000; 0.000 -1.000 0.000 128.000; 0.000 0.000 0.000 1.000; INFO: loading movable Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/example_func.nii.gz Tmov: -------------------- -3.438 0.000 0.000 110.000; 0.000 0.000 3.000 -60.000; 0.000 -3.438 0.000 110.000; 0.000 0.000 0.000 1.000; mkheaderreg = 1, float2int = 0 Computing reg from header (and possibly input matrix) ---- Input registration matrix (computed) -------- 1.000 0.000 0.000 -110.507; 0.000 0.000 1.000 -41.354; 0.000 -1.000 0.000 93.531; 0.000 0.000 0.000 1.000;
---- Input registration matrix -------- 1.000 0.000 0.000 -110.507; 0.000 0.000 1.000 -41.354; 0.000 -1.000 0.000 93.531; 0.000 0.000 0.000 1.000; subject = subject-unknown RegMat --------------------------- 1.000 0.000 0.000 -110.507; 0.000 0.000 1.000 -41.354; 0.000 -1.000 0.000 93.531; 0.000 0.000 0.000 1.000; register: file
Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/reg/freesurfer/anat2exf.register.dat.init
written FSLOUTPUTTYPE NIFTI_GZ tkreg2FSL: mov det = -35.4492, ref det = -1 fslregister --s PT01 --reg
Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/reg/freesurfer/anat2exf.register.dat
--mov
Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/example_func.nii.gz
--dof 6 --bins 256 --cost corratio --fslmat
Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/reg/freesurfer/exf2anat.fsl.mat
--fsvol brainmask --maxangle 90 --initfslmat
Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/reg/freesurfer/exf2anat.init.fsl.mat
--tmp
Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/reg/freesurfer/tmp
--out
Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/reg/freesurfer/tmp/exf-in-anat.nii.gz
Log file is
Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/reg/freesurfer/anat2exf.register.dat.fslregister.log
Mon Jan 12 14:08:56 CST 2009 --s PT01 --reg
Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/reg/freesurfer/anat2exf.register.dat
--mov
Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/example_func.nii.gz
--dof 6 --bins 256 --cost corratio --fslmat
Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/reg/freesurfer/exf2anat.fsl.mat
--fsvol brainmask --maxangle 90 --initfslmat
Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/reg/freesurfer/exf2anat.init.fsl.mat
--tmp
Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/reg/freesurfer/tmp
--out
Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/reg/freesurfer/tmp/exf-in-anat.nii.gz
$Id: fslregister,v 1.16 2007/07/26 16:52:54 greve Exp $ schizo.psych.med.umich.edu Linux schizo.psych.med.umich.edu 2.6.18-92.1.1.el5 #1 SMP Thu
May 22
09:01:47 EDT 2008 x86_64 x86_64 x86_64 GNU/Linux
/net/data4/CCMB08/PTSD/Thickness mri_convert
/net/data4/CCMB08/PTSD/Thickness/PT01/mri/brainmask.mgz
Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/reg/freesurfer/tmp/refvol.fslregister.nii
mri_convert
/net/data4/CCMB08/PTSD/Thickness/PT01/mri/brainmask.mgz
Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/reg/freesurfer/tmp/refvol.fslregister.nii
$Id: mri_convert.c,v 1.146.2.2 2008/03/02 02:53:20 nicks Exp $ reading from
/net/data4/CCMB08/PTSD/Thickness/PT01/mri/brainmask.mgz...
TR=12.27, TE=5.21, TI=650.00, flip angle=9.00 i_ras = (-1, 0, 0) j_ras = (0, 0, -1) k_ras = (0, 1, 0) writing to
Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/reg/freesurfer/tmp/refvol.fslregister.nii...
/net/data4/CCMB08/PTSD/Thickness mri_convert Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/example_func.nii.gz
Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/reg/freesurfer/tmp/movvol.fslregister.nii
--frame 0 mri_convert Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/example_func.nii.gz
Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/reg/freesurfer/tmp/movvol.fslregister.nii
--frame 0 $Id: mri_convert.c,v 1.146.2.2 2008/03/02 02:53:20 nicks Exp $ reading from Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/example_func.nii.gz... TR=1000.00, TE=0.00, TI=0.00, flip angle=0.00 i_ras = (-1, 0, 0) j_ras = (0, 1, 0) k_ras = (0, 0, 1) keeping frame 0 writing to
Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/reg/freesurfer/tmp/movvol.fslregister.nii...
Mon Jan 12 14:09:00 CST 2009 /net/data4/CCMB08/PTSD/Thickness bet
Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/reg/freesurfer/tmp/movvol.fslregister.nii
Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/reg/freesurfer/tmp/movvol.fslregister.nii
-f 0.1 Mon Jan 12 14:09:04 CST 2009 /net/data4/CCMB08/PTSD/Thickness flirt -ref
Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/reg/freesurfer/tmp/refvol.fslregister.nii
-in
Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/reg/freesurfer/tmp/movvol.fslregister.nii
-omat
Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/reg/freesurfer/exf2anat.fsl.mat
-cost corratio -dof 6 -searchrx -90 90 -searchry -90 90
-searchrz -90
90 -init
Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/reg/freesurfer/exf2anat.init.fsl.mat
-o
Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/reg/freesurfer/tmp/outvol.fslregister.nii
Mon Jan 12 14:09:20 CST 2009 /net/data4/CCMB08/PTSD/Thickness mri_convert
Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/reg/freesurfer/tmp/outvol.fslregister.nii
Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/reg/freesurfer/tmp/exf-in-anat.nii.gz
mri_convert
Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/reg/freesurfer/tmp/outvol.fslregister.nii
Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/reg/freesurfer/tmp/exf-in-anat.nii.gz
$Id: mri_convert.c,v 1.146.2.2 2008/03/02 02:53:20 nicks Exp $ reading from
Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/reg/freesurfer/tmp/outvol.fslregister.nii...
TR=1000.00, TE=0.00, TI=0.00, flip angle=0.00 i_ras = (-1, 0, 0) j_ras = (0, 0, -1) k_ras = (0, 1, 0) writing to
Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/reg/freesurfer/tmp/exf-in-anat.nii.gz...
tkregister2_cmdl --s PT01 --mov Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/example_func.nii.gz
--reg
Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/reg/freesurfer/anat2exf.register.dat
--fslreg
Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/reg/freesurfer/exf2anat.fsl.mat
--noedit ---- FSL registration matrix -------- 0.995 -0.035 0.091 20.763; 0.093 0.066 -0.994 171.738; 0.028 0.997 0.069 11.368; 0.000 0.000 0.000 1.000;
INFO: no target volume specified, assuming FreeSurfer orig volume. target volume orig movable volume Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/example_func.nii.gz reg file
Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/reg/freesurfer/anat2exf.register.dat
LoadVol 0 ZeroCRAS 0 $Id: tkregister2.c,v 1.86.2.2 2008/03/10 14:02:35 nicks Exp $ Diagnostic Level -1 INFO: loading target
/net/data4/CCMB08/PTSD/Thickness/PT01/mri/orig.mgz
Ttarg: -------------------- -1.000 0.000 0.000 128.000; 0.000 0.000 1.000 -128.000; 0.000 -1.000 0.000 128.000; 0.000 0.000 0.000 1.000; INFO: loading movable Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/example_func.nii.gz Tmov: -------------------- -3.438 0.000 0.000 110.000; 0.000 0.000 3.000 -60.000; 0.000 -3.438 0.000 110.000; 0.000 0.000 0.000 1.000; mkheaderreg = 0, float2int = 0 FSLOUTPUTTYPE NIFTI fsl2TkReg: mov det = -35.4492, ref det = -1 ---- Input registration matrix (computed) -------- 0.995 -0.028 0.093 4.032; -0.091 0.069 0.994 1.194; -0.035 -0.997 0.066 0.264; 0.000 0.000 0.000 1.000;
---- Input registration matrix -------- 0.995 -0.028 0.093 4.032; -0.091 0.069 0.994 1.194; -0.035 -0.997 0.066 0.264; 0.000 0.000 0.000 1.000; subject = PT01 RegMat --------------------------- 0.995 -0.028 0.093 4.032; -0.091 0.069 0.994 1.194; -0.035 -0.997 0.066 0.264; 0.000 0.000 0.000 1.000; register: file
Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/reg/freesurfer/anat2exf.register.dat
written (standard_in) 1: illegal character: P (standard_in) 1: illegal character: T
Started at Mon Jan 12 14:08:57 CST 2009 Ended at Mon Jan 12 14:09:24 CST 2009
fslregister Done
To check results, run: tkregister2 --mov Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/example_func.nii.gz
--reg
Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/reg/freesurfer/anat2exf.register.dat
--surf orig
tkregister2_cmdl --targ /net/misc/fsl/etc/standard/avg152T1_brain.nii.gz --mov Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/example_func.nii.gz --fslreg
Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/reg/example_func2standard.mat
--reg
Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/reg/freesurfer/std2exf.register.dat
--noedit ---- FSL registration matrix -------- 1.122 0.029 0.006 -34.123; -0.061 1.041 0.128 -3.403; -0.010 -0.119 1.224 25.844; 0.000 0.000 0.000 1.000;
target volume /net/misc/fsl/etc/standard/avg152T1_brain.nii.gz movable volume Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/example_func.nii.gz reg file
Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/reg/freesurfer/std2exf.register.dat
LoadVol 0 ZeroCRAS 0 $Id: tkregister2.c,v 1.86.2.2 2008/03/10 14:02:35 nicks Exp $ Diagnostic Level -1 INFO: loading target
/net/misc/fsl/etc/standard/avg152T1_brain.nii.gz
INFO: target does not conform to COR format, so I'm going to reslice to COR. This will not affect the final registration. Ttarg: -------------------- -1.000 0.000 0.000 128.000; 0.000 0.000 1.000 -128.000; 0.000 -1.000 0.000 128.000; 0.000 0.000 0.000 1.000; INFO: loading movable Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/example_func.nii.gz Tmov: -------------------- -3.438 0.000 0.000 110.000; 0.000 0.000 3.000 -60.000; 0.000 -3.438 0.000 110.000; 0.000 0.000 0.000 1.000; mkheaderreg = 0, float2int = 0 FSLOUTPUTTYPE NIFTI_GZ fsl2TkReg: mov det = -35.4492, ref det = -8 ---- Input registration matrix (computed) -------- 0.890 0.025 0.002 1.505; -0.013 0.092 0.807 4.483; 0.051 -0.947 0.099 3.604; 0.000 0.000 0.000 1.000;
---- Input registration matrix -------- 0.890 0.025 0.002 1.505; -0.013 0.092 0.807 4.483; 0.051 -0.947 0.099 3.604; 0.000 0.000 0.000 1.000; subject = subject-unknown RegMat --------------------------- 0.890 0.002 -0.025 1.505; -0.013 0.807 -0.092 4.483; 0.051 0.099 0.947 3.604; 0.000 0.000 0.000 1.000; register: file
Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/reg/freesurfer/std2exf.register.dat
written mri_matrix_multiply -fsl -im
Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/reg/freesurfer/exf2anat.fsl.mat
-iim
Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/reg/example_func2standard.mat
-om
Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/reg/freesurfer/std2anat.fsl.mat
tkregister2_cmdl --s PT01 --targ /net/data4/CCMB08/PTSD/Thickness/PT01/mri/brainmask.mgz --mov /net/misc/fsl/etc/standard/avg152T1_brain.nii.gz --reg
Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/reg/freesurfer/anat2std.register.dat
--fslreg
Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/reg/freesurfer/std2anat.fsl.mat
--noedit ---- FSL registration matrix -------- 0.885 -0.049 0.075 48.869; 0.074 -0.031 -0.809 195.041; 0.077 0.950 -0.044 18.348; 0.000 0.000 0.000 1.000;
target volume
/net/data4/CCMB08/PTSD/Thickness/PT01/mri/brainmask.mgz
movable volume /net/misc/fsl/etc/standard/avg152T1_brain.nii.gz reg file
Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/reg/freesurfer/anat2std.register.dat
LoadVol 0 ZeroCRAS 0 $Id: tkregister2.c,v 1.86.2.2 2008/03/10 14:02:35 nicks Exp $ Diagnostic Level -1 INFO: loading target /net/data4/CCMB08/PTSD/Thickness/PT01/mri/brainmask.mgz Ttarg: -------------------- -1.000 0.000 0.000 128.000; 0.000 0.000 1.000 -128.000; 0.000 -1.000 0.000 128.000; 0.000 0.000 0.000 1.000; INFO: loading movable
/net/misc/fsl/etc/standard/avg152T1_brain.nii.gz
Tmov: -------------------- -2.000 0.000 0.000 91.000; 0.000 0.000 2.000 -91.000; 0.000 -2.000 0.000 109.000; 0.000 0.000 0.000 1.000; mkheaderreg = 0, float2int = 0 FSLOUTPUTTYPE NIFTI_GZ fsl2TkReg: mov det = -8, ref det = -1 ---- Input registration matrix (computed) -------- 1.116 -0.061 0.100 2.757; -0.105 -0.035 1.225 -4.398; -0.085 -1.046 -0.064 -3.213; 0.000 0.000 0.000 1.000;
---- Input registration matrix -------- 1.116 -0.061 0.100 2.757; -0.105 -0.035 1.225 -4.398; -0.085 -1.046 -0.064 -3.213; 0.000 0.000 0.000 1.000; subject = PT01 RegMat --------------------------- 1.116 -0.061 0.100 2.757; -0.105 -0.035 1.225 -4.398; -0.085 -1.046 -0.064 -3.213; 0.000 0.000 0.000 1.000; register: file
Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/reg/freesurfer/anat2std.register.dat
written
To visually check your results, run: tkregister2_cmdl --s PT01 --mov Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/example_func.nii.gz
--reg
Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/reg/freesurfer/anat2exf.register.dat
--fslreg
Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/reg/freesurfer/exf2anat.fsl.mat
--surf orig --tag
Mon Jan 12 14:09:25 CST 2009 reg-feat2anat done
COMMAND 3: []$ tkregister2 --s PT01 --mov Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/example_func.nii.gz
--reg
call/SM5DF1212/PT01_Recall_SM5DF1212.feat/reg/freesurfer/anat2exf.register.dat
--fslreg
Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/reg/freesurfer/exf2anat.fsl.mat
--surf orig --tag --fsl
Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/reg/example_func2highres.mat
---- FSL registration matrix -------- 0.995 -0.035 0.091 20.763; 0.093 0.066 -0.994 171.738; 0.028 0.997 0.069 11.368; 0.000 0.000 0.000 1.000;
surfname set to orig ---- FSL registration matrix -------- 1.043 -0.038 0.042 18.354; 0.020 1.036 0.060 9.413; -0.046 -0.068 1.128 70.358; 0.000 0.000 0.000 1.000;
INFO: no target volume specified, assuming FreeSurfer orig volume. target volume orig movable volume Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/example_func.nii.gz reg file
Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/reg/freesurfer/anat2exf.register.dat
LoadVol 1 ZeroCRAS 0 $Id: tkregister2.c,v 1.86.2.2 2008/03/10 14:02:35 nicks Exp $ Diagnostic Level -1 INFO: loading target
/net/data4/CCMB08/PTSD/Thickness/PT01/mri/orig.mgz
INFO: changing target type to float Ttarg: -------------------- -1.000 0.000 0.000 128.000; 0.000 0.000 1.000 -128.000; 0.000 -1.000 0.000 128.000; 0.000 0.000 0.000 1.000; INFO: loading movable Recall/SM5DF1212/PT01_Recall_SM5DF1212.feat/example_func.nii.gz Tmov: -------------------- -3.438 0.000 0.000 110.000; 0.000 0.000 3.000 -60.000; 0.000 -3.438 0.000 110.000; 0.000 0.000 0.000 1.000; mkheaderreg = 0, float2int = 0 FSLOUTPUTTYPE NIFTI_GZ fsl2TkReg: mov det = -35.4492, ref det = -1 ---- Input registration matrix (computed) -------- 0.956 0.037 0.032 3.452; -0.038 0.882 -0.059 1.969; -0.021 0.050 0.961 1.159; 0.000 0.000 0.000 1.000;
---- Input registration matrix -------- 0.956 0.037 0.032 3.452; -0.038 0.882 -0.059 1.969; -0.021 0.050 0.961 1.159; 0.000 0.000 0.000 1.000; subject = PT01 Reading lh surface orig Done reading surface Reading rh surface orig Done reading surface Zoom Factor = 2 FOV = 256 Opening window PT01 Setting scale Tagging mov volume tkregister2: interface:
/net/dysthymia/freesurfer/lib/tcl/tkregister2.tcl
tkregister.tcl: startup done tkregister.tcl: default macro interface (to change:
macro,mini,micro)
tkregister.tcl: or: setenv tkregisterinterface {macro,mini,micro} %
-----Original Message----- From: Doug Greve [mailto:greve@nmr.mgh.harvard.edu] Sent: Tue 1/6/2009 12:30 PM To: Wang, Xin Cc: freesurfer Subject: Re: [Freesurfer] reg-feat2anat failed
Hi Xin, I'm not sure what's going wrong here, I think I need more
info.
Can you follow the steps in surfer.nmr.mgh.harvard.edu/fswiki/BugReporting?
doug
Wang, Xin wrote:
hello, Freesurfer group,
I try to use registration.mat files of FSL/FEAT routines,
instead of
anat2example_func.bat, to register the orig.mgz and example_func.nii.gz with the thinking that using
initial_highres in
FEAT may improve the reg. However, tkregister2 catastrophically
failed
to register two brains when I use --fsl flag. Modifying
header of
example_func with mri_convert does not solve the problem. I will appreciate any suggestion on using .mat file for registration.
second question: what is the best registration of
example_func and
orig.mgz? In another word, what should I check besides things
in the
FS tutorial when I check/modify the registration in
tkregister2? Is
there golden reference points I can use?
Thank you in advance for any suggestion.
Xin wang
-----Original Message----- From: Doug Greve [mailto:greve@nmr.mgh.harvard.edu] Sent: Mon 12/15/2008 1:26 PM To: Wang, Xin Cc: freesurfer Subject: Re: [Freesurfer] reg-feat2anat failed
There is a flag in the nifti header saying whether the qform or
sform
are valid, and it looks like your scripts do not set this flag.
You can
use mri_convert (the freesurfer program). This will create valid
files
that can be read into any program.
doug
Wang, Xin wrote:
> Thank you very much for your quick reply, Dr. Greve. I
asked our
> institutional MRI center for correct information on qform and
sform
> since the header of nifti files are created by their scripts.
But I do
> not know what they can do to fix this problem. I will
appreciate if
> you have any further suggestion on how to do the registration
without
> the information on qform and sform in case I can not get them.
Does
> convert to ANALYZE format or manually correct in tkregister2
work?
> > Thank you again, > > Xin > > > -----Original Message----- > From: Doug Greve [mailto:greve@nmr.mgh.harvard.edu] > Sent: Mon 12/15/2008 1:09 PM > To: Wang, Xin > Cc: freesurfer > Subject: Re: [Freesurfer] reg-feat2anat failed > > Hi Xin, the problem is that neither the sform nor the
qform are
valid.
> Did the message below not appear when you ran tkmedit? > > WARNING: neither NIfTI-1 qform or sform are valid > WARNING: your volume will probably be incorrectly oriented > > Fix that, and everything will work fine. > > doug > > > Wang, Xin wrote: > > > Hello, Freesurfer group, > > I am registering the fMRI results to surface using
reg-feat2anat (v
> > 1.26.2.1 2007/08/15) and tkregister2. The automatic
registration
> > catastrophically failed in the subject space, but worked
fine
in the
> > standard space. The orientations of example_func.nii.gz
in the
subject
> > space are flipped in any view. I use the tkmedit to view the > > example_func, and found the orientation is wrong. It seems
Freesurfer
> > can not read nifti header correctly. I attach the
reg-feat2anat log
> > file and an example_func.nii.gz. > > I found the following message from Freesurfer email list. > >
http://www.mail-archive.com/freesurfer@nmr.mgh.harvard.edu/msg05573.html
> > Does Freesurfer have a solution for this type of problems by
now?
> > > > Thank you in advance, > > > > Xin > > >
> > > >_______________________________________________ > >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 > Phone Number: 617-724-2358 > Fax: 617-726-7422 > > In order to help us help you, please follow the steps in: > surfer.nmr.mgh.harvard.edu/fswiki/BugReporting > > >
-- Douglas N. Greve, Ph.D. MGH-NMR Center greve@nmr.mgh.harvard.edu Phone Number: 617-724-2358 Fax: 617-726-7422
In order to help us help you, please follow the steps in: surfer.nmr.mgh.harvard.edu/fswiki/BugReporting
-- Douglas N. Greve, Ph.D. MGH-NMR Center greve@nmr.mgh.harvard.edu Phone Number: 617-724-2358 Fax: 617-726-7422
In order to help us help you, please follow the steps in: surfer.nmr.mgh.harvard.edu/fswiki/BugReporting