Alright that should work. So do I use the pet.anat.nii files (from mri_vol2vol) as the input pet file to the mri_compute_volume_fractions? Also, has anyone published using this pv correction method that you know of?
An alternative method that I was thinking was to do something similar to the cortical thickness correction Salat used and sample on either side of the gray/white boundary and the gray/csf boundary and calculate the ratio. Do you think this could account for the effects of pv in pet?
Corinna
Actually, I think it only will affect the volume measurement for each structure, so that's no good. There is not a flag to mri_segstats that will do this. If you have a registration between your native pet data and the FS anat, then you can use
mri_compute_volume_fractions regfile.dat pet.nii pvf
This will create pvf.{gm,wm,csf}.mgz
You can then correct the pet with something like
fscalc pet.nii div pvf.gm.mgz pet.pvf.nii
You can map the aparc+aseg into the pet space using mri_label2vol (and the reg from above), then use mri_segstats with the files in the native pet space.
doug
corinna bauer wrote:
Hi Doug,
I will give that a try. How will this account for the p.v. in pet data? Does it remove the voxels that are mixed tissue classes? Does it account for partial volume effects between neighbouring grey matter structures?
Corinna
On Wed, Jun 8, 2011 at 1:06 PM, Douglas N Greve <greve@nmr.mgh.harvard.edu <mailto:greve@nmr.mgh.harvard.edu>> wrote:<mailto:Freesurfer@nmr.mgh.harvard.edu>
Hi Corinna,
try adding --pv mri/norm.mgz
doug
corinna bauer wrote:
Hi all,
I am wondering if partial volume is taken into account when
calculating average uptake values of pet data. I know that
aseg.stats takes partial volume into account and I am
wondering if there is something similar in place for pet data?
From what I gather, most people when using FreeSurfer correct
for partial volume effects by regressing out cortical
thickness, but no basis for this correction method has been
provided. Any insight would be appreciated.
n.b. to obtain the pet uptake values I used mri_segstats --seg
mri/aparc+aseg.mgz --sum pet.segstats.dat --i $petanat --ctab
rois.txt
Thank you,
Corinna
------------------------------------------------------------------------
_______________________________________________
Freesurfer mailing list
Freesurfer@nmr.mgh.harvard.edugreve@nmr.mgh.harvard.edu <mailto:greve@nmr.mgh.harvard.edu>
https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
-- Douglas N. Greve, Ph.D.
MGH-NMR Center
Phone Number: 617-724-2358 <tel:617-724-2358> Fax: 617-726-7422
<tel:617-726-7422> <http://surfer.nmr.mgh.harvard.edu/fswiki/BugReporting> <http://www.nmr.mgh.harvard.edu/facility/filedrop/index.html>
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.
--
Douglas N. Greve, Ph.D.
MGH-NMR Center
greve@nmr.mgh.harvard.edu
Phone Number: 617-724-2358 Fax: 617-726-7422
Bugs: surfer.nmr.mgh.harvard.edu/fswiki/BugReporting
FileDrop: www.nmr.mgh.harvard.edu/facility/filedrop/index.html