Hi Feng
the area of the triangls is oriented to detected changes in the jacobian.
That said, they should be positive definite (you can verify
thismris_remove_intersection on the ?h.sphere.reg to make sure that this is
the case. Perhaps Anderson or Doug can comment?
cheers
Bruce
On Tue, 29 Sep 2015, chenhf_uestc wrote:
Dear Freesurfer experts,
I read the raw area file and the area resampled to the fsaverage space and
found several interesting things.
raw58l = read_curv('E:\fs_result\NC058\surf\lh.area'); % the range of the
area value is [-0.0471 6.2118]
raw58r = read_curv('E:\fs_result\NC058\surf\rh.area'); % the range of the
area value is [0 6.8963]
raw59l = read_curv('E:\fs_result\NC059\surf\lh.area'); % the range of the
area value is [-0.3347 5.4892]
raw59r = read_curv('E:\fs_result\NC059\surf\rh.area'); % the range of the
area value is [-1818 5.4269]
The first question, what is the unit of the area file in the surf folder
(such as lh.area), mm2? Why the minimal value of the area can be less than
zero?
Subsequently, I perform further experiments, reading the resampled area
file.
resamp_58l = load_mgh('E:\fs_result\NC058\surf\lh.area.fsaverage.mgh'); %
the range of the area value is [-0.0118 7.4073]
resamp_58r = load_mgh('E:\fs_result\NC058\surf\rh.area.fsaverage.mgh'); %
the range of the area value is [0 19.0191]
resamp_59l = load_mgh('E:\fs_result\NC059\surf\lh.area.fsaverage.mgh'); %
the range of the area value is [-0.3937 9.9084]
resamp_59r = load_mgh('E:\fs_result\NC059\surf\rh.area.fsaverage.mgh'); %
the range of the area value is -0.0909 8.2306[]
The second question, why the value of some area file change a lot? For
example, the maximum value of raw58r is 6.8963, but 19.0191 of resampled
one.
Next,
resamp_smooth5_58l =
load_mgh('E:\fs_result\NC058\surf\lh.area.fwhm5.fsaverage.mgh'); % the
range of the area value is [0 1.7293]
resamp_smooth5_58r =
load_mgh('E:\fs_result\NC058\surf\rh.area.fwhm5.fsaverage.mgh'); % the
range of the area value is [0 2.1877]
resamp_smooth5_59l =
load_mgh('E:\fs_result\NC059\surf\lh.area.fwhm5.fsaverage.mgh'); % the
range of the area value is [0 1.2875]
resamp_smooth5_59r =
load_mgh('E:\fs_result\NC059\surf\rh.area.fwhm5.fsaverage.mgh'); % the
range of the area value is [0 1.3542]
The third question, why all the minimal values after smooth was zero? Why
all the maximal values after smooth reduce a lot?
Finally,
resamp_smooth25_58l =
load_mgh('E:\fs_result\NC058\surf\lh.area.fwhm25.fsaverage.mgh'); % the
range of the area value is [0 1.2014]
resamp_smooth25_58r =
load_mgh('E:\fs_result\NC058\surf\rh.area.fwhm25.fsaverage.mgh'); % the
range of the area value is [0 1.2611]
resamp_smooth25_59l =
load_mgh('E:\fs_result\NC059\surf\lh.area.fwhm25.fsaverage.mgh'); % the
range of the area value is [0 0.9730]
resamp_smooth25_59r =
load_mgh('E:\fs_result\NC059\surf\rh.area.fwhm25.fsaverage.mgh'); % the
range of the area value is [0 1.0119]
The final question, why all the maximal values of fwhm25 smooth are less
than the values of fwhm5 smooth?
How can I get the area value in the recon-all analyses? The cortical
thickness is the distance betwwen the white matter and the pial, but what is
the procedures the area value obtained? If there any reference? My
freesurfer version is stable v5.3.0.
Any help would be greatly appreciately!
Best,
Feng
___
Freesurfer mailing list
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.