Hello FreeSurfer Developers,
I have run mri_gtmpvc command to perform PET PVC (partial volume correction) where input data has 4 frames (dynamic). The output volume as gtm.nii.gz has 4 frames as well which means that PVC was performed for each frame separately. However, gtm.stats.dat file has one value for each ROI. At first I thought it was the average of all frames but after manually comparing some ROIs it seems that the values in the file are taken only from the first frame of dynamic data. Is it the case or those values have another interpretation? Maybe it is the highest value across all frames? In addition, is there an easy way to get the same file where the values would be the average across all frames?
It is the first frame. You'll have to get the average from the gtm.nii.gz file. You can do it easily enough in matlab with
gtm = fast_vol2mat(MRIread('gtm.nii.gz'));
% gtm will be a 4xNrois matrix
mean(gtm)
On 12/09/2016 03:21 AM, Donatas Sederevicius wrote:
Hello FreeSurfer Developers,
I have run mri_gtmpvc command to perform PET PVC (partial volume correction) where input data has 4 frames (dynamic). The output volume as gtm.nii.gz has 4 frames as well which means that PVC was performed for each frame separately. However, gtm.stats.dat file has one value for each ROI. At first I thought it was the average of all frames but after manually comparing some ROIs it seems that the values in the file are taken only from the first frame of dynamic data. Is it the case or those values have another interpretation? Maybe it is the highest value across all frames? In addition, is there an easy way to get the same file where the values would be the average across all frames?
Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
freesurfer@nmr.mgh.harvard.edu
-
Donatas Sederevicius -
Douglas N Greve