Hello FreeSurfer Developers,
I'm attempting to create partial volume (pve) corrected images for some of the subcortical segmentations. I have searched the list and found discussions regarding how the volumes for pve corrected data can be obtained by default with the asegstats2table command, or when using the mri_segstats command with the --pv flag. However, I am wondering if it is possible to save the mixed coefficient map of the pve correction process or save a aseg.mgz file that has been pve corrected?
Sincerely, Andrea
Hi Andrea
sorry, I don't think we do, although could add such a thing I guess, although the format isn't clear to me. I guess you would have to save multiple labels at each voxel with the mixing coefficient of each one.
cheers Bruce
On Tue, 25 Oct 2016, Andrea Shafer wrote:
Hello FreeSurfer Developers,
I'm attempting to create partial volume (pve) corrected images for some of the subcortical segmentations. I have searched the list and found discussions regarding how the volumes for pve corrected data can be obtained by default with the asegstats2table command, or when using the mri_segstats command with the --pv flag. However, I am wondering if it is possible to save the mixed coefficient map of the pve correction process or save a aseg.mgz file that has been pve corrected?
Sincerely, Andrea
Hi Bruce, This would be most helpful. And yes, I think multiple labels for each voxel that are the coefficients that explain the observed intensity of that voxel would be a good format. I am thinking that anyone could then correct for pve of structures taken from the aseg.mgz file in matlab (or python) by keeping or removing voxels on the border of the structure based on which tissue class as the highest coefficient. Does this seem reasonable to you? Is that how the --pv flag for mri_segstats works (i.e., by making the decision to keep or discard a voxel based on the tissue class with the highest coefficient)?
Many thanks, Andrea
Hi Andrea
sorry, I don't think we do, although could add such a thing I guess, although the format isn't clear to me. I guess you would have to save multiple labels at each voxel with the mixing coefficient of each one.
cheers Bruce
On Tue, Oct 25, 2016 at 8:57 PM, Andrea Shafer andrea.t.shafer@gmail.com wrote:
Hello FreeSurfer Developers,
I'm attempting to create partial volume (pve) corrected images for some of the subcortical segmentations. I have searched the list and found discussions regarding how the volumes for pve corrected data can be obtained by default with the asegstats2table command, or when using the mri_segstats command with the --pv flag. However, I am wondering if it is possible to save the mixed coefficient map of the pve correction process or save a aseg.mgz file that has been pve corrected?
Sincerely, Andrea
freesurfer@nmr.mgh.harvard.edu