on reflection, you're probably better off saving in mgh format:
matlab
help read_curv
[curv, fnum] = read_curv(fname) reads a binary curvature file into a vector
[thickness, f] = read_curv('lh.thickness') ; [area, f] = read_curv('lh.area') ; volume = thickness .* area; save_mgh(volume, 'lh.volume.mgz', eye(4)); more on thickness
thickness = 1.5533 1.8863 2.184 . . .
area
area = 0.27881 0.56441 0.8965 0.36047
volume
volume = 0.43309 1.0647 1.958 0.77892 0.81443 0.37083 cheers, Bruce
On Thu, 30 Aug 2007, Juergen Haenggi wrote:
Dear Freesurfer users
We would like to analyze cortical volumes vertex-wise. We recognized that ?h.volume are not produced by the recon-all command. I have read the questions posted some days ago by Rafael Huber concerning the same problem. One proposed solution was to multiply ?h.area times ?h.thickness at each vertices. To do this read_curv.m and the write_curv.m functions in MATLAB were suggested. We tried it but failed maybe due to our insufficient knowledge in MATLAB programming.
The first step (read_curv.m) works and produced a vector of curvature (small negative values) instead a vector of thickness (positive values between 1.5 and 3.5 mm) as we expected. What happens here?
The command write_curv.m does not work at all.
Is there any short documentation about multiplying images in MATLAB?
We appreciate any suggestions about that issue.
Thanks in advance Best regards
Juergen
Juergen Haenggi Ph.D. (Dr. des.) Division of Neuropsychology Institute of Psychology University of Zurich Binzmuehlestrasse 14, PO Box 25 8050 Zurich, Switzerland 0041 44 635 73 97 (phone office) 0041 76 445 86 84 (phone mobile) 0041 44 635 74 09 (fax office) BIN 4.D.04 (office room number) j.haenggi@psychologie.uzh.ch (email) http://www.psychologie.uzh.ch/neuropsy/ (website) http://www.juergenhaenggi.ch (private website)
Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer