Dear Freesurfer experts:
My project objective is to build 3D reconstruction of brain with inserted electrodes. I want register the Brian volume of CT and MR via freesurfer. The registration result looks great when I used th "mri_coreg" instruction, but when I tried to use the information in lta data to register by myself, some trouble happens.
The LTA data -------------------------------------------------------------------------------------------------------------------------------------------------------------- # transform file p19_movCT2refMR.lta # created by caca on Tue Jan 24 17:45:08 2017
type = 0 # LINEAR_VOX_TO_VOX nxforms = 1 mean = 0.0000 0.0000 0.0000 sigma = 10000.0000 1 4 4 4.795561134815216e-01 -1.318652089685202e-02 -2.057363465428352e-02 1.383309745788574e+01 1.247252989560366e-02 4.785549342632294e-01 -6.737518310546875e-02 1.315915203094482e+01 1.107159163802862e-02 3.306149691343307e-02 9.974335432052612e-01 -6.133582305908203e+01 0.000000000000000e+00 0.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00 src volume info valid = 1 # volume info valid filename = patient19CT.nii volume = 512 512 228 voxelsize = 4.873275756835938e-01 4.873275756835938e-01 9.999971985816956e-01 xras = -9.993908405303955e-01 -0.000000000000000e+00 3.489949554204941e-02 yras = -0.000000000000000e+00 1.000000000000000e+00 -0.000000000000000e+00 zras = 3.489949554204941e-02 0.000000000000000e+00 9.993908405303955e-01 cras = 7.397460937500000e-02 -1.105121917724609e+02 1.181088867187500e+02 dst volume info valid = 1 # volume info valid filename = patient19MR.nii volume = 256 256 192 voxelsize = 1.015599966049194e+00 1.015599966049194e+00 1.000001072883606e+00 xras = -1.000000000000000e+00 -0.000000000000000e+00 0.000000000000000e+00 yras = -0.000000000000000e+00 1.000000000000000e+00 -0.000000000000000e+00 zras = -0.000000000000000e+00 -0.000000000000000e+00 1.000000000000000e+00 cras = 3.934204101562500e+00 2.904579925537109e+01 1.511990356445312e+01 subject unknown -------------------------------------------------------------------------------------------------------------------------------------------------------------------
First I thought the 4x4 matrix is the affine matrix (type = VOX_TO_VOX)
4.795561134815216e-01 -1.318652089685202e-02 -2.057363465428352e-02 1.383309745788574e+01 1.247252989560366e-02 4.785549342632294e-01 -6.737518310546875e-02 1.315915203094482e+01 1.107159163802862e-02 3.306149691343307e-02 9.974335432052612e-01 -6.133582305908203e+01 0.000000000000000e+00 0.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00
After I studied the document "FS Coordinates (PDF)" at the following URL : http://freesurfer.net/fswiki/CoordinateSystems there are some question :
I thought in my case correspond to the page 3 condition "Two Fields of View, Two Coordinate Systems" But I am not sure that that does the affine matrix above is the vox_to_vox matrix or not. If it is, could I use it directly to compute the corresponding voxel location between 2 volume? or I need additional coordinate transfer?
Sincerely, Hsintsung