Hi. We have a question about the raw cortical thickness (CT) data exported from FS with read_curv. As an example, I loaded in the “bert” sample dataset, read the left hand cortical
thicknesses and then dumped to a text file (attached).
This comes straight from the FS Wiki (https://surfer.nmr.mgh.harvard.edu/fswiki/UserContributions/FAQ)
Q. How can I make a histogram of cortical thickness?
A: In Matlab you can use the read_curv() function, as in the example below.
- thick = read_curv ('/usr/local/freesurfer/subjects/bert/surfer/lh.thickness');
hist (thick,100);
Notice that there are many zero values that refer to areas where there's no cortical surface. The best procedure is to create a new array without the zero values and then make the histogram.
The matlab code are just in these three lines.
thick = read_curv(bert);
fileID2 = fopen(outfile,'w');
fprintf(fileID2,’%s \n’,thick);
The first image below shows the first 500 thicknesses in the array while the second image shows the first 5,000. There is a distinct oscillation/periodicity in the data that
is written out. We were wondering if this has to do with the way that FS rasters across the brain going from sulci to gyri, etc? If the data is placed in a histogram, in a normal control subject we do get and expect a Gaussian distribution.
Please let us know!
Thanks,
Jon
Jonathan P. Dyke, Ph.D., DABMP (MRI)
Professor of Physics Research in Radiology
Weill Cornell Medicine
Citigroup Biomedical Imaging Center
516 East 72nd Street,
New York, NY 10021
T: (212) 746-5781
https://cbic.weill.cornell.edu
https://radiology.weill.cornell.edu/research/mri-research-institute/jonathan-dyke-laboratory
https://vivo.weill.cornell.edu/display/cwid-jpd2001