On 01/28/2016 06:50 PM, Pradeep wrote:
> Hello Doug,
>
> I have used the gtmseg with --keep-cc flag and the corresponding ctab
> files showed the labels but the mri_gtmpvc step failed.
> ****
> Loading seg for gtm gtmseg.mgz
> Loading seg ctab gtmseg.ctab
> Reading gtmseg.lta
> Replacing 18
> ERROR: CheckSegTissueType() no entry for seg 192
> Failed tissue type check
> ****
What is your mri_gtmpvc command line? What is the rest of the terminal
output?
> My objective is to use the combination of all CC's as a reference
> region and obtain the PVC results, which would be listed in gtm.stats.dat
It will be best to combine them when running mri_gtmpvc using --replace,
eg, --replace 252 251 --replace 253 251 --replace 254 251 --replace 255 251
this will cause all segments of the CC to appear to be a single segment
(251).
>
> Also, I read in the previous email discussions that the default
> ref-region Pons is 'PVC'ed'. So if I want to calculate the SUVR with
> another ROI as a reference region,
> would it be OK take a ratio of the ROI's in gtm.stats.dat table.
Yes, or you can spec the new region, eg --rescale 251
>
> Thanks,
> Pradeep
>
> On Mon, Jan 11, 2016 at 9:25 AM, Douglas N Greve
> <mailto:Freesurfer@nmr.mgh.harvard.edu>> <greve@nmr.mgh.harvard.edu <mailto:greve@nmr.mgh.harvard.edu>> wrote:
>
>
> If you want to use partial volume correction, then you are better off
> using mri_gtmpvc with the bbr registration, something like
>
> 1. To start, run
>
> gtmseg --s subject
>
> This will take a couple of hours and produces some files needed
> for GTM
> PVC (which is used for GTM, MG, RBV).
>
> 2. You'd then register the PET to the anatomical with bbregister
> (probably with --t2 weighting). Make sure to save the output as an LTA
> (--lta). I usually use the mean TAC as the input. You can do this in
> parallel with #1.
>
> 3. You'd then run mri_gtmpvc, something like
>
> mri_gtmpvc --i pet.nii.gz --psf PSF --auto-mask PSF+2 .01 --seg
> gtmseg.mgz
> --reg reg.lta --default-seg-merge --o gtmpvc.output
>
> PSF is the point-spread FWHM of the scanner; reg.lta is the
> registration from #2. By default, this will scale by pons. The output
> will be gtm.stats.dat and gtm.nii.gz. They both basically have the
> same information. gtm.stats.dat is an easy to read text file. Where
> each row is an ROI, something like:
>
> 9 17 Left-Hippocampus subcort_gm 473
> 174.083 1.406 0.1216
>
> 9 = nineth row
> 17 = index for RO
> Left-Hippocampus = name of ROI
> subcort_gm = tissue class
> 473 = number of PET voxels in the ROI
> 174 = variance reduction factor for ROI (based on GLM/SGTM)
> 1.406 = PVC uptake of ROI relative to Pons
> 0.1216 = resdiual varaince across voxels in the ROI
>
> gtm.nii.gz is a nifti file with each "voxel" being an ROI. The value
> is the PVC uptake of ROI relative to Pons. These can easily be
> concatenated together (mri_concat) and used as input to mri_glmfit
> for group analysis.
>
>
>
>
> On 01/08/2016 04:22 AM, Benjamin Spurny wrote:
> > Dear Freesurfer experts!
> >
> > I am currently working on PET analysis using FS
> >
> > I coregistered my PET with the processed MR using bbregister,
> > transfered it to a surface using mri_vol2surf
> > and now createt an overlay in freeview with the lh.inflated and
> used the
> > labels from the lh.aparc.a2009s.annot file.
> >
> > In freeview i get the corresponding BP value for each vertex now but
> > is there a way to get a list of vertices with the corresponding
> BP value
> > and the corresponding ROI this vertex belongs to?
> > Or is there a better to do this analyis?
> >
> > Many thanks in advance!
> >
> > Benjamin
> >
> > _______________________________________________
> > Freesurfer mailing list
> > Freesurfer@nmr.mgh.harvard.edu
> > https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
> >
> >
>
> --
> Douglas N. Greve, Ph.D.
> MGH-NMR Center
> greve@nmr.mgh.harvard.edu <mailto:greve@nmr.mgh.harvard.edu>
> Phone Number: 617-724-2358 <tel:617-724-2358>
> Fax: 617-726-7422 <tel:617-726-7422>
>
> Bugs: surfer.nmr.mgh.harvard.edu/fswiki/BugReporting
> <http://surfer.nmr.mgh.harvard.edu/fswiki/BugReporting>
> FileDrop: https://gate.nmr.mgh.harvard.edu/filedrop2
> www.nmr.mgh.harvard.edu/facility/filedrop/index.html
> <http://www.nmr.mgh.harvard.edu/facility/filedrop/index.html>
> Outgoing:
> ftp://surfer.nmr.mgh.harvard.edu/transfer/outgoing/flat/greve/
>
> _______________________________________________
> Freesurfer mailing list
> Freesurfer@nmr.mgh.harvard.edu <mailto:Freesurfer@nmr.mgh.harvard.edu>
> https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
>
>
> 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.
>
>
>
>
> _______________________________________________
> Freesurfer mailing list
> Freesurfer@nmr.mgh.harvard.edu
> https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
--
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: https://gate.nmr.mgh.harvard.edu/filedrop2
www.nmr.mgh.harvard.edu/facility/filedrop/index.html
Outgoing: ftp://surfer.nmr.mgh.harvard.edu/transfer/outgoing/flat/greve/
_______________________________________________
Freesurfer mailing list
Freesurfer@nmr.mgh.harvard.edu
https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer