Hi Doug,
The geometry of the input is correct in both freeview and tkmedit, and the outlines of the white and pial surfaces align well with with volume in freeview. Looking at the headers, I'm wondering if my initial transform of the image into the proper coordinates might be the problem. The DICOM wasn't oriented correctly at first, so I converted to nifti, and ran "mri_convert -c -iid 0 -1 0 -ijd 0 0 -1 -ikd 1 0 0 mprage.nii.gz mprage.mgz", which produced a conformed image with the correct orientation in freeview. Is there a chance that a problem here could produce the problems later down the pipeline? For example, I see now that mprage.nii.gz's header correctly lists the primary slice direction as axial, whereas the reoriented mprage.mgz incorrectly lists it as coronal. Also, orientation was originally listed as RAS, but changed to LIA after conversion. Is there a chance that this is affecting the later surface registration?
Thanks so much for any help you can provide. I've attached mri_info output for the nifti and mgz files below.
Best regards,
John
mri_info mprage.nii.gz
Volume information for mprage.nii.gz
type: nii
dimensions: 199 x 152 x 183
voxel sizes: 1.0000, 1.0000, 1.0000
type: FLOAT (3)
fov: 199.000
dof: 0
xstart: -99.5, xend: 99.5
ystart: -76.0, yend: 76.0
zstart: -91.5, zend: 91.5
TR: 8.42 msec, TE: 0.00 msec, TI: 0.00 msec, flip angle: 0.00 degrees
nframes: 1
PhEncDir: UNKNOWN
ras xform present
xform info: x_r = 0.7324, y_r = 0.0361, z_r = 0.6799, c_r = 1.7748
: x_a = -0.0157, y_a = 0.9992, z_a = -0.0361, c_a = 25.8853
: x_s = -0.6807, y_s = 0.0157, z_s = 0.7324, c_s = -8.5852
Orientation : RAS
Primary Slice Direction: axial
voxel to ras transform:
0.7324 0.0361 0.6799 -136.0543
-0.0157 0.9992 -0.0361 -45.1875
-0.6807 0.0157 0.7324 -9.0720
0.0000 0.0000 0.0000 1.0000
voxel-to-ras determinant 0.999997
ras to voxel transform:
0.7324 -0.0157 -0.6807 92.7637
0.0361 0.9992 0.0157 50.2057
0.6799 -0.0361 0.7324 97.5155
0.0000 0.0000 0.0000 1.0000
______
mri_info mprage.mgz
Volume information for mprage.mgz
type: MGH
dimensions: 256 x 256 x 256
voxel sizes: 1.0000, 1.0000, 1.0000
type: UCHAR (0)
fov: 256.000
dof: 0
xstart: -128.0, xend: 128.0
ystart: -128.0, yend: 128.0
zstart: -128.0, zend: 128.0
TR: 8.42 msec, TE: 0.00 msec, TI: 0.00 msec, flip angle: 0.00 degrees
nframes: 1
PhEncDir: UNKNOWN
ras xform present
xform info: x_r = -1.0000, y_r = 0.0000, z_r = 0.0000, c_r = 1.7748
: x_a = 0.0000, y_a = 0.0000, z_a = 1.0000, c_a = 25.8853
: x_s = 0.0000, y_s = -1.0000, z_s = 0.0000, c_s = -8.5852
talairach xfm :
Orientation : LIA
Primary Slice Direction: coronal
voxel to ras transform:
-1.0000 0.0000 0.0000 129.7748
0.0000 0.0000 1.0000 -102.1147
0.0000 -1.0000 0.0000 119.4148
0.0000 0.0000 0.0000 1.0000
voxel-to-ras determinant -1
ras to voxel transform:
-1.0000 0.0000 0.0000 129.7748
-0.0000 -0.0000 -1.0000 119.4148
-0.0000 1.0000 -0.0000 102.1147
0.0000 0.0000 0.0000 1.0000