There is no one command that will do it. You can do it in matlab, something like
a = MRIread('spike.mgh');
a.vol = a.vol/a.vol(yourindex);
MRIwrite(a,'spike.rescale.mgh');
Or you can get the vertices and values with
mri_convert spike.mgh --ascii+crsf
to get all the vertices and their values
On 1/4/17 10:27 AM, Song, Daniel wrote:
Hi all,
I'm trying to create a label by making a circle with a fixed radius around a defined vertex. I'm trying to use the method described here:
https://www.mail-archive.com/freesurfer@nmr.mgh.harvard.edu/msg44098.html
I made the initial surface overlay with my vertex of interest having a value of 1 and 0 everywhere else with:
mri_volsynth --pdf delta --delta-crsf 108555 0 0 0 --o pBECTS007_spike.mgh --template rh.thickness
Then I smoothed with fwhm:
mri_surf2surf --s pBECTS007 --hemi rh --sval pBECTS007_spike.mgh --fwhm 10 --tval ./pBECTS007_spike_fwhm.mgh
My question is, is there an easy way to divide all vertices by the value of my vertex of interest after smoothing in tksurfer? If not, is there a command that will allow me to extract the vertex values and indices to do so?
Thank you!
Dan Song Clinical Research Coordinator Massachusetts General Hospital Neurology Department
Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer