Re: [Freesurfer] Question about the negative value of surface area

2015-10-01 Thread Bruce Fischl

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.


[Freesurfer] Question about the negative value of surface area

2015-09-29 Thread chenhf_uestc
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.