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