Yes, this is not unusual. The GTM (which RBV is based on) is just a GLM
where it uses a linear model to estimate the ROI intensities. There are
no constraints on positivity. When the model is not exactly correct
(which it never will be) it is possible to get negative values. Usually
this only happens in small ROIs and/or low uptake regions. Why might the
model not be perfect? It could be a number of things:
1. the MRI-PET registration is inaccurate
2. the MRI segmentation is inaccurate
3. the ideal segmentation is not reflected by the MRI segmentation (eg,
one of the segmentations has a hot spot)
It could also be the case that the OSEM positivity constraint causes
bias in low uptake regions (the GLM expects two-sided Gaussian noise)
On 8/28/20 7:27 PM, Albrecht, Daniel S. wrote:
Hello,
I’ve been running RBV partial volume correction on some amyloid PET
data recently using mri_gtmpvc, and for ~90% of subjects the resulting
rbv.nii.gz looks as expected with the exception of negative values in
the bilateral accumbens (see attached image, thresholded from 0-4000
in scanner units). Any ideas why this might be happening?
This is the command I’ve been using:
mri_gtmpvc --i $PETimg \
--reg $identity_mtrx \
--seg $FSdir/${subject}/mri/gtmseg.mgz \
--psf 6 --default-seg-merge --auto-mask PSF .01 --rbv \
--no-rescale --threads 8 \
--o $outname
Thanks very much,
Dan
_______________________________________________
Freesurfer mailing list
Freesurfer@nmr.mgh.harvard.edu
https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
_______________________________________________
Freesurfer mailing list
Freesurfer@nmr.mgh.harvard.edu
https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer